Analysis of Variance in an Unbalanced Two-Way Mixed Effect Interactive Model ()
Received 9 February 2016; accepted 23 April 2016; published 26 April 2016
1. Introduction
The problem with unbalanced fixed effect interactive model is associated with the appropriate F-test for testing for the main effects when interactions are present. The paper by [1] worked on application of mixed-effects model for exposure assessment by re-analyzing three data sets from published surveys with repeated exposure measurements. The relative contributions of particular characteristics affecting exposure levels were assessed as in a multiple regression model, while controlling for the correlation between repeated measurements.
In [2] , they studied a mixed-effects nonlinear regression for unbalanced repeated measures by estimating and comparing the parameters of a generalized mixed-effects nonlinear regression model. The results are applied to in vitro data on the water transport kinetics of hemodialyzers used in the treatment of patients with chronic renal failure.
Similarly, [3] developed a method for deriving exact tests for variance components in some unbalanced mixed linear models. The derivation was based on a new kind of preliminary orthogonal transformation and a subsequent resampling procedure. The resulting tests are based on mutually independent sums of squares which, under null hypothesis, are distributed as scalar multiples of chi-square variates.
Also [4] derived exact tests for testing hypotheses concerning the variance components of the main effects in an unbalanced random two-way crossed classification with interaction model. The tests are based on four sums of squares that are distributed independently as scalar multiple chi-square variates. These sums of squares can be used to find an exact test concerning the interaction variance components.
However, [5] considered the Two-Way ANOVA model with unequal cell frequencies without the assumption of equal error variances. They used generalized approach to finding p-values, classical F-tests for no interaction effects and equal main effects are extended under heteroscedasticity. The generalized F-test they developed in their article can be utilized in significance testing or in fixed level testing under the Neyman-Pearson theory. The problem in their work is that, the assumption of ANOVA was violated.
Analysis of variance is straightforward when an experimental design is balanced, but unequal cell sizes affect the computation of means, hypotheses tested and F-statistics [6] .
Several solutions have been proposed for the analysis of unbalanced data. Solutions have focused on forcing the unbalanced data to be balanced. Suggestions include imputing cell means as additional data points into the smaller cells.
Given the model
(1)
is the kth observation in ijth cell;
is the overall mean effects;
is the average effects of factor A;
is the average effects of factor B;
is the effects of the interaction between factor A and factor B;
is a random error components and;
is the number of observations per cell.
To derive the expected mean squares for Equation (1), [7] derived the expected mean square for factor A when factor A is fixed and factor B is random as
(2)
According to him, in mixed models, the expected values of the sums of squares contain functions of the fixed effects that cannot be eliminated by considering linear combinations of the sums of squares. He suggested two obvious ways of overcoming the difficulties associated with unbalanced mixed effect data. The first is to ignore the fixed effects and eliminate them from the model. What remains is a random model for which the F-test can be determined. The second possibility is to assume the fixed effects as random and therefore assume the entire model as random effect models. These suggestions are in fact unsatisfactory.
2. Expected Mean Squares
2.1. Two-Way Unbalanced Random Effect Model
From Equation (1) above, [8] derived the expected mean squares for unbalanced two-way interactive random model.
They derive the expected mean squares for Equation (1) as shown in Table 1.
Where
(3)
(4)
Table 1. ANOVA table for unbalanced two-way interactive random model.
(5)
(6)
(7)
(8)
and
(9)
(10)
(11)
From Table 1, they found a linear combination of the mean squares with the expected mean squares and derived an expression for testing for
(12)
With the corresponding F-ratios as
(13)
(14)
(15)
where are degrees of for factor A, factor B, the interaction between factor A and factor B and error components respectively.
And
(16)
(17)
(18)
(19)
The sums of squares for factor A, factor B, the interaction between factor A and factor and the error terms are given by
(20)
(21)
(22)
(23)
The expression for testing for the presence of interaction from Table 1 is
(24)
2.2. Two-Way Unbalanced Mixed Effect Model
From Equation (1), if factor A is fixed and factor B is random
(25)
(26)
and
(27)
Similarly, if factor A is random and factor B is fixed
(28)
(29)
and
(30)
Using Brute Force Method, the expected mean squares of Equation (1) when factor A is fixed and factor B is random the expected mean square are shown in ANOVA Table 2.
From Table 2, if we interested to test for the expression there are no obvious denominator for testing for the factor A.
However if we can obtain the expression
(31)
we would have the F-test as
(32)
where and are the numerator and denominator degrees of freedom respectively
(33)
(34)
(35)
Using Welch Satterthwaite Equation is the degree of freedom for the denominator and is given by
(36)
can be shown to be
(37)
where
Table 2. ANOVA table two way unbalanced mixed interactive model when factor A is fixed and factor B is random.
Statement 1: Equation (37) is an unbiased estimate of Equation (31).
Proof:
We take expectation on Equation (33) to have
But
Similarly, if we are interested to test for factor B we have
there would be no obvious denominator to test for the above hypothesis. However, if we can obtain the expression
(38)
We would have the F-test as
(39)
where and are the numerator and denominator degrees of freedom respectively.
is the degree of freedom for the numerator and is the degree of freedom for the denominator and is given by
(40)
where and are the degrees of freedom for the error components and the interactions respectively.
can be shown to be
(41)
Equation (41) is also an unbiased estimate of Equation (38).
Similarly, if factor A is random, and factor B is fixed, the expected mean square are shown in the ANOVA Table 3.
Where
(42)
(43)
Table 3. ANOVA table two way unbalanced mixed interactive model when factor A is random and factor B is fixed.
(44)
When factor A is random, the hypothesis is given by
and there would be no obvious denominator to test for the above hypothesis. If we can obtain the expression
(45)
the F-test can be shown to be
(46)
(47)
where and are the numerator and denominator degrees of freedom respectively.
(48)
where and are the degrees of freedom for the error components and the interactions respectively.
Equation (47) is also an unbiased estimate of Equation (45).
Similarly, when factor B is fixed, the hypothesis is given by
with no obvious denominator to test for the hypothesis. However, if we can obtain the expression
(49)
the F-test can be shown to be
(50)
(51)
(52)
where and are the numerator and denominator degrees of freedom respectively.
Similarly Equation (51) is also an unbiased estimate of Equation (49).
Finally, the hypothesis for testing for the presence of the interaction is given by
(53)
where and are the numerator and denominator degrees of freedom respectively.
3. Conclusions
Equation (2) contains the functions of the fixed effect which is. If we ignore the fixed ef-
fects and eliminate them from the model, what remains is a random model for which the F-test can be determined. The second possibility is to assume the fixed effects as random and therefore assume the entire model as random effect models. This is completely unreasonable.
From Table 2 and Table 3, when one factor is fixed, we equate the functions of the fixed effect to zero and obtain an expression to determine the denominator for the F-ratio when the hypothesis is specified. Similarly, when the other factor is random, we equate the functions of the random effect to zero and obtain an expression to determine the denominator for the F-ratio when the hypothesis is specified.
To test for the interaction effect for the mixed effect model, we have
This does not involve obtaining any expression and the degrees of freedom for both the numerator and denominator are integer valued whereas the denominator degrees of freedom for the testing for the main effects are non integer valued.
Instead of assuming both effects to be fixed or both effects to be random to enable researchers on mixed effect unbalanced interactive model analyze their data, we highly recommend our method.
This paper is limited to only an unbalanced two-way mixed effect interactive model and cannot be applied to random or fixed effect model.
4. Illustrative Example
Synthetic growth hormone was administered at a clinical research center to growth hormone deficient 18 short children who had not yet reached puberty. The investigator was interested in the effects of a child’s gender (factor A) and bone development (factor B) on the rate of growth induced by hormone administration. A child’s bone development was classified into one of the three categories: severely depressed, moderately depressed and mildly depressed. Three children were randomly selected for each gender-bone development group. The response variable (Y) of interest was the difference between the growth rate during hormone treatment and the normal growth rate prior to the treatment, expressed in centimeters per month. Four of the 18 children were unable to complete the study leading to unequal treatment sample sizes shown below.
Since factor A is random and factor B is assumed to fixed, we shall make use of the information in Table 4.
Source: Netal et al. (1996) Applied Linear Statistical Models [9] .
Our hypothesis for factor A shall be
Using Equations (20), (22) and (23) we have
Similarly, using Equations (43) and (44)
From Equations (47) and (48)
Our conclusion is that we do not reject the null hypothesis.
Similarly, our hypothesis for factor B shall be
Using Equation (21)
From Equations (51) and (52)
Our conclusion is that we do not reject the null hypothesis.
Finally, to test for the interaction we have
Our conclusion is therefore that interaction is present.