Anonymous posted on Thursday, November 03, 2005 - 11:35 am
I'm interested in simulating a multilevel SEM. I want to control my Intraclass correlations. How would I do that in Mplus? In other words, what syntax would I use? If I missed an example in the manual, please tell me the number.
Also, is there any way to have the "results" option of Monte Carlo provide the ICC's?
I am going to use Example 11.4 to describe this. An intraclass correlation is the ratio of between variability to total variability. So you increase the intraclass correlation by increasing between level variability. For example, consider y1. It's between-level variation is the same as the variability of ib because centering is done at time one. The within-level variance of y1 is the same as the iw variance plus the y1 residual variance. Using this example, the intraclass correlation for y1 is 0.2/(0.2 + 1 + 0.5) = .12.
I want to simulate a multilevel CFA model. In within model, I have 6 indicators (y1-y6). Half of them(y1-y3) were predicted by a latent variable, namely, fw1; half(y4-y6) were predicted by fw2. In between model, I have the same factor structure. Half of them(y1-y3) were predicted by a latent variable, namely, fb1; half(y1-y3) were predicted by fb2.
Q1. What is the syntax to set ICC to be equal to .25? Q2. Are the reliabilities of each indicator in between model necessarily larger than the reliabilities in within model? Is there any suggestion/reference when I set the residuals of each indicator in between model?
Q1. ICC is the between variance divided by the between plus within variances. So adjusting the between variance changes the ICC. See the input from the Mplus CD for Example 11.4. There are comments related to ICC's.
Q2. I think this is usually the case. I don't know of any reference.
I just wanted to follow up on these questions and responses, which have been most helpful. If I'm simulating multilevel data and want to manipulate the ICC's for the latent variable, could I just do that by specifying the between and within variances for it only, or would I also need to do so for the individual indicators as well? I'm thinking of a case where I want the ICC=.25, for example. It seems to me that I could set the within variance=.25 and between=.75 for the latent variable and generate what I need. Does that logic seem sound? Thanks.
I'm interested in simulating multilevel data("disaggregated approach", per Bengt's earlier work) and then analyzing these data such that grouping is nuisance ("aggregated approach" per Bengt's earlier work....robust ML, group = cluster).
Specifically, I would like to simulate a simple 2 level model in which person level data (level 1; Sw) are predicted by a single group level (level 2; Sb) predictor (e.g,. child outcomes predicted by treatment that is assigned at classroom level).
1. Can you point me to example code where this has been done?
2. When defining population values in the monte carlo study, should the group level effect be included in within (Sw), between (Sb), or both levels models?
I'm interested in simulating a multilevel, multigroup IRT model in which I control the intraclass correlations for one or more of the observed indicator variables. I've looked through several examples and can't quite see how to put that together. Could you steer me in the right direction? Thanks very much.
You can combine the Monte Carlo counterparts of Example 9.11 which is a two-level multiple group CFA with continuous factor indicators and Example 5.5 which is an IRT model. In Example 9.11, you will need to use the KNOWNCLASS option and TYPE=MIXTURE instead of the GROUPING option.
Thanks very much. I have one other question. Is it possible to assign the size of each of the latent classes analogous to using the ngroups and nobservations options for observed group sizes? Thanks very much.
You need to do that using the means of the latent class variable and then an observed variable that is the same as the latent variable. See the Monte Carlo counterpart of Example 8.8 which has a KNOWNCLASS variable.
QianLi Xue posted on Monday, October 11, 2010 - 4:26 pm
Hi, if I want to simulate a dataset based on a latent class model, say, with 4 classes. How can I specify latent class prevalence parameters in the Monte Carlo step? Could you please elaborate on your answer to the similar question above?
ICC is directly related to the model parameters only for normally distributed variables in principle and not for survival variables. So it is not possible to enter ICC directly in the simulation.
The ICC will increase as you increase the total variance of the random intercept. Consider the model below similar to example 12.10
%BETWEEN% w@1 (var1); t ON w*.2 (beta); t*1 (var2);
This gives a random intercept variance of
You can increase any of the 3 parameters to get higher random intercept variance and higher ICC for the observed T. For each set of parameters compute the ICC of the generated data until you get to the right level (guess and check approach).
QianLi Xue posted on Wednesday, September 21, 2011 - 12:49 pm