taojiea1014的博客

factor.analy1<-function(S, m){p<-nrow(S); diag_S<-diag(S); sum_rank<-sum(diag_S)rowname<-paste("X", 1:p, sep="")colname<-paste("Factor", 1:m, sep="")A<-matrix(0, nrow=

test<-scale(test)ca<-cancor(test[,1:3],test[,4:6])U<-as.matrix(test[, 1:3])%*% ca$xcoefV<-as.matrix(test[, 4:6])%*% ca$ycoefplot(U[,1], V[,1], xlab="U1", ylab="V1")plot(U[,3], V[,3], xlab=...