[R-sig-eco] MANOVA

[Allan Edelsparre]

It was suggested that I post the code I used to construct the model

This is how it was constructed:


time <- factor(rep(c("rep1", "rep2"), c(36, 36)),levels=c("rep1", "rep2"))
morph <- ordered(rep(1:36, 2))
idata <- data.frame(time, morph)

model2<- lm(cbind(a1,b1,c1,d1,e1,f1,g1,h1,i1,j1,k1,l1,m1,n1,o1,p1,q1,r1,s1,t1,u1,v1,w1,x1,y1,z1,aa1,ab1,ac1,ad1,ae1,af1,ag1,ah1,ai1,aj1,a2,b2,c2,d2,e2,f2,g2,h2,i2,j2,k2,l2,m2,n2,o2,p2,q2,r2,s2,t2,u2,v2,w2,x2,y2,z2,aa2,ab2,ac2,ad2,ae2,af2,ag2,ah2,ai2,aj2)~diet*ecomorph*family/ecomorph,data=bothlakes)


Thanks

Allan