This is close to what you want. I created the list by hand, but you
can create it in your processing loop. Once you have the list
created, you can create your own print routine.
> x <- runif(100)
> z <- runif(100)
> y <- runif(100)
>
>
> # I am doing this by hand, but you could easily automate it in a loop or
> function call
> whatIwant <- list() # where results are stored
>
> lm.x <- lm(y~x)
> whatIwant[[1]] <- list(summary(lm.x)$call, summary(lm.x)$coefficients)
>
> lm.z <- lm(y~z)
> whatIwant[[2]] <- list(summary(lm.z)$call, summary(lm.z)$coefficients)
>
> lm.xz <- lm(y~x+z)
> whatIwant[[3]] <- list(summary(lm.xz)$call, summary(lm.xz)$coefficients)
>
> # this is close
> whatIwant
[[1]]
[[1]][[1]]
lm(formula = y ~ x)
[[1]][[2]]
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.46211645 0.05442237 8.4912958 2.244835e-13
x 0.02387787 0.09709648 0.2459190 8.062592e-01
[[2]]
[[2]][[1]]
lm(formula = y ~ z)
[[2]][[2]]
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.49612568 0.05390759 9.2032627 6.488440e-15
z -0.04600174 0.09298692 -0.4947119 6.219108e-01
[[3]]
[[3]][[1]]
lm(formula = y ~ x + z)
[[3]][[2]]
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.48617771 0.07502099 6.4805555 3.804443e-09
x 0.01880206 0.09808697 0.1916877 8.483876e-01
z -0.04400913 0.09402369 -0.4680643 6.407886e-01
>
> #now create a function to print your list
> printMyList <-
+ function(myList){
+ lapply(myList, function(.ele){
+ print(.ele[[1]])
+ print(.ele[[2]])
+ })
+ invisible(NULL)
+ }
>
> printMyList(whatIwant)
lm(formula = y ~ x)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.46211645 0.05442237 8.4912958 2.244835e-13
x 0.02387787 0.09709648 0.2459190 8.062592e-01
lm(formula = y ~ z)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.49612568 0.05390759 9.2032627 6.488440e-15
z -0.04600174 0.09298692 -0.4947119 6.219108e-01
lm(formula = y ~ x + z)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.48617771 0.07502099 6.4805555 3.804443e-09
x 0.01880206 0.09808697 0.1916877 8.483876e-01
z -0.04400913 0.09402369 -0.4680643 6.407886e-01
>
>
On 9/30/07, John Sorkin <[EMAIL PROTECTED]> wrote:
> Jim,
> You are indeed trying to help, again my thanks.
> What I want to do is make a single structure - a table is an apt description
> that will summarize all the regressions, something like:
>
> Estimate Std. Error t value Pr(>|t|)
> (Intercept) lm(formula = y ~ x) 4.927791 2.62115494 1.880007 6.307727e-02
> x 1.887634 0.04372724 43.168382 1.410167e-65
> (Intercept) lm(formula = z ~ x) 6.927791 2.62115494 1.880007 6.307727e-02
> x 1.887634 0.04372724 43.168382 1.410167e-65
> (Intercept) lm(formula = z~x+z) 6.927791 2.62115494 1.880007 6.307727e-02
> x 1.887634 0.04372724 43.168382 1.410167e-65
> z 1.887634 0.04372724 43.168382 1.410167e-65
>
> If you use a non-proportional spaced font (e.g. Courier on a windows system)
> the material above is a single "table" that contains all the data from my
> numerous regressions, including the call and the coefficients.(If I can get
> this to work I will add the R squared values) n.b. In the example above I have
> copied the estimates and changed the label of the lines just to demonstrate
> the
> kind of output I desire. Of course when used on real data each line will have
> different values.
> Again thanks,
> John
>
>
> John Sorkin M.D., Ph.D.
> Chief, Biostatistics and Informatics
> University of Maryland School of Medicine Division of Gerontology
> Baltimore VA Medical Center
> 10 North Greene Street
> GRECC (BT/18/GR)
> Baltimore, MD 21201-1524
> (Phone) 410-605-7119
> (Fax) 410-605-7913 (Please call phone number above prior to faxing)
>
> >>> "jim holtman" <[EMAIL PROTECTED]> 9/30/2007 11:24 PM >>>
> The easiest thing is to same the result from 'lm' and then you can use
> that to extract any of the information that you want. If you are
> going to do thousands of regressions, then you can make a list of 'lm'
> results, or if that is too memory intensive (depends on the size of
> the regressions), you could same them to a file with 'save'. It all
> depends on what you want to do with them. One of my favorite
> questions is "tell me what you want to do, not how you want to do it".
> What is the intended purpose?
>
> On 9/30/07, John Sorkin <[EMAIL PROTECTED]> wrote:
> > Jim,
> > Again thank you for your quick reply. Your suggestion does not give me
> > exactly what I want:
> >
> > > whatIwant<-list(,summary(fitdelete)$call,summary(fitdelete)$coefficients)
> > > whatIwant
> > [[1]]
> > lm(formula = y ~ x)
> >
> > [[2]]
> > Estimate Std. Error t value Pr(>|t|)
> > (Intercept) 4.927791 2.62115494 1.880007 6.307727e-02
> > x 1.887634 0.04372724 43.168382 1.410167e-65
> >
> > What I want is something that looks more like a table:
> >
> > lm(formula = y ~ x) Estimate Std. Error t value Pr(>|t|)
> > (Intercept) 4.927791 2.62115494 1.880007 6.307727e-02
> > x 1.887634 0.04372724 43.168382 1.410167e-65
> >
> > Thanks,
> > John
> >
> >
> >
> > John Sorkin M.D., Ph.D.
> > Chief, Biostatistics and Informatics
> > University of Maryland School of Medicine Division of Gerontology
> > Baltimore VA Medical Center
> > 10 North Greene Street
> > GRECC (BT/18/GR)
> > Baltimore, MD 21201-1524
> > (Phone) 410-605-7119
> > (Fax) 410-605-7913 (Please call phone number above prior to faxing)
> >
> > >>> "jim holtman" <[EMAIL PROTECTED]> 9/30/2007 11:00 PM >>>
> > try using a 'list':
> >
> > whatIwant<-list(call=summary(myreg)$call, coef=summary(myreg)$coefficients)
> >
> > On 9/30/07, John Sorkin <[EMAIL PROTECTED]> wrote:
> > > Widows XP
> > > R 2.3.1
> > >
> > > I have been trying to make a data structure that will contain both the
> > > coefficients from a linear regression along with column and row titles
> > > AND the call, i.e.
> > > myreg<-lm(y~x+y+z)
> > > whatIwant<-cbind(c(summary(myreg)$call,"",""),summary(myreg)$coefficients)
> > >
> > > Neither the statement above, nor any one of twenty variations I have
> > > tried work. I would appreciate any advice.
> > > Thanks,
> > > John
> > >
> > >
> > >
> > >
> > > John Sorkin M.D., Ph.D.
> > > Chief, Biostatistics and Informatics
> > > University of Maryland School of Medicine Division of Gerontology
> > > Baltimore VA Medical Center
> > > 10 North Greene Street
> > > GRECC (BT/18/GR)
> > > Baltimore, MD 21201-1524
> > > (Phone) 410-605-7119
> > > (Fax) 410-605-7913 (Please call phone number above prior to faxing)
> > >
> > > Confidentiality Statement:
> > > This email message, including any attachments, is for the so...{{dropped}}
> > >
> > > ______________________________________________
> > > 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
> > > and provide commented, minimal, self-contained, reproducible code.
> > >
> >
> >
> > --
> > Jim Holtman
> > Cincinnati, OH
> > +1 513 646 9390
> >
> > What is the problem you are trying to solve?
> >
> > Confidentiality Statement:
> > This email message, including any attachments, is for the sole use of the
> > intended recipient(s) and may contain confidential and privileged
> > information. Any unauthorized use, disclosure or distribution is
> > prohibited. If you are not the intended recipient, please contact the
> > sender by reply email and destroy all copies of the original message.
> >
>
>
> --
> Jim Holtman
> Cincinnati, OH
> +1 513 646 9390
>
> What is the problem you are trying to solve?
>
> Confidentiality Statement:
> This email message, including any attachments, is for the ...{{dropped}}
______________________________________________
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
and provide commented, minimal, self-contained, reproducible code.
