Email updates

Keep up to date with the latest news and content from BMC Medical Research Methodology and BioMed Central.

Open Access Research article

Non-linear mixed models in the analysis of mediated longitudinal data with binary outcomes

Emily A Blood12* and Debbie M Cheng1

Author Affiliations

1 Department of Biostatistics Boston University School of Public Health 801 Massachusetts Avenue 3rd Floor Boston, MA 02118 USA

2 Clinical Research Program, Children's Hospital Boston and Harvard Medical School, 300 Longwood Avenue Boston, MA 02115 USA

For all author emails, please log on.

BMC Medical Research Methodology 2012, 12:5  doi:10.1186/1471-2288-12-5

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2288/12/5


Received:15 July 2011
Accepted:24 January 2012
Published:24 January 2012

© 2012 Blood and Cheng; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Structural equation models (SEMs) provide a general framework for analyzing mediated longitudinal data. However when interest is in the total effect (i.e. direct plus indirect) of a predictor on the binary outcome, alternative statistical techniques such as non-linear mixed models (NLMM) may be preferable, particularly if specific causal pathways are not hypothesized or specialized SEM software is not readily available. The purpose of this paper is to evaluate the performance of the NLMM in a setting where the SEM is presumed optimal.

Methods

We performed a simulation study to assess the performance of NLMMs relative to SEMs with respect to bias, coverage probability, and power in the analysis of mediated binary longitudinal outcomes. Both logistic and probit models were evaluated. Models were also applied to data from a longitudinal study assessing the impact of alcohol consumption on HIV disease progression.

Results

For the logistic model, the NLMM adequately estimated the total effect of a repeated predictor on the repeated binary outcome and were similar to the SEM across a variety of scenarios evaluating sample size, effect size, and distributions of direct vs. indirect effects. For the probit model, the NLMM adequately estimated the total effect of the repeated predictor, however, the probit SEM overestimated effects.

Conclusions

Both logistic and probit NLMMs performed well relative to corresponding SEMs with respect to bias, coverage probability and power. In addition, in the probit setting, the NLMM may produce better estimates of the total effect than the probit SEM, which appeared to overestimate effects.

Background

SEMs are a general modeling framework often used in the social sciences to analyze complex relationships between variables, such as mediated relationships between variables. A mediator is a variable in the causal pathway between a predictor and the outcome of interest. SEMs are becoming more common in the clinical research setting and can be used to model hypothesized causal pathways between variables of interest. Extensions of SEMs have been developed to allow for more general types of dependent variables, including binary outcomes [1]. Common statistical techniques for non-mediated longitudinal binary data include non-linear mixed models (NLMM) [2] and generalized estimating equations (GEE) [3]. When interest is primarily in the total effect of a predictor on an outcome, even if mediation may be present, these commonly used techniques may be preferred over SEMs as they specify straightforward predictor-outcome variable relationships and do not require specialized software, as the SEM often does. It is therefore of interest to determine, in a setting conducive to using SEMs, whether a method such as NLMMs adequately models the total effect of a predictor on binary outcomes without directly modeling mediation. We focus on NLMM rather than GEE in this paper as it is more similar to the non-linear SEMs for longitudinal data available in SEM software-both are conditional rather than marginal models.

Comparisons have been made between SEM and other statistical models in different contexts [4-13]. Mixed effect models have been evaluated against SEMs with continuous data [14,15], and found to adequately model mediated predictor-outcome relationships. MacKinnon et al. [16] examined the calculation of mediated effects in cross-sectional binary data with non-SEM techniques using two different methods (difference of coefficients and product of coefficients). While, Palta and Lin [17] compared structural equation models to various marginal models in longitudinal binary data without mediation. To our knowledge, evaluation of NLMMs relative to SEMs has not yet been performed in the context of mediated longitudinal binary data.

Linear and non-linear mixed models differ both in terms of the distributional assumptions and the estimation techniques used for inference. In addition, the parameter estimates in non-linear mixed models using a logit or probit link are inherently scaled to the predictors (and mediators) included in the model. Therefore, comparisons of parameter estimates between NLMMs with different sets of predictors must first be re-scaled in order to make them comparable [16].

In this paper, we evaluate the performance of NLMMs relative to SEMs for the modeling of mediated, binary longitudinal data in a setting where the SEM is presumed to be optimal. The purpose is to assess whether there is an impact of direct modeling of causal pathways in terms of bias, power, and coverage probability when the goal is to determine the total effect of the main independent variable. A simulation study is performed to assess these two classes of models across a variety of settings. We also describe, in an appendix, two different approaches for rescaling estimates when analyzing real world data in order to allow direct comparisons between NLMMs and SEMs or to compute mediated effects via NLMMs only.

Methods

In the current study, we consider a longitudinal data setting with binary outcomes, a repeated binary predictor, a repeated continuous mediator, and a continuous covariate measured at baseline. An example of such a clinical setting would be a prospective cohort study evaluating the impact of heavy alcohol consumption on HIV disease progression, defined as low CD4 cell count (e.g. <350 cells/μL). Heavy alcohol consumption may influence progression of HIV, while also influencing adherence to anti-retroviral therapy (ART). Level of adherence to ART is also a predictor of HIV disease progression. In this setting there is a repeated binary independent variable of primary interest, heavy alcohol consumption (zj), and a longitudinal binary outcome, low CD4 cell count (Yj) signifying HIV progression. In addition, ART adherence (Mj), a continuous mediating variable, is measured repeatedly, and age (w) is a continuous covariate assessed at baseline. ART adherence is said to be a mediator because the primary independent variable, heavy alcohol use, may affect CD4 count directly as well as indirectly through ART adherence. We arbitrarily assume six time-points at which the predictor, outcome and mediator are measured. Time is represented by tj with j = 1, 2,..., 6. In this setting, we considered measurement times to be equally spaced and the same for all individuals. We generated data with a mediated non-linear relationship between the predictor (heavy alcohol consumption) and outcome (low CD4 cell count), i.e. we allowed the mediator (ART adherence) to be directly affected by the predictor and the outcome to be directly affected by both the predictor and mediator. Both the probit and logit links were assessed. We also describe the application of these models to data from a prospective cohort study evaluating the impact of heavy alcohol use on HIV disease progression.

As described by Fitzmaurice, Laird, and Ware [18] and others, binary outcome models can be described equivalently in two ways. The first approach would be to define a linear function of an underlying latent continuous variable (Y*) that when dichotomized represents the observed binary outcome (Y). For example, we could define a continuous latent (unobserved) variable Y* such that observed Y is 1 if Y* > 0 and 0 otherwise and write:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M1">View MathML</a>

(1)

A second approach would be to define a non-linear model of the probability of a binary response. If we consider ϵ ~ N(0,1), the model in Equation 1 defines a univariate probit model that can be equivalently represented using the following non-linear link format:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M2">View MathML</a>

Likewise, if we could consider the errors in Equation 1 to have standard logistic distribution (mean of 0 and variance of <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M3">View MathML</a>) the model defines a univariate logistic model that can be represented as:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M4">View MathML</a>

In more complex situations, such as the longitudinal data we are studying, the same equivalence between model descriptions exist and we use both model formulations for the NLMMs and SEMs that follow. The convention for binary or categorical outcomes in SEMs has been to describe binary regression models with the latent variable format while the NLMMs are often defined using the non-linear link format.

SEM

To evaluate the performance of NLMMs in a setting conducive to the use of SEMs, we generated mediated longitudinal binary outcomes using a non-linear SEM. We then fit the data with a NLMM as well as the non-linear SEM to evaluate the performance of the NLMM relative to the SEM. The non-linear SEM used to generate the data and subsequently fit to the generated data is described below.

Following the notation from above, xj is the independent variable of primary interest, Mj is the continuous mediating variable, w is a continuous time-invariant covariate, and tj represents time-point. Using the latent variable notation, we define a continuous unobserved outcome <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M5">View MathML</a> that takes a value of 1 only if <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M6">View MathML</a> for j = 1 to 6. This model can be expressed as follows (dropping the subject index i for simplicity), where:

Measurement model

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M7">View MathML</a>

(2)

Just as in the simpler models above, if we assume ϵj ~ N(0,1) this defines a probit model and if we assume ϵj ~ Logistic (<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M8">View MathML</a>) this defines a logit model.

Structural model

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M9">View MathML</a>

(3)

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M10">View MathML</a>

(4)

For j = 1 to 6,

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M11">View MathML</a>

(5)

where Ui1 represents a latent intercept, Ui2 represents a latent slope, zij represents the repeated binary predictor and Mij represents the repeated continuous mediator. The errors in the structural model are normally distributed with cov(ζ1, ζ2) = Ψ, cov(ζ3 : ζ8) = Θ and ζ(2+j) ~ N(0, θ). This model can be represented in a path diagram (Figure 1), a visual display of the interrelationships between variables typically presented along with SEMs.

thumbnailFigure 1. Path diagram. The non-linear structural equation model defined in Equations 2 to 5.

The parameters of the SEM defined in Equations 2 - 5 include: λ, which represents the effect of the repeated mediator on the repeated outcome; γ1, which represents the effect the repeated primary independent variables on the repeated mediator; γ2, which represents the effect of the continuous covariate on the repeated outcomes; and κ, which represents the effect of the repeated independent variable on the repeated outcome.

In this simulation study we focused on the total effect of the repeated binary predictor and the repeated binary outcome, which is represented by λγ1 + κ. The interpretation of the parameters of this model is subject-specific since it represents the effect of a predictor on the outcome when the individual intercept, individual slope and mediator value are held constant.

When the structural model Equations (3-5) are substituted into the measurement model Equation 2, the full model can be rewritten as:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M12">View MathML</a>

where ωj = (α1 + λ α3) + γ2w + α2 tj + (κ + λγ1) zj. The following presents the non-linear link formats for the probit and logit SEMs where the structural equations have been substituted into the measurement equation (the subject index i has again been dropped for simplicity):

Probit SEM

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M13">View MathML</a>

(6)

Logit SEM

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M14">View MathML</a>

(7)

To fit these models, Mplus uses maximum likelihood estimation when a logit link is used and weighted least squares estimation with a robust estimation of standard errors (WLSMV) when a probit link is used [19].

Non-linear mixed effects model

The following NLMM was evaluated in comparison to the SEM:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M15">View MathML</a>

(8)

where b1 is a random individual intercept and b2 is a random individual slope. Since the objective is to evaluate the total effect of the main independent variable, the mediator is excluded from this model [20]. The regression coefficient associated with the primary predictor (β3) therefore represents its total (i.e. direct plus indirect) effect on the outcome [14].

Probit NLMM

The probit model assumes that ej ~ N(0,1) and can be written as:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M16">View MathML</a>

(9)

where νj = β0 + β1 w + β2 tj + β3 zj.

Logit NLMM

the logit model assumes that ej ~ Logistic (<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M17">View MathML</a>) and can be written as:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M18">View MathML</a>

(10)

These models can be fit with SAS PROC NLMIXED which estimates parameters via maximum likelihood [21]. We note that the regression coefficients of the NLMM are interpreted conditional on the random individual intercept and random individual slope, but marginal on the residual error of the mediator (since the mediator is not included in the model).

Comparing NLMMs to SEMs

As noted previously, the SEM and NLMM condition differently on the mediating variable. Specifically, the SEM conditions on the random intercept and slope as well as on the residual variance of the mediating variable, while the NLMM conditions only on the random intercept and random individual slope. Thus estimates from the two types of models are not directly incomparable. Instead, to compare parameters from the NLMM to that of the SEM, we must first re-scale the regression coefficient from the NLMM so that it represents the effect of the primary predictor variable zj conditional on the mediator. To determine the scaling factor, we rewrite the SEM (for both the probit and logistic models) conditional only on the variance of the random intercept and slope to mimic the conditioning in the NLMM.

Comparing probit models

For the probit SEM, we generated the data according to the model described in Equations 6 and 7. Conditioning only on the variances of the random intercept and slope (ζ1 and ζ2), but not on the variance of the mediator (ζ2+j), it can be shown that:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M19">View MathML</a>

(11)

The sum of terms on the left-hand side of the inequality do not have a standard normal distribution since λζ2+j is added to ϵj which itself has a standard normal distribution. In order to express the probability in Equation 11 using the standard normal cumulative probability function, we re-scale the terms on either side of the inequality by the standard deviation of ϵj +λζ2+j to create a standard normal random variable:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M20">View MathML</a>

Conditioning on only the variance of the random individual intercept and slope, all regression coefficients are divided by the factor <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M21">View MathML</a>. For example, the regression coefficient associated with zj, which was κ + λ γ1, is now <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M22">View MathML</a>. Thus, the model parameters from the SEM are scaled to the variance of ϵj + λζ2+j which is 1 + λ2 θ and the model parameters from the NLMM are scaled to the variance of ϵj which is 1, resulting in a scaling factor of <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M23">View MathML</a> Parameter estimates from the SEM and NLMM must be on the same scale before making direct comparisons. For example, the total effect of the main independent variable from the probit NLMM, β3, which is also conditioned only on the random individual intercept and slope (Equation 9) should be multiplied by a factor of <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M21">View MathML</a> before it is compared to the total effect from the probit SEM, κ + λγ71. Direct comparisons of parameter estimates from the NLMM to those from the SEM without first re-scaling would underestimate effects by a factor of <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M21">View MathML</a>. In the current study we present the conditional total effect estimates from the SEM and compare them to scaled and unscaled NLMM estimates. Note that in the analysis of real (i.e. non-simulated) data, true parameter values are unknown and therefore must be estimated. We describe in the appendix two approaches for rescaling estimates in practice to allow direct comparisons between NLMMs and SEMs or to compute mediated effects via NLMMs only.

Comparing logistic models

Unlike the probit model, when the logit SEM is conditioned on only the random intercept and slope, the true relationship between the predictor and the outcome no longer follows a logistic model. That is, the distribution of the terms on the left-hand side of Equation 11 in a logit SEM does not follow a logistic distribution since the sum of a normal random variable (ϵj) and logistic random variable (ζ2+j) does not follow a logistic distribution. The result of this is that the scaled coefficients from the logit NLMM only approximate the mediated relationship described in a logit SEM. A similar situation occurs, for example, when comparing a non-linear mixed model to a non-linear generalized estimating equation as noted by Fitzmaurice, Laird, and Ware [18].

The scale factor for the logit model is created in the same was as it was for the probit model. The regression coefficient representing the total effect of the main independent variable (β3) from the logit NLMM, can be multiplied by the standard deviation of ϵj + λζ2+j and divided by the standard deviation of ϵj. The scaling factor for the logit model is therefore: <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M24">View MathML</a>.

Simulation plan

Data generation and model fitting

Because the goal of this study was to evaluate the performance of NLMM relative to SEMs in the setting where the SEM is presumed to be optimal, the SEM framework was used to generate the mediated binary data for the simulation studies. Data were generated according to Equations 2-5. For the probit model, errors in Equation 2 were assumed to be independent standard normal random variates. For the logistic model, the errors in Equation 2 were assumed to be independent standard logistic random variates. Data generation was repeated to create 1000 datasets. NLMM were fit with SAS (Version 9.2) PROC NLMIXED and SEM were fit with Mplus (Version 5.2).

Simulated data scenarios

We evaluated the performance of NLMM against SEM across several scenarios by examining the following:

• Sample size: ranging from 100 to 1000. The range of sample sizes was chosen to evaluate sample sizes that achieved adequate power with a moderate effect size.

• Effect size: ranging from 0.2 to 0.5. The range of effect sizes represent small to moderate effect sizes as defined by Cohen [22].

• Distribution of effects: three cases were evaluated: equally distributed direct and indirect effects, primarily direct effects of the main independent variable, and primarily indirect effects of the main independent variable. The total effect sizes (0.3 for probit models and 0.4 for logit models) were chosen such that adequate power was obtained when direct and indirect effects were equally distributed.

Model performance was assessed based on the following: 1) Bias- the difference between the true parameter value and the mean observed parameter value divided by the true parameter value; 2) Coverage probability- the percentage of the 1000 95% confidence intervals that contained the true parameter value; 3) Power - the percentage of the 1000 datasets in which the null hypothesis that the total effect of the main independent variable is equal to zero was statistically significant.

Results

Logistic link

The results evaluating the conditional total effect from simulations evaluating sample size are displayed in Table 1. For the logistic model, a sample size of 700 was required to obtain adequate (83%) power to detect a moderate effect size (0.3) for both the SEM and NLMM. Across all sample sizes, the estimates of power for the scaled NLMM were comparable to the SEM. As expected, for the unscaled NLMM, effects were underestimated in all cases. For example with a sample size of 600, the bias of the SEM was -0.6% while the bias for the unscaled NLMM was -3.8% and the bias of the scaled NLMM was -1.6%. Once the scale factor was applied to the NLMM, however, the estimated bias decreased and the power and coverage probability estimates were similar to that of the SEM. The results for the unscaled NLMM are therefore not described in subsequent tables as they are not directly comparable to the SEM results.

Table 1. Impact of sample size. Based on 1000 simulated datasets with moderate effect size (0.3) equally distributed between direct and indirect effects. Impact of sample size on model performance in evaluating the total effect of the repeated independent variable on the repeated outcome.

The results of the simulations in which the effect size (Table 2) and effect distribution (Table 3) were varied also demonstrated that the appropriately scaled NLMM produced comparable estimates of the total effect of the exposure relative to the SEM. The scaled NLMM results were similar to the comparison SEM with low bias at all effect sizes (ranging from -2.0% to 0.4%) and effect distributions (ranging from -1.4% to -1.0%). As expected, the models showed increasing power with increasing effect sizes (e.g., from 37% power for an effect size of 0.2 to 98% power for an effect size of 0.5) and power was highest (91%) when the effect distribution was primarily direct.

Table 2. Impact of effect size. Based on 1000 simulated datasets with sample size of 500 equally distributed between direct and indirect effects. Impact of effect size on model performance in evaluating the total effect of the repeated independent variable on the repeated outcome.

Table 3. Impact of effect distribution. Based on 1000 simulated datasets with sample size of 500 and effect size of 0.4 for the logit link and effect size of 0.3 for the probit link. Impact of effect distribution on model performance in evaluating the total effect of the repeated independent variable on the repeated outcome.

Probit models

For the probit model, the SEM showed consistent positive bias (i.e., no simulation scenario with the probit model resulted in a negative bias for the probit SEM). For example, the estimated bias for a sample size of 100 was 51.5% and decreased to 7.5% with a sample size of 1000. In comparison, the scaled NLMM had bias ranging from 2.7% to -2.1%. Notably, the estimated power for the NLMM was consistently higher than that of the comparison SEM. In the effect size simulation scenarios, bias in the SEM appeared to increase with effect size (effect sizes of 0.2, 0.3, and 0.5 resulted in biases of 7.6%, 8.8%, and 11.7%, respectively). In contrast, the probit NLMM showed relatively small bias (-0.7% to 1.7%) for all effect sizes. For the probit model, the scaled NLMM generally performed better than the SEM, across a range of sample sizes, effect sizes, and effect distributions, with higher estimated power and lower bias.

Positive bias in probit models

We explored possible explanations of the unexpected positive bias observed for the probit SEM models. We repeated simulations by first eliminating mediation from the simulated scenarios to evaluate whether the positive bias in the SEM persisted in unmediated settings. In addition, we assessed different estimation methods, including WLSMV (implemented in MPlus software) and also maximum likelihood using iteratively re-weighted least squares (implemented in Splus and SAS). We found that the mean regression coefficient estimates using WLSMV had larger bias in nearly all cases compared to estimates using maximum likelihood and iteratively reweighted least squares (Table 4). The positive bias for the WLSMV method was avoided only with a sample size of 5000. Our results suggest that the WLSMV, which is used by Mplus, may produce positively biased results. This bias was also present, however to a smaller degree, in the NLMM fit with maximum likelihood.

Table 4. Univariate Probit Model Results

Real data example: alcohol and HIV disease progression

To demonstrate the application of both the logit and probit NLMMs and SEMs evaluated in the simulation study, we analyzed data from a prospective cohort study evaluating the effect of alcohol use on HIV disease progression. Samet et al. have previously reported the analyses from this longitudinal cohort study [23]. The original analyses combined data from two cohorts (the HIV-ALC and HIV-LIVE cohorts), however, to illustrate the models evaluated in this paper, we have used data from the HIV-LIVE study only. For clarity of presentation, we limited the analyses to subjects who reported any ART use during follow-up, had complete data on the first four time-points (as Mplus and SAS have different methods for handling missing data in these models), and examined only the following key variables: heavy alcohol consumption (yes vs. no), the main independent variable; ART adherence (percentage of pills taken in the last three days), the mediator; age, a potential confounder; and low CD4 cell count (dichotomized at <350 cellsL), the primary outcome. Each variable was assessed every six months for up to four years, however for the current example only the first four time-points were analyzed in order to maximize the number of subjects with complete data. The resulting dataset was composed of 98 individuals contributing 392 observations. The total effect of heavy alcohol consumption on low CD4 cell count was not significant for any of the SEMs or NLMMs fit to the data. For the logit SEM, the total effect estimate (SE) was 0.554(1.246) with an associated p-value of 0.66. The scaled result from the logit NLMM was similar: estimated total effect (SE) = 0.5107(0.701), p = 0.47. However for the probit link, the estimated total effects (SE) from the SEM and NLMM appeared to differ substantially (probit SEM: 6.287(52.661),p = 0.91; scaled probit NLMM: 0.303(0.391), p = 0.44). Thus consistent with the results from the simulation study, the logit SEM and NLMM produced similar estimates in the real data example, whereas the probit SEM produced estimated effects that appeared much larger in magnitude in comparison to the probit NLMM.

Discussion

The purpose of this study was to evaluate the performance of NLMMs relative to SEMs in the analysis of mediated longitudinal binary outcomes in a setting where the SEM is presumed to be optimal. We found model performance differed based on the link function that was used in the non-linear portion of the models. Based on simulations performed across a variety of settings, the logistic NLMM performed well with respect to bias, coverage probability and power relative to the logistic SEM. The results were similar for the SEM and scaled NLMM in the logistic model setting, with both accurately estimating the effect of the time-dependent predictor on the longitudinal binary outcome. Application of these techniques to a real-date example from a prospective cohort assessing the effect of heavy alcohol consumption on low CD4 cell count also illustrate the similarity of results from the logit SEM and NLMM.

For the probit model, however, the SEM consistently overestimated the total effects of the predictor and generally had larger bias and lower power compared to the NLMM in both mediated and non-mediated data. The larger bias may be due to the weighted least squares estimation method used for the probit SEMs (fit with Mplus), which differs from the maximum likelihood estimation method used for the NLMM (fit with SAS) and for the logistic SEM (fit with Mplus). In contrast to the SEM, the scaled probit NLMM had good performance (low bias and high power and coverage probability) with adequate sample sizes. Similar results were observed in the real data example where estimates from the probit SEM appeared larger than those from the probit NLMM.

The results showing similar estimated effects for the SEM and NLMM in the logistic model setting are similar to results seen in the non-mediated case where SEM was compared to the generalized estimating equations (a non-SEM) technique. Palta and Lin [17] compared probit models for SEM and generalized estimating equations in the analysis of data from a cohort study and found that when appropriately scaled, the two models yielded similar results. They noted, however, that the SEM allowed for more flexible specification of variance structure and therefore allowed coefficients to be scaled to provide marginal or cluster-specific interpretation.

To obtain scaled NLMM estimates in practice, it may be preferable to model the mediation by fitting separate equations, one for each pathway, using maximum likelihood rather than weighted least squares (the only estimation method currently available for probit link models in MPlus). The potential burden of fitting multiple equations separately rather than simultaneously using SEMs may be outweighed by the benefit of using maximum likelihood estimation which, in the probit model simulations, appeared to produce less biased results. In addition, estimating scale parameters and using the product of coefficients method appears to produce acceptable estimates of the total effect of the exposure. Our study demonstrated that results using this approach were similar to those obtained when NLMM results were scaled using true parameter values. If indirect effects are of interest and the NLMM is used to analyze the mediated longitudinal binary data, scaling will also be necessary. Unlike the case with linear models for continuous outcome data, the product of coefficients method is not equivalent to the "difference of coefficients" method of determining the indirect or mediated effect [24] in the case of binary outcomes. Using the difference of coefficients approach in linear models, the total effect is obtained by fitting a model that excludes the mediating variable and the direct effect is obtained by fitting a model including the mediating variable. The indirect effect is then determined by taking the difference between the total effect and direct effect. However, in the binary case, the scale of the direct effect obtained from a model that includes the mediating variable is different from the scale of the total effect obtained from a model that excludes the mediating variable [16]. As demonstrated by MacKinnon et al., to obtain comparable estimates of the indirect effect in binary outcome models, the total effect must be appropriately scaled before the difference is taken.

This study presents results based on simulated data from a single-mediator model. Conclusions from these results may not be generalizable to scenarios with different data characteristics. For example, in scenarios with multiple mediators and pathways, the advantages and disadvantages of NLMMs relative to SEMs may differ. The performance of NLMMs and SEMs in other scenarios, such as the analysis of nominal and ordinal outcomes as well as the case of multiple mediators, should be evaluated in future studies.

Conclusions

Overall, we found the NLMM performed sufficiently well in the analysis of mediated longitudinal binary outcomes with respect to bias, coverage probability, and power. Under the logistic model, both the NLMM and SEM had acceptable performance and the results for the two types of models were similar. The NLMM requires scaling of the regression parameters and this scaling requires fitting additional models to separately estimate direct effects of the predictor, and effects of the primary predictor on the mediator. An advantage of the SEM is that it can fit all of the linear and non-linear models simultaneously, avoiding the burden of fitting multiple models. For the probit model, however, the SEM estimated using weighted least squares may overestimate effects. In contrast, the NLMM appears to perform adequately across a range of settings and therefore is preferred over the SEM for probit models.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

EAB was involved in the conception of the study, designing and performing the simulation analysis and drafting the manuscript. DMC was involved in conception of the study, designing the simulation analysis and drafting the final manuscript. Both authors read and approved the final manuscript.

Appendix

Estimating the Scaling Factor for Total Effect Estimates

Scaling effect estimates is necessary to obtain total effect estimates that represent the total effect of the predictor on the outcome, conditional on the mediation. This is the case when one wishes to compare estimates from SEMs to estimates from NLMMs, but also necessary if one wishes to compute indirect (i.e. mediated) effects using results from NLMMs only. The difference in scale for models with differing sets of predictors is due to the variance of the residual error term being fixed. This is not an issue in linear models where this variance is estimated rather than fixed [16-18]. In practice, the unknown parameters necessary for scaling, i.e. θ (the variance of ϵj and λ (the effect of the mediator on the outcome), can be estimated by fitting the following additional models:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M25">View MathML</a>

(12)

<a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M26">View MathML</a>

(13)

The first model (Equation 12) is similar to part of the structural model in the SEM (Equation 5), but instead is fit as a general linear model for longitudinal data with Mj as the outcome and zj as the predictor that allows correlation between the repeated observations. Using data from the alcohol and HIV example described earlier, this would be a model with longitudinal ART adherence measures as the outcome and measures of heavy alcohol consumption as time-varying predictors. The second model is an NLMM modeling Yj as a function of a random intercept and slope, and fixed effects for the continuous covariate (w), time (tj), the repeated binary predictor (zj), and the repeated mediator (Mj). For the HIV example, this would be a model with low CD4 count as the outcome and include both heavy alcohol consumption and ART adherence as predictors.

Using Equations 12 and 13, the estimated variance of ϵj provides an estimate of θ and the estimated coefficient β4 associated with Mj provides an estimate of λ. The original, unscaled NLMM estimate of the total effect of the primary predictor, β3 from Equation 8), can then be rescaled by multiplying by the factor: <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M27">View MathML</a> for a probit model (or <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M28">View MathML</a> for a logit model).

An alternative approach to using scaled regression coefficients would be to model indirect and direct pathways separately and obtain the total effect by summing the indirect and direct effects. That is, estimates of the coefficients β3 and β4 associated with zj and Mj, respectively (from Equation 13) can be used along with estimates of the coefficient α1 associated with zj (from Equation 12) to obtain the estimated total effect of the main independent variable <a onClick="popup('http://www.biomedcentral.com/1471-2288/12/5/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/12/5/mathml/M29">View MathML</a>.

In the current simulation study, we calculated the total effects of the main independent variable using both of the approaches described above. For the probit model, both methods yielded results comparable to those where the true parameters were known. The rescaled NLMM resulted in a parameter estimate (standard error) of 0.407 (0.105) which is a bias of 1.7%. The product of coefficients method yielded an estimate of 0.400 (0.102), which is a bias of 1.7%. Both estimates were very similar to those obtained using true values for the scaling factor, parameter estimate of 0.407 (0.096) and bias of 1.7%.

Acknowledgements

The authors wish to thank Howard Cabral and Timothy Heeren for thoughtful discussions of the problem. We are also grateful to Jeffrey H. Samet, PI of the HIV-LIVE Study (NIAAA: R01-AA13216), for use of the data.

References

  1. Muthen BO: A General Structural Equation Model With Dichotomous, Ordered Categorical, and Continuous Latent Variable Indicators.

    Psychometrika 1984, 49:115-132. Publisher Full Text OpenURL

  2. Davidian M, Giltinan D: Nonlinear Models for Repeated Measurement Data. New York, New York: Chapman and Hall; 1995. OpenURL

  3. Diggle PJ, Liang KY, Zeger SL: Analysis of Longitudinal Data. New York, New York: Oxford University Press; 1994. OpenURL

  4. McArdle JJ, Epstein D: Latent growth curves within developmental structural equation models.

    Child Development 1987, 58:110-133. PubMed Abstract | Publisher Full Text OpenURL

  5. Bollen KA: Structural Equations With Latent Variables. New York, New York: Wiley; 1989. OpenURL

  6. McArdle JJ, Aber MS: Statistical Methods in Longitudinal Research. New York: Academic Press; 1990:151-224.

    chap. Patterns of change within latent variable structural equations models

    OpenURL

  7. Muthen B: Best Methods for the Analysis of Change: Recent Advances, Unanswered Questions, Future Directions. Washington, D.C.: APA; 1991:1-17.

    chap. Analysis of Longitudinal Data Using Latent Variable Models with Varying Parameters

    OpenURL

  8. Muthen BO, Curran PJ: General Longitudinal Modeling of Individual Differences in Experimental Designs: A Latent Variable Framework for Analysis and Power Estimation.

    Psychological Methods 1997, 2(4):371-402. OpenURL

  9. Rovine MJ, Molenaar PCM: A Structural Modeling Approach to a Multilevel Random Coefficients Model.

    Multivariate Behavioral Research 2000, 35:51-88. Publisher Full Text OpenURL

  10. Rovine MJ, Molenaar PCM: New methods for the analysis of change. Washington, D.C.: APA; 2001:67-96.

    chap. A Structural Equations Modeling Approach to the General Linear Mixed Model

    OpenURL

  11. Stoel RD, van Den Wittenboer G, Hox J: Estimating Multilevel Linear Models as Structural Equation Models.

    Metodologia de las Ciencias del Comportamiento 2003, 5:21-42. OpenURL

  12. Curran PJ: Have Multilevel Models Been Structural Equation Models All Along?

    Multivariate Behavioral Research 2003, 38(4):529-569. Publisher Full Text OpenURL

  13. Bauer DJ: Estimating Multilevel Linear Models as Structural Equation Models.

    Journal of Educational and Behavioral Statistics 2003, 28(2):135-167. Publisher Full Text OpenURL

  14. Blood EA, Cabral H, Heeren T, Cheng DM: Performance of mixed effects models in the analysis of mediated longitudinal data.

    BMC Medical Research Methodology 2010., 16(10) OpenURL

  15. Blood EA, Cheng DM: The Use of Mixed Models For the Analysis of Mediated Data with Time-Dependent Predictors.

    Journal of Environmental and Public Health 2011, 2011. OpenURL

  16. MacKinnon D, Lockwood C, Brown C, Wang W, Hoffman J: The intermediate endpoint effect in logistic and probit regression.

    Clinical Trials 2007, 4:499-513. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Palta M, Lin CY: Latent Variables, Measurement Error and Methods for Analysing Longitudinal Binary and Ordinal Data.

    Statistics in Medicine 1999, 18:385-396. PubMed Abstract | Publisher Full Text OpenURL

  18. Fitzmaurice GM, Laird NM, Ware JH: Applied Longitudinal Analysis. Hoboken, New Jersey: John Wiley and Sons, Inc; 2004. OpenURL

  19. Muthen L, Muthen B: Mplus User's Guide. Fourth edition. Los Angeles, CA: Muthen & Muthen; 1998.

    chap. 3

    PubMed Abstract | Publisher Full Text OpenURL

  20. Rosner B: Fundamentals of Biostatistics. Pacific Grove, California: Dubxbury Press; 2005. OpenURL

  21. SAS Institute Inc: SAS 9.2 Help and Documentation. Cary, NC: SAS Institute Inc; 2002. OpenURL

  22. Cohen J: Statistical Power Analysis for the Behavioral Sciences. Hillsdale, New Jersey: Lawrence Erlbaum Associates; 1988. OpenURL

  23. Samet J, Cheng D, Libman H, Nunes D, Alperen J, Saitz R: Alcohol consumption and HIV disease progression.

    Journal of Acquired Immune Deficiency Syndrome 2007, 46(2):194-199. Publisher Full Text OpenURL

  24. MacKinnon DP, Warsi G, Dwyer JH: A Simulation Study of Mediated Effect Measures.

    Multivariate Behavioral Research 1995, 30:41-62. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2288/12/5/prepub