有关主成分分析的简单思路
代码
#first,input the data and transform the chinese into english
mydata<-read.table(“biao5.1.csv”,sep = “,”,header = TRUE)
mydataKaTeX parse error: Expected 'EOF', got '#' at position 77: …matrix(mydata)
#̲standarlize the…values
Characteristicvector<-special$vectors
which(Characteristicroots==0)
#select principal component about 85% variance
n<-length(Characteristicroots)
retain<-Characteristicroots[1]
var<-sum(Characteristicroots)
cnt<-1
while((retain/var)<0.85){cnt<-cnt+1;
retain<-(retain+Characteristicroots[cnt])}
#retain the principal component about 85% variance into a matrix
out <-Characteristicvector[,1:cnt]
课本上对于主成分分析的流程
后期学会作图再来添加,还有线性多重响应未做,因为这里特征根刚好不为零