[R] Test for equality of coefficients in multivariate multiple regression
Andrew Robinson
A.Robinson at ms.unimelb.edu.au
Wed Jul 19 03:51:44 CEST 2006
Hi Uli,
I suggest that you try to rewrite the model system into a single
mixed-effects model, which would allow direct parameterization of the
tests that you're interested in. A useful article, which may be
overkill for your needs, is:
Hall, D.B. and Clutter, M. (2004). Multivariate multilevel nonlinear
mixed effects models for timber yield predictions, Biometrics, 60:
16-24.
See the publications link on Daniel Hall's website:
http://www.stat.uga.edu/~dhall/
Cheers
Andrew
On Tue, Jul 18, 2006 at 08:15:12PM +0200, Ulrich Keller wrote:
> Hello,
>
> suppose I have a multivariate multiple regression model such as the
> following:
>
> > DF<-data.frame(x1=rep(c(0,1),each=50),x2=rep(c(0,1),50))
> > tmp<-rnorm(100)
> > DF$y1<-tmp+DF$x1*.5+DF$x2*.3+rnorm(100,0,.5)
> > DF$y2<-tmp+DF$x1*.5+DF$x2*.7+rnorm(100,0,.5)
> > x.mlm<-lm(cbind(y1,y2)~x1+x2,data=DF)
> > coef(x.mlm)
> y1 y2
> (Intercept) 0.07800993 0.2303557
> x1 0.52936947 0.3728513
> x2 0.13853332 0.4604842
>
> How can I test whether x1 and x2 respectively have the same effect on y1
> and y2? In other words, how can I test if coef(x.mlm)[2,1] is
> statistically equal to coef(x.mlm)[2,2] and coef(x.mlm)[3,1] to
> coef(x.mlm)[3,2]? I looked at linear.hypothesis {car} and glh.test
> {gmodels}, but these do not seem the apply to multivariate models.
> Thank you in advance,
>
> Uli Keller
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
--
Andrew Robinson
Department of Mathematics and Statistics Tel: +61-3-8344-9763
University of Melbourne, VIC 3010 Australia Fax: +61-3-8344-4599
Email: a.robinson at ms.unimelb.edu.au http://www.ms.unimelb.edu.au
More information about the R-help
mailing list