[R] biplot drawing conc ellipses

Pierre BADY pierre.bady at univ-lyon1.fr
Fri Nov 5 00:15:02 CET 2004


Hi all,

you can see in the library(ade4).
there is  a nice graphical function called ‘s.class’.

Example for a principal components analysis:
?s.class

data(meaudret)
pca1 <- dudi.pca(meaudret$mil, scan = FALSE, nf = 4)
# plot row coordinates 
s.class(pca1$li, meaudret$plan$dat, sub = "Principal Component Analysis",add.
p=F,col=c(1,2,3,4),cstar=0)
# add principal axes
s.arrow(pca1$c1,clabel=0.7,add.p=T)


# &#8220;&#8230;  'scatter.dudi' is a factorial map of individuals and the
# projection of the vectors of the canonical basis multiplied by a
# constante of rescaling. &#8230;&#8221;

?scatter.dudi
scatter.dudi(pca1,clab.row=0)
s.class(pca1$li, meaudret$plan$dat, sub = "Principal Component Analysis",add.
p=T,col=c(1,2,3,4), cstar=0)

hope this helps,


P.BADY

Quoting "T. Murlidharan Nair" <nair at sdsc.edu>:

> Is there an option to draw concentration ellipses in biplots ? It
> seems
> really nice to summarize large number of points of each group.
> 
> Cheers../ Murli
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
> 
> 



---------------------------------------------------------
Pierre BADY     <°)))><
Université Claude Bernard Lyon 1
UMR CNRS 5023, LEHF
bat Alphonse Forel
43 boulevard du 11 novembre 1918 
F-69622 VILLEURBANNE CEDEX 
FRANCE
TEL : +33 (0)4 72 44 62 34 
FAX : +33 (0)4 72 43 28 92 
MEL : pierre.bady at univ-lyon1.fr
http://limnologie.univ-lyon1.fr




More information about the R-help mailing list