# Create an Interval-Data object containing the intervals for 899 observations
# on the temperatures by quarter in 60 Chinese meteorological stations.
ChinaT <- IData(ChinaTemp[1:8],VarNames=c("T1","T2","T3","T4"))
#MANOVA tests, assuming that one of C3, C4 or C5 restricted configurations holds
#Classical (homocedastic) MANOVA tests
ManvChina <- MANOVA(ChinaT,ChinaTemp$GeoReg,Config=3:5)
cat("China, MANOVA by geografical regions results =\n")
print(ManvChina)
#Heterocedastic MANOVA tests
HetManvChina <- MANOVA(ChinaT,ChinaTemp$GeoReg,Mxt="Het",Config=3:5)
cat("China, heterocedastic MANOVA by geografical regions results =\n")
print(HetManvChina)
Run the code above in your browser using DataLab