[R] how can I save the estimates of a regression model in a file?

From: pilar schneider <pilar.schneider_at_gmail.com>
Date: Wed, 05 Nov 2008 12:12:08 -0500


Dear all
I need some help with R.
How can I save the estimates of a regression model in a file?

here is what I did:

  1. this is my regression model: fit1 <- lm(logmilk ~ logdays + days, data=data2)
  2. however, I want to get the parameters estimates for each individual (by group): so i did the following: by(data2, list(data2$V2),function(.data2) lm(logmilk ~ logdays + days, data= .data2))
  3. Then to keep the estimates in a file I did: res1 <- by(data2, list(data2$V2),function(.data2) lm(logmilk ~ logdays + V5, data= .data2))
  4. and this is what I get for each individual: : 1 Call: lm(formula = logmilk ~ logdays + days, data = .data2) Coefficients:
    (Intercept) logdays days
    3.414105 0.069387 -0.001732

: 2
Call:
lm(formula = logmilk ~ logdays + days, data = .data2)

Coefficients:
(Intercept) logdays days
  3.2276114 0.1223412 -0.0006836

and so on:

5) THE QUESTION is:

there is a way in R to get an output file as you get in SAS when you use: prog reg data=xxx outest=estimate;
I would need an output that looks like this:

individual     intercept     logdays   days           etc
1           3.414105        0.069387   -0.0006836
2           3.2276114       0.1223412  -0.0006836

n           ........               ..........       .........
THANKS FOR YOUR HELP Maria
-- 
-----------------------------
pilar.schneider_at_gmail.com

	[[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 Wed 05 Nov 2008 - 17:17:09 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 05 Nov 2008 - 20:30:23 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