Statalist The Stata Listserver


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

Re: st: GLLAMM Model convergence


From   [email protected]
To   [email protected]
Subject   Re: st: GLLAMM Model convergence
Date   Mon, 30 Apr 2007 13:59:15 -0600

I had the same experience about the lenght of processing time.

My model was much simpler than yours with much smaller sample size (about 9,000) with probably around 5 explanatory varibales in each level.

Once I incorporated one random coefficient in the model, the program run about a day but nothing came out. So, I terminated the program.

I hear people running random coefficient models (although using different programs) and let the program run for several days (or even a week). One paper mentioned that it took them about 40 hours or so to run a model. So I guess a long processing time seems not to be uncommon.

Keiko







Quoting "Richardson, Kelly K." <[email protected]>:


Hi Stata Listers,

I am trying to run a GLLAMM model with a random intercept for HRR and a
random coefficient for black race. I have a sample size of 798,565 with
98 HRRs.  There are 30 independent variables. Age and income are
continuously measured but all other variables are dichotomous including
the dependent variable. I thought about collapsing the data but I didn't
have any cases in which multiple patients had the same values on the
independent variables. I used the following code:


gen cons=1
eq hrr_cons: cons
eq hrr_black: race_black
gllamm died_admitp30 age female race_black zip_median~e elix_chf
elix_coagu~t elix_hyper~n elix_pvd elix_fluid~s elix_diabe~s elix_liver
elix_renal~l elix_copd elix_valve elix_obesity hannan97_cvd spec_prima~i
spec_cc_cr~h spec_ami_s~l spec_prevc~g spec_ami_i~t adtype_ele~e
adtype_urg~t adsour_tra~n adsour_eme~m yr9596 yr9798 yr9900 yr0102
yr0304, link(logit) fam(binom) i(hrrcode) nrf(2) eqs(hrr_cons hrr_black)
nip(4) adapt

The program ran for several days and finally the power went out in my
building and, well, I got nothin'. I'm wondering, before I start this
process again, am I running too large of a model? Is my sample size too
big? Is there something wrong with my code? I know the documentation
says these models take awhile to converge but should it be several days?
If my model is too large, what are the limits for GLLAMM models? Should
I do this model in HLM instead? Is it my machine?

I am working with:
Microsoft Windows XP
Professional
Version 2002
Service Pack 2

Dell Precision PWS670
  Intel(R)
Xeon(TM) CPU 3.80 GHz
2.77 GHz, 3.25 GB of RAM


Any advice you can give would be greatly appreciated.

Thanks,
Kelly
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Kelly K. Richardson, PhD

Data Analyst

Center for Research in the Implementation of Innovative Strategies in
Practice (CRIISP)

Iowa City Department of Veterans Affairs (VA) Medical Center (152)

601 Highway 6 West

Iowa City, IA  52246-2208

Phone:: (319) 338-0581 ext. 97641

Fax: (319) 887-4932

E-mail: [email protected]



*
*   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/




*
*   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/



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