iris<-matrix(scan(file='iris.txt'),ncol=4,byrow=T)
type<-factor(rep(LETTERS[1:3],c(50,50,50)))
ir.df<-data.frame(iris,type)
#
postscript(horizo=F,file='ir1.ps')
par(mfrow=c(2,2),bty='l',las=1,oma=c(0,0,4,0))
plot.design(ir.df)
mtext(' Factor means for iris data',
line=1,side=3,cex=1.5,outer=T)
dev.off()
#
postscript(horizo=F,file='ir1.ps')
par(mfrow=c(2,2),bty='l',las=1,oma=c(0,0,4,0))
boxplot(ir.df$X1~ir.df$type,main='X1')
boxplot(ir.df$X2~ir.df$type,main='X2')
boxplot(ir.df$X3~ir.df$type,main='X3')
boxplot(ir.df$X4~ir.df$type,main='X4')
mtext(' Factor boxplots for iris data',
line=1,side=3,cex=1.5,outer=T)
dev.off()
#
m1<-manova(iris~type,ir.df)
summary(m1,test='Pillai')
summary(m1,test='Wilks')
summary(m1,test='Hotelling-Lawley')
summary(m1,test='Roy')