[R-sig-ME] qqnorm question

Doran, Harold HDoran at air.org
Fri Mar 30 20:53:18 CEST 2007


Here is another way to get the qqplot to evaluate the random effects

data(star, package='mlmRev')
fm <- lmer2(math ~ yrs + (yrs|sch), star)

qqmath(~ranef(fm)$sch$'(Intercept)')
qqmath(~ranef(fm)$sch$yrs)

Don't forget about mcmcsamp() and the ancillary functions in coda for
also evaluating distributional assumptions.

Harold


> -----Original Message-----
> From: r-sig-mixed-models-bounces at r-project.org 
> [mailto:r-sig-mixed-models-bounces at r-project.org] On Behalf 
> Of Roberts, J. Kyle
> Sent: Friday, March 30, 2007 1:29 PM
> To: Chaudhari, Monica; r-sig-mixed-models at r-project.org
> Subject: Re: [R-sig-ME] qqnorm question
> 
> That's what I want. Thanks!
> 
> So the attached picture would be bad, right?!?!? ;)
> 
> Kyle
> 
> 
> ***************************************
> J. Kyle Roberts, Ph.D.
> Baylor College of Medicine
> Center for Educational Outreach
> One Baylor Plaza, MS:  BCM411
> Houston, TX   77030-3411
> 713-798-6672 - 713-798-8201 Fax
> jkrobert at bcm.edu
> ***************************************
> 
> -----Original Message-----
> From: Chaudhari, Monica [mailto:mchaudhari at deltadentalwa.com]
> Sent: Friday, March 30, 2007 11:31 AM
> To: Roberts, J. Kyle; r-sig-mixed-models at r-project.org
> Subject: RE: [R-sig-ME] qqnorm question
> 
> Hello Kyle,
> 
> Are you trying to test the normality of random effects? If 
> yes, then below is the way to go:
> str(re<-ranef(modl,postVar = TRUE))
> qqmath(re)
> or qqmath(~re[[1]]); qqmath(~re[[2]]);
> 
> Hope this helps.
> 
> Thanks,
> Monica
> 
> 
> -----Original Message-----
> From: r-sig-mixed-models-bounces at r-project.org
> [mailto:r-sig-mixed-models-bounces at r-project.org] On Behalf 
> Of Roberts, J. Kyle
> Sent: Friday, March 30, 2007 8:56 AM
> To: r-sig-mixed-models at r-project.org
> Subject: [R-sig-ME] qqnorm question
> 
> Dear All,
>  
> Forgive my lack of knowledge, but I can't get qqnorm to work 
> with lmer.
> I get the following error (dep.lmer is a lmer class object):
>  
> > qqnorm(dep.lmer)
> Error in min(..., na.rm = na.rm) : invalid 'type' (list) of 
> argument In addition: Warning message:
> is.na() applied to non-(list or vector) in: is.na(y)
>  
> Any ideas?
>  
> Kyle
>  
> ***************************************
> J. Kyle Roberts, Ph.D.
> Baylor College of Medicine
> Center for Educational Outreach
> One Baylor Plaza, MS:  BCM411
> Houston, TX   77030-3411
> 713-798-6672 - 713-798-8201 Fax
> jkrobert at bcm.edu
> ***************************************
>  
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list 
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
> 
> 
> #########################################################
> The information contained in this e-mail and subsequent 
> attachments may be privileged, confidential and protected 
> from disclosure.  This transmission is intended for the sole 
> use of the individual and entity to whom it is addressed.  If 
> you are not the intended recipient, any dissemination, 
> distribution or copying is strictly prohibited.  If you think 
> that you have received this message in error, please e-mail 
> the sender at the above e-mail address.
> #########################################################
> 
>




More information about the R-sig-mixed-models mailing list