For the purpose of power analysis, I am attempting to use monte carlo simulation to generate a four category observed nominal variable. I will regress a continuous(observed) independent variable on this nominal variable. I am having trouble locating examples that show how to generate a multinomial independent variable for such a monte carlo simulation. Any suggestions on how to approach this?
So you want to regress y on c, where y is continuous observed and c is nominal observed.
You can do this in Mplus by working with a latent class variable that is identical to your nominal variable. How to do this is shown in a FAQ on our web site. The influence of c on y is portrayed by the y means changing as a function of the classes. You don't say "y ON c".
I have two follow up questions. Q1: Regarding the document that you referred me to, "Making an observed nominal variable u equivalent to a latent class variable c", I suppose that you are able to manipulate the threshold values that you specified to change the proportion of individuals that fall into each class. Would you give guidance or refer me to a reference that describes how to specify thresholds to create the nominal variable appropriately? Q2: What model setup would I use to portray the y means varying as a function of the classes? I hope to use monte carlo simulation to determine the power to detect significant differences in y across classes. Thank you.
1. You need to give the program the intercept values in a logit scale. Let's say you want the proportions of a three-category nominal variable to be .20, .30, and .50 where .50 is the reference class. You need to give logit values for the first two classes. You compute them as the log of the ratio of the probability of being in one class to the probability of being in the reference class, for example,
log (.20/.50) and log (.30/.50)
2. Add the y means to the MODEL and MODEL POPULATION commands. Example 8.6 has a distal outcome. Look at the Monte Carlo counterpart to that example.
mpduser1 posted on Thursday, June 02, 2011 - 8:14 pm
I am trying to perform a Monte Carlo power analysis for a multinomial logistic regression model in Mplus. I wish to have my dependent variable (Y) be a 3-level nominal categorical variable, while my predictor (X) is dichotomous. I've been experimenting with the Mplus syntax and have found something like the following to work:
NAMES ARE Y X; NOBSERVATIONS = 500; GENERATE = Y (n 3); NOMINAL = Y; CUTPOINTS = X(0);
MODEL POPULATION: X*.50 ; [X*.50] ; Y#1 on X*.3; Y#2 on X*.4; Y#3 on X*.5;
ANALYSIS: TYPE = GENERAL; ESTIMATOR = ML;
MODEL: Y#1 on X; Y#2 on X; Y#3 on X;
My question is why does Mplus require me to specify a regression model for Y#3 -- isn't Y=3 the baseline for my multinomial logistic regression model?