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)
T2<-T2subI(estat,20,10)
dados<-datasubI(20,10,mu,Sigma,2)
# estat is a list with the auxiliary statistics. T2 is a matrix with the values of the T2 statistic.
graphsubI(estat,T2,20,10,2)
Run the code above in your browser using DataLab