主成分分析，自编函数

最新推荐文章于 2023-06-04 18:40:06 发布

蔡定诚

最新推荐文章于 2023-06-04 18:40:06 发布

阅读量331

点赞数

本文链接：https://blog.csdn.net/weixin_43207304/article/details/89735141

版权

cost<-as.matrix(scale(read.table(‘D:/main.txt’,sep = ‘\t’,header = FALSE)))
covariancematrix <-function(data=cost){
averagex<-colMeans(data)
I<-matrix(1,nrow =nrow(data),ncol = 1)
Aminus<-cost-I%%t(as.matrix(averagex))
E<-t(Aminus)%%Aminus
return(E)
}
mainmaterial<-function(data=cost,n1=ncol(cost)){
EMT<-covariancematrix(data)
values<-eigen(EMT) $v a l u e s v e c t o r s < - e i g e n (E M T)$ vectors
mainmtrix<-matrix(data = NA,ncol = n1,nrow =ncol(vectors) )
for (i in 1:n1) {
a<-sort(values,decreasing=T)[i]
b<-which(values==a)
mainmtrix[,i]<-vectors[,b]
}
return(mainmtrix)
}