Notice: On March 31, it was **announced** that Statalist is moving from an email list to a **forum**. The old list will shut down on April 23, and its replacement, **statalist.org** is already up and running.

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

From |
"Michael N. Mitchell" <Michael.Norman.Mitchell@gmail.com> |

To |
statalist@hsphsun2.harvard.edu |

Subject |
Re: st: exact logistic regression-further details |

Date |
Fri, 25 Jun 2010 22:46:55 -0700 |

Dear Nick

I hope this helps, Michael N. Mitchell Data Management Using Stata - http://www.stata.com/bookstore/dmus.html A Visual Guide to Stata Graphics - http://www.stata.com/bookstore/vgsg.html Stata tidbit of the week - http://www.MichaelNormanMitchell.com On 2010-06-25 10.52 AM, Nikolaos Pandis wrote:

Dear Michael and all, I tried exlogistic because I could not fit the logistic model with the variables I was interested in doing so. I was getting the message that var predicts the outcome perfectly. Perhaps using exlogistic is not the appropriate approach, but anyways I gave it a try. I am giving below details of dataset and intended analysis: Variables id:patinet id missence: binary variable missing or not missing any teeth jaw: 4 categories corresponding to jaw location of missing teeth (0=no jaw-always 0 when var missinence is 0, 1=upper jaw, 2=lower jaw, 3=missing teeth in both jaws) sex-binary rightleft: 4 levels indicating if no side is missing any teeth, the right, the left or bot sides have missing teeth missingpattern: 5 categories and corresponds to different combinations of missing teeth. For example category 0 is combination of missing 2 particular teeth in the upper jaw. The idea is to find how the odds ratio of missing teeth vs no missing teeth vary after adjusting for the variables such as sex, jaw rightleft ect I fitted: .logistic missence sex // to get the odds ratio between male/female The above model did produce reasonable odds ratios and CIs I tried .xi:logistic missence sex i.jaw // to get the OR of missing in jaw category 1,2,3 vs jaw==0 which is equivelant to missence==0. This model does not converge and states that we have perfect data prediction. Is it because jaw==0 corresponds exactly to missence ==0 and jaw=1,2,3 corresponds to missence==1? I was wondering if there is way to produce the ORs that give the odds of jaw 1, 2 or 3 vs jaw==0 (adjusted for sex) or is it not possible? I also fitted: .mlogit missingpattern sex which run fine if I run mlogit missingpattern sex i.jaw I get z values close to 0 and zero values for CIs undicating again perfect prediction. I think I can see the problem and I am wondering if I should expand my data where for each id we have as many rows as the number of missing teeth? Any suggestions would be appreciated. Many thanks, Nick --- On Mon, 6/21/10, Michael N. Mitchell<Michael.Norman.Mitchell@gmail.com> wrote:From: Michael N. Mitchell<Michael.Norman.Mitchell@gmail.com> Subject: Re: st: exact logistic regression To: statalist@hsphsun2.harvard.edu Date: Monday, June 21, 2010, 11:29 PM Dear Nick I created an example based on the -auto- dataset like this... clear all sysuse auto replace rep78 = rep78 - 1 if rep78 !=1 tab rep78 expand 28 keep in 1/2061 xi: exlogistic foreign i.rep78 , memory(1400m) and I was not able to complete this either, getting the message observation 415: enumerations = 9333408 observation 416: enumerations = 9470664 the operating system refuses to provide the requested memory r(909); Do you have a very small or empty cell when you crosstabulate -a- by -b-? Perhaps you could share the crosstab of -a- and -b- with us to illustrate the need for exact logistic regression? Best regards, Michael N. Mitchell Data Management Using Stata - http://www.stata.com/bookstore/dmus.html A Visual Guide to Stata Graphics - http://www.stata.com/bookstore/vgsg.html Stata tidbit of the week - http://www.MichaelNormanMitchell.com On 2010-06-21 12.29 PM, Nikolaos Pandis wrote:Hi to all. 1. I try to run in Stata 11 IC: xi:exlogistic a i.b a= binary (0,1) b=4 level categorical (0-3) 2061 observations I get the following error message: exceeded memory limit of 10.0M bytes; use the memory()option to increase the memory limiterror(499) I use set memory to values larger than 10m and stillget the same message.exceeded memory limit of 10.0M bytes; use the memory()option to increase the memory limiterror(499) 2. I was wondering if it is possible to runmultinomial logistic regression using exact methods.Any comments would be appreciated. Many thanks, Nick * * 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/

**References**:**Re: st: exact logistic regression-further details***From:*Nikolaos Pandis <npandis@yahoo.com>

- Prev by Date:
**st: Summarizing Kappa Output** - Next by Date:
**st: kdensity with binned data** - Previous by thread:
**Re: st: exact logistic regression-further details** - Next by thread:
**st: estimation of heteroskedasticity 3SLS using cmp (conditional mixed process) estimator** - Index(es):