*** Demonstration *** MANOVA in R
> help(manova)
> help(summary.manova)
> help(gl)  % generates factors
> da=read.table("t6-9.dat")
> dim(da)
[1] 48  4
> y=cbind(da[,1],da[,2],da[,3])
> gen=factor(gl(2,24))
> gen
 [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2
[39] 2 2 2 2 2 2 2 2 2 2
Levels: 1 2
> m1=manova(y~gen)
> m1
Call:
   manova(y ~ gen)

Terms:
                      gen Residuals
resp 1           6165.333 13576.583
resp 2           2451.021  5100.792
resp 3           1541.333  1747.917
Deg. of Freedom         1        46

Residual standard error: 17.17973 10.53028 6.164267 
Estimated effects may be unbalanced

> summary(m1,test="Wilks")
          Df   Wilks approx F num Df den Df    Pr(>F)    
gen        1  0.3886  23.0782      3     44 3.967e-09 ***
Residuals 46                                             
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

> da=read.table("T6-4.dat")
> da
   V1 V2  V3   V4  V5
1   0  0 6.5  9.5 4.4
2   0  0 6.2  9.9 6.4
3   0  0 5.8  9.6 3.0
4   0  0 6.5  9.6 4.1
5   0  0 6.5  9.2 0.8
6   0  1 6.9  9.1 5.7
7   0  1 7.2 10.0 2.0
8   0  1 6.9  9.9 3.9
9   0  1 6.1  9.5 1.9
10  0  1 6.3  9.4 5.7
11  1  0 6.7  9.1 2.8
12  1  0 6.6  9.3 4.1
13  1  0 7.2  8.3 3.8
14  1  0 7.1  8.4 1.6
15  1  0 6.8  8.5 3.4
16  1  1 7.1  9.2 8.4
17  1  1 7.0  8.8 5.2
18  1  1 7.2  9.7 6.9
19  1  1 7.5 10.1 2.7
20  1  1 7.6  9.2 1.9
> y=cbind(da[,3],da[,4],da[,5])
> fac1=factor(da[,1])
> fac1
 [1] 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1
Levels: 0 1
> fac2=factor(da[,2])
> m2=manova(y~fac1+fac2+fac1*fac2)
> m2
Call:
   manova(y ~ fac1 + fac2 + fac1 * fac2)

Terms:
                   fac1    fac2 fac1:fac2 Residuals
resp 1           1.7405  0.7605    0.0005    1.7640
resp 2           1.3005  0.6125    0.5445    2.6280
resp 3           0.4205  4.9005    3.9605   64.9240
Deg. of Freedom       1       1         1        16

Residual standard error: 0.3320392 0.4052777 2.014386 
Estimated effects may be unbalanced
> summary(m2,test="Wilks")
          Df  Wilks approx F num Df den Df   Pr(>F)   
fac1       1 0.3819   7.5543      3     14 0.003034 **
fac2       1 0.5230   4.2556      3     14 0.024745 * 
fac1:fac2  1 0.7771   1.3385      3     14 0.301782   
Residuals 16                                          
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

> summary(m2,test="Pillai")
          Df Pillai approx F num Df den Df   Pr(>F)   
fac1       1 0.6181   7.5543      3     14 0.003034 **
fac2       1 0.4770   4.2556      3     14 0.024745 * 
fac1:fac2  1 0.2229   1.3385      3     14 0.301782   
Residuals 16                                          
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1