[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)
# “… 'scatter.dudi' is a factorial map of individuals and the
# projection of the vectors of the canonical basis multiplied by a
# constante of rescaling. …”
?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