APIM with count data and proportion data PreviousNext
Mplus Discussion > Structural Equation Modeling >
 Biru Zhou posted on Thursday, October 10, 2013 - 2:21 pm
Hello, I am very new to SEM and Mplus. I have some questions regarding anlayzing proportion data or count data in actor partner interdependence models (APIM). My questions are as follows:
1) What are the appropriate estimators for count and proprotion data for APIM in the context of SEM?

2) APIM with count and proportion data can also be analyzed in the MLM approach using link functions = log and logit repectively (Loeys et al., 2013). What would the interpretation be like in SEM with the appropriate estimators? Would it be the same as in MLM, the results can be interpreted in terms of odds ratios?

Many thanks!
 Linda K. Muthen posted on Friday, October 11, 2013 - 5:27 pm
1. Maximum likelihood estimation is used for count variables. See the COUNT option in the user's guide to see the available models.

Proportion=Count / N ... you can then use offset modeling with a count variable see the following post.


2. The results are in the log rate scale. When it is exponentiated, it is a mean.
 Biru Zhou posted on Wednesday, October 16, 2013 - 4:24 pm
Thank you for your reply, Linda. If I have two variables that are proportions (i.e., S2 and P2) and I want to run a simple APIM with S1 and P1 as IVs and S2, P2 as DVs, nS1 as the number of trails for S1 and nP1 as the number trials for P1, would the following syntax be correct:

USEVARIABLES ARE S1 S2 P1 P2 offset1 offset2;
COUNT = S2 P2;
DEFINE: offset1=log(nS1);
ANALYSIS: Type = General;
Estimator = mlr;
MODEL: S2 ON offset1@1 S1;
S2 ON offset1@1 P1;
P2 ON offset2@1 S1;
P2 ON offset2@1 P1;

Many thanks!


 Biru Zhou posted on Thursday, October 17, 2013 - 5:11 am
Hello, in my previous post, I realized that I missed one statement in the model:
P2 ON offset2@1 P1;

Since S2 and P2 are proportions, should I model is as:

offset1@1 WITH offset2@1 ?


 Bengt O. Muthen posted on Thursday, October 17, 2013 - 6:20 am
WITH statements are not allowed with count DVs. Instead, a factor can be defined by the 2 DVs; this lets them correlate beyond what their predictors explain. You should not use WITH for the offsets.
 Catheryn Koss posted on Wednesday, February 03, 2016 - 7:54 am
I am analyzing dyadic data using the APIM model to predict 2 dichotomous outcome variables (1 for husbands, 1 for wives). The data are dyadically structured with one row per dyad. To account for correlations between the DVs, I defined a factor as suggested in an earlier post which seems to work well. However, several sources on APIM suggest to also estimate covariances among the predictor variables using WITH statements. When I do this, my model fails to converge. Are these WITH statements necessary, or does the Mplus default of allowing predictors to covary already account for the nonindependence of the IVs? Thanks in advance for any advice you can give.
 Bengt O. Muthen posted on Thursday, February 04, 2016 - 6:51 pm
Predictors (independent, x vbles) are correlated by default even though these correlations are not estimated for predictors.
 Catheryn Koss posted on Wednesday, February 10, 2016 - 6:53 am
Thank you for the confirmation. I have one more question about modeling an APIM with categorical variables. In order to account for the covariance of the two dichotomous outcome variables, I am following advice given earlier and creating one factor loading on my two outcome variables (adH and adW). I have been defining the factor as follows:
f1 BY adH*(1) adW(1);

Does this seem correct? Thanks in advance.
 Bengt O. Muthen posted on Wednesday, February 10, 2016 - 1:17 pm
That's one way of doing it. It is ok because you achieve the goal of having a single parameter represent the residual covariance.
 Sarah Arpin posted on Saturday, February 24, 2018 - 10:41 am

I am trying to estimate a path model wherein my dependent variable is a proportion (# of insomnia symptoms/# days of survey data). Values on this variable are between 0 and 1, and are zero-inflated. As this DV is not count, and includes non-integers, what is the appropriate way of modeling this dependent variable? I've read a bit about censored-inflated estimators. Is this an appropriate way to handle this data, or is there an alternative estimator/approach you would suggest?
 Bengt O. Muthen posted on Saturday, February 24, 2018 - 2:01 pm
You can do censored, censored-inflated, two-part etc. This is described in our Regression and Mediation book in Chapter 7. See also our YouTube video and handout for Short Course Topic 11 on our website.
 Dora posted on Friday, May 11, 2018 - 12:47 pm
I am very new to SEM in Mplus.I would like to run the APIM with dummy coded categorical predictor variables. I followed the basic saturated APIM:
VARIABLE: names are tanmal tanfem konmal konfem;
usevariables are tanmal tanfem konmal konfem;
missing is all (999);
ANALYSIS: estimator = ml;
bootstrap = 5000;
MODEL: konmal on tanmal;
konmal on tanfem;
konfem on tanmal;
konfem on tanfem;
tanmal with tanfem;
konmal with konfem;
OUTPUT: stand;
I am doubting whether I have used the correct estimation for dummy coded predictors.What is the correct syntax that I should use to run the APIM?
 Bengt O. Muthen posted on Friday, May 11, 2018 - 1:35 pm
Which are your dummy coded predictors? tanmal and tanfem? How are they obtained/scored?
 Dora posted on Saturday, May 12, 2018 - 3:58 am
The categorical predictor variable has 4 categories (1=practical goals, 2=learning goals, 3=work goals, 4=health goals). Females and males were classified into these four categories. To be able to run regression we created dummy variables. In our example using the variable "goals", the first new variable called "tan" had a value of one for each observation in which goal is learning, and zero for all other observations. (tanmal refers learning goal for males, tanfem refers learning goal for females). Likewise, we created "mun" to be 1 when the person has work goals, and 0 otherwise, and "eg" is 1 when the person has health goals, and 0 otherwise. The level of the categorical variable that is coded as zero in all of the new variables is the reference level, or the level to which all of the other levels are compared. In our example, "practical goals" is the reference level.
In total we have 6 dummy coded predictor variables:
tanmal (learning goals for males)
tanfem (learning goals for females)
munmal (work goals for males)
munfem (work goals for females)
egmal (health goals for males)
egfem (health goals for females)
The dependent variable is relationship conflict:
konmal: relationship conflict for males
konfem: relationship conflict for females.
 Bengt O. Muthen posted on Monday, May 14, 2018 - 4:28 pm
With 2*4=8 categories, you should be able to handle 7 dummy variables. But this is a question better suited for SEMNET.
 Dora posted on Tuesday, May 15, 2018 - 11:29 am
Thank you for your help.
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