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