Hmm, after reading one of your other posts, I am thinking you may *just* want all pairwise combinations. This worked for me:
# Create a sample data frame with 60 named columns x <- data.frame(matrix(rnorm(1200), ncol = 60, dimnames = list(NULL, paste("Col", 1:60, sep='')))) # calculate the correlation matrix and the R^2 mycorrmatrix <- cor(x) myr2s <- mycorrmatrix^2 ###To get all the linear regressions # get the column names mynames <- colnames(x) #Use sapply() and paste() to create formulae for all pairwise combinations myformula <- sapply(mynames, function(x) {paste(x, "~", mynames)}) #Remove the diagonals since these regress a variable on itself diag(myformula) <- NA myformula <- myformula[which(!is.na(myformula))] #check the structure of myformula str(myformula) #initialize a list to store all the models mylms <- vector("list", 3540) # calculate all the models for(i in 1:3540) { mylms[[i]] <- lm(formula = myformula[i], data = x) } # As an example summary(mylms[[1]]) On Wed, Aug 18, 2010 at 7:21 AM, Joshua Wiley <jwiley.ps...@gmail.com> wrote: > On Wed, Aug 18, 2010 at 6:09 AM, ashz <a...@walla.co.il> wrote: >> >> Hi, >> >> I have 20*60 data matrix (with some NAs) and I wish to perfom a Pearson >> correlation coefficient matrix as well as simple linear regression equation > > The correlation matrix can be readily obtained by calling cor() on the > entire matrix. > >> and coefficient of determination (R2) for every possible combination. Any >> tip/idea/library/script how do to so. > > So you have 60 variables, and you want every possible combination? I > may be mistaken, but isn't that 60! if you ignore any interactions? > If so, this strikes me as a job for > > library(fortunes) > fortune("hundreds") > > >> >> Thanks, >> As hz >> >> -- >> View this message in context: >> http://r.789695.n4.nabble.com/Linear-regression-equation-and-coefficient-matrix-tp2329804p2329804.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. >> > > > > -- > Joshua Wiley > Ph.D. Student, Health Psychology > University of California, Los Angeles > http://www.joshuawiley.com/ > -- Joshua Wiley Ph.D. Student, Health Psychology University of California, Los Angeles http://www.joshuawiley.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.