mu<-c(5.682,88.22)
Sigma<-symMatrix(c(3.770,-5.495,13.53),2)
datum<-dataindI(50,mu,Sigma,2)
# estat is the list with the values of the auxiliary statistics.
estat<-estatauxind(datum,50,2)
T2indI(estat,50)
Run the code above in your browser using DataLab