[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

From |
"Nick Cox" <n.j.cox@durham.ac.uk> |

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

Subject |
st: RE: KcKelvey & Zavoina's R2 in -fitstat- after -gologit2- |

Date |
Wed, 5 Mar 2008 21:41:34 -0000 |

Richard Williams gave a good answer. A detail remains over output shown as a stop or period (.). Just as with data, a numeric result shown as . is Stata's way of saying missing, or don't know. Just from missing, however, it is difficult to reverse engineer what lies behind it. There are two main possibilities. 1. Stata did the calculation, but can't show a number because the result is indeterminate. . di 7/0 . is a simple example. 7/0 is indeterminate. 2. Stata doesn't know what you are referring to. . summarize foobar . di r(SUM) . You ask for a display of r(SUM). But so such result is known to Stata. It shrugs its shoulders and puts missing. (You might have meant r(sum). Such typos are quite a common reason for a display of missing.) Nick n.j.cox@durham.ac.uk Sara Mottram I have fitted a partial proportional odds model in Stata 9.2 using -gologit2-. In their book, Regression Models for Categorical Dependent Data using Stata, Long and Freese suggest that McKelvey & Zavoina's R2 most closely approximates the R2 from linear regression models (which I assume makes this the most suitable R2 to use). Following the advice on Richard William's website to use the v1 option, I have been able to use -fitstat- with -gologit2-, however, it does not produce the M&Z R2. Please could someone tell me if this is because it does not make sense with a partial proportional odds model, in which case, is there a more appropriate measure of fit that I can use? Or do I need to specifically ask for M&Z R2 via an option? There was a similar question posed last year, but the solution suggested to that post (see below) gives me an estimate of "." I'm not sure what that means. sysuse auto logit foreign price fitstat estadd scalar r2_mz = r(r2_mz) and then tabulate using: estout ., eform cells(n se(par)) stats(r2_mz) style(fixed) or: esttab ., eform se scalars(r2_mz) * * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

**Follow-Ups**:**Re: st: RE: KcKelvey & Zavoina's R2 in -fitstat- after -gologit2-***From:*Richard Williams <Richard.A.Williams.5@ND.edu>

**References**:**st: KcKelvey & Zavoina's R2 in -fitstat- after -gologit2-***From:*Sara Mottram <s.mottram@cphc.keele.ac.uk>

- Prev by Date:
**st: interpret treatment effect coefficient** - Next by Date:
**Re: st: RE: RE: RE: RE: Re: Compositional data** - Previous by thread:
**Re: st: KcKelvey & Zavoina's R2 in -fitstat- after -gologit2-** - Next by thread:
**Re: st: RE: KcKelvey & Zavoina's R2 in -fitstat- after -gologit2-** - Index(es):

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