Next Article in Journal
Bayesian Semiparametric Regression Analysis of Multivariate Panel Count Data
Previous Article in Journal
Opening the Black Box: Bootstrapping Sensitivity Measures in Neural Networks for Interpretable Machine Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Repeated-Measures Analysis in the Context of Heteroscedastic Error Terms with Factors Having Both Fixed and Random Levels

Department of Statistics, University of South Africa, Johannesburg 1709, South Africa
*
Author to whom correspondence should be addressed.
Stats 2022, 5(2), 458-476; https://doi.org/10.3390/stats5020027
Submission received: 12 March 2022 / Revised: 10 April 2022 / Accepted: 15 April 2022 / Published: 6 May 2022

Abstract

:
The design and analysis of experiments which involve factors each consisting of both fixed and random levels fit into linear mixed models. The assumed linear mixed-model design matrix takes either a full-rank or less-than-full-rank form. The complexity of the data structures of such experiments falls in the model-selection and parameter-estimation process. The fundamental consideration in the estimation process of linear models is the special case in which elements of the error vector are assumed equal and uncorrelated. However, different assumptions on the structure of the variance–covariance matrix of error vector in the estimation of parameters of a linear mixed model may be considered. We conceptualise a repeated-measures design with multiple between-subjects factors, in which each of these factors has both fixed and random levels. We focus on the construction of linear mixed-effects models, the estimation of variance components, and hypothesis testing in which the default covariance structure of homoscedastic error terms is not appropriate. We illustrate the proposed approach using longitudinal data fitted to a three-factor linear mixed-effects model. The novelty of this approach lies in the exploration of the fixed and random levels of the same factor and in the subsequent interaction effects of the fixed levels. In addition, we assess the differences between levels of the same factor and determine the proportion of the total variation accounted for by the random levels of the same factor.

1. Introduction

Linear mixed-effects models can handle both fixed and random effects simultaneously. Not only are the models convenient for modelling the means of the data, but they are also convenient for modelling covariances [1,2]. Situations arise in fields such as psychology and medicine in which longitudinal or correlated data display notorious heterogeneity of responses to stimuli and treatment [3]. When the classical linear mixed-effects model is inappropriate, alternative approaches, which identify and consider subgroups in fixed and random effects, are necessary. The authors of [4] propose a heterogeneity model as an extension to that proposed by the authors of [5,6] by replacing the normality assumption for random effects and by assuming a non-zero mean vector and a common error vector.
Repeated-measures data are often dependent, a property that does not conform to the generality of a mixed-effects model [1,7,8]. Linear mixed-effects models are some of the most convenient statistical approaches, which account for this dependency [9]. However, setting them up for data analysis requires some care, especially in choosing the most appropriate covariance structure to keep the type I error down [10]. Selecting a suitable model for the covariance is important because the precision of the confidence intervals and the tests of hypotheses concerning model parameters depend on the correct model [11]. According to [12], selecting an appropriate mixed-effects model and construction approach, such as partitioning of fixed and random effects, allows for analysis of variance in correlated data. Numerous diagonal and non-diagonal covariance structures for correlated data, which cover a range of assumptions about the associations between responses from the same cluster data, are available [13,14]. Modern statistical software, such as SAS Studio in Linux Environment, provides options for selecting candidate covariance structures through the PROC MIXED procedure [15] (see Appendix A). The residual maximum likelihood (REML) is one of the most famous methods for estimating covariance parameters associated with linear mixed models [16,17], among other alternatives [18].
We present the construction and analysis of a three-factor linear mixed model for repeated-measures designs when the between-subjects factors consist of both fixed and random levels, and the structure of the variance–covariance matrix of the error terms is not the identity. We consider a repeated-measures design setting under a linear mixed-effects model, with factors sharing both fixed and random components of the model. For experimental designs consisting of factors with a unique composition of both fixed and random levels, we propose a partitioning approach based on factor levels in model construction, estimation, and hypothesis testing. In such a situation, the fixed levels allow for the comparison of specific levels of interest within the factor, whereas the random levels allow for the assessment of variation within the same factor. We assess the effect of introducing heterogeneity of error terms in the selection of the most appropriate covariance structure, in the assessment of the changes that occur in the estimation, and when drawing inferences. To simplify the proposed approach, we focus on the diagonal covariance structure for the repeated-measures linear mixed-effects model.
Section 2 presents the approach to model construction, for a general linear mixed-effects model in a completely randomised design (CRD) and in a repeated-measures design (RMD) under the default assumption of error terms. Section 3 presents results from a numerical example. Discussion of the results is in Section 4, followed by the conclusions in Section 5.

2. Materials and Methods

2.1. An Illustrative Data Structure

We motivate the approach using data collected from a study that investigated the impact of combining carbon tetrachloride ( C C l 4 ) with four levels (0, 1.0, 2.5, and 5.0 mM) and chloroform ( C H C l 3 ) with four levels (0, 5, 10, and 25 mM) on the toxicity of cells on in vitro toxicity of isolated hepatocyte suspensions [19]. Four flasks were assigned to each of the 16 treatments. Cell toxicity is measured by the amount of lactic dehydrogenase (LHD) enzyme percentage leakage from each of the 64 flasks after 0.01, 0.25, 0.5, 1, 2, and 3 h since the application of the treatment. For illustration purposes, we consider the between-subjects factor C C l 4 levels 2.5 and 5.0 as fixed (new technology), and levels 0 and 1.0 are considered as existing random levels (old technology). Similarly, we consider the between-subjects factor C H C l 3 levels 10 and 25 as fixed (new technology), and levels 0 and 5 are taken as random levels (old technology). Of interest in the analysis is the percentage leakage observed at times 1, 2, and 3. We demonstrate the model construction procedure under certain assumptions in a completely randomised design (CRD) and in a repeated-measures design (RMD).

2.2. Construction of a Linear Mixed-Effects Model in CRD

Consider a three-way treatment structure in a balanced, completely randomised design (CRD) with full interaction of factors A, B, and C, each consisting of f fixed and r random levels. Assume we have f A ,     f B , and f C fixed levels and r A ,   r B , and r C random levels of factor A, respectively. We partition the dataset based on the combinations of factor levels and construct a partitioned model in each partition. For example, the FRF partitioned model is built from the f A , r B , and f C levels. Similarly, other possible partitions are FFF, FFR, RFF, RRF, RFR, FRR, and RRR. We illustrate the model construction using the FRF linear mixed-effects model in CRD, having at least one replication per treatment combination and expressed as
y FRF i j k l = μ FRF + φ A i + φ B j + φ C k + π 1 + + π t + ϵ FRF i j k l ,
where y FRF i j k l is the lth observation in the (ijk)th treatment cell of the FRF partition; l = 1 , ,   r h   are the replicates (where all r h = r for balanced data); μ FRF is the overall mean; φ A i ,   φ B j ,   φ C k are the main effects of the three factors; π 1 ,   ,   π t are the interaction effects; and φ A i   i = 1 ,   2 , , f A ,   f A + 1 , f A + 2 , ,   a   a = f A + r A , φ B j   j = r B + 1 , r B + 2 , ,   b   b = f B + r B , and φ C k k   k = 1 ,   2 , , f C ,   f C + 1 , f C + 2 , ,   c   c = f C + r C are unknown parameters corresponding to fixed factor A, random factor B, and fixed factor C, respectively. Defining the random main effect as φ R and the random interaction effect as π R in (1), the random effects and the random error term ϵ i j k l , are commonly assumed to have zero mean and variance, i.e., φ R   ~   N 0 , σ φ R 2 , π R   ~   N 0 , σ π R 2 , and ϵ i j k l   ~   N 0 , σ ϵ 2 .
For a balanced data scenario with r replications per cell, for example, the general linear mixed model Equation (1) is normally expressed in matrix form as
y FRF = X FRF β + Z FRF u + ϵ FRF ,
where y FRF : N × 1 is a vector of response observations in the FRF partition; matrix X FRF : N × p is a known incidence matrix associated with the vector of p fixed-effects β : p × 1 in the model; matrix Z FRF : N × q is a known incidence matrix associated with the vector of q random-effects u : q × 1 in the model; and ϵ FRF : N × 1 is a vector of random errors. The usual assumption under this model is that the random effects are u   ~   N 0 , G , and the random residuals are ϵ FRF   ~   N 0 , R , where R = σ 2 N and G is a diagonal matrix of variance components (i.e., different variances, and all zero covariances).
G = σ 1 2 I r 0 0 σ t 2 I r   and   R = σ ϵ 2 0 0 σ ϵ 2
The total variance–covariance is the structured matrix V = Z G Z + R , a structure that guarantees independence and homogeneity of residual errors. This implies that the variance of y is modelled through Z , G , and R . The simple total variance–covariance, V , has a block-diagonal structure given by the matrix
V = σ 1 2 I r + σ ϵ 2 I r 0 0 σ t 2 I r + σ ϵ 2 I r

2.3. Linear Mixed-Effects Model in RMD

Traditionally, between-subject and within-subject factors in repeated-measures experiments are designated as either fully fixed effects or random effects. Depending on the objectives of the experiment, some factors in linear mixed-effects models may exist with both fixed and random levels [12,20,21]. The same scenario is common with a repeated-measures experiment, in which either the between-subjects factor or the repeated measures consist of both fixed and random levels. For instance, for improved results, a researcher may decide to consider additional levels of a between-subjects factor, in addition to the old and existing levels. In that case, the improved analysis needs to consider the new factor levels as fixed levels, and the old and existing levels are considered random. The approach creates an opportunity to compare and evaluate the effectiveness of new factor levels (fixed) against existing ones (random) and/or to compile a combined analysis of both. In addition, random levels allow for the assessment of variation between and within the factor levels in the entire population.
We consider a three-factor repeated-measures experiment with n experimental units (EU) that are randomly assigned to each of the a levels of the between-subjects factor A (with f a fixed levels and r a random levels, a = f a + r a , for example), the between-subjects factor B (with f b fixed levels and r b be random levels, b = f b + r b , for example), and the within-subjects factor C with t measurements (all considered to be fixed in this case) taken on each of the experimental units (EU).
y p i j k l = μ p + α i + β j + α β i j + γ k i j + τ l + α τ i l + β τ j l + α β τ i j l + ϵ p i j k l ,
where the subscript p in y p i j k l and ϵ p i j k l denotes the partition p = 1 ,   2 , ,   4 and where y p i j k l is the lth measurement ( l = 1 ,   2 , , t ) of the kth experimental unit ( k = 1 ,   2 , n ) in the ith level ( i = 1 ,   2 , a ) of factor A and the jth level ( j = 1 ,   2 , b ) of factor B in the pth partition. Depending on whether the effects are fixed or random, the model parameters are as defined in (1).

2.4. Model Assumptions

In experimental research that involves analysis of variance (ANOVA) as a technique for comparing different treatment means, a set of assumptions that includes the usual normality and homogeneity of variance must be checked before the analysis of data [22]. These traditional normality tests, such as Q–Q plots or the Shapiro–Wilk test, and outlier detection approaches (e.g., box plots), are appropriate for diagnosing violations of assumptions.

2.4.1. Sphericity (Circularity) Assumption

Similar to the homogeneity of variances in a between-subjects analysis of variance, sphericity holds in ANOVA for repeated-measures designs when the variances of the differences among all possible pairs of the within-subject factor group means are equal [23,24]. The assumption is usually unrealistic in repeated-measures designs in which observations are correlated. Univariate tests for within-subjects effects apply when sphericity holds. If the sphericity assumption is violated, several alternatives for adjusting the numerator and denominator degrees of freedom have been suggested [25]. The most popular degrees of freedom corrections are the Greenhouse–Geisser [26] or the Huynh–Feldt [27] adjustments, which are extensions of the Box correction factor [28]. These adjustments provide a more accurate adjusted p-value [29,30]. Failure to address the problem of sphericity when conducting an analysis of variance often leads to inflated F-ratios, type I errors, biased post hoc tests for group mean differences, and inaccurate conclusions [24]. Tests for sphericity are offered in most statistical computing packages (SAS, R, SPSS, etc.).

2.4.2. Compound Symmetry Assumption

An overly restrictive assumption closely related to sphericity is the compound symmetry assumption. It states that there is a constant variance and correlation between observations from the same subject, which is not always realistic in many repeated-measures applications [31]. The Greenhouse–Geisser [26] or the Huynh–Feldt [27] adjustments are used to circumvent restrictive compound symmetry assumptions and to accommodate more general covariance structures for the repeated measures [1]. Compound symmetry implies sphericity, but not vice versa. Univariate tests for within-subject effects apply when compound symmetry holds. If compound symmetry (and hence sphericity) holds, then a split-plot analysis is used as an appropriate approximation to the repeated-measures experiment, since it provides relatively more accurate p-values for testing treatment effects [31]. Compound symmetry is simplified in the Huynh–Feldt [27] condition for sphericity. Hence, both are tested in many software packages by Mauchly’s test [32]. However, Mauchly’s test for sphericity has been criticised for its over sensitivity and tendency to reject compound symmetry [33].

2.5. Estimation Techniques

The general linear mixed-effects model (2) is used to describe data from partitioned repeated measurements, wherein the fixed-effects component, X p β , consists of the design matrix X p , and the fixed-effects coefficients β are as defined in (3). The random-effects component, Z p u , contains the block-diagonal random-effects design matrix Z p with design matrices for the individual subjects ( Z i ); u is a vector of random coefficients (the between-subject variance–covariance components), and ϵ p denotes the within-subject errors from the pth partition. The random effects and the residuals follow the distribution u   ~   N 0 , G and ϵ p   ~   N 0 , R , respectively, where G is a block-diagonal covariance matrix of the random effects and where R is a diagonal covariance matrix with partitions corresponding to each subject (within-subject errors) in the analysis. The covariance matrix for the repeated-measures data is composed of matrices Z p , G , and R , which is a block-diagonal Σ p = V a r y p = Z p G Z P + R . The non-singular components G and R are usually estimated by two principal likelihood methods for estimating variance components [34], i.e., the maximum likelihood (ML) method and restricted maximum likelihood (REML). These procedures are available in various mixed-model statistical software, such as SAS (PROC MIXED procedure) and R, with the REML estimates generally preferred unless the data sets are quite large [34].
Assuming that both the random effects and the error terms are normally distributed, the likelihood function for the repeated-measures mixed model is given by [35]:
l = log L y p = N 2 log 2 π 1 2 log Σ p 1 2 y p X p β Σ p 1 y p X p β = C 1 2 log Σ p 1 2 y p X p β Σ p 1 y p X p β
where y p and V = Σ p are as defined in (2). Similarly, a modification of the ML procedure through factorisation of the likelihood function is proposed as an alternative method of estimating covariance parameters, which is the restricted maximum likelihood function [36]:
l r e = log L y p = N 2 log 2 π 1 2 log Σ p 1 2 y p X p β Σ p 1 y p X p β = C 1 2 log X p Σ p 1 X p 1 2 log Σ p 1 2 y p X p β ^ Σ p 1 y p X p β ^  
where the available covariance matrix Σ p is used to estimate the fixed-effects parameters,
β ^ = X p Σ p 1 X p 1 X p Σ p 1 y p .  
The main challenge in repeated-measures analysis of variance is to determine the adequate correlation structure, because the constant variance assumption for the distribution of the error terms is likely not to be reasonable for the distribution of error terms within subjects. There are various possible choices of covariance structures for repeated measures within each subject depending on the chosen parameterisation for G and R . The choices are usually guided by limitations of the software and by insight from the researcher. The most common covariance structures include variance components, compound symmetry (common covariance plus diagonal), unstructured (general covariance), and autoregressive [37]. With the PROC MIXED statement in SAS, one can specify any repeated-measurements covariance structure for G by using the Random statement and by specifying the form of R with the Repeated statement, in conjunction with the Type option [38]. Excluding the Repeated statement specifies the classical R , which is assumed to be equal to σ 2 I .
There are numerous ways of identifying the most appropriate covariance structure amongst set candidate structures [37]. The most recommended approach is to select the structure that gives the smallest Akaike’s Information Criterion (AIC) [39], a statistic that is defined by the model and the maximum likelihood estimates of the parameters from specifying the variance–covariance as
AIC = 2 L β ^ , Σ ^ p + 2 g ,  
where g is the effective number of independently adjusted parameters in the covariance matrix from the pth partition, and L β ^ , Σ ^ p = l o g M L is the value of the likelihood function evaluated at β ^ , Σ ^ p . A better model is the one with the smallest AIC value. Different forms of R can be compared for adequacy using the likelihood ratio test statistic [38]. The hypotheses involved are the following: H 0 :   R 1 is as adequate as R 2 ; and H 1 :   R 1 is not as adequate as R 2 , where R 1 is a special case of R 2 . Suppose R 1 and R 2 have g 1 and g 2 parameters, respectively, with ( g 1 < g 2 ) , and the test statistic is Q = 2 L β ^ 1 , Σ ^ 1 L β ^ 2 , Σ ^ 2 , which is distributed as χ 2 g 2 g 1 . We reject H 0 when Q χ α 2 2 g 2 g 1 .

2.6. Methods of Inference

We present the algorithm for obtaining expected mean squares using the ANOVA approach for the FRF model.

2.6.1. Process for Deriving Expected Mean Squares

(a)
Based on the model involved, construct a two-way table with column headings corresponding to the source of variation, effect labels, each of the subscripts included in the model, and row headings corresponding to each source of variation in the ANOVA table.
(b)
Above each subscript, write the associated number of factor levels and insert on top either an “F” if the factor levels are fixed or an “R” if the factor levels are random.
(c)
Create an extra column on the extreme right for the variance components corresponding to the source of variation, and insert the appropriate random variance component ( σ . 2 ) or fixed variance component ( θ . ) for each source of variation.
(d)
Compare the column subscript and the factor effect in each row, and write the number of levels corresponding to that subscript if the column subscript is not included in the factor effect label. Otherwise, leave it blank.
(e)
For rows that have effects that contain bracketed subscripts, write a “1” under the column if the subscript is included in the bracket.
(f)
For each row that has a fixed variance component ( θ . ), put a zero in the cell headed by an “F” when the subscript is included in the effect label.
(g)
Enter a “1” in all remaining blank cells.
(h)
To obtain the expected mean squares for each effect, identify all the variance components associated with that effect label. Cover the column(s) headed by the effect subscript(s) in that effect, and obtain the coefficient of each of the identified components from the product of the entries in the column(s) headed by the uncovered subscript(s). Include the variance component σ ϵ 2 with the coefficient of 1 in the list.
Figure 1 illustrates the key steps to be taken when constructing the coefficients of variance components for a partitioned linear mixed-effects model.
Once the coefficients of variance components are established in step (m) of Figure 1, the expected mean squares are found by Step (h) of the Process in Section 2.6.1. Table 1 summarises the variance components obtained for a three-factor repeated-measures design, with two between-subjects factors ( A and B ) on one within-subject factor ( C ) using the Process in Section 2.6.1. We use this process to build an FRF model, which assumes factors A and C as fixed, factor B as random, and experimental units (EU) as random. The F, R, F, and Rep denote the effect; a ,   b ,   t , and n denote levels of factors A and B, treatment number, and sample size, respectively. The subscripts i ,   j ,   l , and k denote the indices of the corresponding levels.
For example, E M S A , with effect α i , is composed of the variance components θ α , σ α β 2 , σ γ α β 2 , θ α τ , σ α β τ 2 , and σ ϵ 2 as follows:
E M S A = σ ϵ 2 + t σ γ α β 2 + t n σ α β 2 + n σ α β τ 2 + b t n θ α .
Table 2 displays the ANOVA layout and the expected mean squares for a three-factor repeated-measures design when one of the factors is a within-subject factor. The FRF model is considered for illustration, with experimental units assumed to be random.

2.6.2. Hypothesis Testing for Fixed Effects

We are interested in testing the main and interaction effects of the between-subjects factor and of the within-subjects factor in both the partitioned and combined repeated-measures linear mixed-effects model. In addition to checking model assumptions, the following hypotheses are of interest for each partitioned model:
Hypothesis 1 (H1).
Between-subjects main and interaction effects (i.e., H 0 : α i = 0 ).
Hypothesis 2 (H2).
Within-subjects main and interaction effects (i.e., H 0 : τ i = 0 ).
The test statistic for H1 is given by F = M S A M S U A × B   ~   F a 1 , a b n 1 α , and the test statistic for H2 is given by F = M S P M S E   ~   F t 1 , a b t 1 n 1 α . The interaction effects among the between-subjects and the interaction effects among within-subjects are tested by the M S U A × B and the M S E on the denominator, respectively.

2.6.3. Hypothesis Testing for Random Effects

Variance components are estimated by equating the mean square to the expected mean squares derived in Table 2. Where there are no valid F tests, approximate F tests are constructed for the sources of variability in random effects [40]. For random factor B, the hypothesis of interest may be
Hypothesis 3 (H3).
Random effects (e.g., H 0 : σ β 2 = 0 , against H 1 : σ β 2 > 0 ).

2.6.4. Combined Analysis

The individual partitioned models provide pieces of information which are needed for an integrated analysis. The combined model is built by combining the degrees of freedom and sum of squares associated with each source of variation for each appropriate hypothesis test. For example, the combined effect of the within-subjects factor ( P e r i o d ) in the FC model is obtained from the partitions in which the factor P e r i o d is fixed, i.e., the pieces of information are supplied by the partitioned models FFF, FRF, RFF, and RRF. Similarly, the other main and interaction effects for the combined model are obtained by summing up the associated degrees of freedom and sums of squares.

3. Results

3.1. Checking Model Assumptions

A three-factor ( C C l 4 , C H C l 3 , and T i m e ) repeated-measures experiment on one of the factors ( T i m e ) is proposed. Table 3 shows the multivariate data (wide format) layout for a three-factor repeated-measures experiment.

3.1.1. Normality and Outliers

Q–Q plots are used to check normality and outlier assumptions simultaneously. The FFF, FRF, RFF, and RRF data sets do not show any serious deviations from normality. Figure 2 shows the normal Q–Q plots for the four data subsets.
Furthermore, the Q–Q plots do not show any influential point (outlier) that warrants exclusion, since all plots are not very far from the diagonal.

3.1.2. Sphericity and/or Compound Symmetry

In order to test the suitability of using a repeated-measures design in the experiment, the sphericity (or compound symmetry) assumption is tested for each partitioned data subset. Sphericity test results are produced by fitting two models using the PROC MIXED procedure: one specifying the unrestricted covariance structure and the other with a less conservative Huynh–Feldt (H–F) adjustment in the Type option [41]. The null hypothesis for the test is:
Hypothesis 4 (H4).
The covariance structure fits the sphericity structure.
The difference ( D ) between the two 2 log-likelihoods of the two competing models follows a Chi-square distribution with degrees of freedom equal to the difference in the numbers of parameters in them.
  • FFF model: χ 35 2 0.05 55.76 , d f = 35 , D = 1722.6 , significant;
  • FRF:   χ 38 2 0.05 55.76 , d f = 38 , D = 184.54 , significant;
  • RFF: χ 37 2 0.05 55.76 , d f = 37 , D = 288.8 , significant;
  • RRF: χ 36 2 0.05 55.76 , d f = 35 , D = 384.1 , significant.
The sphericity assumption H4 fails in all the four (4) partitioned data sets.

3.2. Analysis of Results

We use the restricted maximum likelihood estimation (REML) approach to estimate these variance components. Table 4 contains the estimated AIC from each of the covariance structures as well as the number of covariance parameters estimated with a non-zero value.
The covariance structure ARH(1) has the minimum AIC values in models FFF, FRF, and RRF, and covariance structure AR(1) has the smallest AIC value in model RFF. The covariance structure ARH(1) is chosen as the most adequate covariance structure for the partitioned models based on AIC.
When the appropriate covariance structure is incorporated into the model, the relationships among the errors and variance components are specified. The PROC MIXED procedure does not require further assumptions on them [42], making it a robust and flexible procedure [15] that produces valid F-tests regardless of whether the sphericity assumption is satisfied or not. Table 5 summarises the PROC MIXED F-tests for the main and interaction effects of the between-subjects and the within-subjects factors for the four partitioned models.
The C H C l 3 factor has a significant effect on the leakage percentage (p-value < 0.001) over time in the FFF partition. Based on the fixed-effects (FFF) repeated-measures study, it may be concluded that chloroform ( C H C l 3 ) has a significant impact on the amount of lactic dehydrogenase (LHD) enzyme percentage leakage (toxicity of cells) over time, whereas neither carbon tetrachloride ( C C l 4 ) in isolation nor the interaction thereof has significant influence. Furthermore, the T i m e factor plays an important role in determining the amount of leakage as well. However, the interaction of old (random) and new (fixed) levels of the between-subjects factor levels has non-significant effects in FRF, RFF, and RRF models when C H C l 3 fixed levels are involved.
Table 6 gives a summary of the estimated covariance parameters and the proportion (in percentage) of variation they contribute in each of the partitioned models. A zero variance for the random levels of chloroform ( C H C l 3 ) and very small estimates in other factors are obtained, which result in very low estimates of covariance parameters in the models.
The CHCl3 and CCl4 random levels and their interaction have a noticeable contribution to the proportion of variation in the amount of lactic dehydrogenase (LHD) enzyme percentage leakage in the RRF partition. Generally, Time has very little interaction effect with the between-subjects factors in FRF and RRF models in determining the proportion of variation of the toxicity of cells.
Since the mixed-model methodology directly computes neither sums of squares nor F-statistics from the ratio of mean squares, and since the R CRAN package l m e 4 for mixed models currently does not have options for other covariance structures to cater for correlated error variances, generating a combined analysis is not a straightforward exercise. We scrap the targeted data subset based on the effects of interest before fitting the models. Analogous to the partitioned analyses, a comparison of model fit via the AIC approach is conducted. For convenience purposes, let the factors C C l 4 , C H C l 3 , and T i m e be labelled as factor A , B , and C , respectively. The PROC MIXED procedure is used to fit the repeated-measures linear mixed models for intended narrow inference spaces [12,43]. Table 7 shows the Type III tests for the combined models.
Of the possible candidate covariance structures (CS, CSH, AR(1), and ARH(1)), structure ARH(1) is selected as the most appropriate covariance structure for the combined fixed-effects model FA, and AR(1) is appropriate for FB and FA × FB. The factors B (CHCl3) and C (Time) have significant effects (p- value < 0.05 ) in the combined models. The broad inference scope results for the combined models (assuming random factor A or B effects) are similarly analysed.

4. Discussion

Based on the illustrative example results, the approach manages to isolate the effects of new and old factor levels over time. The combined analysis confirmed the results of the partitioned analysis on the percentage leakage in cells. The proposed approach conforms to the model construction and the analysis procedures in a repeated-measures design. It can be used as a planning tool in which factor combination and time are of interest in designing experiments that involve repeated measures. In such experiments, blindly adopting the assumption of homogeneous error terms without exploring possible candidate covariance structures may compromise the ability of an experiment to detect sufficient variation in the response variable. In addition, our approach enhances the accuracy of inferences by providing partitioned analysis of heterogeneous variances and covariance structures, which sometimes are not identical in the data subsets.
Given the increased complexity of research data in various research fields, the application of a linear mixed-model methodology must be in line with data covariance structures for accurate results to be achieved. One of the approaches that has proved to be a reliable tool for managing big data complexity issues is the partitioning approach [12,20,21], in which the traditional homogeneous error variance structure is assumed. The current study extends the new approach to a three-factor treatment structure in a repeated-measures design in which linear mixed-effects models are applicable. In essence, the approach can be extended to cater for repeated-measures experiments in which any number of between-subjects and within-subjects factors are involved. In most cases, repeated-measures experiments do not assume equal and uncorrelated error vectors, since regularly timed measurements taken on the same subject over time are usually correlated [34].
For the fixed-effects partitions, the linear mixed-effects models for a repeated-measures design are fit by the PROC GLM procedure, and a combined analysis of these can be obtained by syncretising the sum of squares and degrees of freedom from the fit models. However, obtaining a combined analysis using the SAS PROC MIXED procedure is impossible using the sum of squares approach, since the PROC MIXED procedure uses a likelihood-based estimation scheme instead of the least-squares method. A comparable alternative to reduction in the sum of squares for the fit model in PROC MIXED is to consider the amount of information retained by the fit model when compared to the null model.

5. Conclusions

The proposed approach allows for construction and hypothesis testing in repeated-measures data if a heterogeneous error structure is assumed. Although the MLE method used in SAS PROC MIXED does not estimate sums of squares, data scrapping based on targeted factor levels work equally well as an alternative approach to obtain the combined analysis. The proposed approach can be adopted as an essential tool for the comparison of new inventions against existing strategies and equipment. It leads to an exploration of the fixed and random levels of the same factor and the subsequent interaction of levels of factors of interest. We can assess the differences between levels of the same factor and understand variation within the same factor. In addition, modelling allows for the assessment of various covariance structures. Although this paper focuses on a limited scope of covariance structures, we pose an open research problem on the application of the proposed approach to other designs of experiments that incorporate more complex and non-diagonal covariance structures.

Author Contributions

Conceptualisation, L.C. and P.N.; methodology, L.C.; software, L.C.; validation, L.C. and P.N.; formal analysis, L.C.; writing—original draft preparation, L.C.; writing—review and editing, P.N.; supervision, P.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received publication funding from Sol Plaatje University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data for this research are available in [19,31].

Acknowledgments

The authors would like to acknowledge the support from the University of South Africa, and from the Sol Plaatje University for the publication costs provided.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AICAkaike’s Information Criteria
ANOVAAnalysis of variance
AR(1)First-order autoregressive
ARH(1)Heterogeneous first-order autoregressive
CCl4Carbon tetrachloride
CHCl3Chloroform
CRDCompletely randomised design
CSCompound symmetry
CSHHeterogeneous compound symmetry
EUExperimental unit
FFFFixed-Fixed-Fixed
FFRFixed-Fixed-Random
FRFFixed-Random-Fixed
FRRFixed-Random-Random
LHDLactic dehydrogenase
REMLRestricted maximum likelihood
RFFRandom-Fixed-Fixed
RFRRandom-Fixed-Random
RRFRandom-Random-Fixed
RRRRandom-Random-Random
SASStatistical Analysis Systems
VCVariance components

Appendix A

SAS Code for Mixed-Model Analysis
/* Fitting FFR Model in Proc Mixed */
FILENAME REFFILE ‘/home/u35581214/LDH Leakage Data.sav’;
PROC IMPORT DATAFILE = REFFILE
   DBMS = SAV
   OUT = LDH;
RUN;
/* View Repeated-Measures Data in Multivariate Form */
Proc print data = LDH;
run;
/* Set Repeated-Measures Data to Univariate form */
Data LDH_mult(keep = CCl4 CHCl3 Flask Time4 Time5 Time6)
LDH_univ(keep = CCl4 CHCl3 Flask Time Leakage);
set LDH;
output LDH_mult;
Leakage = Time4;Time = 1; output LDH_univ;
Leakage = Time5;Time = 2; output LDH_univ;
Leakage = Time6;Time = 3; output LDH_univ;
run;
/* View Data in Univariate and Multivariate Form */
Proc print data = LDH_univ;
run;
Proc print data = LDH_mult;
run;
/* Subset or partition FRF from LHD univariate original Data */
Data FRF;
set LDH_univ;
if (CCl4 = 2.5 AND CHCl3 = 0) then output;
if (CCl4 = 2.5 AND CHCl3 = 0) then output;
if (CCl4 = 2.5 AND CHCl3 = 0) then output;
if (CCl4 = 2.5 AND CHCl3 = 5) then output;
if (CCl4 = 2.5 AND CHCl3 = 5) then output;
if (CCl4 = 2.5 AND CHCl3 = 5) then output;
if (CCl4 = 5 AND CHCl3 = 0) then output;
if (CCl4 = 5 AND CHCl3 = 0) then output;
if (CCl4 = 5 AND CHCl3 = 0) then output;
if (CCl4 = 5 AND CHCl3 = 5) then output;
if (CCl4 = 5 AND CHCl3 = 5) then output;
if (CCl4 = 5 AND CHCl3 = 5) then output;
run;
Proc print data = FRF;
run;
/* Plot differences in leakages contributed by predictors */
proc means noprint data = FRF nway;
var Leakage;
class CCl4 CHCl3 Flask Time;
output out = avgFRF mean = avgLeakage;
run;
proc print data = avgFRF;
run;
/* New data set called avg created*/
/* Plot differences in leakage by predictor CHCl3 and Time */
Proc gplot data = avgFRF;
plot avgLeakage*Time = CHCl3/haxis = 0 to 8 by 1 hminor = 0 vminor = 0;
symbol1 v = star c = blue i = join l = 1;
symbol2 v = plus c = red i = join l = 2;
titlePercentage leakage per time per CHCl3”;
run; Quit;
/* Partitioning FRF multivariate data for covariance analysis */
Data FRF_mult;
set LDH_mult;
if (CCl4 = 2.5 AND CHCl3 = 0) then output;
if (CCl4 = 2.5 AND CHCl3 = 0) then output;
if (CCl4 = 2.5 AND CHCl3 = 0) then output;
if (CCl4 = 2.5 AND CHCl3 = 5) then output;
if (CCl4 = 2.5 AND CHCl3 = 5) then output;
if (CCl4 = 2.5 AND CHCl3 = 5) then output;
if (CCl4 = 5 AND CHCl3 = 0) then output;
if (CCl4 = 5 AND CHCl3 = 0) then output;
if (CCl4 = 5 AND CHCl3 = 0) then output;
if (CCl4 = 5 AND CHCl3 = 5) then output;
if (CCl4 = 5 AND CHCl3 = 5) then output;
if (CCl4 = 5 AND CHCl3 = 5) then output;
run;
Proc print data = FRF_mult;
run;
/* Sphericity Test using PROC MIXED */
/* Sphericity test H0: Sphericity holds */
proc mixed data = FRF method = reml cl ic covtest;
class CCl4 CHCl3 Time Flask;
model Leakage = CCl4|CHCl3;
random CHCl3 CCl4*CHCl3 CHCl3*Time CCl4*CHCl3*Time /s;
repeated / subject = Flask(CCl4*CHCl3) type = un;
run;
proc mixed data = FRF method = reml cl ic covtest;
class CCl4 CHCl3 Time Flask;
model Leakage = CCl4|CHCl3;
random CHCl3 CCl4*CHCl3 CHCl3*Time CCl4*CHCl3*Time /s;
repeated / subject = Flask(CCl4*CHCl3) type = HF;
run;
/* Normality Q–Q plots */
ods graphics on;
proc mixed data = FRF plots = influenceestplot;
class CCl4 CHCl3 Time Flask;
model Leakage = CCl4 Time CCl4*Time/residual;
random CHCl3 CCl4*CHCl3 CHCl3*Time CCl4*CHCl3*Time;
repeated/subject = Flask(CCl4*CHCl3) type = cs r;
run;
ods graphics off;
/* Checking Covariance Structure */
proc corr data = FRF_mult cov;
var Time4 Time5 Time6;
run;
/* Fit the model by PROC MIXED and compare covariance structures */
proc mixed data = FRF method = reml cl ic covtest;
class CCl4 CHCl3 Time Flask;
model Leakage = CCl4 Time CCl4*Time/s;
random CHCl3 CCl4*CHCl3 CHCl3*Time CCl4*CHCl3*Time;
repeated/subject = Flask(CCl4*CHCl3) type = cs r;
run;
proc mixed data = FRF method = reml cl ic covtest;
class CCl4 CHCl3 Time Flask;
model Leakage = CCl4 Time CCl4*Time;
random CHCl3 CCl4*CHCl3 CHCl3*Time CCl4*CHCl3*Time/s;
repeated/subject = Flask(CCl4*CHCl3) type = arh(1) r;
lsmeans CCl4/pdiff cl adjust = tukey;
run;
proc mixed data = FRF method = reml cl ic covtest;
class CCl4 CHCl3 Time Flask;
model Leakage = CCl4 Time CCl4*Time;
random CHCl3 CCl4*CHCl3 CHCl3*Time CCl4*CHCl3*Time/s;
repeated/subject = Flask(CCl4*CHCl3) type = ar(1) r;
lsmeans CCl4/pdiff cl adjust = tukey;
run;
/* Scrapping Data for the Combined Model FB */
Data FB_univ;
set LDH_univ;
if (CHCl3 = 10) then output;
if (CHCl3 = 25) then output;
run;
/* Scrapping Data for the FA x FB Combined Model */
Data FAFB_univ;
set LDH_univ;
if (CCl4 = 2.5 and CHCl3 = 10) then output;
if (CCl4 = 2.5 and CHCl3 = 25) then output;
if (CCl4 = 5 and CHCl3 = 10) then output;
if (CCl4 = 5 and CHCl3 = 25) then output;
run;
Proc print data = FB_univ;
run;
/* Fitting the combined model FB (for narrow inferential scope) */
proc mixed data = FB_univ method = reml cl ic covtest;
class CCl4 CHCl3 Time Flask;
model Leakage = CCl4 CHCl3 CCl4*CHCl3 Time CCl4*Time CHCl3*Time CCl4*CHCl3*Time/s;
repeated/subject = Flask(CCl4*CHCl3) type = arh(1) r;
run;
/* Fitting the combined model FB (for broad inferential scope) */
proc mixed data = FB_univ method = reml cl ic covtest;
class CCl4 CHCl3 Time Flask;
model Leakage = CHCl3 Time CHCl3*Time/s;
random CCl4 CCl4*CHCl3 CCl4*Time CCl4*CHCl3*Time;
repeated/subject = Flask(CCl4*CHCl3) type = ar(1) r;
run;

References

  1. Fitzmaurice, G.M.; Davidian, M.; Verbeke, G.; Molenberghs, G. (Eds.) Longitudinal Data Analysis; Chapman and Hall/CRC Handbooks of Modern Statistical Methods; Chapman & Hall/CRC: Boca Raton, FL, USA, 2009. [Google Scholar]
  2. Pan, J.; Shang, J.A. Simultaneous Variable Selection Methodology for Linear Mixed Models. J. Stat. Comput. Simul. 2018, 88, 3323–3337. [Google Scholar] [CrossRef]
  3. Demidenko, E. Mixed Models: Theory and Applications with R, 2nd ed.; Wiley Series in Probability and Statistics; Wiley: Hoboken, NJ, USA, 2013. [Google Scholar]
  4. Verbeke, G.; Molenberghs, G. Linear Mixed Models for Longitudinal Data; Springer Series in Statistics; Springer: New York, NJ, USA; Berlin, Heidelberg, 2000. [Google Scholar]
  5. Verbeke, G. The Linear Mixed Model. A Critical Investigation in the Context of Longitudinal Data Analysis. Ph.D. Dissertation, Catholic University of Leuven, Faculty of Sciences, Department of Mathematics,, Leuven, Belgium, 1995. [Google Scholar]
  6. Verbeke, G.; Lesaffre, E. The Effect of Misspecifying the Random-Effects Distribution in Linear Mixed Models for Longitudinal Data. Comput. Stat. Data Anal. 1997, 23, 541–556. [Google Scholar] [CrossRef]
  7. Davis, C.S. Statistical Methods for the Analysis of Repeated Measurements; Springer Texts in Statistics; Springer: New York, NY, USA, 2002. [Google Scholar]
  8. Muller, K.E.; Edwards, L.J.; Simpson, S.L.; Taylor, D.J. Statistical Tests with Accurate Size and Power for Balanced Linear Mixed Models. Statist. Med. 2007, 26, 3639–3660. [Google Scholar] [CrossRef] [PubMed]
  9. Hickey, G.L.; Mokhles, M.M.; Chambers, D.J.; Kolamunnage-Dona, R. Statistical Primer: Performing Repeated-Measures Analysis. Interact. Cardiovasc. Thorac. Surg. 2018, 26, 539–544. [Google Scholar] [CrossRef] [PubMed]
  10. Matuschek, H.; Kliegl, R.; Vasishth, S.; Baayen, H.; Bates, D. Balancing Type I Error and Power in Linear Mixed Models. J. Mem. Lang. 2017, 94, 305–315. [Google Scholar] [CrossRef]
  11. Fitzmaurice, G.M.; Laird, N.M.; Ware, J.H. Applied Longitudinal Analysis, 2nd ed.; Wiley Series in Probability and Statistics; Wiley: Hoboken, NJ, USA, 2011. [Google Scholar]
  12. Chaka, L.; Njuho, P. Construction of a Linear Mixed Model with Each Factor Having Both Fixed and Random Levels: A Case of Split-Split-Plot Structure in a RCBD. Int. J. Agric. Stat. Sci. 2021, 17, 501–518. Available online: https://connectjournals.com/03899.2021.17.501 (accessed on 20 January 2022).
  13. Crowder, M.J.; Hand, D.J. Analysis of Repeated Measures, 1st ed.; Monographs on Statistics and Applied Probability; Chapman and Hall: London, UK; New York, NY, USA, 1990. [Google Scholar]
  14. Barnett, A.G.; Koper, N.; Dobson, A.J.; Schmiegelow, F.; Manseau, M. Using Information Criteria to Select the Correct Variance-Covariance Structure for Longitudinal Data in Ecology: Selecting the Correct Variance-Covariance. Methods Ecol. Evol. 2010, 1, 15–24. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, Z.; Goonewardene, L.A. The Use of MIXED Models in the Analysis of Animal Experiments with Repeated Measures Data. Can. J. Anim. Sci. 2004, 84, 1–11. [Google Scholar] [CrossRef]
  16. Patterson, H.D.; Thompson, R. Recovery of Inter-Block Information When Block Sizes Are Unequal. Biometrika 1971, 58, 545–554. [Google Scholar] [CrossRef]
  17. Diffey, S.M.; Smith, A.B.; Welsh, A.H.; Cullis, B.R. A New REML (Parameter Expanded) EM Algorithm for Linear Mixed Models. Aust. N. Z. J. Stat. 2017, 59, 433–448. [Google Scholar] [CrossRef] [Green Version]
  18. Dempster, A.P.; Laird, N.M.; Rubin, D.B. Maximum Likelihood from Incomplete Data via the EM Algorithm. J. R. Stat. Society. Ser. B (Methodol.) 1977, 39, 1–22. [Google Scholar]
  19. Gennings, C.; Chinchilli, V.M.; Carter, W.H. Response Surface Analysis with Correlated Data: A Nonlinear Model Approach. J. Am. Stat. Assoc. 1989, 84, 805–809. [Google Scholar] [CrossRef]
  20. Njuho, P.M.; Milliken, G.A. Analysis of Linear Models with One Factor Having Both Fixed and Random Levels. Commun. Stat. Theory Methods 2005, 34, 1979–1989. [Google Scholar] [CrossRef]
  21. Njuho, P.M.; Milliken, G.A. Analysis of Linear Models with Two Factors Having Both Fixed and Random Levels. Commun. Stat. Theory Methods 2009, 38, 2348–2365. [Google Scholar] [CrossRef]
  22. Kotchaporn, S.; Araveeporn, A. Modifications of Levene’s and O’Brien’s Tests for Testing the Homogeneity of Variance Based on Median and Trimmed Mean. Thail. Stat. 2018, 16, 106–128. [Google Scholar]
  23. Sullivan, L.M. Repeated Measures. Circulation 2008, 117, 1238–1243. [Google Scholar] [CrossRef] [Green Version]
  24. Armstrong, R.A. Recommendations for Analysis of Repeated-Measures Designs: Testing and Correcting for Sphericity and Use of MANOVA and Mixed Model Analysis. Ophthalmic Physiol. Opt. 2017, 37, 585–593. [Google Scholar] [CrossRef] [Green Version]
  25. Freund, R.J.; Wilson, W.J.; Mohr, D.L. Design of Experiments. In Statistical Methods; Elsevier: Amsterdam, The Netherlands, 2010; pp. 521–576. [Google Scholar] [CrossRef]
  26. Geisser, S.; Greenhouse, S.W. An Extension of Box’s Results on the Use of the F Distribution in Multivariate Analysis. Ann. Math. Statist. 1958, 29, 885–891. [Google Scholar] [CrossRef]
  27. Huynh, H.; Feldt, L.S. Estimation of the Box Correction for Degrees of Freedom from Sample Data in Randomized Block and Split-Plot Designs. J. Educ. Stat. 1976, 1, 69. [Google Scholar] [CrossRef]
  28. Box, G.E.P. Some Theorems on Quadratic Forms Applied in the Study of Analysis of Variance Problems, I. Effect of Inequality of Variance in the One-Way Classification. Ann. Math. Statist. 1954, 25, 290–302. [Google Scholar] [CrossRef]
  29. Verma, J.P. Repeated Measures Design for Empirical Researchers; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2015. [Google Scholar]
  30. Conover, W.J.; Guerrero-Serrano, A.J.; Tercero-Gómez, V.G. An Update on ‘a Comparative Study of Tests for Homogeneity of Variance’. J. Stat. Comput. Simul. 2018, 88, 1454–1469. [Google Scholar] [CrossRef]
  31. Ott, L.; Longnecker, M. An Introduction to Statistical Methods & Data Analysis, 7th ed.; Cengage Learning: Boston, MA, USA, 2016. [Google Scholar]
  32. Mauchly, J.W. Significance Test for Sphericity of a Normal N-Variate Distribution. Ann. Math. Statist. 1940, 11, 204–209. [Google Scholar] [CrossRef]
  33. Statistical Methods and Data Analytics. UCLA: Statistical Consulting Group. Available online: https://stats.oarc.ucla.edu/sas/seminars/sas-repeatedmeasures/ (accessed on 6 February 2022).
  34. Moskowitz, D.S.; Hershberger, S.L.; American Psychological Association (Eds.) Modeling Intraindividual Variability with Repeated Measures Data: Methods and Applications; Multivariate Applications Book Series; L. Erlbaum Associates: Mahwah, NJ, USA, 2002. [Google Scholar]
  35. Hocking, R.R. The Analysis of Linear Models; Brooks/Cole Pub. Co: Monterey, CA, USA, 1985. [Google Scholar]
  36. Harville, D.A. Maximum Likelihood Approaches to Variance Component Estimation and to Related Problems. J. Am. Stat. Assoc. 1977, 72, 320–338. [Google Scholar] [CrossRef]
  37. SAS Institute Inc. SAS/STAT® 14.3 User’s Guide; SAS Institute Inc.: Cary, NC, USA, 2017. [Google Scholar]
  38. Milliken, G.A.; Johnson, D.E. Analysis of Messy Data. 3: Analysis of Covariance; Chapman & Hall/CRC: Boca Raton, FL, USA, 2002. [Google Scholar]
  39. Akaike, H. A New Look at the Statistical Model Identification. IEEE Trans. Automat. Contr. 1974, 19, 716–723. [Google Scholar] [CrossRef]
  40. Kuehl, R.O. Design of Experiments: Statistical Principles of Research Design and Analysis, 2nd ed.; Brooks/Cole, Cengage Learning: Belmont, CA, USA, 2000. [Google Scholar]
  41. Wolfinger, R.D.; Chang, M. Comparing the SAS GLM and Mixed Procedures for Repeated Measures. In Proceedings of the Twentieth Annual SAS Users Groups Conference; SAS Institute Inc: Cary, NC, USA, 1995. [Google Scholar]
  42. Hamer, R.M. Mixed-Up Mixed Models: Things That Look Like They Should Work but Don’t, and Things That Look Like They Shouldn’t Work but Do. In Proceedings of the Twenty-Fifth Annual SAS® Users Group International Conference, Indianapolis, Indiana, 9–12 April 2000. [Google Scholar]
  43. McLean, R.A.; Sanders, W.L.; Stroup, W.W. A Unified Approach to Mixed Linear Models. Am. Stat. 1991, 45, 54–64. [Google Scholar] [CrossRef]
Figure 1. Process of constructing coefficients of variance components.
Figure 1. Process of constructing coefficients of variance components.
Stats 05 00027 g001
Figure 2. Q–Q plots for (a) FFF, (b) FRF, (c) RFF, and (d) RRF data sets.
Figure 2. Q–Q plots for (a) FFF, (b) FRF, (c) RFF, and (d) RRF data sets.
Stats 05 00027 g002
Table 1. Variance components for a three-factor repeated-measures model.
Table 1. Variance components for a three-factor repeated-measures model.
Source of
Variation
Factor
Effect
FRFRepComponents
of Variance
a b t n
i j l k
A α i 0btN θ α
B β j a1tN σ β 2
AB α β i j 11tN σ α β 2
EU γ k i j 11t1 σ γ α β 2
P τ l ab0N θ τ
PA α τ i l 0b0N θ α τ
PB β τ j l a11N σ β τ 2
PAB α β τ i j l 111N σ α β τ 2
Error ε l i j k 1111 σ ϵ 2
Table 2. Expected mean squares for a three-factor repeated-measures design.
Table 2. Expected mean squares for a three-factor repeated-measures design.
Source of
Variation
Sum of
Squares
Degrees of
Freedom
E(MS)
A SSA a 1 σ ϵ 2 + t σ γ α β 2 + t n σ α β 2 + n σ α β τ 2 + b t n θ α
B SSB b 1 σ ϵ 2 + n σ α β τ 2 + a n σ β τ 2 + t σ γ α β 2 + t n σ α β 2 + a n t σ β 2
A × B SS A × B a 1 b 1 σ ϵ 2 + t σ γ α β 2 + n σ α β τ 2 + t n σ α β 2
Unit A × B SSU A × B a b n 1 σ ϵ 2 + t σ γ α β 2
P e r i o d SSP t 1 σ ϵ 2 + n σ α β τ 2 + a n σ β τ 2 + a b n θ τ
P e r i o d × A SSP × A a 1 t 1 σ ϵ 2 + n σ α β τ 2 + b n θ α τ
P e r i o d × B SSP × A b 1 t 1 σ ϵ 2 + n σ α β τ 2 + a n σ β τ 2
P e r i o d × A × B SSP × A × B a 1 b 1 t 1 σ ϵ 2 + n σ α β τ 2
Residual SSE a b t 1 n 1 σ ϵ 2
Total SST a b t n 1
Note that: θ α = 1 a 1 i = 1 a α i 2 ; θ α τ = 1 a 1 t 1 i = 1 a l = 1 t α τ i l 2 ; and θ τ = 1 t 1 l = 1 t τ l 2 .
Table 3. Data layout for a three-factor repeated-measures experiment.
Table 3. Data layout for a three-factor repeated-measures experiment.
FactorsTime Period
CCl4CHCl3Flask12t
111 y 1111 y 1112 y 111 t
n y 11 n 1 y 11 n 2 y 11 n t
21 y 1211 y 1212 y 121 t
n y 12 n 1 y 12 n 2 y 12 n t
b 1 y 1 b 11 y 1 b 12 y 1 b 1 t
n y 1 b n 1 y 1 b n 2 y 1 b n t
211 y 2111 y 2112 y 211 t
n y 21 n 1 y 21 n 2 y 21 n t
21 y 2211 y 2212 y 221 t
n y 22 n 1 y 22 n 2 y 22 n t
b 1 y 2 b 11 y 2 b 12 y 2 b 1 t
n y 2 b n 1 y 2 b n 2 y 2 b n t
a 11 y a 111 y a 112 y a 11 t
n y a 1 n 1 y a 1 n 2 y a 1 n t
21 y a 211 y a 212 y a 21 t
n y a 2 n 1 y a 2 n 2 y a 2 n t
b 1 y a b 11 y a b 12 y a b 1 t
Table 4. Akaike’s Information Criteria (AIC) for the covariance structures in partitioned models.
Table 4. Akaike’s Information Criteria (AIC) for the covariance structures in partitioned models.
Model FFFModel FRF
Covariance StructureNumber of ParametersAICCovariance StructureNumber of ParametersAIC
CS2−282.6CS4−361.3
AR(1)2−443.1AR(1)3−492.4
ARH(1)10−525.6ARH(1)13−506.3
CSH10−356.7
Model RFFModel RRF
Covariance StructureNumber of ParametersAICCovariance StructureNumber of ParametersAIC
CS4−450.1CS4−665.8
AR(1)4−522.4AR(1)4−782.4
ARH(1)13−509.3ARH(1)13−785.3
AR(1): First-order autoregressive; ARH(1): Heterogeneous first-order autoregressive; CS: Compound symmetry; VC: Variance components; CSH: Heterogeneous compound symmetry.
Table 5. Type III tests of fixed effects of FFF, FRF, FFR, and FRR models.
Table 5. Type III tests of fixed effects of FFF, FRF, FFR, and FRR models.
ModelEffectNumerator
Degrees of
Freedom
Denominator
Degrees of
Freedom
Fp-Value
FFFCCl41120.240.6363
CHCl311217.330.0013 **
C C l 4 × CHCl32121.350.2678
Time22495.99<0.0001 **
T i m e × CCl42242.660.0908
T i m e × CHCl322420.91<0.0001 **
T i m e × C C l 4 × CHCl32241.260.3023
FRFCHCl3110.50.6079
Time222.170.3156
T i m e × CCl4222.830.2610
RFFCHCl3117.070.2290
Time225.420.1557
T i m e × CHCl3223.260.2348
RRFTime221.010.4978
‘**’ significant at α = 0.05 .
Table 6. Covariance parameter estimates of FRF, RFF, and RRF, models based on ARH(1) covariance structure.
Table 6. Covariance parameter estimates of FRF, RFF, and RRF, models based on ARH(1) covariance structure.
ModelCovariance
Parameter
EstimateStandard ErrorProportion of Variation Accounted for
FRFCHCl3000
C C l 4 × CHCl30.042490.0447022.5
T i m e × CHCl3000
T i m e × C C l 4 × CHCl30.0006210.0007590.3
RFFCCl40.0006070.0030560.05
C C l 4 × CHCl3000
T i m e × CCl40.0001490.0008090.02
T i m e × C C l 4 × CHCl30.0008060.0009840.08
RRFCCl40.00110.004012.6
CHCl30.003060.006777.1
C C l 4 × CHCl30.002320.003965.4
T i m e × CCl40.0000020.000020.0
T i m e × CHCl30.0000650.0000960.1
T i m e × C C l 4 × CHCl30.0000130.0000250.0
Table 7. Fixed-effects F-tests for the combined models in narrow inference scope.
Table 7. Fixed-effects F-tests for the combined models in narrow inference scope.
Type III Tests of Fixed Effects in Combined Models
ModelAIC [CS]EffectNum DFDen DFFPr > F
FA−127.6 [ARH(1)]A1243.360.0794
B3245.810.0039 **
A × B 3249.100.0003 **
C24812.330.0001 **
A × C 2482.520.0908
B × C 6482.330.0468 **
A × B × C 6480.760.6027
FB−127.6 [AR(1)]A3240.530.6631
B12416.30.0005 **
A × B 3240.220.8801
C24824.82<0.0001 **
A × C 6481.090.3795
B × C 2487.820.0012 **
A × B × C 6480.500.8064
FA × FB−69.1 [AR(1)]A1120.140.7106
B11210.630.0068 **
A × B 1120.830.3807
C22419.75<0.0001 **
A × C 2240.580.5667
B × C 2244.650.0196 **
A × B × C 2240.260.7725
‘**’ significant at α = 0.05 .
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chaka, L.; Njuho, P. Repeated-Measures Analysis in the Context of Heteroscedastic Error Terms with Factors Having Both Fixed and Random Levels. Stats 2022, 5, 458-476. https://doi.org/10.3390/stats5020027

AMA Style

Chaka L, Njuho P. Repeated-Measures Analysis in the Context of Heteroscedastic Error Terms with Factors Having Both Fixed and Random Levels. Stats. 2022; 5(2):458-476. https://doi.org/10.3390/stats5020027

Chicago/Turabian Style

Chaka, Lyson, and Peter Njuho. 2022. "Repeated-Measures Analysis in the Context of Heteroscedastic Error Terms with Factors Having Both Fixed and Random Levels" Stats 5, no. 2: 458-476. https://doi.org/10.3390/stats5020027

Article Metrics

Back to TopTop