###MANOVA

> data(iris, package="datasets")

> summary (iris)
  Sepal.Length    Sepal.Width     Petal.Length    Petal.Width   
 Min.   :4.300   Min.   :2.000   Min.   :1.000   Min.   :0.100  
 1st Qu.:5.100   1st Qu.:2.800   1st Qu.:1.600   1st Qu.:0.300  
 Median :5.800   Median :3.000   Median :4.350   Median :1.300  
 Mean   :5.843   Mean   :3.057   Mean   :3.758   Mean   :1.199  
 3rd Qu.:6.400   3rd Qu.:3.300   3rd Qu.:5.100   3rd Qu.:1.800  
 Max.   :7.900   Max.   :4.400   Max.   :6.900   Max.   :2.500  
       Species  
 setosa    :50  
 versicolor:50  
 virginica :50  
                
                
                

> iris1 <- manova(cbind(Sepal.Length, Sepal.Width) ~ Species, data = iris)

> summary (iris1)
           Df  Pillai approx F num Df den Df    Pr(>F)    
Species     2 0.94531   65.878      4    294 < 2.2e-16 ***
Residuals 147                                             
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

> iris2 <- manova(cbind(Sepal.Length, Sepal.Width, Petal.Length, Petal.Width) ~ Species, data = iris)

> summary (iris2)
           Df Pillai approx F num Df den Df    Pr(>F)    
Species     2 1.1919   53.466      8    290 < 2.2e-16 ***
Residuals 147                                            
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

> summary.aov (iris1)
 Response Sepal.Length :
             Df Sum Sq Mean Sq F value    Pr(>F)    
Species       2 63.212  31.606  119.26 < 2.2e-16 ***
Residuals   147 38.956   0.265                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

 Response Sepal.Width :
             Df Sum Sq Mean Sq F value    Pr(>F)    
Species       2 11.345  5.6725   49.16 < 2.2e-16 ***
Residuals   147 16.962  0.1154                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1


> summary (iris1, test = "Wilks")
           Df   Wilks approx F num Df den Df    Pr(>F)    
Species     2 0.16654   105.88      4    292 < 2.2e-16 ***
Residuals 147                                             
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

> scatterplotMatrix(~Petal.Length+Petal.Width+Sepal.Length+Sepal.Width | Species, 
+                        regLine=FALSE, smooth=FALSE, diagonal=list(method="density"), by.groups=TRUE, data=iris)

> scatterplotMatrix(~Petal.Length+Petal.Width+Sepal.Length+Sepal.Width | Species, 
+                        regLine=FALSE, smooth=FALSE, diagonal=list(method="boxplot"), by.groups=TRUE, data=iris)

###########################################################
################SCRIPT##########################

###MANOVA
data(iris, package="datasets")
summary (iris)
iris1 <- manova(cbind(Sepal.Length, Sepal.Width) ~ Species, data = iris)
summary (iris1)
iris2 <- manova(cbind(Sepal.Length, Sepal.Width, Petal.Length, Petal.Width) ~ Species, data = iris)
summary (iris2)
summary.aov (iris1)
summary (iris1, test = "Wilks")