1 2 3 4 5 6 7 8 | pc.biplot(SV, x = 1, y = 2, obs.opt = list(cex = 2, pch = 20),
obs.names = FALSE, obs.col.palette = NULL, obs.pch.palette = NULL,
var.opt = NULL, var.names = NULL, var.col.palette = NULL,
lwd = 2, filename = NULL, addPercEV = TRUE, var.suppl = NULL,
var.suppl.color = NULL, var.text.size = 3, var.text.col = "gray30",
var.text.relpos = NULL, var.text.cex = 1, var.arrow.size = 0.2,
main = "Biplot", asp = 1, alpha = 0.5, xlim = NULL,
ylim = NULL, legend = NULL, xylabs = "PC", percToAdd = NULL, ...)
|
SV |
the result of a call to princomp, prcomp, svd function. It can also be a metaMDS object (from vegan library). |
x |
dimension on the x axis (number or name) |
y |
dimension on the y axis (number or name) |
obs.opt |
\item obs.names \itemobs.col.palette \itemvar.opt \itemvar.names \itemvar.col.palette \itemlwd \itemfilename \itemaddPercEV \itemvar.suppl \itemvar.suppl.color \itemvar.text.size \itemvar.arrow.size \itemmain \itemasp \itemalpha \itemxlim \itemylim \itemlegend \itemxylabsstring or vector of labels. percToAddusually ... \itemtitlea string. "Biplot" by default. \itemobs.color.title \itemobs.pch.title SV e' il risultato di prcomp, di princomp o di svd. assi calcolati come in http://en.wikipedia.org/wiki/Biplot si veda anche http://www.multivariatestatistics.org/biplots.html data(pal.uno) par(mar=c(1,1,1,1)) Y=matrix(rnorm(30),10,3) rownames(Y)=paste("obs",1:nrow(Y)) sv=svd(Y) pc.biplot(sv) #i nomi vengono persi con svd pc.biplot(sv,obs.names = TRUE) # soluzione a mano rownames(sv$u)=rownames(Y) pc.biplot(sv,obs.names = TRUE) ########### sv=svd(scale(Y,center=TRUE,scale=FALSE)) pc.biplot(sv,obs.names = TRUE) pc=princomp(Y) pc.biplot(pc,obs.names = TRUE) pc=prcomp(Y) pc.biplot(pc,obs.names = TRUE) pc.biplot(sv,asp=1,obs.col.palette = pal.uno) pc.biplot(sv,obs.opt = list(col=rep(1:2,5))) ~kwd1 ~kwd2 |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.