1. Repeated Measures from Daily Diary Studies

Data from daily diary studies often display fluctuations in scores within individuals. An example is a set of daily positive affect scores for 9 individuals shown in Fig. 1 that come from the Midlife in the United States (MIDUS 2): Daily Stress Project (Ryff and Almeida Reference Ryff and Almeida2004-2009). The data are part of a larger longitudinal study that began in 1995 (MIDUS 1) and included a nationally representative random-digit-dial sample of noninstitutionalized, English-speaking adults residing in the contiguous USA. Participants of MIDUS 1 were invited to participate in a second longitudinal wave (MIDUS 2) about 10 years later, along with a newly added adult cohort. From MIDUS 2, a subset of participants were selected at random to participate in the Daily Stress Project. Participants were interviewed by telephone to obtain daily self-reports for 8 consecutive days. Multiple scales were administered daily to measure different psychological constructs. Here, responses for 435 women with complete data on 5 survey items designed to reflect daily positive affect are considered. For instance, the first item asked "Did you feel in good spirits?" Items were measured using a 5-point ordinal response scale: 0=none of the time to 4=all of the time. In Fig. 1 are daily scores for 9 women, with scores computed as the daily mean of the 5 scale items. Calculating an average or summed score is a common way to handle responses obtained for a set of items designed to measure a common construct (e.g., Hedeker et al. Reference Hedeker, Mermelstein and Demirtas2012; Rast et al. Reference Rast, Hofer and Sparks2012). Although a measurement model for these items is provided later in Section 4, it is useful to first consider the usual approach to handling responses to multiple items. As shown in Fig. 1, the solid horizontal line that is common to all panels of the display is the mean score of 2.8 calculated across days and women. The dashed horizontal line unique to each plot is the individual’s response averaged across days. From the display, it can be seen that within individuals, scores fluctuate from day to day relative to the individual’s mean, and the extent of the daily variation within person varies between people. The mean responses of individuals differ from the overall mean and from each other. It is also interesting to note that for some individuals, scores reach the maximum possible score of 4 (Case 10186) while also displaying vary little variation across days. Although standard mixed-effects model might be effective in characterizing the differences between individuals with regard to their mean affect levels, the model would not address the between-subject heterogeneity of the intraindividual variation. Indeed, it is the between-subject differences in the within-subject variation that is a notable feature of the affect scores. Thus, a mixed-effects location scale model that permits the within-subject variance to vary by subject would seem to be quite suitable in characterizing the scores. This, in particular, would help to address possible ceiling effects in measurement (Hedeker et al. Reference Hedeker, Mermelstein and Demirtas2012).

Figure 1.

Daily positive affect scores \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$(y_i)$$\end{document} and daily positive affect scores averaged across days \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$(ybar_i)$$\end{document} for 9 study participants, plotted with the mean score (2.8) across days and all participants.

Common to many studies of a psychological variable is an implicit assumption that the variable is composed of trait- and state-level components (Nesselroade Reference Nesselroade, Baltes, Featherman and Lerner1987). For example, latent state trait theory suggests that a psychological measure is comprised of three independent sources: a latent trait, a latent state residual and measurement error (Steyer et al. Reference Steyer, Ferring and Schmitt1992). The latent trait component reflects the part of the score that is not specific to the measurement occasion, but rather, is unique to the person. Conversely, the latent state residual component is specific to the measurement occasion and the individual and, as such, is expected to fluctuate about the trait-level component according to the measurement occasions. Importantly, both the latent trait and state components are assumed to be free of measurement error. To permit study of the two latent components of a psychological variable, a second-order model can be applied to account for measurement error in the observed measures (Steyer et al. Reference Steyer, Ferring and Schmitt1992). An interesting application of a mixed-effects location scale model is to problems such as the study of daily positive affect (Hedeker et al. Reference Hedeker, Mermelstein and Demirtas2012; Rast et al. Reference Rast, Hofer and Sparks2012). The model provides a framework for separating the trait and state components of a measured variable (albeit assumed to be measured without error) while allowing the state residual to be truly unique to the individual by inclusion of a random subject effect in a model for the occasion-level residual variance (i.e., state residual variance). However, the standard formulation of the model does not address measurement error. To address the issue of measurement error in applications of latent state trait theory, a common factor model is applied to repeated measures data. The variance of a general common factor is assumed to be due to the latent trait factor and the variance of the occasion-specific factor is assumed to be due to the latent state residual factor. Thus, a potentially useful model is a latent variable mixed-effects location scale model that would permit partitioning of the observed item variance into the distinct components, and additionally, provide a truly person- and occasion-specific measure of the state residual variance.

2. The Mixed-Effects Location Scale Model

A two-level linear mixed-effects location scale model is described here (for a 3-level model, see Lin et al. Reference Lin, Mermelstein and Hedeker2018; for a generalized linear version of the model, see Hedeker et al. Reference Hedeker, Mermelstein, Demirtas and Berbaum2016). Let \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$y_{t i}$$\end{document} be an observed response at time \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$t =1 ,\ldots ,n_{i}$$\end{document} , for individual i, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$i =1 ,\ldots ,N$$\end{document} , where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$n_{i}$$\end{document} denotes that the number of observations can vary by the individual. Missing data are assumed to be missing at random. A linear mixed-effects location scale model for \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$y_{t i}$$\end{document} can be given as (cf: Hedeker et al. Reference Hedeker, Mermelstein and Demirtas2008):

(1) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} y_{t i} =\mathbf {x}_{t i}^{ \prime } (\mathbf {\beta }_{x} +\mathbf {\upsilon }_{i}) +\mathbf {w}_{i}^{ \prime } \mathbf {\beta }_{w} +\varepsilon _{t i} , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {x}_{t i}$$\end{document} is a p x 1 vector of covariates that vary according to the measurement occasion t and individual i and would usually include a "1" for the intercept of model, and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} is a q x 1 vector of covariates that vary according to the individual i. The coefficient \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\beta }_{x}$$\end{document} is a set of p fixed regression coefficients linking \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {x}_{t i}$$\end{document} to \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$y_{t i}$$\end{document} , and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\upsilon }_{i}$$\end{document} is a set of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p^{ *} \le p$$\end{document} random subject effects that correspond to one or more of the fixed effects in \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\beta }_{x}$$\end{document} . The coefficient \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\beta }_{w}$$\end{document} is a set of q fixed regression coefficients linking \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} to \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$y_{t i}$$\end{document} . The random effect \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\upsilon }_{i}$$\end{document} is assumed to be independent and normally distributed across individuals as

(2) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\upsilon }_{i} \sim N \left( \mathbf {0} ,\mathbf {\Phi }_{\upsilon } \left( \mathbf {w}_{i} ,\mathbf {\alpha } ,\mathbf {\rho }_{\upsilon }\right) \right) , \end{aligned}$$\end{document}

where the between-subject covariance matrix \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Phi }_{\upsilon }$$\end{document} may be a function of between-subject covariates ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} ) and corresponding coefficient vector \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\alpha }$$\end{document} , with \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\alpha } =(\alpha _{0} ,\ldots ,\alpha _{q}$$\end{document} ) denoting the set of coefficients that link the covariates to elements of the matrix \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Phi }_{\upsilon }$$\end{document} . The set of coefficients \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\rho }_{\upsilon }$$\end{document} \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$ =\left( \rho _{1} ,\ldots \right) $$\end{document} are the correlations between the random effects conditional on \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} . For instance, a common application of a mixed-effects location scale model is one in which only the intercept is random (e.g., Hedeker et al. Reference Hedeker, Mermelstein and Demirtas2008), and so \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Phi }_{\upsilon } =\phi _{\upsilon _{i}}^{2}$$\end{document} with the variance of the random intercept modeled as

(3) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \phi _{\upsilon _{i}}^{2} =\exp \left\{ \alpha _{0} +\alpha _{1} w_{1 i} +\ldots +\alpha _{q} w_{q i}\right\} , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\alpha _{0}$$\end{document} , when exponentiated, is the variance of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\upsilon _{i}$$\end{document} conditional on \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} . The remaining coefficients are the between-subject covariate effects on the between-subject variance. A positive effect of a coefficient indicates that a unit increase in a covariate corresponds to an increase in the variance of the random effect; a negative effect indicates that a unit increase in a covariate corresponds to a decrease in the variance. As will be shown later using an empirical example, allowing a between-subject variance to be a function of covariates makes it possible address between-subject heterogeneity of variance in a random effect.

At the occasion level, the residual \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon _{t i}$$\end{document} in (1) is assumed to be independent and normally distributed as:

(4) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \varepsilon _{t i} \sim N \left( 0 ,\sigma _{\varepsilon }^{2} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}, \mathbf {\tau }\text {,}a_{i}\right) \right) , \end{aligned}$$\end{document}

where the variance \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\sigma _{\varepsilon }^{2}$$\end{document} may be a function of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {x}_{t i}$$\end{document} (that would usually include a "1" for the intercept) and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} , with corresponding coefficient vector \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\tau } =(\mathbf {\tau }_{x} ,\mathbf {\tau }_{w})$$\end{document} , where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\tau }_{x} =(\mathbf {\tau }_{x 1} ,\ldots ,\mathbf {\tau }_{x p})$$\end{document} and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\tau }_{w} =(\tau _{w 1} ,\ldots ,\tau _{w q}$$\end{document} ) denoting sets of coefficients that link the within-subject and between-subject covariates, respectively, to the within-subject variance \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\sigma _{\varepsilon }^{2}$$\end{document} . In (4), the coefficient \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} denotes a between-subject random effect. Specifically, the variance of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon _{t i}$$\end{document} is modeled as

(5) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \sigma _{\varepsilon }^{2} =\exp \left\{ \tau _{0} +\tau _{x 1} +\ldots +\tau _{x p} +\ldots +\tau _{w 1} +\ldots +\tau _{w q} +a_{i}\right\} , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\tau _{0}$$\end{document} , when exponentiated, is the residual variance conditional on the covariates and for a subject for whom \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i} =0$$\end{document} . Similar to the interpretation of the covariate effects in (2) that serve to account for between-subject heterogeneity of variance in a random effect, a positive value for the effect of a covariate in the model for the within-subject variance indicates that a unit increase in a covariate corresponds to an increase in the residual variance, and a negative value indicates that a unit increase in a covariate corresponds to a decrease in the variance. As will be shown later using an empirical example, allowing the within-subject variance to be a function of covariates makes it possible address heterogeneity of the individual and occasion-specific residual variance.

The model for the variance of the occasion-specific residual \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon _{t i}$$\end{document} in (4) includes \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} to address heterogeneity of variance due to unobserved effects. That is, the variance of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon _{t i}$$\end{document} may differ between individuals, even after conditioning on covariates, and so including the random effect \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} allows for heterogeneity that may remain. This may be important, for instance, in problems for which there are between-subject differences in the variability in scores across occasions, such as if a ceiling effect in the measurement of the response is present and this results in reduced variation in scores for individuals whose average scores are at the extremes of the response scale. The random effect is assumed to be independent and normally distributed as

(6) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} a_{i} \sim N \left( 0 ,\phi _{a}^{2}\right) . \end{aligned}$$\end{document}

At the subject level of the model, the random effect \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\upsilon }_{i}$$\end{document} in (1) and the random effect \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} in (4) may covary. To allow this, the random effects are stacked, with the set defined here as \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {c}_{i} \left( r\text {x}1\right) =\left( \mathbf {\upsilon }_{i}^{ \prime } ,a_{i}\right) ^{ \prime }$$\end{document} to include all \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$r =p^{ *} +1$$\end{document} random effects. It is assumed that the joint distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {c}_{i}$$\end{document} is

(7) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {c}_{i} \sim N \left( \mathbf {0} ,\mathbf {\Phi }_{c}\right) , \end{aligned}$$\end{document}

where the covariance matrix is

(8) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\Phi }_{c} =\left( \begin{array}{cc}\mathbf {\Phi }_{\upsilon } \left( \mathbf {w}_{i} ,\mathbf {\alpha } ,\mathbf {\rho }_{\upsilon }\right) &{} \mathbf {\phi }_{\upsilon a}^{{}} \\ \mathbf {\phi }_{\upsilon a}^{ \prime } &{} \phi _{a}^{2}\end{array}\right) . \end{aligned}$$\end{document}

It follows that the probability distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {c}_{i}$$\end{document} is

(9) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} h_{c} \left( \mathbf {c}_{i}\right) =2 \pi ^{r/2} \left| \mathbf {\Phi }_{c}\right| ^{ -1/2} \exp \left( \left( \mathbf {c}_{i} -\mathbf {\mu }_{c}\right) ^{ \prime } \mathbf {\Phi }_{c}^{ -1} \left( \mathbf {c}_{i} -\mathbf {\mu }_{c}\right) \right) . \end{aligned}$$\end{document}

Letting \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i} =\left( y_{1 i} ,\ldots ,y_{n_{i}}\right) ^{ \prime }$$\end{document} and conditional on \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {c}_{i}$$\end{document} , the distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} is

(10) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {y}_{i}\vert \mathbf {c}_{i} \sim N \left( \mathbf {\mu }_{y\vert c} ,\mathbf {\Sigma }_{y\vert c}\right) , \end{aligned}$$\end{document}

where

(11) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\mu }_{y\vert c} =\mathbf {\Gamma }_{i} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}\right) \mathbf {\beta } \end{aligned}$$\end{document}

and

(12) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\Sigma }_{y\vert c} =\mathbf {\Gamma }_{i} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}\right) \mathbf {\Phi }_{\upsilon } \left( \mathbf {w}_{i} ,\mathbf {\alpha } ,\mathbf {\rho }_{\upsilon }\right) \mathbf {\Gamma }_{i} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}\right) ^{ \prime } +\mathbf {\Theta }_{\varepsilon } \left( \mathbf {x}_{t i} ,\mathbf {w}_{i} , \mathbf {\tau },\,a_{i}\right) , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Gamma }_{i} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}\right) $$\end{document} is a design matrix containing the within- and between-subject covariates, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\beta } =\left( \mathbf {\beta }_{x}^{ \prime } ,\mathbf {\beta }_{w}^{ \prime }\right) ^{ \prime }$$\end{document} is a stacked set of fixed coefficients linking the covariates to \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} , and assuming independence of the occasion-specific residuals between occasions and subjects, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Theta }_{\varepsilon } \left( \mathbf {x}_{t i} ,\mathbf {w}_{i} ,\mathbf {\tau }\text {,}a_{i}\right) =\mathbf {I}_{n_{i}} \sigma _{t i}^{2} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i} ,\mathbf {\tau }\text {,}a_{i}\right) $$\end{document} , with \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {I}_{n_{i}}$$\end{document} being an identity matrix of order \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$n_{i}$$\end{document} . Then, the conditional probability density function can be defined as

(13) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} h_{y\vert c} \left( \mathbf {y}_{i}\vert \mathbf {c}_{i}\right) =2 \pi ^{ -p/2} \left| \mathbf {\Sigma }_{y\vert c}\right| ^{ -1/2} \exp \left( \left( \mathbf {y}_{i} -\mathbf {\mu }_{y\vert c}\right) ^{ \prime } \mathbf {\Sigma }_{y\vert c}^{ -1} \left( \mathbf {y}_{i} -\mathbf {\mu }_{y\vert c}\right) \right) , \end{aligned}$$\end{document}

and the marginal distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} is then

(14) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} g \left( \mathbf {y}_{i}\right) ={\displaystyle \int _{\mathbf {c}_{i}}}h_{y} \left( \mathbf {y}_{i}\vert \mathbf {c}_{i}\right) h_{c} \left( \mathbf {c}_{i}\right) d \mathbf {c}_{i} . \end{aligned}$$\end{document}

The parameters of the model are \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\gamma } =\left( \mathbf {\beta } ,\mathbf {\alpha } ,\mathbf {\tau } ,\mathbf {\phi }_{\upsilon } ,\mathbf {\phi }_{\upsilon a} ,\phi _{a}^{2}\right) $$\end{document} . Given a random sample of observations, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{1} ,\ldots ,\mathbf {y}_{N}$$\end{document} , the marginal log-likelihood of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\gamma }$$\end{document} is

(15) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \text {ln}L \left( \mathbf {\gamma }\right) ={\displaystyle \sum _{i =1}^{N}}\ln g \left( \mathbf {y}_{i}\right) . \end{aligned}$$\end{document}

The model for \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} is nonlinear in the scale random effect, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} , and consequently, the integral in (14) does not have a closed-form expression. Estimation of the integral generally would rely on an approximation to the integral (Pinheiro and Bates Reference Pinheiro and Bates1995), such as by implementing SAS \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$^{1}$$\end{document} PROC NLMIXED that permits random effects to be incorporated into the within-subject covariance matrix in (4) (see Hedeker et al. Reference Hedeker, Mermelstein and Demirtas2012).

2.1. Example: Manifest Measures of Daily Mean Positive Affect

The daily mean positive affect scores described earlier for the 435 women are analyzed using a mixed-effects location scale model (Hedeker et al. Reference Hedeker, Mermelstein and Demirtas2012; Rast et al. Reference Rast, Hofer and Sparks2012). Within individuals, daily measures of positive affect may not be expected to change systematically with the passage of days, but rather, vary across days. Additionally, scores averaged across days within individuals would be expected to naturally vary between individuals. Accordingly, the manifest measures of mean positive affect are assumed to follow a random intercept model:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} y_{t i} =\beta _{0 i} +\varepsilon _{t i}\text {,} \end{aligned}$$\end{document}

where, for participant i, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$y_{t i}$$\end{document} is the observed mean positive affect score on day t, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\beta _{0 i} =\beta _{0} +\upsilon _{i}$$\end{document} is the expected response across days, and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon _{t i}$$\end{document} is the discrepancy between the observed and fitted responses on day t. The random intercept \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\beta _{0 i}$$\end{document} represents the trait component that varies by individual, and the residual \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon _{t i}$$\end{document} represents the state component that varies by both individual and occasion. The random intercept is assumed to be normally distributed: \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\beta _{0 i} \sim N \left( \beta _{0} ,\phi _{\upsilon }^{2}\right) $$\end{document} , where

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned}\phi _{\upsilon }^{2} =\exp \left\{ \alpha _{0}\right\} \text {.} \end{aligned}$$\end{document}

This expression of the variance of the random intercept becomes useful later when covariates are added to the model to address between-subject heterogeneity in the random intercept variance. The residual is assumed to be normally distributed: \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon _{t i} \sim N \left( 0 ,\mathbf {I} \sigma _{\varepsilon }^{2}\right) $$\end{document} , where

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned}\sigma _{\varepsilon }^{2} =\exp \left\{ \tau _{0} +a_{i}\right\} \text {,} \end{aligned}$$\end{document}

and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} is a random subject effect that is assumed to be normally distributed: \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i} \sim N \left( 0 ,\phi _{a}^{2}\right) \text {.}$$\end{document} At the subject level, the random intercept and the random effect of the within-subject variance may covary: \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$c o v \left( \upsilon _{i} ,a_{i}\right) =\phi _{\upsilon a}^{{}}$$\end{document} . The model has five parameters: \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\gamma } =\left( \beta _{0} ,\alpha _{0} ,\tau _{0} ,\phi _{\upsilon }^{2} ,\phi _{\upsilon a}^{{}}\right) \text {.}$$\end{document}

To assess the need to include the random effect in the model for \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\sigma _{\varepsilon }^{2}$$\end{document} , estimates were obtained for a reduced form of the model, henceforth called Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{1}$$\end{document} , that assumed homogeneity of the residual variance, and so \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\phi _{\upsilon }^{2} =0$$\end{document} (1st set of estimates in Table 1 corresponding to Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{1}$$\end{document} ). Thus, Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{1}$$\end{document} is a random intercept model with a fixed intercept \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\beta _{0}$$\end{document} and includes parameter \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\alpha _{0}$$\end{document} of the variance model for the random intercept ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\phi _{\upsilon }^{2} =\exp \left\{ \alpha _{0}\right\} $$\end{document} ), and parameter \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\tau _{0}$$\end{document} of the variance model for the residual ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\sigma _{\varepsilon }^{2} =\exp \left\{ \tau _{0}\right\} $$\end{document} ). Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{1}$$\end{document} and the mixed-effects location scale model, henceforth called Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} , were fit using SAS PROC NLMIXED version 9.4 using nonadaptive Gaussian quadrature with 7 quadrature points. Increasing the number of quadrature points beyond 7 or using a Laplacian approximation (Vonesh Reference Vonesh1996) did not result in a converged solution for Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} . Maximum likelihood (ML) estimates of Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} are given in Table 1 (second set of estimates). Starting values for the estimation of these models involved fitting a fixed-effects model first, followed by adding a random intercept (for Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{1}$$\end{document} ), then a random effect in the within-subject variance model and the covariance between it and the random intercept (for Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} ). ML estimation carried out using SAS PROC NLMIXED with a Laplacian approximation). Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} that assumed \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\phi _{\upsilon }^{2} >0$$\end{document} yields an improvement in model fit in terms of the AIC, from 6, 087 (Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{1}$$\end{document} ) down to 5, 547 (Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} ) suggesting between-subject heterogeneity of the within-person variance. A likelihood ratio test (LRT) of the goodness-of-fit between the two models using 2 df is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\chi _{2}^{2} =6,081 -5,537 =544\text {,}$$\end{document} \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p <.001\text {,}$$\end{document} supporting the need to include the random effect in the occasion-level residual variance model. From Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} , the estimated common variance of the residual is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\exp ( -1.84) =0.16$$\end{document} , and the estimated standard deviation of the random effect in the model for \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\sigma _{\varepsilon }^{2}$$\end{document} is 1.16. These results suggest between-subject differences in the extent to which daily positive affect scores vary from day to day. The estimated correlation between the random intercept and the random subject effect for the within-subject variance is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$r = -.36$$\end{document} , suggesting a tendency for those with a relatively low day-to-day variation in positive affect scores to also have a relatively high positive affect score averaged across days. Thus, either those with relatively high levels of positive affect across days are also relatively stable in their affect, or the result is an indication of a ceiling effect in the measurement. Finally, from the second of the two models, the estimated mean positive affect score across days and participants is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\beta }_{0} =3.00$$\end{document} \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\left( s e =0.006\right) $$\end{document} , indicating a rather high average across days and participants.

Table 1.

ML Estimates of a mixed-effects location scale model of manifest measures of daily mean positive affect ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$n=435$$\end{document} ).

Notes: \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\ ^{a}$$\end{document} ML estimation carried out using SAS PROC NLMIXED with a Laplacian approximation to the integral. \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$^{b}$$\end{document} ML estimation carried out using SAS PROC NLMIXED with nonadaptive Gaussian quadrature with 7 quadrature points

2.1.1. Daily Stressors

Extensive research has demonstrated a clear relationship between negative affect and daily stress, but positive reactivity to daily stress is not as well understood. Several studies report mixed findings in terms of the direction of the relationship, suggesting that the relationship is complex and highly individualistic (Schilling and Diehl Reference Schilling and Diehl2014). Here, the mixed-effects location scale model is extended to include the number of daily stressors reported for each interview day, denoted by daily stressors \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{t i}$$\end{document} . Scores were based on a sum of responses to 7 questions scored as 1=yes and 0=no; the first in the set of questions, for instance, was “Did you have an argument or disagreement with anyone since (this time/we spoke) yesterday?” The possible minimum and maximum scores were 0 and 7, respectively. Given repeated measures of the number of stressors, both the within- and between-subject effects of daily stressors on daily mean positive affect scores were studied. To do this, the variable was centered about the individual’s mean number of stressors across days: DS \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{t i}^{ *}$$\end{document} = daily stressors \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{t i}- \textit{DS}_{ .i}$$\end{document} ), where DS \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{t i}^{ *}$$\end{document} denotes an individual’s daily stress count centered about their mean count across days, with the latter denoted by DS \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{ .i}$$\end{document} . Daily stressor counts for the sample ranged from 0 to 5; individual-level mean counts ranged from 0 to 2.38, with an overall mean of 0.53 and standard deviation of 0.42. The person-mean centered and mean values of the daily stressor measures were included in the model to predict the observed daily mean positive affect score:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned}y_{t i} =\beta _{0 i} +\beta _{1} \text {\textit{DS}}_{t i}^{ *} +\beta _{2} \text {\textit{DS}}_{ .i} +\varepsilon _{t i}\text {.} \end{aligned}$$\end{document}

To evaluate the effect of daily stressors on the variation of positive affect scores within and across days, after adjusting for the within- and between-person effects of daily stressors, both measures of daily stressors were included in the model of the occasion-level residual variance:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned}\sigma _{\varepsilon }^{2} =\exp \left\{ \tau _{0} +\tau _{1} \text {\textit{DS}}_{t i}^{ *} +\tau _{2} \text {\textit{DS}}_{ .i} +a_{i}\right\} \text {.} \end{aligned}$$\end{document}

Thus, this model was used to evaluate whether higher daily counts of stressors or an overall high level of stressors across days was related to the daily variation in mean positive affect scores across days. The variance of the random intercept was also evaluated for between-subject heterogeneity by including the mean daily stressor count in the model of the random intercept variance:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \phi _{\upsilon }^{2} =\exp \left\{ \alpha _{0} +\alpha _{1} \text {\textit{DS}}_{ .i}\right\} \text {.} \end{aligned}$$\end{document}

Starting values for the effects of covariates were set to 0, and the remaining parameters were set equal to the values obtained from fitting the unconditional model described in the previous section. The model was fit using SAS PROC NLMIXED version 9.4 using nonadaptive Gaussian quadrature with 7 quadrature points. Estimates of the conditional model, referred to as Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{3}$$\end{document} , are in the last column of Table 1. The mean positive affect score for those with no daily stressors is estimated to be 3.05 ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =.012$$\end{document} ), a value that is slightly greater than that obtained under Model A \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} . The estimates further suggest an overall negative within-person effect ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\beta }_{1} = -0.075$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.012$$\end{document} ), as well as a negative between-person effect ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\beta }_{2} = -0.36$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.030$$\end{document} ), of the number of daily stressors on daily mean positive affect, with the magnitude of the between-person effect being nearly five times the size of the within-person effect (H \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{0} :\beta _{1} =\beta _{2}$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\chi ^{2} \left( 1\text {df}\right) =164$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p <.001$$\end{document} ). The within-person effect indicates that, on average, a higher number of daily stressors correspond to a relatively low daily mean positive affect score, and the between-person effect indicates that a higher mean number of daily stressors corresponds to a relatively low daily mean positive affect score. Conditional on these effects, greater daily variation in affect scores corresponds to a higher number of daily stressors ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\tau }_{1} =0.30$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.049$$\end{document} ), as well as a higher mean number of daily stressors ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\tau }_{2} =0.78$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.15$$\end{document} ), with the latter of the two effects being the largest (H \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{0} :\tau _{1} =\tau _{2}$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\chi ^{2} \left( 1\text {df}\right) =7.5$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p =.006$$\end{document} ). Additionally, the variance of the random intercept is positively related to the mean number of daily stressors ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\alpha }_{1} =0.36$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.070$$\end{document} ), such that the between-subject variation in mean positive affect scores increases with an increase in the mean number of daily stressors. That is, individuals increasingly differ from each other with an increase in the mean number of daily stressors. Overall, it seems that between-subject differences, relative to within-subject variation, in the number of stressors relates most strongly to the different aspects of daily positive affect.

3. The Latent Variable Mixed-Effects Location Scale Model

The mixed-effects location scale model is a two-level model for observations nested within individuals (Hedeker et al. Reference Hedeker, Mermelstein and Demirtas2008). The responses are assumed to be free of measurement error. To address measurement error, a second-order model is developed here. In a second-order version of the model, the repeated measure is a latent variable. That is, at each occasion, a set of manifest variables is assumed to reflect a latent variable that varies by time and person. The subject-specific model is then applied to the repeated measures of the latent variable. Let \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{t i}$$\end{document} be a set of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$m_{t i}$$\end{document} indicators of the latent variable \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\eta _{t i}$$\end{document} observed at occasion t for individual i. Across \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$T_{i}$$\end{document} occasions, let \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i} =\left( \mathbf {y}_{1 i} ,\ldots ,\mathbf {y}_{T_{i}}\right) ^{ \prime }$$\end{document} and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$M_{i}$$\end{document} denote the number of observations in \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} across indicators and occasions. It is not necessary that each measurement occasion have the same set of indicator variables to represent the latent variable, as long as all share a common indicator that anchors the latent variable over time (Bollen and Curran Reference Bollen and Curran2006), and data across occasions need not be complete. Missing data are assumed to be missing at random. A longitudinal measurement model for \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} is

(16) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {y}_{i} =\mathbf {\tau }_{i} +\mathbf {\Lambda }_{i} (\mathbf {\lambda }) \mathbf {\eta }_{i} +\mathbf {\delta }_{i} , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\tau }_{i}$$\end{document} is an \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$M_{i}$$\end{document} x 1 intercept vector, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Lambda }_{i} (\mathbf {\lambda })$$\end{document} is an \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$M_{i}$$\end{document} x Q matrix of factor loadings, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\eta }_{i}$$\end{document} is a Q x 1 latent variable vector observed according to \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$T_{i}$$\end{document} , and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\delta }_{i}$$\end{document} is an \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$M_{i}$$\end{document} x 1 unique item factor vector. The subscript i used for both \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\tau }_{i}$$\end{document} and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Lambda }_{i}$$\end{document} indicates that both may vary by i with regard to their dimensions but not otherwise, whereas both the values and dimensions of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\eta }_{i}$$\end{document} and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\delta }_{i}$$\end{document} may vary by i. The latent variable \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\eta _{t i}$$\end{document} is assumed to follow a mixed-effects location scale model:

(17) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \eta _{t i} =\mathbf {x}_{t i}^{ \prime } (\mathbf {\beta }_{x} +\mathbf {\upsilon }_{i}) +\mathbf {w}_{i}^{ \prime } \mathbf {\beta }_{w} +\varepsilon _{t i} , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {x}_{t i}$$\end{document} is a p x 1 vector of within-person covariates that would usually include a "1" for the intercept of the model, and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} is a q x 1 vector of between-subject covariates. The coefficient \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\beta }_{x}$$\end{document} contains p fixed coefficients linking \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {x}_{t i}$$\end{document} to \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\eta _{t i}$$\end{document} ; \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\upsilon }_{i}$$\end{document} is a set of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p^{ *} \le p$$\end{document} random subject effects that correspond to one or more of the fixed effects in \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\beta }_{x}$$\end{document} . The coefficient \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\beta }_{w}$$\end{document} contains q fixed coefficients linking \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} to \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\eta _{t i}$$\end{document} . The random effect \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\upsilon }_{i}$$\end{document} is assumed to be independent and normally distributed across individuals as

(18) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\upsilon }_{i} \sim N \left( \mathbf {0} ,\mathbf {\Phi }_{\upsilon } \left( \mathbf {w}_{i} ,\mathbf {\alpha } ,\mathbf {\rho }_{\upsilon }\right) \right) , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Phi }_{\upsilon }$$\end{document} is a function of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\alpha } =\left( \alpha _{0} ,\ldots ,\alpha _{q}\right) \text {,}$$\end{document} and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\rho }_{\upsilon } =\left( \rho _{1} ,\ldots \right) $$\end{document} . The residual \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\varepsilon }_{i}$$\end{document} , specific to the individual with elements specific to the occasions, is assumed to be normally distributed:

(19) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\varepsilon }_{i} \sim N \left( \mathbf {0} ,\mathbf {\Theta }_{\varepsilon } \left( \mathbf {x}_{t i} ,\mathbf {w}_{i} ,\mathbf {\tau }\text {,}a_{i}\right) \right) , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Theta }_{\varepsilon }$$\end{document} is a function of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {x}_{t i}$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\tau }$$\end{document} , and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} , with \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\tau } =\left( \mathbf {\tau }_{x} ,\mathbf {\tau }_{w}\right) $$\end{document} similar to (4) in which the coefficients are the effects of within- and between-subject covariates, respectively. Assuming the residuals are independent between occasions, then \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Theta }_{\varepsilon } =\sigma _{\varepsilon }^{2} \mathbf {I}_{n_{i}}$$\end{document} , where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\sigma _{t i}^{2}$$\end{document} is a function of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {x}_{t i}$$\end{document} and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {w}_{i}$$\end{document} with random effect \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} , and where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} is assumed to be normally distributed: \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i} \sim N \left( 0 ,\phi _{a}^{2}\right) \text {.}$$\end{document} Thus, the model for the observed response \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} is

(20) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {y}_{i} =\mathbf {\tau }_{i} +\mathbf {\Lambda }_{i} (\mathbf {\lambda }) \left( \mathbf {\Gamma }_{i} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}\right) \mathbf {\beta }_{i} +\mathbf {\varepsilon }_{i}\right) +\mathbf {\delta }_{i} . \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Gamma }_{i} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}\right) $$\end{document} is a design matrix containing the within- and between-subject covariates, and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\beta }_{i} =\left( \left( \mathbf {\beta }_{x} +\mathbf {\upsilon }_{i}\right) ^{ \prime } ,\mathbf {\beta }_{w}^{ \prime }\right) ^{ \prime }$$\end{document} is a stacked set of coefficients linking the within- and between-subject covariates, respectively, to \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} . The residual \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\delta }_{i}$$\end{document} is assumed to be normally distributed:

(21) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\delta }_{i} \sim N \left( \mathbf {0} ,\mathbf {\Psi }_{i}\right) , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Psi }_{i}$$\end{document} is a diagonal covariance matrix of the unique item factors.

Similar to the model for manifest variables, the random effect vector \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\upsilon }_{i}$$\end{document} of (18) and the random effect \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} in (19) may covary. To allow this, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\upsilon }_{i}$$\end{document} and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} are stacked in a vector defined by \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {c}_{i} \left( r\text {x}1\right) =\left( \mathbf {\upsilon }_{i}^{ \prime } ,a_{i}\right) ^{ \prime }$$\end{document} that includes all \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$r =p^{ *} +1$$\end{document} random effects. The joint distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {c}_{i}$$\end{document} is assumed to be

(22) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {c}_{i} \sim N \left( \mathbf {\mu }_{c} ,\mathbf {\Phi }_{c}\right) , \end{aligned}$$\end{document}

where the mean vector and covariance matrix are

(23) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\mu }_{c} =\left( \begin{array}{c}\mathbf {\beta }_{x} \\ 0\end{array}\right) \end{aligned}$$\end{document}

and

(24) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\Phi }_{c} =\left( \begin{array}{cc}\mathbf {\Phi }_{\upsilon } \left( \mathbf {w}_{i} ,\mathbf {\alpha } ,\mathbf {\rho }_{\upsilon }\right) &{} \mathbf {\phi }_{\upsilon a} \\ \mathbf {\phi }_{\upsilon a}^{ \prime } &{} \phi _{a}^{2}\end{array}\right) , \end{aligned}$$\end{document}

respectively. Corresponding to this, the probability distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {c}_{i}$$\end{document} is

(25) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} h_{c} \left( \mathbf {c}_{i}\right) =2 \pi ^{r/2} \left| \mathbf {\Phi }_{c}\right| ^{ -1/2} \exp \left( \left( \mathbf {c}_{i} -\mathbf {\mu }_{c}\right) ^{ \prime } \mathbf {\Phi }_{c}^{ -1} \left( \mathbf {c}_{i} -\mathbf {\mu }_{c}\right) \right) . \end{aligned}$$\end{document}

Conditional on the random effects in \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {c}_{i}$$\end{document} , the distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} is

(26) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {y}_{i}\vert \mathbf {c}_{i} \sim N \left( \mathbf {\mu }_{y\vert c} ,\mathbf {\Sigma }_{y\vert c}\right) , \end{aligned}$$\end{document}

where

(27) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\mu }_{y\vert c} =\mathbf {\tau }_{i} +\mathbf {\Lambda }_{i} \mathbf {\Gamma }_{i} \mathbf {\beta } \end{aligned}$$\end{document}

and

(28) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\Sigma }_{y\vert c} =\mathbf {\Lambda }_{i} \left( \mathbf {\Gamma }_{i} \mathbf {\Phi }_{\upsilon } \mathbf {\Gamma }_{i}^{ \prime } +\mathbf {\Theta }_{\mathbf {\varepsilon }}\right) \mathbf {\Lambda }_{i}^{ \prime } +\mathbf {\Psi }_{i} , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\beta } =\left( \mathbf {\beta }_{x}^{ \prime } ,\mathbf {\beta }_{w}^{ \prime }\right) ^{ \prime }$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Gamma }_{i} =\mathbf {\Gamma }_{i} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}\right) $$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Phi }_{\upsilon } =\mathbf {\Phi }_{\upsilon } \left( \mathbf {w}_{i} ,\mathbf {\alpha } ,\mathbf {\rho }_{\upsilon }\right) $$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Theta }_{\mathbf {\varepsilon }} =\mathbf {\Theta }_{\mathbf {\varepsilon }} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i} ,\mathbf {\tau } ,a_{i}\right) $$\end{document} , and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\Psi }_{i} =\mathbf {\Psi }_{i} \left( \mathbf {\psi }\right) $$\end{document} . The conditional probability density function is defined as

(29) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} h_{y\vert c} \left( \mathbf {y}_{i}\vert a_{i}\right) =2 \pi ^{ -p/2} \left| \mathbf {\Sigma }_{y\vert c}\right| ^{ -1/2} \exp \left( \left( \mathbf {y}_{i} -\mathbf {\mu }_{y\vert c}\right) ^{ \prime } \mathbf {\Sigma }_{y\vert c}^{ -1} \left( \mathbf {y}_{i} -\mathbf {\mu }_{y\vert c}\right) \right) . \end{aligned}$$\end{document}

The marginal distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} is then

(30) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} g \left( \mathbf {y}_{i}\right) ={\displaystyle \int _{\mathbf {c}_{i}}}h_{y} \left( \mathbf {y}_{i}\vert \mathbf {c}_{i}\right) h_{c} \left( \mathbf {c}_{i}\right) d \mathbf {c}_{i} . \end{aligned}$$\end{document}

The parameters of the model are \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\gamma } =\left( \mathbf {\beta } ,\mathbf {\alpha } ,\mathbf {\rho }_{\upsilon } ,\mathbf {\tau } ,\mathbf {\phi }_{\upsilon } ,\mathbf {\phi }_{\upsilon a}^{{}} ,\phi _{a}^{2} ,\mathbf {\psi }\right) $$\end{document} . Given a sample of observations, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{1} ,\ldots ,\mathbf {y}_{N}$$\end{document} , the marginal log-likelihood of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\gamma }$$\end{document} is

(31) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \text {ln}L \left( \mathbf {\gamma }\right) ={\displaystyle \sum _{i =1}^{N}}\ln g \left( \mathbf {y}_{i}\right) . \end{aligned}$$\end{document}

Due to the random effect \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} , the model for \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} is nonlinear in one of the random effects, and consequently, the integral in (30) does not have a closed-form expression. Steps taken to estimate the integral are described in the next section.

3.1. Simplifying the Likelihood Function

The latent variable mixed-effects location scale model is more complicated than the version of the model for manifest variables. This is due to the measurement model needed for the set of response variables at multiple occasions, in addition to the random effects (one of which enters the model in a nonlinear manner) that are used in the model for the repeated measures of the latent variable. Interestingly, it is possible to simplify the calculation of the \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$i^{t h}$$\end{document} individual’s contribution to the likelihood function by analytically eliminating the random effect \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\upsilon }_{i}$$\end{document} from the marginal distribution in (30). This can also be done for the marginal distribution in (14) based on the model for manifest variables (see Appendix A). Based on the fact that, conditional on \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} , the model is linear in \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\upsilon }_{i}$$\end{document} , only the conditional distribution given \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} is needed (cf: du Toit and Cudeck Reference du Toit and Cudeck2009). Because the random effects contained in \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {c}_{i}$$\end{document} are jointly normally distributed, the conditional distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\upsilon }_{i}$$\end{document} given \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} is also normal,

(32) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\upsilon }_{i}\vert a_{i} \sim N \left( \mathbf {\mu }_{\beta .a} ,\mathbf {\Phi }_{\beta .a}\right) , \end{aligned}$$\end{document}

with mean vector and covariance matrix (Morrison Reference Morrison1990, 3rd ed., p. 92)

(33) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\mu }_{\beta .a} =\mathbf {\beta } +\frac{\mathbf {\phi }_{\upsilon a}}{\phi _{a}^{2}} a_{i} \end{aligned}$$\end{document}

and

(34) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\Phi }_{\upsilon .a} =\mathbf {\Phi }_{\upsilon } \left( \mathbf {w}_{i} ,\mathbf {\alpha } ,\mathbf {\rho }_{\upsilon }\right) -\mathbf {\phi }_{\upsilon a} (\phi _{a}^{2})^{ -1} \mathbf {\phi }_{\upsilon a}^{ \prime } , \end{aligned}$$\end{document}

respectively. It follows that the distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} conditional on \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} is normal:

(35) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {y}_{i}\vert a_{i} \sim N \left( \mathbf {\mu }_{y .a} ,\mathbf {\Sigma }_{y .a}\right) , \end{aligned}$$\end{document}

with

(36) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\mu }_{y .a} =\mathbf {\tau }_{i} +\mathbf {\Lambda }_{i} \mathbf {\Gamma }_{i} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}\right) \mathbf {\mu }_{\beta .a} \end{aligned}$$\end{document}

and

(37) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\Sigma }_{y .a} =\mathbf {\Lambda }_{i} \left( \mathbf {\Gamma }_{i} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}\right) \mathbf {\Phi }_{\upsilon .a} \mathbf {\Gamma }_{i} \left( \mathbf {x}_{t i} ,\mathbf {w}_{i}\right) ^{ \prime } +\mathbf {\Theta }_{\varepsilon }^{a}\right) \mathbf {\Lambda }_{i}^{ \prime } +\mathbf {\Psi }_{i}\text {.} \end{aligned}$$\end{document}

Assuming independence between the residuals relating to the latent variable between occasions,

(38) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \mathbf {\Theta }_{\varepsilon }^{a} =\mathbf {I}_{n_{i}}\exp \left\{ \tau _{0} +\tau _{x 1} +\ldots +\tau _{x p} +\ldots +\tau _{w 1} +\ldots +\tau _{w q} +a_{i}\right\} . \end{aligned}$$\end{document}

The probability density function of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} given \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} is

(39) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} h_{y\vert a} \left( \mathbf {y}_{i}\vert a_{i}\right) =2 \pi ^{ -p/2} \left| \mathbf {\Sigma }_{y .a}\right| ^{ -1/2} \exp \left( \left( \mathbf {y}_{i} -\mathbf {\mu }_{y .a}\right) ^{ \prime } \mathbf {\Sigma }_{y .a}^{ -1} \left( \mathbf {y}_{i} -\mathbf {\mu }_{y .a}\right) \right) . \end{aligned}$$\end{document}

The marginal distribution of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} given \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} is then

(40) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} g \left( \mathbf {y}_{i}\right) ={\displaystyle \int _{a_{i}}}h_{y\vert a} \left( \mathbf {y}_{i}\vert a_{i}\right) h_{a} \left( a_{i}\right) d a_{i} , \end{aligned}$$\end{document}

thus requiring a one-dimensional integral. The log-likelihood function is then

(41) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \text {ln}L \left( \mathbf {\gamma }\right) ={\displaystyle \sum _{i =1}^{N}}\ln g \left( \mathbf {y}_{i}\right) , \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\gamma } =\left( \mathbf {\beta } ,\mathbf {\alpha } ,\mathbf {\rho }_{\upsilon } ,\mathbf {\tau } ,\mathbf {\phi }_{\upsilon } ,\mathbf {\phi }_{\upsilon a}^{{}} ,\phi _{a}^{2} ,\mathbf {\psi }\right) $$\end{document} . Again, as the model for \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {y}_{i}$$\end{document} is nonlinear in \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} , the integral in (40) does not have a closed-form expression, and so estimation relies on an approximation. In the empirical example presented in the next section, Gauss–Hermite quadrature is used to evaluate the integral. A computer program, such as one written with SAS IML, may be used to carry out estimation.

3.2. Example: Daily Latent Positive Affect

One shortcoming of the previous analysis of the daily mean positive affect scores was the fact that the measurement error of the manifest variables was not addressed. Here, a latent variable mixed-effects location scale model is applied to the item-level response data. Measurement invariance is generally required for a meaningful interpretation of a latent variable that is studied over repeated occasions because it is necessary that an instrument perform in the same manner regardless of the measurement occasion (Chan, 1998; Meredith Reference Meredith1993; Millsap Reference Millsap2011). Assuming a single dimension for the set of item responses, measurement invariance across the 8 days was evaluated. First, equality of factor loadings was tested, followed by similar tests about the item intercepts and then the uniqueness. Relative to a model without restrictions placed on the factor loadings, item intercepts or uniquenesses (Deviance = 31,009, BIC = 32,759), LRTs suggested equal factor loadings by item across days ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\chi ^{2} \left( 28 d f\right) =35.4\text {,}$$\end{document} \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p =.16$$\end{document} ; Deviance = 31,045, BIC = 32,625), equal item intercepts by item across days (with the exception of item 4 on day 8) ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\chi ^{2} \left( 27 d f\right) =36.8\text {,}$$\end{document} \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p =.10$$\end{document} , Deviance = 31,046, BIC = 32,632), and both equal factor loadings and item intercepts by item across days (with the exception of allowing for unique intercepts for item 10 on days 1 and 8 and item 11 day 8) ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\chi ^{2} \left( 53 d f\right) =62.47\text {,}$$\end{document} \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p =.18$$\end{document} , Deviance = 31,072, BIC = 32,500). Although the BIC was lowest for a model that additionally constrained the uniquenesses to be equal by item and across days (Deviance = 31,226, BIC = 32,441), a LRT did not support this ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\chi ^{2} \left( 35 d f\right) =154.03\text {,}$$\end{document} \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p <.001$$\end{document} ). The analysis proceeded assuming equal uniquenesses, along with equal factor loadings and item intercepts, by item across days, with the results evaluated for the tenability of this assumption.

A latent variable model based on the daily responses to the 5-item set assumed that the intercept of the first item equation was equal to 0 and that the factor loading was equal to 1:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned}\left[ \begin{array}{c}y_{1 1 i} \\ y_{2 1 i} \\ y_{3 1 i} \\ y_{4 1 i} \\ y_{5 1 i} \\ \vdots \\ y_{1 8 i} \\ y_{2 8 i} \\ y_{3 8 i} \\ y_{4 8 i} \\ y_{5 8 i}\end{array}\right] =\left[ \begin{array}{c}0 \\ \tau _{2} \\ \tau _{3} \\ \tau _{4} \\ \tau _{5} \\ \vdots \\ 0 \\ \tau _{2} \\ \tau _{3} \\ \tau _{4} \\ \tau _{5}\end{array}\right] +\left[ \begin{array}{cccc}1 &{} 0 &{} \cdots &{} 0 \\ \lambda _{2} &{} 0 &{} \cdots &{} 0 \\ \lambda _{3} &{} 0 &{} \cdots &{} 0 \\ \lambda _{4} &{} 0 &{} \cdots &{} 0 \\ \lambda _{5} &{} 0 &{} \cdots &{} 0 \\ \vdots &{} \vdots &{} \ddots &{} \vdots \\ 0 &{} 0 &{} \cdots &{} 1 \\ 0 &{} 0 &{} \cdots &{} \lambda _{2} \\ 0 &{} 0 &{} \cdots &{} \lambda _{3} \\ 0 &{} 0 &{} \cdots &{} \lambda _{4} \\ 0 &{} 0 &{} \cdots &{} \lambda _{5}\end{array}\right] \left[ \begin{array}{c}\eta _{1 i} \\ \vdots \\ \eta _{8 i}\end{array}\right] +\left[ \begin{array}{c}\delta _{1 1 i} \\ \delta _{2 1 i} \\ \delta _{3 1 i} \\ \delta _{4 1 i} \\ \delta _{5 1 i} \\ \vdots \\ \delta _{1 8 i} \\ \delta _{2 8 i} \\ \delta _{3 8 i} \\ \delta _{4 8 i} \\ \delta _{5 8 i}\end{array}\right] \text {.} \end{aligned}$$\end{document}

Elements of the unique factor score vector \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\delta }_{i} =\left( \delta _{1 1 i} ,\ldots ,\delta _{5 8 i}\right) $$\end{document} were assumed to be independently and normally distributed across individuals as \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\delta }_{i} \sim N \left( \mathbf {0} ,\mathbf {\Psi }\right) $$\end{document} , with scores assumed to be independent between days with constant variance across days. Let \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\tau } =\left( 0 ,\tau _{2} ,\tau _{3} ,\tau _{4} ,\tau _{5}\right) $$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\lambda } =\left( 1 ,\lambda _{2} ,\lambda _{3} ,\lambda _{4} ,\lambda _{5}\right) $$\end{document} , and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\psi } =\left( \psi _{1} ,\psi _{2} ,\psi _{3} ,\psi _{4} ,\psi _{5}\right) $$\end{document} , where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\psi }$$\end{document} is the vector of the unique factor score variances relating to the five scale items. To then fit a latent variable mixed-effects location scale model, the latent measure of positive affect was assumed to follow a mixed-effects location scale model. Similar to the daily scores based on averages of the item responses, the latent measure of daily positive affect was expected to fluctuate across days at the individual level and so was assumed to follow a random intercept model:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned}\eta _{t i} =\beta _{0 i} +\varepsilon _{t i}\text {,} \end{aligned}$$\end{document}

where, for participant i, \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\beta _{0 i}$$\end{document} is the latent measure of positive affect on day t, and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon _{t i}$$\end{document} is the discrepancy between the latent score and the fitted value on day t. The random intercept is assumed to be normally distributed:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned}\beta _{0 i} \sim N \left( \beta _{0} ,\phi _{\upsilon }^{2}\right) \text {,} \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\beta _{0}$$\end{document} is the expected latent response across days, and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\phi _{\upsilon }^{2} =\exp \left\{ \alpha _{0}\right\} \text {.}$$\end{document} At the occasion level, the residual is assumed to be normally distributed:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \varepsilon _{t i} \sim N \left( 0 ,\mathbf {I} \sigma _{\varepsilon }^{2}\right) \text {,} \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\sigma _{\varepsilon }^{2} =\exp \left\{ \tau _{0} +a_{i}\right\} \text {,}$$\end{document} and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} is a random subject effect that is assumed to be normally distributed:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned}a_{i} \sim N \left( 0 ,\phi _{a}^{2}\right) \text {.} \end{aligned}$$\end{document}

At the subject level, the random intercept and the random effect of the within-subject residual variance may covary: \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\phi _{\upsilon a} =c o v \left( \upsilon _{i} ,a_{i}\right) $$\end{document} . The model, henceforth referred to as Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} , has five parameters of central interest, along with the parameters relating to the measurement model of y: \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathbf {\gamma } =\left( \mathbf {\tau } ,\mathbf {\lambda } ,\mathbf {\psi } ,\beta _{0} ,\alpha _{0} ,\tau _{0} ,\phi _{a}^{2} ,\phi _{\upsilon a}\right) \text {.}$$\end{document} Code developed for use with SAS IML for estimation of the model is provided in Appendix B.

Table 2.

ML estimates of a latent variable mixed-effects location scale model for daily positive affect ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$n=435$$\end{document} ).

Notes: ML estimation was carried out using Gauss–Hermite quadrature with 10 quadrature points to evaluate the integral. \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$^{a}$$\end{document} ML estimates are given in Appendix C

ML estimates of Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} are given in Table 2. The values presented in the table are based on an analysis carried out using a trust-region optimization method (Moré and Sorensen Reference Moré and Sorensen1983) available through SAS IML and used to optimize the likelihood function. The solution was evaluated using quadrature points that ranged from 5 to 30 with increments of 5 points, resulting in no appreciable differences in the parameter estimates or standard errors when using the different number of quadrature points. The estimates provided in Table 2 are based 10 quadrature points. To obtain starting values for estimation, results from the factor analysis model supplied values for the measurement model. A fixed-effects model was then fit, followed by adding a random intercept, then a random effect in the within-subject residual variance model with its covariance with the random intercept. To assess whether there was appreciable heterogeneity of variance in the occasion-level residuals, estimates were obtained for a model (Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{1}$$\end{document} ) that assumed \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\phi _{a}^{2} =0$$\end{document} and compared to those under the model that assumed \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\phi _{a}^{2} >0$$\end{document} (Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} ). The model assuming \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\phi _{a}^{2} >0$$\end{document} results in an improvement in model fit in terms of the AIC, from 33,710 down to 33,070, suggesting heterogeneity of the within-person residual variance. A LRT of the goodness-of-fit between the two models using 2 df is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\chi _{2}^{2} =33,678 -33,034 =644\text {,}$$\end{document} \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p <.001\text {,}$$\end{document} supporting the need for the random effect in the residual variance model. As noted earlier, the assumption that the uniquenesses for each item were equal across the days of measurement was not tenable according to a LRT. Thus, Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} was fit relaxing the assumption that the uniquenesses were equal, with this model henceforth called Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{3}$$\end{document} . Estimates from this model are given in Table 2 (with additional estimates provided in Appendix C). A LRT comparing Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{2}$$\end{document} and Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{3}$$\end{document} supports the need to allow for the unique values. From the estimates based on Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{3}$$\end{document} , the estimated correlation between the random intercept and the random subject effect for the within-subject residual variance is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$r = -.43$$\end{document} , an increase from \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$r = -.36$$\end{document} in the model that relied on mean scores, but again suggesting a tendency for those with a relatively low daily variation to also have a relatively high mean positive affect score across days. Similar to the analysis that relied on the means of the item response variables, the negative correlation suggests a tendency for those who have a relatively high mean level of positive affect across days to also have a relatively low degree of day-to-day variability. Thus, individuals who had higher estimated averages in affect across days also showed greater stability in their affect across days. Again, two possible conclusions from this are that those with greater overall positive affect are more stable in their daily positive affect or there is a ceiling effect of measurement. Available in supplementary materials are empirical Bayes estimates of the location and scale values based on Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{3}$$\end{document} , along with estimated conditional residuals between the observed item scores and the fitted values based on the ML estimates from fitting Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{3}$$\end{document} .

In the previous analysis that relied on sample mean scores, it was not possible to estimate the state residual variance of positive affect scores without the measurement error naturally inherent in manifest scores. In the latent variable model that accounts for measurement error, the estimated state residual variance, assuming homogeneity of the within-person residual variance of daily latent scores, is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\exp \left\{ -1.93\right\} =0.15$$\end{document} , and the estimated trait variance is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\exp \left\{ -2.41\right\} =0.090$$\end{document} , with the state variance almost double the size of the trait variance. With a relatively high residual state variance, variation in positive affect seems to be due more to daily variability rather than between-person differences. Relaxing the assumption of homogeneity of the within-person residual variance, however, gives a different description of the data. First, given the improvement in model fit, it does not seem reasonable to conclude that the residual state variance is common across individuals, but rather, that it varies according to person. The estimated residual state variance for the typical individual is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\exp \left\{ \tau _{0} +.5 \phi _{a}^{2}\right\} $$\end{document} (Hedeker et al. Reference Hedeker, Mermelstein and Demirtas2008) \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$ =\exp \left\{ -2.56 +.5 \left( 1.32\right) ^{2}\right\} =0.18$$\end{document} and the estimated trait variance is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\exp \left\{ -1.21\right\} =0.30$$\end{document} , with the trait variance now being greater than the magnitude of the residual state variance, thus suggesting that, for the typical individual, variation in positive affect scores is more highly attributable to between-person differences than within.

As discussed earlier, applications of a mixed-effects location scale model generally do not account for measurement error in response data. If scores include measurement error, the residuals reflect a combination of within-subject variation both about the fitted model and that due to measurement error. To evaluate the impact of partitioning the variance of the manifest scores as can be done using a latent variable model, a version of Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{3}$$\end{document} was fit in which a single variance was used to represent these combined sources. The factor loadings, item intercepts, factor mean, standard deviation of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} and correlation between \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$a_{i}$$\end{document} and \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\upsilon _{i}$$\end{document} were set equal to the estimates resulting from fitting Models B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{3}$$\end{document} , resulting in an estimated combined common variance of \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\exp \left\{ -0.77\right\} =0.46$$\end{document} . From Model B \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{3}$$\end{document} that partitions the variance due to measurement error and variation due to the fitted model, the common variance at the occasion level is estimated as \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\exp \left\{ -2.56\right\} =0.08$$\end{document} , a substantial decrease from that in which the two sources of variation are combined.

3.2.1. Daily Stressors

To study the within- and between-subject effects of daily stressors on latent positive affect, the person-mean centered and mean number of daily stressors were included in the latent variable location scale model, similar to what had been done when evaluating positive affect using the manifest scores directly:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned}\eta _{t i} =\beta _{0 i} +\beta _{1} \text {\textit{DS}}_{t i}^{ *} +\beta _{2} \text {\textit{DS}}_{ .i} +\varepsilon _{t i}\text {.} \end{aligned}$$\end{document}

Additionally, the two measures of daily stressors were included in the model of the occasion-level residual variance relating to the latent affect measures:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned} \sigma _{\varepsilon }^{2} =\exp \left\{ \tau _{0} +\tau _{1} \text {\textit{DS}}_{t i}^{ *} +\tau _{2} \text {\textit{DS}}_{ .i} +a_{i}\right\} \text {,} \end{aligned}$$\end{document}

and the individual mean number of stressors was included in the variance model of the random intercept:

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\begin{aligned}\phi _{\upsilon }^{2} =\exp \left\{ \alpha _{0} +\alpha _{1} \text {\textit{DS}}_{ .i}\right\} \text {.} \end{aligned}$$\end{document}

Starting values for the effects of covariates were set to 0. Estimates of this model are in the last column of Table 2. For those with no reported daily stressors, the point estimate of the overall positive affect score across days and individuals is \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\beta }_{0} =3.20$$\end{document} ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.13$$\end{document} ), a value that is higher than the point estimate of 3.05 obtained from the model that relied on the mean of the manifest variables. For daily positive affect, the estimates suggest an overall negative within-person ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\beta }_{1} = -0.10$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.011$$\end{document} ), as well as between-person ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\beta }_{2} = -0.38$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.060$$\end{document} ), effect of the number of daily stressors on daily mean positive affect. Thus, on average, the within-person effect indicates that a higher number of daily stressors tends to correspond to a lower daily positive affect level, and the between-person effect indicates that a higher mean number of daily stressors across days tends to correspond to a lower daily mean positive affect, with the latter being greater than the former (H \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{0} :\beta _{1} =\beta _{2}$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\chi ^{2} \left( 1\text {df}\right) =22$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p <.001$$\end{document} ). In the previous analysis based on sample mean scores to represent positive affect, the relative magnitudes of these two effects were somewhat different relative to the estimates from the latent variable model, but inference from both models in terms of the direction of the effects is the same with the between-subject differences in the number of daily stressors having the greater impact on positive affect. Conditional on these effects, greater daily variation in affect scores corresponds to a high number of daily stressors ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\tau }_{1} =0.30$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.071$$\end{document} ), as well as a higher number of mean daily stressors ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\tau }_{2} =0.89$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.17)\text {,}$$\end{document} , with the latter being greater than the former (H \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$_{0} :\tau _{1} =\tau _{2}$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\chi ^{2} \left( 1\text {df}\right) =9.7$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p =.002$$\end{document} ). These estimates are comparable to those obtained under the model that relied on manifest scores. Also under the latent variable model, the variance of the random intercept does not have a clear relationship with the mean number of daily stressors ( \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hat{\alpha }_{1} = -0.20$$\end{document} , \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$s e =0.17$$\end{document} ), unlike the result under the standard model where it seemed that individuals differed to a greater extent from each other as their mean number of daily stressors increased. From this, one might conclude that there is no appreciable association between the between-subject variation in positive affect levels and the mean number of stressors under the latent variable model but that the association is instead strong and positive when studied under the standard version of the model. Finally, a comparison of the estimated standard errors of the two models overall indicates that the precision of the estimates is generally better under the model that used the sample means scores. As the latent variable model is preferred for the sake of accounting for measurement error in the manifest measures, it is possible that the precision of the estimates is negatively biased in the standard version of the model when scores do contain measurement error and are best reflected by those values obtain under the latent variable model.