Many thanks, works great! Charilaos Skiadas-3 wrote: > > Try this: > > lapply( 1:2, function(i) lm( y~x, data=list(x=xdat[,i], y=ydat[,i]) ) ) > > Haris Skiadas > -- View this message in context: http://www.nabble.com/Using-lm-with-a-matrix--tp17708207p17829661.html Sent from the R help mailing list archive at Nabble.com.