[R] How to test combined effects?

Dieter Menne dieter.menne at menne-biomed.de
Tue Oct 30 14:08:35 CET 2007


Gang Chen <gangchen <at> mail.nih.gov> writes:

> 
> Suppose I have a mixed-effects model where yij is the jth sample for  
> the ith subject:
> 
> yij= beta0 + beta1(age) + beta2(age^2) + beta3(age^3) + beta4(IQ) +   
> beta5(IQ^2) + beta6(age*IQ)  + beta7(age^2*IQ)  +  beta8(age^3 *IQ)  
> +random intercepti + eij
> 
> In R how can I get an F test against the null hypothesis of  
> beta6=beta7=beta8=0? In SAS I can run something like contrast  age*IQ  
> 1, age^2*IQ 1, age^3 *IQ 1, but is there anything similar in R?

Check packages multcomp and gmodels for contrast tests that work with lme.

Dieter



More information about the R-help mailing list