So far we were able to solve the repeated measures design only within the frameworks of a repeated measures t-test, or by using the Microsoft Excel spreadsheet for the double classification analysis of variance design. The multiple regression solution of the repeated measures design approach is quite exacting. A separate coding vector has to be constructed for all subjects involved in the analysis save the last one. These coding vectors are not orthogonal. They are filled with 0s save few 1s marking the location of the subject within the different conditions of the experiment. To compare these different approaches to the repeated measures design, let us first review its idealized model.
This model makes an assumption that, initially, the means and variances of all components of the model are the same.
|
|
Y0 |
Y1 |
Y |
|
Allen |
1 |
|
1 |
|
Becky |
2 |
|
2 |
|
Cathy |
3 |
|
3 |
|
Debra |
|
1 |
1 |
|
Edgar |
|
2 |
2 |
|
Francis |
|
3 |
3 |
|
M |
2 |
2 |
2 |
|
σ2 |
.67 |
.67 |
.67 |
Following some type of intervention,
|
|
Y0 |
Y1 |
Y |
|
Allen |
1 |
1+2=3 |
1 |
|
Becky |
2 |
2+0=2 |
2 |
|
Cathy |
3 |
3+1=4 |
3 |
|
Allen |
|
|
3 |
|
Becky |
|
|
2 |
|
Cathy |
|
|
4 |
|
M |
2 |
3 |
2.5 |
|
σ2 |
.67 |
.67 |
.92 |
the total variances of the data will change. This change of the total variance component is due to the changed variance between the column
|
|
M |
|
Y0 |
2 |
|
Y1 |
3 |
|
M |
2.5 |
|
σ2 |
.25 |
and row
|
|
Y0 |
Y1 |
Y0 + Y1 |
M |
|
Allen |
1 |
3 |
4 |
2 |
|
Becky |
2 |
2 |
4 |
2 |
|
Cathy |
3 |
4 |
7 |
3.5 |
|
M |
2 |
3 |
5 |
2.5 |
|
σ2 |
.67 |
.67 |
2 |
.50 |
means. The residual variance component is due to the mean differences between the changed row values. For the example,
|
|
Y0 |
Y1 |
Y0 - Y1 |
M |
|
Allen |
1 |
3 |
-2 |
-1 |
|
Becky |
2 |
2 |
0 |
0 |
|
Cathy |
3 |
4 |
-1 |
-.5 |
|
M |
2 |
3 |
-1 |
-.5 |
|
σ2 |
.67 |
.67 |
.67 |
.17 |
For the example, .92 = .25 + .50 + .17.
The column component of variance can be obtained by using the coded regression analysis as
The column variance component equals .25.
To obtain the row variance component, one has to use the multiple regression analysis where the group membership is coded by codes specifying which scores in different conditions of the experiment were obtained from the same subject, as shown below.
Note that the code for the last subject which would have been [0 0 1 0 0 1] have not been included, as it is redundant. You might wish to experiment with including a code for the last subject only to observe that inclusion of this code does not make any difference, since the information about subject-group memberships is fully conveyed by k-1 coding vectors.
Since the predictor variables are not orthogonal, the variance contributions of the predictor variables are not additive. Select the variable X as the initial variable and partial out the redundant variance components of the X1 and X2 variables, as
Multiple regression analysis
partitions variance of the variable Y (.92) into the row variance component (.13 + .37 = .50), into the column variance component (.25) and the residual component of variance (.17). Thus, .92 = .25 + .50 + .17, identical to values obtained by the conceptual analysis of the same data.
For the example, the variance components obtained by the coded multiple regression can be summarized as
|
Source of Variance |
Degrees of Freedom |
Variance Components |
Unbiased Variance Components |
F |
Probability |
|
Columns |
1 |
.25 |
.25 |
3.0 |
.1440 |
|
Rows |
2 |
.50 |
.25 |
3.0 |
|
|
Residual |
2 |
.17 |
.08 |
|
|
|
Total |
5 |
.92 |
|
|
|
For the example, the variance components obtained by the coded multiple regression can be summarized as
|
Source of Variance |
Degrees of Freedom |
Standard Variance Components |
Unbiased Variance Components |
F |
Probability |
|
Columns |
1 |
.27 |
.27 |
3.0 |
.1440 |
|
Rows |
2 |
.55 |
.27 |
3.0 |
|
|
Residual |
2 |
.18 |
.09 |
|
|
|
Total |
5 |
1.00 |
|
|
|