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.>
Can someone please help?
--
David Winsemius, MD
West Hartford, CT
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
David Winsemius, MD
West Hartford, CT
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
