####################################### # SMP Componenti principali. # Matrice delle omponenti principali. # Si possono scegliere colonne con j. Smp = function (X,j) {XE=Smg.cen(X)%*%Smp.autovettori(X) if (missing(j)) XE else XE[,j]} # Alternativa ad Smp. Smp.R = function (X) {p=princomp(X)$scores dimnames(p)=NULL; p} # Autovalori di X^{COM}. Smp.autovalori = function (X) {Com=Sm.com(X) Mm.autovalori(Com,simm=T)} # Autovettori di X^{COM}. Smp.autovettori = function (X) {Com=Sm.com(X) Mm.autovettori(Com,simm=T)}