主成分分析(R软件).doc

D=read.csv(file=file.choose(),head=T) a=array(0,dim=c(38,6)) for(i in 1:38){for(j in 1:6)a[i,j]=(D[i+1,j]-D[i,j])/D[i,j]} m=cor(a) eigen(m) 74.5 P-princomp(a, cor=TRUE) summary(P, loadings=TRUE) screeplot(P) screeplot(P,type=lines) biplot(P) 3. X=read.csv(file=file.choose(),head=T) d - dist(scale