mu<-c(5.682,88.22)
Sigma<-symMatrix(c(3.770,-5.495,13.53),2)
datum<-datasubI(20,10,mu,Sigma,2)
estat<-estatauxsub(datum,20,10,2)
datum2<-datasubII(estat,10,p=2)
T2II<-T2subII(datum2,estat,10)
# for the first sample j=1. T2II is a vector with the value of the firts T2 statistic.
graphsubII(T2II,20,10,1,25,2)
Run the code above in your browser using DataLab