[R] Contrasts in manova

paola aikaterini07 at hotmail.com
Sat Nov 3 15:49:49 CET 2012


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.



More information about the R-help mailing list