taojiea1014的博客

student.pr <- princomp(student, cor = TRUE) summary(student.pr, loadings=TRUE)predict(student.pr)screeplot(student.pr,type="lines")R<-matrix(0, nrow=16, ncol=16, dimnames=list(names, names))for ...

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=...