Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
From | Milena Przheska <milena.przheska@cosmicdevelopment.com> |
To | statalist@hsphsun2.harvard.edu |
Subject | Re: st: bayesmlogit |
Date | Wed, 21 Aug 2013 19:10:04 +0200 |
Hi Matthew, First of all, I like to express appreciation for the package you created. I was able to get the estimates. As for the r(1400) that I got, I will be more precise - the error message was: 1927 (group size) take 674 (# positives) combinations results in numeric overflow; computations cannot proceed This happened when I tried to run the code on a data set that had 7708 rows for 4 levels of the group variable and 1927 unique values for the identifier variable. When I limited the scope of the data set to the first 2000 rows the command was executed. The first part of the output was: Bayesian Mixed Logit Model Acceptance rates: Fixed coefs = 0.234 Random coefs(ave,min,max)= 0.070, 0.005, 0.267 Observations = 2000 Groups = 500 Choices = 4 Total draws = 1000 Burn-in draws = 0 I did not see any limitations in the help file for the command, so I didn't know it could be a problem. Best regards, Milena On Wed, Aug 21, 2013 at 6:35 PM, Matthew Baker <matthew.baker@hunter.cuny.edu> wrote: > Well, that was easy! Nothing like reading through email, finding that > there was a problem with your package, and that it was solved so > efficiently. But let me ask: is the "variable overflow" issue related > to bayesmlogit? It doesn't seem to be but I just wanted to be sure. > > And thanks Nick! > > Matt Baker > > On Wed, Aug 21, 2013 at 12:13 PM, Milena Przheska > <milena.przheska@cosmicdevelopment.com> wrote: >> Hi Nick, >> >> Thank you very much for your reply. >> The command did not return r(2000) after I encoded the variables, so >> your guess was spot on. >> >> It did return an error for variable overflow r(1400), but I know that >> is a totally different issue. >> >> I am relatively new to Stata, so I apologize if I violated some code >> of conduct or failed to supply appropriate reference in my previous >> message. >> >> In that regard, I would like to supplement my previous message with >> the reference to the package I used to run the command >> >> package bayesmlogit from http://fmwww.bc.edu/RePEc/bocode/b >> Author: Matthew Baker, Hunter College >> >> Once again, thank you very much. >> >> Best, >> Milena >> >> >> >> On Wed, Aug 21, 2013 at 5:16 PM, Nick Cox <njcoxstata@gmail.com> wrote: >>> -bayesmlogit- (SSC) is user-written. Please remember to explain where >>> user-written programs come from. >>> >>> At a guess, your -rand- is a string variable but -bayesmlogit- wants >>> it to be numeric. If so, use -encode- to produce a numeric variable >>> with value labels. >>> Nick >>> njcoxstata@gmail.com >>> >>> >>> On 21 August 2013 16:09, Milena Przheska >>> <milena.przheska@cosmicdevelopment.com> wrote: >>>> Dear Statalist members, >>>> >>>> I got error "no observations" r(2000) after running -bayesmlogit- >>>> >>>> I have a data set that contains choices of subjects at 4 separate >>>> occasions, and a dependent variable - a dichotomous variable. >>>> >>>> The data set looks like this: >>>> >>>> id time prior y >>>> 1 1 a 0 >>>> 1 2 a 0 >>>> 1 3 b 1 >>>> 1 4 c 0 >>>> 2 1 c 1 >>>> 2 2 c 0 >>>> 2 3 b 1 >>>> 2 4 a 0 >>>> 3 1 d 0 >>>> 3 2 d 0 >>>> 3 3 a 1 >>>> 3 4 c 1 >>>> >>>> The commands I ran were >>>> >>>> bayesmlogit y, group(time) identifier(id) rand(prior) >>>> >>>> >>>> and >>>> >>>> >>>> bayesmlogit y, group(id) identifier(id) rand(prior) >>>> >>>> >>>> Got the same error that there are no observations. >>>> >>>> >>>> Can anyone please explain why this error occurs? >>>> >>>> >>>> I appreciate your help and I thank you all in advance. >>>> >>>> >>>> Regards, >>>> Milena Przheska >>>> * >>>> * For searches and help try: >>>> * http://www.stata.com/help.cgi?search >>>> * http://www.stata.com/support/faqs/resources/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/faqs/resources/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/faqs/resources/statalist-faq/ >> * http://www.ats.ucla.edu/stat/stata/ > > > > -- > Dr. Matthew J. Baker > Department of Economics > Hunter College and the Graduate Center, CUNY > * > * For searches and help try: > * http://www.stata.com/help.cgi?search > * http://www.stata.com/support/faqs/resources/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/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/