Re: [R] using mle2 for multinomial model optimization

2010-02-12 Thread Ben Bolker
Ravi Varadhan jhmi.edu> writes: > > > Use "Nelder-Mead" as the method in optim. This will give you the correct solution. > > m1 <- mle2(mfun, start=svec, vecpar=TRUE, fixed=svec[1:(l-1)], > method="Nelder") > > Hope this is helpful, > Ravi. > {Resending because it doesn't seem to have g

Re: [R] using mle2 for multinomial model optimization

2010-02-12 Thread Ben Bolker
Ravi Varadhan jhmi.edu> writes: > > > Use "Nelder-Mead" as the method in optim. This will give you the correct solution. > > m1 <- mle2(mfun, start=svec, vecpar=TRUE, fixed=svec[1:(l-1)], > method="Nelder") > > Hope this is helpful, > Ravi. Well-meaning though it is, this solution doesn'

Re: [R] using mle2 for multinomial model optimization

2010-02-12 Thread Ravi Varadhan
radhan, Ph.D. Assistant Professor, Division of Geriatric Medicine and Gerontology School of Medicine Johns Hopkins University Ph. (410) 502-2619 email: rvarad...@jhmi.edu - Original Message - From: Benedikt Gehr Date: Friday, February 12, 2010 8:26 am Subject: [R] using mle2 for multino

[R] using mle2 for multinomial model optimization

2010-02-12 Thread Benedikt Gehr
Hi there I'm trying to find the mle fo a multinomial model ->*L(N,h,S¦x)*. There is only *N* I want to estimate, which is used in the number of successes for the last cell probability. These successes are given by: p^(N-x1-x2-...xi) All the other parameters (i.e. h and S) I know from somewhere