Ok, I am sorry, I was building a reproducible example to post here and I made a mistake. Th second loop (in the variable k) updates ycorr. Therefore, it would be
ycorr <- rnorm(10) Nindiv <- 10 x <-matrix(sample(c(0:2),50,re = T),10) gtemp <- rnorm(10) for (i in 1:Nindiv) { for (k in 1:ncol(x)) { ycorr[i] <- ycorr[i] - x[i,k]*gtemp[k] } } -- View this message in context: http://r.789695.n4.nabble.com/Help-with-apply-tp3518773p3519270.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.