Dear Rui and Ozgur Thanks to both of you for your swift and helpful responses!
The idea with the loop is great and the cbind() function is also what I am looking for. However, is it also possible to summarize the estimators and p-values for each data set y in one table only? For example, something like estimate p.value (Intercept) -0.8698107 1.445023e-01 x1 0.9991003 1.539969e-108 x2 1.1200786 3.691776e-14 estimate p.value (Intercept) -1.6826360 7.523808e-01 x1 0.6564982 1.632763e-14 x2 1.7081632 1.358624e-01 Is it even possible to call the R^2 of each regression? In the end I need to provide a table like this (for ~50 data sets): Intercept P-Value Beta 1 P-Value R2 data set1 ... ... ... ... ... data set2 ... ... ... ... ... Is an output like this even possible? Thank you once more for your kind help. Best, Chris -- View this message in context: http://r.789695.n4.nabble.com/Customized-R-Regression-Output-tp4631497p4631505.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.