
Hello everybody

 

After doing a MANOVA on a bunch of data, I want to be able to make some comment on the amount of variation in the data that is explained by the factor of interest. I want to say this in the following way: XX% of the data is explained by A.

 

I can acheive something like what I want by doing the following: 

 

 

X <- structure(c(9, 6, 9, 3, 2, 7), .Dim = as.integer(c(3, 2)))
Y <- structure(c(0, 2, 4, 0), .Dim = as.integer(c(2, 2)))
Z <- structure(c(3, 1, 2, 8, 9, 7), .Dim = as.integer(c(3, 2)))

U <- rbind(X,Y,Z)
m <- manova(U~as.factor(rep(1:3, c(3, 2, 3))))

summary(m,test="Wilks")

SS<-summary(m)$SS

(a<-mean(SS[[1]]/(SS[[1]]+SS[[2]])))

 

and concluding that 94% of variation is explained.

 

Is my desire misguided? If it is a worthy aim, is this a valid way of acheiving it?

 

Thanks a lot!

 

Sam

 

 

 

Samuel Brown

Research assistant

Bio-Protection Research Centre

PO Box 84

Lincoln University

Lincoln 7647

Canterbury

New Zealand

sam.brown@lincolnuni.ac.nz

http://www.the-praise-of-insects.blogspot.com

 

 

 

 

 

 
 		 	   		  
_________________________________________________________________

ws Live 

	[[alternative HTML version deleted]]

