Re: [R] lm() of one matrix against another

From: Henrique Dallazuanna <wwwhsd_at_gmail.com>
Date: Sat, 12 Apr 2008 16:08:49 -0300

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é <ba208_at_exeter.ac.uk>:
> 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_at_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_at_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.
Received on Sat 12 Apr 2008 - 19:13:51 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sun 13 Apr 2008 - 13:30:28 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive