Thank you very much David for the advice ! I haven't checked the plyr and data.table packages but here is the script to replace the for loop ::
dfrlm1<-as.data.frame(do.call("rbind", lapply(split(df1, df1$Name), function(X){coef(rlm(col2 ~ col3, data=X))}))) names(dfrlm1)<-c("Intercept", "Slope") Have a nice end of day, Ptit Bleu. -- View this message in context: http://r.789695.n4.nabble.com/Add-values-of-rlm-coefficients-to-xyplot-tp3032166p3034373.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.