Message/Author 

Anonymous posted on Sunday, June 27, 2004  1:07 pm



Is there anyway way to model multilevel data in discretetime survival analysis? Our sample consists of siblings coming from the same families. Thanks. 

bmuthen posted on Sunday, June 27, 2004  4:15 pm



Yes you can do this  in two ways. One way is to do a singlelevel analysis where you model the variables for all of the siblings (the sample size is the number of families). The other way is to do a twolevel analysis with siblings nested within families, using random effects (e.g. a random intercept) that vary across families 


Hi Bengt, as a followup to Anonymous's question, I assume that you would have to use the "f by event indicator notation" instead of the LCA parameterization. Would "f" then be allowed to vary across clustering units? Also, is there an example output file available? Thanks. Best, Hanno 

bmuthen posted on Wednesday, July 28, 2004  1:37 pm



Yes. Ex6.19 in the Version 3 User's Guide would have to be combined with the examples of the multilevel chapter 9, say ex9.6, but perhaps with only the fb factor on between. 


This is a somewhat naive question but will appreciate feedback and guidance. We are estimating twolevel survival analyses and need to get estimates of the random effects  variance, interquartile hazard ratio (HR) and median HR. This is discussed in the article: Chaix & Merlo. Am J Epidemiol 2005;162:171–182. My collaborators used a Bayesian approach in SAS to estimate the variance parameter. How do I do specify it in MPLUS to get estimates of the random effect? This is the current form of the model. I had also specified it as a discrete time model. WITHIN = ... ; BETWEEN = x2x5; CLUSTER = tractid; SURVIVAL = p_yralld (ALL); TIMECENSORED = deadcnsr (1=NOT 0=RIGHT); ANALYSIS: TYPE = TWOLEVEL; BASEHAZARD = OFF; ALGORITHM = INTEGRATION; INTEGRATION = MONTECARLO; ESTIMATOR=ML; PROCESSORS=2(STARTS); MODEL: %WITHIN% p_yralld ON entr_age sex bmicur marriage ltcoll coll raceblk raceoth diabetes heart stroke smoke logcal logfiber fairpoor; %BETWEEN% p_yralld ON x2x5; 


In your setup you estimate on Between the residual variance of the random intercept for the p_yralld survival variable. If you delete x2x5 as covariates on Between, you will estimate the variance of the random intercept. 


Thank you so much. I modified the model as follows and used MLR: BETWEEN = x2x5; ...... ANALYSIS: .... ESTIMATOR=MLR; MODEL: %WITHIN% p_yralld ON entr_age ; %BETWEEN% p_yralld ON ; 1.) The standardized variances were 1.0 and se=0. Should I expect that? 2.) Also, I modeled quintiles of the BETWEEN variable. What would be your advice for deriving the IHR, and Median HR? 3.) The formula in Chaix & Merlo's article seems to be based on continuous rather than dummies. Would I need to use the continuous BETWEEN variable rather than the "quintile" variables? 


Regarding number 1, please send the output and your license number to support@statmodel.com. Regarding 2, 3, and 4 I am not familiar with the paper. 


I think I figured it out. I was looking to derive the random effects (frailty) parameter with Cox frailty models. This is not possible in the multilevel framework  is it? 


Yes, you can have frailties in the Cox model, and on both levels. 


Hello, I am working on a multilevel survival analysis using cox regression (continuous time survival). The researcher I am working with have found meaningful person level predictors of returning to hospitalization (only the first return to treatment). The researcher would like to test to see if there is a random effect of hospitals, in particular to determine if there are different survival rates between hospitals. I have received some recommendations that it may be best to dummy code each hospital (N=30) and include them in the analyses as a person level predictor, but I have not seen any literature to necessarily support that approach. I am wondering if there is a way to test if there are differences between hospitals (between level variables) on survival. An added difficulty is that there are no other hospital level predictors in the model. Here is the syntax for the model I have proposed to test, but I realize that the variance estimation at the between level does not answer my question if survival differs as a function of hospitals: ANALYSIS: TYPE = TWOLEVEL; BASEHAZARD = off; ALGORITHM=Integration; MODEL: %WITHIN% Days_Ret ON age race11 race3 R_arr3 ExtBeh Intern CareIss1 los_627; %BETWEEN% Days_Ret; Thank you. 


Either fixed effect model with dummy variables or random effect model that you have above should illuminate the issue. 


Hi Bengt and Tihomir, as a followup to this conversation I have another question. I've never used Mplus and it can be my sofware for the future. I'm working on a multilevel frailty model in order to study survival of children. The structure implies a household (second) and regional (third) level. The measurement model is a 2level model with a latent variables at the household and regional level. These is then included as predictors in the 3level structural model for mortality I wonder if there is the possibility to include the variance of the householdlevel latent variable as a predictor in the structural model. In the case it could be done, is it feasible to include it in both the two and threelevel models? Any advice would be extremely useful! Thank you! 


You will need to use type=twolevel and the household variables have to be setup as a multivariate vector (wide format). Currently only twolevel modeling is available in Mplus with survival variables. To use the variance variable on the higher level, see section 4.1 in https://www.statmodel.com/download/NCME12.pdf or see slides 24 and 28 (two different approaches) in https://www.statmodel.com/download/handouts/MuthenV7Part2.pdf 

Back to top 