[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From |
"Martin Weiss" <martin.weiss1@gmx.de> |

To |
<statalist@hsphsun2.harvard.edu> |

Subject |
st: Re: RE: Re: RE: RE: R2 stats using statsby or parmby??? |

Date |
Sun, 28 Jun 2009 11:23:56 +0200 |

<>

*** use http://www.ats.ucla.edu/stat/stata/dae/logit.dta, clear global tflist "" global modseq=0 foreach X in gre topnotch gpa { global modseq=$modseq+1 tempfile tf$modseq parmby "logit admit `X'", label /// command format(estimate min95 max95 /// %8.2f p %8.1e) idn($modseq) /// saving(`tf$modseq',replace) /// flist(tflist) escal(r2_p) /// rename(es_1 r2_p) } dsconcat $tflist sort idnum parmseq describe bys idnum command:list parm /// label estimate min95 /// max95 p r2_p,noobs *** HTH Martin _______________________

To: <statalist@hsphsun2.harvard.edu> Sent: Sunday, June 28, 2009 12:38 AM Subject: st: RE: Re: RE: RE: R2 stats using statsby or parmby???

Thanks Martin - you've got me moving forward again. You wouldn't happen to know how to get the model pseudo r-square displayed would you? In Roger's post below he suggests that parmby calculates and stores it as r2_p, but I don't see it. -----Original Message----- From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Martin Weiss Sent: Saturday, June 27, 2009 3:44 PM To: statalist@hsphsun2.harvard.edu Subject: st: Re: RE: RE: R2 stats using statsby or parmby??? <> Not sure whether that blank after -flist-- hurt you, but you should use line continuation "///". Also note that -by- likes to be passed a dataset sorted by its -varlist-, so append an "s" to make it "bysort"... *** use http://www.ats.ucla.edu/stat/stata/dae/logit.dta, clear global tflist "" global modseq=0 foreach X in gre topnotch gpa { global modseq=$modseq+1 tempfile tf$modseq parmby "logit admit `X'", label /// command format(estimate min95 max95 /// %8.2f p %8.1e) idn($modseq) /// saving(`tf$modseq',replace) /// flist(tflist) } dsconcat $tflist sort idnum parmseq describe bys idnum command:list parm /// label estimate min95 max95 p,noobs *** HTH Martin _______________________----- Original Message -----From: "Hugh Robinson" <hugh.robinson@umontana.edu>To: <statalist@hsphsun2.harvard.edu> Sent: Saturday, June 27, 2009 11:10 PM Subject: st: RE: RE: R2 stats using statsby or parmby???Roger, I was intrigued by reading your response to this posting, I think the foreach and parmby commands could be very helpful for me in running a number of univariat logits and producing a single table containingeachvariable name, its coefficient, p-value, and r-square value. Towards that goal I thought I would start with the UCLA sample logit data and try to modify the command list from the parmby help file to create something looks like what I'm after. When I run the following code the foreach loop only gets as far the first covariate "gre" then produces an error that reads "variabletflistnot found". Can you see what's wrong? Thanks HR use http://www.ats.ucla.edu/stat/stata/dae/logit.dta, clear global tflist "" global modseq=0 foreach X in gre topnotch gpa { global modseq=$modseq+1 tempfile tf$modseq parmby "logit admit `X'", label command format(estimate min95 max95 %8.2f p %8.1e) idn($modseq) saving(`tf$modseq',replace) flist (tflist) } dsconcat $tflist sort idnum parmseq describe by idnum command:list parm label estimate min95 max95 p,noobs -----Original Message----- From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Newson,RogerB Sent: Friday, June 26, 2009 1:07 PM To: 'statalist@hsphsun2.harvard.edu' Subject: st: RE: R2 stats using statsby or parmby??? You do not specify what goodness-of-fit statistic you want, or which program you are using to do your logistic regression (O prefer -glm-, but many other people prefer -logit- or -logistic-). However, if the goodness-of-fit statistic is saved in e(), then either -statsby- or -parmby- can save it. In the case of -parmby-, you use the -escal()- option to save extra scalar results. As in parmby "glm y x, link(logit) family(bin) eform robust", eformnorestoreby(group) escal(aic bic chi2 dispers_p) which should save -e(aic)-, -e(bic)-, -e(chi2)- and -e(dispers_p)- in scalar variables named -es_1-, -es_2-, -es_3-, and -es_4-,respectively.You can give them more informative names using the -rename()- option. Similarly, the -logit- and -logistic- commands save thepseudo-R-squaredin -e(r2_p)-. I hope this helps. Best wishes Roger Roger B Newson BSc MSc DPhil Lecturer in Medical Statistics Respiratory Epidemiology and Public Health Group National Heart and Lung Institute Imperial College London Royal Brompton Campus Room 33, Emmanuel Kaye Building 1B Manresa Road London SW3 6LR UNITED KINGDOM Tel: +44 (0)20 7352 8121 ext 3381 Fax: +44 (0)20 7351 8322 Email: r.newson@imperial.ac.uk Web page: http://www.imperial.ac.uk/nhli/r.newson/ Departmental Web page:http://www1.imperial.ac.uk/medicine/about/divisions/nhli/respiration/popgenetics/reph/ Opinions expressed are those of the author, not of the institution. -----Original Message----- From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Thomas M Holbrook Sent: 26 June 2009 18:09 To: statalist@hsphsun2.harvard.edu Subject: st: R2 stats using statsby or parmby??? I'm trying to find a way the generate goodness-of-fit stats for logit models run separately over subgroups of my my data (I'm running a vote model using individual-level data and I want to generate fit stats by day of the campaign). I can get the slopes and standard errors using "statsby" or "parmby" but I don't set a way of generating the fitstats.Any ideas??? -Tom Thomas M. Holbrook Wilder Crane Professor of Government Department of Political Science University of Wisconsin-Milwaukee 3210 North Maryland Avenue Milwaukee, WI 53211 www.uwm.edu/~holbroot www.election08data.blogspot.com 414-229-6468 * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/ * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/ * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/* * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/ * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

* * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

**Follow-Ups**:**st: RE: Re: RE: Re: RE: RE: R2 stats using statsby or parmby???***From:*"Hugh Robinson" <hugh.robinson@umontana.edu>

**References**:**st: R2 stats using statsby or parmby???***From:*Thomas M Holbrook <holbroot@uwm.edu>

**st: RE: R2 stats using statsby or parmby???***From:*"Newson, Roger B" <r.newson@imperial.ac.uk>

**st: RE: RE: R2 stats using statsby or parmby???***From:*"Hugh Robinson" <hugh.robinson@umontana.edu>

**st: Re: RE: RE: R2 stats using statsby or parmby???***From:*"Martin Weiss" <martin.weiss1@gmx.de>

**st: RE: Re: RE: RE: R2 stats using statsby or parmby???***From:*"Hugh Robinson" <hugh.robinson@umontana.edu>

- Prev by Date:
**st: Re: RE: creating a text slide** - Next by Date:
**st: lag endog variables as instruments in ivreg2** - Previous by thread:
**RE: st: Stata 11 data format** - Next by thread:
**st: RE: Re: RE: Re: RE: RE: R2 stats using statsby or parmby???** - Index(es):

© Copyright 1996–2015 StataCorp LP | Terms of use | Privacy | Contact us | What's new | Site index |