ACP2.R

# preparación de los datos
x <- read.table("http://pastebin.com/raw.php?i=xeV7h24k",header=T,sep=",")
rownames(x) <- x[,1]
x <- x[-1]

x.pca <- prcomp(x)

plot(x.pca$x,pch="")
text(x.pca$x,labels=rownames(x))
abline(h = 0, v = 0, col = "gray60")

cat("\nValores propios\n")

## 
## Valores propios

print(x.pca$sdev)

##  [1] 7.923089e-01 7.662847e-01 6.504913e-01 5.891933e-01 5.154883e-01
##  [6] 3.943740e-01 2.887579e-01 2.748080e-01 1.767189e-01 8.847078e-02
## [11] 4.013925e-02 3.330074e-17

cat("\nValores propios en % \n")

## 
## Valores propios en %

print(round(100*(x.pca$sdev)/sum(x.pca$sdev),1))

##  [1] 17.3 16.7 14.2 12.9 11.3  8.6  6.3  6.0  3.9  1.9  0.9  0.0

cat("\nValores propios acumulados en % \n")

## 
## Valores propios acumulados en %

print(cumsum(round(100*(x.pca$sdev)/sum(x.pca$sdev),1)))

##  [1]  17.3  34.0  48.2  61.1  72.4  81.0  87.3  93.3  97.2  99.1 100.0
## [12] 100.0