From: pilar schneider <pilar.schneider_at_gmail.com>

Date: Wed, 05 Nov 2008 12:12:08 -0500

: 2

Call:

lm(formula = logmilk ~ logdays + days, data = .data2)**THANKS FOR YOUR HELP
**
Maria

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?

- this is my regression model: fit1 <- lm(logmilk ~ logdays + days, data=data2)
- 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))
- Then to keep the estimates in a file I did: res1 <- by(data2, list(data2$V2),function(.data2) lm(logmilk ~ logdays + V5, data= .data2))
- 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 ........ .......... .........

-- ----------------------------- 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.
*