Power analysis PreviousNext
Mplus Discussion > Latent Variable Mixture Modeling >
 Bonnie J. Taylor posted on Friday, February 11, 2000 - 12:55 pm
I would like to perform a power analysis following the procedure of Satorra and Saris (1985). How can I output the model implied covariance matrix and mean vector?
 Bengt O. Muthen posted on Friday, February 11, 2000 - 4:25 pm
You do a run where you fix all the parameters to the desired values and request RESIDUAL to get the estimated mean vector and covariance matrix. In this run you can use means and covariance input with dummy values, for example using zeros for the means and an identity matrix for the covariance matrix.

We have a couple of pages that describe the steps for Satorra-Saris power calculation in a growth model setting which we discuss in our courses and we will shortly post this information on our web site. This includes how to get the power value from the non-centrality parameter (the "chi-square" value).
 Jichuan Wang posted on Tuesday, November 21, 2000 - 8:58 am
Dear Dr. Muthen:
In Muthen & Curran (1996) and Curran & Muthen (1996), parameter values in the growth model were chosen so that the difference between the treatment and control group means at the last time point scaled by the pooled s.d. at the last time point represented a effect size of 0.20.
Could you please explain how to shoose parameter values (e.g., slope mean and variance, correlation between the intercept and slope, et al.) in the first step for Satorra-Saris power calculation for a given effect size (e.g., 0.20)?
Thank you very much for your help.


Jichuan Wang, Ph.D.
School of Medicine
Wright State University
Dayton, OH 45435
 bmuthen posted on Tuesday, November 21, 2000 - 3:49 pm
You need to know how to calculate means and variances for the outcomes given a certain growth model with a certain choice of population parameter values. Given this, the simplest way is choosing parameter values that give the desired effect size using trial and error.
 Susan E. Collins posted on Wednesday, October 03, 2007 - 6:01 pm
Hi there,
I am new to mplus and am confused about setting the start values in monte carlo for an LTA model. I would like to run an mc for a power analysis and am trying to use example 8.13 as a guide.

I would like to see what power I have at N=300, 400 and 500 to test the possibility that x (continuous variable), is posititvely correlated with membership in high-risk drinking classes c2#1 and c1#1. I want to test the x effects as small, med and large effects over the different sample sizes. I originally wanted to use the classic cohen values, r=.1, .3 and .5 to represent small, med, and large effects, respectively, but realize that this is more like a logistic regression and am thus confused as to what start values I should use to test these different effect size levels.

In the example 8.13, you suggest "c2#1 ON c1#1*.5 x*1;" but I have no idea what x*1 would actually be in terms of an effect size. Can you tell me how best to represent different effect sizes in the mplus code?

I was also curious why "c2#1 on x*.2;" under the specific model "%c1#1%" while it appears differently in the %overall% command (i.e., c2#1 ON c1#1*.5 x*1;)

Thanks so much for your help!
 Bengt O. Muthen posted on Thursday, October 04, 2007 - 12:48 pm
Effect size has a different meaning when the dependent variable is categorical and is not as easily settled. You are interested in an effect of x on c in the multinomial logistic regression of c on x. So instead you have to ask yourself how much the probability of being in a certain c class changes as a function of x changing, say by 1 SD. You have to decide what a small/medium/large probability change is. Note, however, that the probability change is different for a 1 SD x change at low, medium, and high x values. You can use UG Chapter 13 formulas to compute the probabilities.

You may also take a look at logistic regression books to see if any counterpart to effect size appears, but I doubt that.

Regarding your last question, the overall c2#1 on x will apply to the last class whereas the c1#1-specific statement refers to the first class. This is how you represent the model picture's broken line.
 shahadut hossain posted on Friday, April 24, 2009 - 10:45 am

I ran a mixture model in mplus with two latent classes. I am now interested in simulating a 2 class model with 235 observations to see how much power Iíd have. I canít seem to get the thing to run.

Especially, for two classes I get two estimates for threshold (one for each class). My question is, how I can specify twi threshold GENERATE command under MONTECARLO.


Shahadut Hossain
UBC school of nursing
 Linda K. Muthen posted on Friday, April 24, 2009 - 11:16 am
Find the example in Chapter 7 that is like your model. Then use the Monte Carlo counterpart of that example as a starting point. If you run into problems, send your output and license number to support@statmodel.com.
 Owis Eilayyan posted on Monday, February 10, 2014 - 8:07 am

Can we do power analysis for path analysis model? I just saw examples for CFA and growth model.

Thank you
 Linda K. Muthen posted on Monday, February 10, 2014 - 10:51 am
You can do a power analysis for any model.
 Owis Eilayyan posted on Tuesday, February 11, 2014 - 12:13 pm
Thank you for the reply.
I have run power analysis, does "% Sig coff" indicate the power?

Thank you,
 Linda K. Muthen posted on Tuesday, February 11, 2014 - 1:27 pm
The column labeled % Sig Coeff gives the proportion of replications for which the null hypothesis that a parameter is equal to zero is rejected at the .05 level (two-tailed test with a critical value of 1.96). The statistical test is the ratio of the parameter estimate to its standard error, an approximately normally distributed quantity (z-score) in large samples. For parameters with population values different from zero, this value is an estimate of power with respect to a single parameter, that is, the probability of rejecting the null hypothesis when it is false. For parameters with population values equal to zero, this value is an estimate of Type I error, that is, the probability of rejecting the null hypothesis when it is true.
 Owis Eilayyan posted on Tuesday, February 11, 2014 - 6:45 pm
So when I get zero for all parameters means that my null hypothesis is true but I reject it. But I already got p-value<0.05 for some of these parameters!
Does that also mean the model does not have enough power to run with the current sample size?

Thank you,
 Linda K. Muthen posted on Wednesday, February 12, 2014 - 10:51 am
Please send your output and license number to support@statmodel.com. Be explicit about where in the output you are looking.
 Johanna van Rijn posted on Wednesday, May 07, 2014 - 2:07 am

I like to do a poweranalysis for my Latent Profile Analysis model. However, I'm using z-scores as input for the LPA. If I specify the power analysis with the values the LPA gave me, this error occurs:


Is it possible to do the power analysis? And what would be the best way to do this?

Thank you in advance!
 Linda K. Muthen posted on Wednesday, May 07, 2014 - 11:42 am
Please send the output and your license number to support@statmodel.com.
Back to top
Add Your Message Here
Username: Posting Information:
This is a private posting area. Only registered users and moderators may post messages here.
Options: Enable HTML code in message
Automatically activate URLs in message