^{1}

^{2}

^{*}

Background
: Repeatability is a statement on the magnitude of measurement error. When biomarkers are used for disease diagnoses, they should be measured accurately. **Objectives**: We derive an index of repeatability based on the ratio of two variance components. Estimation of the index is derived from the one-way Analysis of Variance table based on the one-way random effects model. We estimate the large sample variance of the estimator and assess its adequacy using bootstrap methods. An important requirement for valid estimation of repeatability is the availability of multiple observations on each subject taken by the same rater and under the same conditions. **Methods**: We use the delta method to derive the large sample variance of the estimate of repeatability index. The question related to the number of required repeats per subjects is answered by two methods. In first methods we estimate the number of repeats that minimizes the variance of the estimated repeatability index, and the second determine the number of repeats needed under cost-constraints. **Results and Novel Contribution**: The situation when the measurements do not follow Gaussian distribution will be dealt with. It is shown that the required sample size is quite sensitive to the relative cost. We illustrate the methodologies on the Serum Alanine-aminotransferase (ALT) available from hospital registry data for samples of males and females. Repeatability is higher among females in comparison to males.

Repeatability and reproducibility are ways of measuring precision, particularly in the fields of biochemistry, radiology, and medical diagnoses. In general, scientists perform the same experiment several times in order to confirm their findings. These findings may show variations. In the context of an experiment, repeatability measures the variation in measurements taken by a single instrument or person under the same conditions, while reproducibility measures whether an entire study or an experiment can be reproduced. There has been confusion in the literature about the way that repeatability and reproducibility are quantified. Both concepts were often reported as either standard deviations or coefficient of variations.

The main focus of this paper is on the concept of repeatability, which was first introduced by Bland and Altman [

What’s known as “the repeatability” is in fact a measurement of precision, which denotes the absolute difference between a pair of repeated test results. We note that when we have more than two readings per subject the idea of pairing produces several repeatability coefficients and the concept becomes unclear.

Repeatability is also known as test-retest reliability indicating the closeness of the agreement between the results of successive measurements of the same measurand carried out under the same conditions of measurement. A less-than-perfect test-retest reliability causes test-retest variability. Such variability can be caused by, for example, intra-individual variability and intra-observer variability. A measurement may be said to be repeatable when this variation is smaller than a pre-determined acceptance criterion. A complete account on the reliability literature can be found in Shoukri [

One of the most important applications of the concept of repeatability is in the construction of the normal range or reference range in clinical medicine, which relies on the availability of a large sample of healthy individuals. Research has shown that the distribution of these measurements is affected by two main sources of variations: the between subjects and the within subject-components of variations.

This paper has three-fold objectives: Firstly; we define a proposed index of repeatability, as the ratio of the within-subjects’ variations to the between subject variation. The within subject variation is expected to be quite small relative to the between subjects-variations. To formalize the presentation, we assume that a single measurement y i from subject i = 1 , 2 , ⋯ , k is written as:

y i = μ + s i + e i (1)

Hence s i represents the sources of between-subjects biological variation, and e i represents sources of within subject variations, while μ denotes the population mean. Note that the assumption of additivity of components is made to simplify the presentation. However, a multiplicative model may be made additive under the logarithmic transformation. Following Harries and De Mets [

The salient point is that θ cannot be estimated unless we have at least two repeated measurements on any subject in the study.

In Section 2 we specify the model generating the observations and discuss a general method of estimating RIP from a sample of k subjects when there is an opportunity to have n repeated samples per subject. In Section 3, we provide two alternatives for the sampling strategies. The first, we assume that the investigator has decided to acquire on total number of measurements N = k n , and the question becomes; what is the best split between ( k , n ) , that maximizes the accuracy of estimating RIP?

One of the biggest obstacles in clinical studies is the cost constraints. Therefore, the second strategy is to find the optimal split of N = k n , so that IRP is estimated with maximum precision under cost restrictions (constrained optimization). The third objective of the study is to address the issue of estimating RIP when the assumption of the Gaussian distribution of observation is not tenable.

We assume that for subject i, n replicates of the same variable of interest y i j are taken by the same instrument at the same time, so that

y i j = μ + s i + e i j (2)

i = 1 , 2 , ⋯ , k

k is the number of subjects

j = 1 , 2 , ⋯ , n

n is the number of replications per subject

We further assume that the components of the model described by (2) are such that:

s i ~ N ( 0 , σ s 2 ) , are independently distributed random variables measuring the subjects effect, are independently distributed of the within subjects variation denoted by e i j ~ N ( 0 , σ e 2 ) .

Under the additivity assumption of the model components, we have:

var ( y i j ) = σ s 2 + σ e 2 = σ s 2 ( 1 + σ e 2 / σ s 2 ) = σ s 2 ( 1 + θ )

The parameter θ = σ e 2 / σ s 2 is the target parameter of interest, named “Repeatability Index Parameter” (RIP). The components of variation of the model set-up can be estimated using the well-known one-way Analysis of Variance (ANOVA) with random effects (

S.O.V = Source of variation, DF = Degrees of freedom associated with the corresponding sum of squares, S.O.S = Corrected sums of squares, MS = Mean square error = S.O.S/DF, EMS = Expected mean square. The sample statistics needed for the ANOVA computations based on the available observations are given as:

S.O.V | DF | S.O.S | M.S | EMS |
---|---|---|---|---|

Between subjects | k − 1 | SSB | MSB | σ e 2 + n σ s 2 |

Within subjects | N − k | SSW | MSW | σ e 2 |

Total | N − 1 |

y ¯ .. = 1 n k ∑ i ∑ j y i j

S S W = ∑ i ∑ j ( y i j − y ¯ i . ) 2

y ¯ i . = 1 n ∑ j = 1 n y i j

S S B = ∑ i = 1 k n ( y ¯ i − y ¯ .. ) 2

The moment estimator of the parameter θ , and hence the maximum likelihood estimator (under balanced design) is given by:

θ ^ = n [ M S W M S B − M S W ] (3)

The parameter estimator θ ^ given in Equation (3) is a nonlinear function of the sample statistics, and therefore an exact expression for its variance is not available. We use the delta method (Kendall and Stuart, 1989) [

var ( θ ^ ) = var ( M S B ) ( ∂ ˙ θ ∂ M S B ) 2 + 2 cov ( M S B , M S W ) ( ∂ ˙ θ ∂ M S B ) ( ∂ ˙ θ ∂ M S W ) + var ( M S W ) ( ∂ ˙ θ ∂ M S W ) 2 (4)

Substituting the required quantities in (4) and simplifying we get the first order approximation of the variance of θ ^ as:

v = var ( θ ^ ) = 2 θ 2 ( n + θ ) 2 N ( n − 1 ) ( 1 + θ ) 8 (5)

An ( 1 − α ) 100 % approximate Wald’s confidence interval on θ may be constructed as:

θ ^ ± z v (6)

where z in Equation (6) is the ( 1 − α / 2 ) 100 % cut-off point from the standard normal table.

Our first approach to estimate the optimal number of replications is to assume that N = k n is fixed a priori, and one needs to determine the number of replicates that minimizes the variance v, as given in (5).

Minimizing v with respect to the number of replications per subjects and solving for n we get:

∂ v ∂ n = 0 ⇒ n = 2 + θ (7)

∂ 2 v ∂ n 2 = 1 2 ( 1 + θ ) 2 > 0

This means that var ( θ ^ ) is minimized (i.e. maximizing precision) when at least 2 repeats are attained from each subject as shown in Equation (7). When θ = 0 (no within subject-variation) then n = 2 precisely.

We may also estimate the number of repeats for fixed width confidence interval as follows:

Suppose that we have decided on the number of subjects k. The question now is how many repeats per-subject are needed to estimate θ with 95% confidence such that the width of the confidence interval has a maximum given length w.

Since the length of the Wald’s confidence interval is given as:

w = 2 ( 1.96 ) var (θ^)

w 2 = 4 ( 1.96 ) 2 ⋅ 2 θ 2 ( n + θ ) 2 ( 1 + θ ) 8 k n ( n − 1 )

Let A = k w 2 ( 1 + θ ) 8 8 ( 1.96 ) 2 θ 2 > 1

Solving for n we have:

n = ( A + 2 θ ) + [ A ( 4 θ 2 + 4 θ + A ) ] 1 / 2 2 ( A − 1 ) (8)

This closed form expression is quite simple, and the computation of n from Equation (8), is straight forward. Substituting θ = 0.25 , k = 100 , and w = 1 in (8), then n = 156 .

It is an extremely expensive, and in some circumstances, it is a difficult task to obtain repeated samples from each subject. Some of these difficulties are related to cost and time (which may be translated into cost). Clearly too small a sample may lead to a study that produces many false negatives, too large a sample may result in many false positives and additional cost. Thus, a critical decision in constructing accurate estimate of normal range is to balance the cost of recruiting healthy normal with the need to obtain accurate estimate of RIP. In this section we shall address the issue of obtaining the combination ( n , k ) that minimizes the variance of θ ^ subject to cost constraints. The sampling cost depends primarily on the size of the sample, and includes the data collection costs, subjects recruiting costs, management and technicians’ costs. On the other hand, overhead costs remain fixed regardless of the sample size. The total cost is assumed this additive formula

T = t 0 + k t 1 + n k t 2 (9)

In Equation (9) t 0 is the fixed cost, t 1 is the cost of recruiting a healthy subject, and t 2 is the cost of taking a single measurement. Denoting the variance of θ ^ by V, the main objective is to determine the number of repeated measurements that minimize the variance of θ ^ subject to cost constraints T. In terms of language of optimization, we construct the objective function

Q = V + λ ( T − t 0 − k t 1 − n k t 2 ) (10)

The parameter λ in Equation (10) is the Lagrange-multiplier. The necessary conditions for minimization of Q are:

∂ Q ∂ n = 0 , ∂ Q ∂ k = 0 , and ∂ Q ∂ λ = 0.

Differentiating Q with respect to n, k, and λ and equating to zero we get:

n 3 − n 2 ( 2 + θ ) − n R ( 1 + 2 θ ) + θ R = 0 (11)

Note that from Equation (9) we have:

k = T − t 0 t 1 + n t 2 , where R = t 1 / t 2

The cubic Equation (11) has an explicit solution given by:

n o p t = 2 + θ 3 + α 1 / 3 ( 1 + θ 3 ) + β / 3 (12)

where

β = 1 ( 1 + θ ) α 1 / 3 { ( 2 + θ ) 2 + 3 R ( 1 + 2 θ ) }

and

α = 3 1 + θ [ 3 R { ( R + 1 ) 2 1 ( 1 + θ ) 4 − ( 6 R 2 + 4 R − 2 ) 1 ( 1 + θ ) 3 + 12 R ( R + 1 ) 1 ( 1 + θ ) 2 − ( 8 R 2 + 10 R + 2 ) 1 1 + θ − R − 1 } ] 1 / 2 + 9 R ( 1 ( 1 + θ ) 3 − 1 ( 1 + θ ) 2 + 1 1 + θ ) + ( 2 + θ 1 + θ ) 3

Equation (12) is the optimum number of replicates per subject that is needed to minimize the variance of the estimated RIP when the total cost of the investigation is held fixed.

Note that, when t 1 = 0 and t 2 = 1 (i.e. R = 0 ), then n o p t = 2 + θ , as given in Equation (7).

This means that a special cost structure is implied by the optimal allocation procedure discussed in the previous section. Note also, when θ = 0 , n o p t = 1 + ( 1 + R ) 1 / 2 ≽ 2 , implying that the ratio R = t 1 / t 2 is as important factor in determining the optimal allocation ( n , k ) .

Examples:

R = 0.1 , θ = 0.1 , then n = 2

R = 0.1 , θ = 0.5 , then n = 3

R = 0.2 , θ = 3 , then n = 6

R = 0.5 , θ = 4 , then n = 7

Remarks:

We set as a bench mark to the value of the estimator of RIP a maximum of 1%. That is if the within subject variation relative to the between-subjects variation is above 1%, then repeatability is low, and visa-versa.

Note also, that the estimator of θ is a non-linear function of the sample data, and hence is potentially biased estimator. Moreover, the derived variance is just a first order approximation of the actual variance. Finally, if the measurements are not normally distributed, then construction of confidence interval on the population parameter using the normal quantiles will not be acceptable unless the sample size is quite large. One way to assess the properties of the proposed estimator is to use the nonparametric-bootstrap sampling techniques. We shall address this issue in the data analysis section.

Not all biological markers that are measured on continuous scale have Gaussian distributions. In this section we drop the assumptions of normality regarding the distributions of b i and e i j , and evaluate the effect of non-normality on the estimation of the RIP. The immediate consequences of dropping the assumption of non-normality of the measurements are:

1) The one-way ANOVA mean squares MSB and MSW will not have chi-square distributions.

2) The mean squares MSB and MSW are no longer independent, and hence the ratio of the mean squares will not have the usual F-distribution.

Relaxing the assumption of normality both the measures of Kurtosis of b i and e i j are needed in the calculation of the asymptotic variance of θ ^ [

Let δ e and δ b denote respectively the coefficients of kurtosis of e i j and b i . These quantities are defined as:

δ b = { E ( b i 4 ) / σ b 4 } − 3

δ e = { E ( e i j 4 ) / σ e 4 } − 3

Using results for the balanced one-way ANOVA [

var ( M S W ) = c 1 σ e 4 , var ( M S B ) = c 2 σ e 4 , and cov ( M S W , M S B ) = c 12 σ e 4 , where

c 1 = { k ( n − 1 ) } − 1 [ 2 + δ e n − 1 n ]

c 2 = 2 ( k ) − 1 [ ( θ ) − 2 n 2 ( 1 + δ b ) + 2 n ( θ ) − 1 + ( 1 + δ e n ) ]

c 12 = ( k n ) − 1 δ e

Using the delta method [

Simplifying we get:

var ( θ ^ ) = 1 n 2 { [ θ ( n + θ ) ] 2 c 1 − 2 θ 3 ( n + θ ) c 12 + c 2 θ 4 } (13)

Comments:

The first question that needs to be answered is: which component of variation has the largest effect on the variance of the RIP estimate, and hence on the number of repeats. We answer this question in a heuristic manner. We note from Equation (13) that δ e is divided by the factor {kn} in c 1 , c 2 , and c 12 . The implication is that, as the number of subjects increase, the kurtosis of the error term has negligible effect on the variance of the estimated RIP.

We may also demonstrate the effect of non-normality using tools of probability and power calculations. This can be illustrated through testing of statistical hypotheses on the RIP. Suppose that we need to determine the number of subjects to detect the departure from the null hypothesis H 0 : θ = θ 0 in the direction of the one-sided alternative H 1 : θ = θ 1 < θ 0 , with type-one error rate α and power 1 − β . For fixed n, we can show that:

k = [ z α υ [ θ 0 ] ] + z β υ [ θ 1 ] 2 ( θ 0 − θ 1 ) 2 (14)

If we set the Type I error rate at 5% and power at 80%, for given values of θ 0 , θ 1 , n , δ b and δ e , the estimated values of k can be easily calculated.

Specifically, for an effect size ( θ 0 − θ 1 ) = ( 0.05 − 0.02 ) , δ b = 0 , and δ e = 6 , and n = 5 , then from Equation (14) we need to recruit 6 subjects, while for the same range of values of the RIP we need to recruit 21 subjects if δ b = 6 , and δ e = 0 . The worst situation is when the two components of variation are far from being normal. For example, for the same values under the null and alternative hypotheses, with δ b = 3 , δ e = 3 , and n = 30 , then k = 67 . However, when δ b = δ e = 0 , we need to recruit k = 18 only. These computations illustrate the impact of the departure from normality of the distribution of between and within subject-variations on the sample size requirements.

In this section we apply the methodology presented in this paper on Serum Alanine-aminotransferase ( ALT ). The ALT is a critical parameter for both the assessment and follow-up of patients with liver disease. Therefore, establishing the repeatability and the precision of ALT measurements as a diagnostic marker are of paramount importance. Regardless of gender or body mass index ( BMI ) [

Furthermore, lately published HBV guidelines suggested that treatment decisions should be based on these new ALT levels [

From a large tertiary hospital-based registry, the available data were grouped into female group with 20 subjects and another male group of 30 subjects. In both groups, each subject’s ALT was evaluated three times according to the rules set in [

Bootstrap results

We used R to bootstrap the data. We set the number of bootstrap samples at 1000 for both data sets.

Bootstrap Statistics for females’ data:

original bias std. error

0.001 0.00117 0.0004

As can be seen from

Bootstrap Statistics for males’ data:

original bias std. error

0.002 - 0.0001 0.0006

In contrast to females’ data, the histogram of the sample statistics as shown in

ALT Mean | Standard Deviation | N | Minimum | Maximum |
---|---|---|---|---|

ALT1 75.05 | 39.16 | 20 | 1 | 141 |

ALT2 76.65 | 39.34 | 20 | 2 | 143 |

ALT3 75.84 | 39.30 | 20 | 1 | 142 |

ALT Mean | Standard Deviation | N | Minimum | Maximum |
---|---|---|---|---|

ALT1 106.97 | 81.91 | 30 | 6 | 317 |

ALT2 108.87 | 81.96 | 30 | 8 | 319 |

ALT3 108.00 | 82.06 | 30 | 7 | 318 |

As can be seen from the histograms and the Q-Q plots, the distribution of the estimated RIP = t1* is far from being normally distributed. But we expect that the distributional properties may be closer to normality when the number of subjects is much larger than the number of have here. When one attempts to establish the population-based reference range of health populations, the number of subjects is typically in the hundreds, and the issue of normality may be irrelevant.

Further investigations for the case of categorical measurements and when the number of replications per subject is not fixed, are needed.

The authors declare no conflicts of interest regarding the publication of this paper.

Al-Eid, M. and Shoukri, M.M. (2019) On the Index of Repeatability: Estimation and Sample Size Requirements. Open Journal of Statistics, 9, 530-541. https://doi.org/10.4236/ojs.2019.94035