[R] PCA: Showing file datalabels on biplot
colin1
colinhindmarch at talktalk.net
Mon Jan 25 17:39:48 CET 2010
The script below successfully produces a biplot of the data but the 'site
names' (rows) and the names of the 'response variables' (columns) are shown
as simple numerals (rather than the column and row names). How might I
'enforce' the use of the row/column names used in the datafile (section of
datafile shown below)?
Can anyone help, please?
Section of datafile
sample a b c d ---n42
HBR-6 1 1 1 1
HBR-7 1 1 1 2
HBR-8 1 1 1 1
HBR-9 1 1 1 2
HBR-10 1 1 1 1
HBR-11 1 1 1 1
HBR-12 1 1 1 1
HBR-13 1 1 1 1
HBR-14 2 1 2 2
HBR-15 2 1 2 1
HBR-16 1 1 2 1
HBR-17 1 1 2 1
HBR-18 1 1 2 1
HBR-19 1 1 2 1
HBR-20 1 1 2 1
HBR-21 1 2 1 1
HBR-30 1 1 2 2
HBR-31 1 2 1 2
HBR-32 2 2 1 2
HBR-33 1 2 1 2
HBR-34 1 2 2 2
HBR-35 1 2 2 1
HBR-36 1 2 1 2
HBR-37 1 2 1 2
cat <-read.table("R855868825matrix.csv", sep=",", header=T)
head (cat, n=24)
catopsis <- (cbind(cat$a, cat$b, cat$c, cat$d, cat$e, cat$f, cat$g, cat$h,
cat$i, cat$j, cat$k, cat$l, cat$m, cat$n, cat$o, cat$p, cat$q, cat$r, cat$s,
cat$t, cat$u, cat$v, cat$w, cat$x, cat$y, cat$z, cat$aa, cat$ab, cat$ac,
cat$ad, cat$ae, cat$af, cat$ag, cat$ah, cat$ai, cat$aj, cat$ak, cat$al,
cat$am, cat$an, cat$ao, cat$ap))
pca <- prcomp (t(catopsis))
plot(pca, catopsis [])
summary(prcomp(catopsis))
biplot(prcomp(catopsis, scale = TRUE))
--
View this message in context: http://n4.nabble.com/PCA-Showing-file-datalabels-on-biplot-tp1289524p1289524.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list