[R] Observation in a confidence ellipse

2011-05-28 Thread Jessica Minkue
Hello everyone I really need some help here. I made a confidence ellipse using the function ellipse from the package ellipse: ellipse(SD, centre=colMeans(pcsref),t=sqrt((p * (n-1)/(n-p))*qf(0.99, p,n-p)) Now, I want to write a function whom return TRUE or FALSE if a given observation is in th

[R] help: Using hotelling for a confidence region for PCA scores

2011-05-16 Thread Jessica Minkue
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