fit is initialized as a vector of integers. How can you assign an mle
fit to an element of an integer vector?
Initialize fit as a list, use lapply, or whatever. Have you read An
Intro to R (ships with R) or other R (e.g. web) tutorial? This looks
like the sort of basic misunderstanding that one who has not made much
effort to learn how R works would make.If not, please do so before
posting further.
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
(650) 467-7374
"Data is not information. Information is not knowledge. And knowledge
is certainly not wisdom."
H. Gilbert Welch
On Sat, Feb 15, 2014 at 5:38 AM, John Hodgson <john at formby.plus.com>
wrote:> I have a 7-parameter model to fit using mle. I would like to generate
> fits for all
> pairs of parameters (with others fixed)
>
> The following code looked like it should work:
>
> library(stats4)
>
> # dummy mll function for sake of example
> mll = function (lnk=1.5,lnhs=-5,lnhi=-5,lnss=-5,lnsi=-5,lnws=-5,lnwi=-5)
> lnk^2 + lnhs^2 + lnhi^2+ lnss^2 + lnsi^2+ lnws^2 + lnwi^2
>
> fit=1:6
> pars = list(lnk=1.5,lnhs=-5,lnhi=-5,lnss=-5,lnsi=-5,lnws=-5,lnwi=-5)
>
> for (i in 2:7) {
> fit[i] = mle(mll,start=pars[c(1,i)],fixed>
pars[-c(1,i)],method="Nelder")
> }
>
>
> but it gives the following error:
> }
> Error in fit[i] = mle(mll, start = pars[c(1, i)], fixed = pars[-c(1,
> i)], :
> incompatible types (from S4 to integer) in subassignment type fix
>
> What have I missed?
>
> Thanks for any suggestions
>
> John Hodgson
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at 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.