# Example 3.1 (Table 2.1) in Davis, 2002

temp <- c(-10,25,37,50,65,80)
m <- matrix(scan("e:\\vhm\\vhm882\\data\\deal.dat"),ncol=7,byrow=T)
subject <- m[,1]
y <- m[,2:7]
diff <- matrix(c(y[,1]-y[,2],y[,2]-y[,3],y[,3]-y[,4],y[,4]-y[,5],y[,5]-y[,6]),nrow=8,ncol=5,byrow=FALSE)
summary(manova(diff~1),intercept=TRUE)