Hi everybody
I am trying to find contrast in MANOVA.
I used next code
contrasts(ffage)<-ctr
contrasts(ffage)
MANOVA.agec<-manova(Y1~ffage,data=vol18.df)
summary(MANOVA.agec, split =list (ffage=list("0-17 v over 18"=0,
"18-25 v
over 26"=1, "26-31 v over 32"=2, "32-42 v over 43"=3,
"43-65 v 66+"=4)))
But the output was only the overall fit
> contrasts(ffage)<-ctr
> contrasts(ffage)
[,1] [,2] [,3] [,4] [,5]
0 5 0 0 0 0
1 -1 4 0 0 0
2 -1 -1 3 0 0
3 -1 -1 -1 2 0
4 -1 -1 -1 -1 1
5 -1 -1 -1 -1 -1> MANOVA.agec<-manova(Y1~ffage,data=vol18.df)
> summary(MANOVA.agec, split =list (ffage=list("0-17 v over 18"=0,
"18-25 v
> over 26"=1, "26-31 v over 32"=2, "32-42 v over
43"=3, "43-65 v 66+"=4)))
Df Pillai approx F num Df den Df Pr(>F)
ffage 5 0.30607 2.1543 20 520 0.002681 **
Residuals 130
---
Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1
I try to use linearHypothesis, but i couldnt find anything, maybe I didnt
use in a corect way, *any advice is welcome*
I always get the error
Error in solve.default(wcrossprod(model.matrix(model), w = wts)) :
Lapack routine dgesv: system is exactly singular>
Thanks in advance
Paola
--
View this message in context:
http://r.789695.n4.nabble.com/Contrasts-in-manova-tp4648306.html
Sent from the R help mailing list archive at Nabble.com.