
Hello everyone.

 In my last post I did not explained my problem quite well. I made a principal component analysis and took the 2 first principal components. I made ​​a chart of my points based on the score of the 2 PC. I would like to add on this graph a  95% confidence region. To do this I used the ellipse function as follows: 


pcsref=PC$score[data[,1]==ref,1:2]  #matrix containing the scores of the two PC
SD=cov(pcsref[,1:2]) #covariance matrice of the two PC scores
lines (ellipse(SD, centre=colMeans(pcsref), level=0.95),type="l", lty=5)

But I want the confidence region calculated by the Hotelling T2 method and I'm not sure that ellipse function does it. Does anyone know how to do this?

 Thank you very much for your help and sorry for my english.
 		 	   		  
	[[alternative HTML version deleted]]

