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.
______________________________________________
R-help@r-project.org 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.
