Optimal Repeated Measurements for Two Treatment Designs with Dependent Observations : The Case of Compound Symmetry

In this paper, we construct optimal repeated measurement designs of two treatments for estimating direct effects, and we examine the case of compound symmetry dependency. We present the model and the design that minimizes the variance of the estimated difference of the two treatments. The optimal designs with dependent observations in a compound symmetry model are the same as in the case of independent observations.


Introduction
In repeated measurement designs, a sequence of treatments is applied to each experimental unit (e.u.).In particular, one treatment is applied in each period.For example, for two treatments, A and B, and three periods, a possible sequence is ABA, which means that the treatments A, B, and A are respectively applied at the beginning of each of the three periods.The direct effect of a treatment is the effect of the treatment which is applied in the period that is examined.The residual effect is the effect of the treatment which is applied in the period preceding the period that is examined.In the case of two treatments, A and B, the direct τ A and τ B can be estimated.In every period, a treatment is applied, so either τ A or τ B is estimable.In this paper, the parameter of interest is the difference of direct effects Most researchers who have investigated repeated measurement designs, such as [1][2][3][4][5][6], have been occupied with universally optimal designs where the observations are independent.However, researchers have also shown interest in designs with dependent observations, as in the cases of [7][8][9][10][11].
The model we use in this paper, and which is presented below, was first introduced by Hedayat and Afsarinejad [12,13].In previous research [14,15] using this model, the author of this article studied two treatment designs under the assumption that consecutive observations were independent.Building on that previous work, in the present article the author examines the case of compound symmetry dependency.The aim is to find a design that corresponds to a minimum variance estimator.

The Model
A compound symmetry model has the following characteristics: (i) For each sequence, the variance matrix is of the form Σ m = aI m + bJ m , where I m is the unit m × m matrix, and J m is the m × m matrix where all elements are equal to 1 (m is the number of periods).(ii) The observations corresponding to different treatment sequences (different e.u.) are independent, and the number of sequences is 2 m .
The goal is to find the design that corresponds to the minimum variance estimator.I show that, in this case, the optimal design regarding the direct effect is the same as in the model of independent observations, and only the variance of the estimator is different.
The model is [12]: j corresponds to the j-th period, j = 1, 2, . . ., m; i corresponds to the i-th sequence, i = 0, 1, . . . 2 m − 1; k corresponds to the unit k = 1, 2, . . ., n; τ A , τ B : are direct effects of treatments A and B; π j : is the effect of the j-th period; δ A , δ B : are the residual effects of A and B; γ i : is the effect of the i-th sequence; and is the effect of the k-th e.u.(subject effect), which is a random variable, independent of the error e i jk .
The errors e ijk are assumed to be independent.However, the quantities ζ k + e ijk are independent only between sequences and not within sequences.
The overparameterized model vector form of the above model is written as: where are 1 × mn vectors; the direct effect vector is 1 if the treatment is A, and zero if it is B. For example, for the sequence ABB . . .
Keeping only the linear independent vector [16], the model ( 2) is transformed to where q = 2 m − 1.In a vector form: where Y is (mn) × 1, the design matrix X is (mn) × s, b is s × 1, e is (mn) × 1, and s is the number of unknown parameters.If we are interested only in some and not in all of the parameters, then we write b = b 1 b 2 , where b 1 is the r parameters of interest, and b 2 is the s-r remaining parameters.We assume only one parameter of interest for the difference of the direct effects, τ = τ A − τ B , which can be considered as the direct effect of A in the case of τ B = 0.In order to guarantee the estimability of the model, we postulate the restrictions τ B = 0, π m = 0, γ 2 m −1 = 0.
The matrix X 1 corresponds to the coefficients of τ, and the matrix X 2 corresponds to the coefficients of the rest of the non-random variables.Let us assume matrix, P the projection matrix of X 2 , P = X 2 (X 2 T X 2 ) −1 X 2 T and Σ are the (mn) × (mn) variance matrix of the observations.
From the ordinal least-squares equations, we derive the following relation for the estimation of the main effect τ:

The Case of Compound Symmetry
The observations are dependent within sequences with variance matrix Σ m .The observations from different sequences are independent, therefore: where q = 2 m − 1 and In order to obtain a sequence enumeration, the binary enumeration system was used, with 0 corresponding to A, and 1 to B. Thus, we obtained the enumerations 0, 1, . . ., 2 m − 1.For example, if we have five periods and the sequence BABBA, then this is the 13th sequence, since For two periods, we have four sequences, that is, AA ↔ 0, BA ↔ 1, AB ↔ 2, BB ↔ 3 .For three periods (two treatments) we have eight sequences: where u i i = 0, 2, 3, 4, 5, 6, 7 is the number of units that received the i-th sequence of treatments.The sequences that we obtain by substituting A for B and vice versa are called dual or reversal designs.
Proposition 1.For a repeated measurement design with m periods, n experimental units, and a variance matrix Σ that consists of n diagonal block matrices of the form where In other words, P is the matrix of the orthogonal projection to R( X 2 ).
For example, for three periods (m = 3), we have the matrices: For the linear space R( X 2j ) and for any sequence (m observations) R( X 2j ) = R(X 2j ), we observe the following: (i) The matrix (aI m + bJ m ) is positive definite, so the matrix (aI m + bJ m ) − 1 2 is also positive definite, and we conclude that: The coefficients of the general mean are 1, so 1 m ∈ R(X 2j ) and.
(iv) If P m is the matrix of the orthogonal projection to the linear space R( X 2j ), then P m j = P m j , where P mj = X 2 j (X 2 j T X 2 j ) −1 X 2 j T is the matrix of the orthogonal projection to R(X 2j ) and P mj • 1 m = 1 m ⇒ P mj • J m = J m .From the above, we conclude that: Corollary 1.The designs that result in the estimators with the minimum variance, i.e., minvar( τ) are exactly the optimal designs of the model with independent observations.In this case, the variance var( τ) is multiplied by α: is the variance of the optimal designs in the model with independent observations).
Proof.From the previous proof, we conclude that the variance of the estimator of the direct effect, which is given by Formula (3), equals to var(τ) = σ 2 a • (Q * ) −1