[R] How to test combined effects?
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.
More information about the R-help