Re: [R] linear regression "group by"

From: Henrique Dallazuanna <wwwhsd_at_gmail.com>
Date: Mon, 14 Apr 2008 14:56:37 -0300

Try:

 by(x, list(x$STNID), function(.x)lm(TEMP ~ DEWP, data = .x))

On Mon, Apr 14, 2008 at 2:03 PM, Ryan Lauritsen <ryanlauritsen_at_gmail.com> wrote:

> Hi all. I'm brand new to R.
>
> My dataset (stored in MySQL) is a list of weather stations in rows by
> year with various weather variables in columns, for example:
> STNID YEAR TEMP DEWP
> station1 1990 54 50
> station1 1991 23 10
> station1 1992 34 18
> station2 1990 45 41
> station2 1991 32 25
> station2 1992 21 11
>
> I'm trying to run linear regression and get the basic output (i.e.
> intercept, slope, and significance) for each station. I'm able to run
> the regression on the entire dataset using:
> lm(TEMP~DEWP, data=select)
> But is there a way to aggregate the data ("group by" in MySQL) by
> STNID during the regression? Ideally I would just have a list of
> stations and their approriate summary output, which I could use for
> further analysis.
>
> I've searched the manual, etc. for solutions, but have been
> unsuccessful. Any assistance is greatly appreciated.
>
> Thank you,
> Ryan
>
> ______________________________________________
> 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

	[[alternative HTML version deleted]]


______________________________________________ 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 Mon 14 Apr 2008 - 21:29:00 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 Mon 14 Apr 2008 - 21:30:29 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