Try mapply(function(x,y)lm(y~x),as.data .frame(data.x),as.data.frame(data.y),SIMPLIFY=F)
2008/4/12, baptiste Auguié <[EMAIL PROTECTED]>: > Hello R list, > > > I have two matrices of identical dimensions, and I want to fit a > straight line for each pair of columns and plot the resulting lines. > I got it to work with a for loop, but there must be a better way, > > > > n<-5 > > N<-10 > > > > data.x<-matrix(1:(n*N),ncol=n) > > data.y<-matrix(1:(n*N) + rnorm(n*N,sd=1),ncol=n) > > > > matplot(data.x,data.y,t="p",pch=1:n,bty="n") > > > > for (ii in 1:n) > > { > > test <- lm(y~x,data=list(x=data.x[,ii],y=data.y[,ii]) ) > > lines(data.x[,ii],test$coefficients[1] + test$coefficients[2] * > > data.x[,ii],lty=2,col=ii) > > } > > > > Thanks, > > baptiste > > > _____________________________ > > Baptiste Auguié > > Physics Department > University of Exeter > Stocker Road, > Exeter, Devon, > EX4 4QL, UK > > Phone: +44 1392 264187 > > http://newton.ex.ac.uk/research/emag > http://projects.ex.ac.uk/atto > > ______________________________________________ > 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. > -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O ______________________________________________ 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.