Hi Jim:

Thanks a lot for helping me with this question.

The "xxxx" are just a model formula.

For example,

lm(y~x, data=data1) #model 1
lm(y~x, data=data2) #model 2
...
lm(y~x, data=data10) #model 10


All I want is extracting the coefficient estimates and p values and saving
them into one data frame, so that I can easily compare these 10 models. I
think your way works. But can you tell me how to write those 10 models into
your list variable "myModels"? Thanks.

FD


On Mon, Sep 15, 2008 at 9:19 PM, jim holtman <jholtman@gmail.com> wrote:

> It is hard to tell since you have not told us what 'xxxxxxxxxxxxx' is,
> but one way is to create a list of the objects you want to send to
> 'lm' and then save them in a list as output:
>
> fit.mod <- list()
> for (i in seq_along(myModels)){
>    fit.mod[[i]] <- lm(myModels[[i]]
>  }
>
> On Mon, Sep 15, 2008 at 4:49 PM, FD <fhduan@gmail.com> wrote:
> > Hi All:
> >
> > I need to run a bunch of models one-by-one and retrieve the outputs (like
> > estimates, p values) accordingly. Some of models may not have outputs due
> to
> > singularity or even no observations. If this happens, assign NA to the
> > outputs.
> >
> > For example,
> >
> > #running models
> > fit.mod <- lm(xxxxxxxxxx)  #model 1
> >  fit.mod <- lm(xxxxxxxxxx)  #model 2
> > .........
> > fit.mod <- lm(xxxxxxxxxx)  #model 10
> >
> > #outputting estimates and p values to a data frame
> > ....
> >
> >
> > Any suggestions? Thank you very much in advance.
> >
> > FD
> >
> >        [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help@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<http://www.r-project.org/posting-guide.html>
> > and provide commented, minimal, self-contained, reproducible code.
> >
>
>
>
> --
> Jim Holtman
> Cincinnati, OH
> +1 513 646 9390
>
> What is the problem that you are trying to solve?
>

	[[alternative HTML version deleted]]

