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@gmail.com

	[[alternative HTML version deleted]]

