On Sep 5, 2011, at 11:41 AM, David Winsemius wrote:

>

> On Sep 5, 2011, at 10:30 AM, colstat wrote:

>

>> Hi,

>> First time using the optim(), can someone please tell me what I am

>> doing

>> wrong? The error looks like this

>>

>> Error in .Internal(pnorm(q, mean, sd, lower.tail, log.p)) :

>> 'sd' is missing

>>

>

> You should be using a textbook. Which one are you consulting?

>

>>

>> An example of the error

>> dat = c(20, 19, 9, 8, 7, 4, 3, 2)

>> dat_mu=mean(dat)

>> dat_s=sd(dat)

>>

>> max.func = function(dat, mu, sd) {

>> pnorm(dat, mu, sd)

>> }

>>

>> optim(fn=max.func, dat=dat, par=c(mu=dat_mu, s=dat_s))

>>

>> I get sd is missing error. If I wrote par=c(s=dat_s, mu=dat_mu) ,

>> then it

>> tells me mu is missing.

>

> You constructed a function with two

< I meant to say that the function had three arguments (none of them

"par") and you were incorrectly assuming the function would be able to

find what was inside par despite not passing it.>

> arguments but optim only passes two argument to the objective

> function. Build your objective function a) so that it accepts two

> arguments

<And you don't actually need to use two args. You could set it up so

that your max.fun gets the "dat" data from the calling environment.>

> and b) so that it returns one value. At the moment it does not do

> either.

<But I remain convinced you need to be reading something didactic. I

suggest Bolker's "Ecological Models and Data in R" ch 6 or Venables

and Ripley's MASS4 ch 16. Bolker's text seems to start the reader out

a bit more gently. V&R's first example is a mixture model. I have a

text I looked at Crawley's "R Book but it has no examples of the use

of 'optim'. It did have a tiny example of using 'optimize' that I

found difficult to follow because he uses attache and the data

variables had been defined pages earlier.>

David Winsemius, MD

West Hartford, CT

______________________________________________

[hidden email] mailing list

https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide

http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.