Re: [R] apply lm() for all the columns of a matrix

From: Costas Douvis <cdouvis_at_geol.uoa.gr>
Date: Wed, 09 Apr 2008 16:44:08 +0300 (EEST)


Thank you all very much for replying. Of course you are absolutely right but unfortunately I really deal with the case of a 4-d matrix so what you said does not apply. I should have specified but being a new R user I hadn't realized the difference between a matrix and an array.

So please tell me if you know a fast way (not using a loop) to perform a linear fit on all the vectors of the 4-th dimension of a 4-d array.

Thanks again
Kostas

> If you have the same design matrix then you can specify a matrix of
> responses in lm(), e.g.,
>
> Y <- matrix(rnorm(100*10), 100, 10)
> x <- rnorm(100)
>
> fit <- lm(Y ~ x)
> fit
> summary(fit)
>
>
> I hope it helps.
>
> Best,
> Dimitris
>
> ----
> Dimitris Rizopoulos
> Biostatistical Centre
> School of Public Health
> Catholic University of Leuven
>
> Address: Kapucijnenvoer 35, Leuven, Belgium
> Tel: +32/(0)16/336899
> Fax: +32/(0)16/337015
> Web: http://med.kuleuven.be/biostat/
> http://www.student.kuleuven.be/~m0390867/dimitris.htm
>
>
> ----- Original Message -----
> From: "Costas Douvis" <cdouvis_at_geol.uoa.gr>
> To: <r-help_at_r-project.org>
> Sent: Wednesday, April 09, 2008 12:55 PM
> Subject: [R] apply lm() for all the columns of a matrix
>
>
>> Hi all,
>>
>> My question is not really urgent. I can write a loop and solve the
>> problem. But I know that I'll be in a similar situation many more
>> times so
>> it would be useful to find out the answer
>>
>> Is there a fast way to perform linear fit to all the columns of a
>> matrix?
>> (or in the one dimension of a multi-dimensional array.) I'm talking
>> about
>> many single linear fits, not about a multiple fit. I thought that a
>> combination of apply and lm would do it but I can't make it work
>>
>> Thank you
>> Kostas
>>
>>
>> --
>> Kostas Douvis
>> PhD Student
>> University of Athens - Department of Geography and Climatology
>> Academy of Athens - Research Centre for Atmospheric Physics and
>> Climatology
>> email: cdouvis_at_geol.uoa.gr
>> tel: +30-210-8832048
>>
>> ______________________________________________
>> 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.
>>
>
>
> Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm
>
>



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 Wed 09 Apr 2008 - 13:50:17 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 Wed 09 Apr 2008 - 15:00: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