Dear Users, I have one problem here, I tried many time and even read a few notes on writing function but still.....
Can anyone help me on how to simplify Part B (please refer the programming below), so that I don't have to repeat the calculation of "num" and "r" ? Thank you very much..every help is very much appreciated... ## Part A n=1400 m=matrix(c(0,0,0,0),4,1) m2=matrix(c(0,2,0,0),4,1) c4<-matrix(c(1.0,0.2,-0.5,0.3,0.2,1,0.2,-0.5,-0.5,0.2,1,0.2,0.3,-0.5,0.2,1),4,4,byrow=T) set.seed(428) X=mvrnorm(n,m2,c4) X.pca<-princomp(X) loadings(X.pca) pc=X.pca$loadings[,1] ## Part B num1=rep(1:n) for(i in 1:n)num1[i]=pc[1]%*%X[i,1] num2=rep(1:n) for(i in 1:n)num2[i]=pc[2]%*%X[i,2] num3=rep(1:n) for(i in 1:n)num3[i]=pc[3]%*%X[i,3] num4=rep(1:n) for(i in 1:n)num4[i]=pc[4]%*%X[i,4] den=rep(1:n) for(i in 1:n)den[i]=pc%*%X[i,] r1=num1/den r2=num2/den r3=num3/den r4=num4/den MLAV=sum(r2>r1 & r2>r3 & r2>r4) MLAV Best regards, ayu -- View this message in context: http://n4.nabble.com/How-to-repeat-for-function-tp1458806p1458806.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
