Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-21 Thread Douglas Bates
.11.1 tcltk_2.11.1  tools_2. > > ) > > > >>> >>> Whoever suggested using the verbose option to see the progress of the >>> iterations has the right idea. >>> >>> >>> On Mon, Sep 20, 2010 at 2:50 PM, array chip wrote: >>>>

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-21 Thread Peter Dalgaard
chip wrote: >>> Thank you Peter and Ben for your comments. >>> >>> John >>> >>> >>> - Original Message ---- >>> From: Peter Dalgaard >>> To: array chip >>> Cc: r-help@r-project.org; r-sig-mixed-mod...@r-project

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-21 Thread Douglas Bates
ts. >> >> John >> >> >> - Original Message >> From: Peter Dalgaard >> To: array chip >> Cc: r-help@r-project.org; r-sig-mixed-mod...@r-project.org >> Sent: Mon, September 20, 2010 12:28:43 PM >> Subject: Re: [R] lmer() vs.

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-21 Thread Douglas Bates
ter Dalgaard > To: array chip > Cc: r-help@r-project.org; r-sig-mixed-mod...@r-project.org > Sent: Mon, September 20, 2010 12:28:43 PM > Subject: Re: [R] lmer() vs. lme() gave different variance component estimates > > On 09/20/2010 08:09 PM, array chip wrote: >> Thank you Pe

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-20 Thread array chip
Thank you Peter and Ben for your comments. John - Original Message From: Peter Dalgaard To: array chip Cc: r-help@r-project.org; r-sig-mixed-mod...@r-project.org Sent: Mon, September 20, 2010 12:28:43 PM Subject: Re: [R] lmer() vs. lme() gave different variance component estimates

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-20 Thread Peter Dalgaard
), and I would expect glht() to give meaningful results. Do check the residuals for at least approximate normality, though. > > Thank you very much! > > John > > > > > > - Original Message > From: Peter Dalgaard > To: array chip > Cc: r-hel

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-20 Thread array chip
org; r-sig-mixed-mod...@r-project.org Sent: Sat, September 18, 2010 1:35:45 AM Subject: Re: [R] lmer() vs. lme() gave different variance component estimates For a nested design, the relation is quite straightforward: The residual MS are the variances of sample means scaled to be comparable with th

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-18 Thread Peter Dalgaard
On 09/17/2010 10:50 PM, array chip wrote: > > > Thank you Peter. Actually 3 people from mixed model mailing list tried my > code > using lmer(). They got the same results as what I got from lme4(). So they > couldn't replicate my lmer() results: > > Random effects: > Groups NameVari

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-17 Thread Kevin Wright
gt; > John > > > > > > > ----- Original Message > From: Peter Dalgaard > To: array chip > Cc: r-help@r-project.org > Sent: Fri, September 17, 2010 1:05:27 PM > Subject: Re: [R] lmer() vs. lme() gave different variance component > estimates > >

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-17 Thread array chip
: Fri, September 17, 2010 1:05:27 PM Subject: Re: [R] lmer() vs. lme() gave different variance component estimates On 09/17/2010 09:14 PM, array chip wrote: > Hi, I asked this on mixed model mailing list, but that list is not very > active, > > so I'd like to try the general R ma

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-17 Thread Peter Dalgaard
On 09/17/2010 09:14 PM, array chip wrote: > Hi, I asked this on mixed model mailing list, but that list is not very > active, > so I'd like to try the general R mailing list. Sorry if anyone receives the > double post. > > > Hi, I have a dataset of animals receiving some eye treatments. There