Re: [R] Help with optim() to maximize log-likelihood

2015-03-10 Thread Prof J C Nash (U30A)
1) It helps to include the require statements for those of us who work outside your particular box. lme4 and (as far as I can guess) fastGHQuad are needed. 2) Most nonlinear functions have domains where they cannot be evaluated. I'd be richer than Warren Buffett if I got $5 for each time someon

Re: [R] Help with optim() to maximize log-likelihood

2015-03-09 Thread Sophia Kyriakou
yes Ben, this works indeed! Thanks a million!! On Mon, Mar 9, 2015 at 7:17 PM, Ben Bolker wrote: > Sophia Kyriakou gmail.com> writes: > > > > > hello, I am using the optim function to maximize the log likelihood of a > > generalized linear mixed model and I am trying to replicate glmer's > > es

Re: [R] Help with optim() to maximize log-likelihood

2015-03-09 Thread Ben Bolker
Sophia Kyriakou gmail.com> writes: > > hello, I am using the optim function to maximize the log likelihood of a > generalized linear mixed model and I am trying to replicate glmer's > estimated components. If I set both the sample and subject size to q=m=100 > I replicate glmer's results for the