library(MVA)
library(lattice)

# dataset sull'inquinamento di alcune cittą USA
X = USairpollution

# dataset standardizzato
Z = scale(X)

round(apply(Z, 2, mean),2)
round(apply(Z, 2, var),2)

# matrice delle distanze
D =as.matrix(dist(Z))
D = round(D, 2)

rownames(D)=colnames(D)=abbreviate(rownames(D),3)
head(D)

# grafico di D
levelplot(D)

# norma di D
norm(D)

# grafico su PCA1,2
windows()
plot(princomp(D)$scores[,1:2], type="n")
text(princomp(D)$scores[,1:2], labels = rownames(X), cex = 0.6, col ="blue")
abline(h=0, lty=2, col="green")
abline(v=0, lty=2, col="green")