Detecting Selection Using Time-Series Data of Allele Frequencies with Multiple Independent Reference Loci

Recently, in 2013 Feder et al. proposed the frequency increment test (FIT), which evaluates natural selection at a single diallelic locus by the use of time-series data of allele frequencies. This test is unbiased under conditions of constant population size and no sampling noise. Here, we expand upon the FIT by introducing a test that explicitly allows for changes in population size by using information from independent reference loci. Various demographic models suggest that our proposed test is unbiased irrespective of fluctuations in population size when sampling noise can be ignored and that it has greater power to detect selection than the FIT if sufficient reference loci are used.

ABSTRACT Recently, in 2013 Feder et al. proposed the frequency increment test (FIT), which evaluates natural selection at a single diallelic locus by the use of time-series data of allele frequencies. This test is unbiased under conditions of constant population size and no sampling noise. Here, we expand upon the FIT by introducing a test that explicitly allows for changes in population size by using information from independent reference loci. Various demographic models suggest that our proposed test is unbiased irrespective of fluctuations in population size when sampling noise can be ignored and that it has greater power to detect selection than the FIT if sufficient reference loci are used. In population genetics, most data are obtained from a single point in time. When genetic time-series data are available, the use of such data to detect and estimate natural selection is an attractive concept. Timeseries data may have direct information about natural selection because they affect allele frequencies in time. Bollback et al. (2008) introduced a statistical framework for estimating and testing natural selection by using time-series data of allele frequencies at a single diallelic locus. These authors applied their framework to an ancient human DNA sequence (Hummel et al. 2005) and a sample from an experimental evolution study of the bacteriophage, MS2 (Bollback and Huelsenbeck 2007). Recent advances in high-throughput sequencing technology, including pooled DNA sequencing, have facilitated the acquisition of time-series data, and Bollback et al.'s (2008) method has been extended to more complicated situations. Mathieson and McVean (2013) applied a lattice model of population subdivision that enabled joint estimation of migration rate and spatially varying selection coefficients. The allele age also is an important parameter in selection inferences because it can provide information regarding the origin of a particular phenotype associated with the allele. Malaspinas et al. (2012) developed a method to estimate the selection coefficient and the allele age simultaneously. In addition, there has been increasing numbers of studies in which researchers focus on specific settings in each evolutionary experiment (e.g., Illingworth and Mustonen 2011;Gallet et al. 2012;Illingworth et al. 2012). Feder et al. (2013) reported recently that Bollback et al.'s (2008) standard x2 -based test for selection is biased for realistic data with few sampled time points. When the number of sampled time points is sufficiently large, the likelihood ratio statistic (LRS) follows a x 2 distribution. However, the actual number of sampled time points rarely exceeds a few dozen. Particularly, when the null hypothesis is composite and the profile likelihood is used, the estimation of nuisance parameters can substantially bias inferences of the parameters of interest (e.g., see Chapter 10 of Pawitan 2001). For the problem described in this report, the nuisance parameter is the population size.
To avoid bias, Feder et al. (2013) proposed two methods that both were modeled under conditions of constant population size and no sampling noise. In the empirical likelihood ratio test (ELRT), the population size is preliminarily estimated under neutrality as a first approximation. The estimated population size then is used to generate the empirical distribution of the LRS by computer simulation. Neutrality then can be evaluated by comparing the observed LRS with the empirical distribution. Although the ELRT was shown to be unbiased, this approach can be computationally intensive. To reduce the computational load, Feder et al. (2013) proposed the frequency increment test (FIT). The statistic used for FIT is defined as the following: Let x 0 ; x 1 ; . . . ; x L be the population frequencies of one allele at a diallelic locus of interest at the sampled time, t 0 ¼ 0; t 1 ; . . . ; t L . The sampling time scales are short compared to the population size. Then the standardized allele frequency increment, is approximately normally distributed with mean 0 under the null model, that is, neutral evolution. The variance of Y i is equal to 1/ (2N) in the Wright2Fisher model with N diploids. However, the variance is unknown because N is unknown. In such a situation, natural selection can be evaluated by letting where Y and S are the sample mean and variance of Y i , respectively. We then perform a t-test using the fact that t FIT follows the Student's t distribution with L 2 1 degrees of freedom under the null model. This is the FIT.
The FIT treats the nuisance parameter, N, as an unknown parameter instead of estimating it. When Y i for any i follows the normal distribution with the same variance under the null model, the FIT is an exact and unbiased test. Feder et al. (2013) verified that actual type I error rates approach the nominal significance level for various parameter settings. These investigators also demonstrated that the power of the FIT is equal to or greater than that of the ELRT. Although the FIT does not account for the sampling process from a population explicitly, the test was shown to work well even when the sampling process exists if the sample size is not small.
The FIT is a simple and bias-controlled method to detect selection. However, it is not clear whether the FIT works well when the population size fluctuates. Theoretically, under the null model with fluctuating population size Y i does not follow the same normal distribution for all i, and therefore, t FIT (Data) does not follow the Student's t distribution.
This study is an extension of Feder et al.'s (2013) FIT that allows for fluctuations in population size by using reference loci. First, the FIT's actual type I error rates in a fluctuating population is investigated. Then, a new test is introduced, the frequency increment test with reference loci (FITR). Given a fluctuating population size, the FITR's actual type I error rates are almost the same as the nominal significance level. Then, the powers of the FITR and the FIT to detect natural selection were evaluated. Finally, the simulation method used in this study was validated and the properties of the FITR in practical situations were investigated. Model descriptions are presented just below and added before introducing the FITR.

Model and simulation methods
Let us consider a population evolving according to the Wright-Fisher model with fluctuating population size. The population size fluctuates as a function of generation time, t, and is denoted by N(t). To investigate the actual type I errors and the powers of Feder et al.'s (2013) FIT and the FITR introduced in this study, we conducted computer simulations under the five demographic models shown in Figure 1. The two alleles at a diallelic locus of interest are denoted by A 0 and a 0 , respectively. At generation times t 0 ¼ 0; t 1 ; . . . ; t L ¼ T, the frequencies of a 0 are denoted by x 0;0 ; x 0;1 ; . . . ; x 0;L . Here, t 0 = 0 and t L = T are the first and the last sampling times, respectively, and the number of sampled times is L + 1. The fitnesses of genotypes A 0 A 0 , A 0 a 0 , and a 0 a 0 are assumed to be 1, 1 þ 0:5s 0 , and 1 þ s 0 , respectively (i.e., no dominance is assumed). The population size, N(t), is independent of the frequency of a 0 . As described in the next section, the FITR also uses neutral reference loci.
In the Wright-Fisher model, the allele frequency can be obtained exactly every generation as a binomial distribution. However, the generation of these data poses an extreme computational burden that is impractical for large populations (Figure 1). To avoid this burden and simulate changes in allele frequencies, we applied the pseudosampling method (Kimura and Takahata 1983), which is an improved version of the methods of Kimura (1980). In this method, to determine the allele frequency every generation, a uniform random number with the same mean and variance as those of the exact binomial distribution is used when the allele frequency is moderate. When the allele frequency is high or low, a Poisson random number with the same mean as that of the exact binomial random number is used. In this study, a frequency of #5 minor alleles in the population was used as the criteria for high or low allele frequencies. In addition, the normal distribution was used instead of the uniform distribution because the normal distribution better approximates the binomial distribution, which next-generation allele frequencies follow under the Wright-Fisher model.
to have less power to detect selection when the population size is fluctuating.

Frequency increment test with reference loci (FITR)
Here we propose a new test, the FITR. Consider R reference loci in addition to the focal locus. It is assumed that these R þ 1 loci are evolving independently and that R reference loci are evolving under neutrality. We denote by x h;0 ; x h;1 ; . . . ; x h;L ðh ¼ 1; 2; :::; RÞ the population frequencies of one allele at the h-th reference diallelic locus at times t 0 ¼ 0; t 1 ; . . . ; t L ¼ T. Recall that x 0,i is the allele frequency of the focal locus at t i . Suppose that N(t) is a step function and let N i be the population size from t i 2 1 to t i such that Nðt i21 , t # t i Þ ¼ N i . Note that although the following FITR discussion assumes NðtÞ is a step function, the same discussion can apply even the case that N(t) is a continuous function ( Figure 1) because N i can be interpreted as the variance effective size over the period from t i21 to t i . For this reason, the FITR is unbiased irrespective of fluctuations in population size.
n Table 1 Actual type I error rates (%) of FIT Values indicate the actual type I error rates obtained by 100,000 simulations under a nominal significance level of 5%. The initial allele frequency, x 0,0 , was assumed to be 0.5. FIT, frequency increment test. a Dt ¼ t i 2 t i 2 1 ði ¼ 1; 2; . . . ; LÞ.
n Values indicate the actual type I error rates obtained by 100,000 simulations under a nominal significance level of 5%. The initial frequencies for all R +1 loci, x h,0 , are assumed to be 0.5. FITR, frequency increment test with reference loci.
The numbers of reference loci, R, are randomly assigned to each demographic model.

Figure 2
Powers of the FITR and the FIT in various demographic models. Powers of the FIT (black line) and the FITR with R reference loci (colored lines) are shown as functions of the selection coefficients in the five demographic models. Each point corresponds to the power obtained by 100,000 simulations at the 5% significance level. The duration of sampling time and the number of sampled points were T = 1000 and (L + 1) = 11, respectively. The intervals between any two adjacent sampled points were the same at Dt ¼ t i 2 t i 2 1 ¼ 100 ði ¼ 1; 2; . . . ; LÞ. The initial frequency for all R+ 1 loci, x h;0 , was assumed to be 0.5.

When x h;i21
is not close to 0 or 1 and where follows the standard normal distribution for h = 0 under the null model and for h ¼ 1; 2; :::; R under the null or alternative models. We then consider whether the allele frequency change from t i 2 1 to t i for the focal locus, Y 0,i , is significant. If N i is known, we can test for neutrality using the fact that Y 0,i follows the standard normal distribution. In this case, however, N i is unknown. Let us then define a statistic, The statistic t FITRðiÞ is independent of N i , as seen in (4), because N i in (3) is canceled out. In addition, t FITRðiÞ is independent of Dt i . In (3), the numerator follows the standard normal distribution, and the denominator is equal to the square root of the x 2 random variable divided by its degrees of freedom, R. Because the numerator and denominator are independent, t FITRðiÞ follows a Student's t distribution with R degrees of freedom (Fisher 1925). Although we determined the form of the test statistic, t FITRðiÞ , intuitively, t FITRðiÞ can be derived as the exact LRS using data, Dx i ¼ ðDx 0;i ; Dx 1;i ; . . . ; Dx R;i Þ, in a plausible setting (see Appendix). In other words, the aforementioned t-test is equivalent to the likelihood ratio test under conditions of normality, as observed in several statistical situations (see, e.g., Lehmann and Romano 2005). Next, let us define a statistic using all the data from t 0 ¼ 0 to t FITR is the standardized sum of t FITRðiÞ overall i. The standardization factor 1= ffiffiffi L p allows for a straightforward interpretation of the statistic because t FITR asymptotically follows the standard normal distribution as R becomes large. The exact distribution of t FITR with infinite R is difficult to express explicitly, but the distribution of t FITR can be obtained empirically by generating L Student's t random variables with R degrees of freedom and summing them. This approach is valid because each t FITRðiÞ follows a Student's t distribution with R degrees of freedom. We obtained t FITR using the statistical package R (http://www.R-project.org) with 100,000 simulations of t FITR for each combination of R and L. The test using t FITRðiÞ is the FITR, an exact significance test assuming Y h,i follows the standard normal distribution. That is, the actual type I error rate of the test is expected to be close to the nominal significance level regardless of fluctuations in population size. Unlike t FITRðiÞ , the t FITR statistic is not the exact LRS, which is very difficult to express explicitly. An ad hoc interpretation of the test statistic, t FITR , is presented in the Appendix.
Type I error rate of FITR and powers of FITR and FIT Table 2 shows the actual type I error rates of the FITR. As expected, for all demographic models, the actual type I error rates are close to the nominal level. Figure 2 shows the powers of the FITR and the FIT as a function of the strength of selection. In all demographic models, including the constant-size model, the FITR had more power than the FIT if five or more reference loci were used. For model 2 (slow-growth model) and model 3 (moderate-bottleneck model), the power of the FIT was acceptable. However, for model 4 (rapid-growth model) and model 5 (severe-bottleneck model), the power of the FIT was relatively small, and the FITR demonstrated much larger power than the FIT. Figure 3A displays the powers of the FITR and the FIT as functions of the number of sampling points, L. In many cases, the powers approached certain asymptotic values with increasing L. In this case, the values of L for which the powers approached their asymptotes were relatively small (e.g., L = 10 or 20). The powers of the FITR for R = 1 or 2 in Model 1 decreased somewhat as the sampling points increased. This trend also was observed for the FITR with R = 1 and for the FIT in Model 5. The powers of the FITR and the FIT as functions of the duration of sampling time, T, are given in Figure 3B. For all cases, the powers increased with increasing T, as expected. For Figure 3, A and B, the FITR had more power than the FIT if 10 or more reference loci were used even in the constant-size model. This difference in power was more obvious for model 5 (severe-bottleneck model).

Applying the simulation method and FITR to practical situations
The FITR was developed and evaluated for type I error rate and power under ideal conditions of "selectively neutral" reference loci evolving independently of the focal locus and of each other, allele frequencies at n Values indicate the rejection rates (%) obtained by 10,000 simulations for binomial sampling or by 100,000 simulations for pseudo-sampling under a nominal significance level of 5%. The number of reference were R = 10 The initial frequencies for all R + 1 loci, x h,0 , are assumed to be 0.5. a Dt ¼ t i 2 t i 2 1 ði ¼ 1; 2; . . . ; LÞ. b r, recombination fraction per generation between two adjacent loci of R + 1 loci. "Free" refers to free recombination. c Binomial, the binomial sampling.

Figure 4
The effects of selection at reference loci on the power of the FITR. The powers of the FITR under various selection strengths, s 0 , at focal loci are shown as functions of the selection coefficient, s h ðh 6 ¼ 0Þ, at reference loci for demographic models 1 and 5. Each point corresponds to the power obtained by 100,000 simulations at the 5% significance level. The number of reference loci, the duration of sampling time, and the number of sampled points were R = 10, T = 1000, and (L + 1) = 11, respectively. The intervals between any two adjacent sampled points were the same at Dt ¼ t i 2 t i 2 1 ¼ 100 ði ¼ 1; 2; . . . ; LÞ. The initial frequency for all R þ 1 loci, x h;0 , was assumed to be 0.5. reference loci 6 ¼ 0 or 1, definable FITR statistics, and exactly known population allele frequencies. In practice, these ideal conditions may be violated. Before this section, the computer simulation method used in this study had not been validated. Here, we discuss the applicability of the simulation and describe cases that violate the aforementioned assumptions.
For the determination of allele frequencies in successive generations at R + 1 loci, the exact binomial sampling is computationally intensive and impractical for realistic population sizes ( Figure 1). Therefore, we used a pseudosampling method to simulate the binomial sampling process (Wright2Fisher model). Even for a small, constant-size (n = 20) Wright2Fisher population, the fixation times for a mutant obtained by the pseudosampling were consistent with those obtained by binomial sampling (Kimura 1980). However, we considered only relatively short time scales, and the performance of the pseudo-sampling method was not obvious.
In Table 3, the rows denoted by r = "free" correspond with an assumption of free recombination (i.e., evolving independently) n Values indicate rejection rates (%) obtained by 100,000 simulations under a nominal significance level of 5%. Values in parentheses correspond to the mean number of reference loci used to calculate the FITR statistics. The number of reference loci at t 0 were R = 10. The initial frequency of the focal locus, x 0,0 , was assumed to be 0.5. a Dt ¼ t i 2 t i 2 1 ði ¼ 1; 2; . . . ; LÞ. b x h,0 , the allele frequencies of the reference loci at t 0 . Figure 5 The effects of sampling error on the power of the FITR. The powers of the FITR under various sampling regimes are shown as functions of the selection coefficients for demographic models 3 and 5. All, or 5000, 1000, 500, or 100 individuals in a population were assumed to be sampled. Sampling was assumed to be conducted by binomial sampling at each (R + 1) locus and at each (L + 1) time point. Each point corresponds to the power obtained by 100,000 simulations at the 5% significance level. The number of reference loci, the duration of sampling time, and the number of sampled points were R = 10, T = 1000, and (L + 1) = 3 (upper graphs) or 11 (lower graphs), respectively. The intervals between any two adjacent sampled points were the same at D t = t i -t i-1 = 500 (top graphs) or 100 (bottom graphs) (i = 1,2,. . .,L). The initial frequency for all R þ 1 loci, x h;0 , were assumed to be 0.5.
among R + 1 loci. Rejection rates simulated by the exact binomial sampling and by the pseudosampling are given. Demographic Models 1 and 5 with reduced population sizes were used (see Table 3, legend). We did not observe any differences in the results generated by the binomial sampling vs. the pseudosampling under neutral or selective cases. These findings support the applicability of pseudo-sampling to our problem of concern. Next, we considered the case in which the reference loci and focal loci were not independent. We limited our analysis to the case in which R + 1 loci were in linkage equilibrium (LE) at t=0. For closely linked loci, linkage disequilibrium (LD) is a distinct possibility. In addition, selection at the focal locus can drastically promote LD (e.g., Sabeti et al. 2002). However, for example, empirical studies of human genomes suggest that LD can be extended, at most, to several megabase pairs from the selective locus (e.g., Saunders et al. 2005). Because the genome is large compared with the megabase pairs scale, we can select R reference loci such that R + 1 loci are in LE. For this reason, our discussion is limited to the case in which R + 1 loci are in LE.
In Table 3, the rows indicated by r = 0.1, 0.01, and 0 describe results corresponding to a case in which the per-generation recombination fraction between any two adjacent loci are 0.1, 0.01, and 0, respectively. At t = 0, R + 1 loci are assumed to be in LE. That is, the alleles at R + 1 loci are randomly combined to form haplotypes. The simulations were conducted by the exact binomial sampling. For the neutral or selective cases, we observed no obvious differences between free recombination and limited recombination (r = 0.1, 0.01, and 0; Table 3). That is, the type I error rates and powers were maintained regardless of recombination fractions.
A case in which the reference loci are under selection is evaluated in Figure 4. The selection model is the same because the focal loci and R loci are under the same degree of selection. That is, for all hð6 ¼ 0Þ loci, the fitnesses of genotypes A h A h , A h a h , and a h a h are assumed to be 1, 1 + 0.5s h , and 1 + s h (s 1 = s 2 = ÁÁÁ = s R ), respectively. The effects of selection at the reference loci are conservative for type I error rates (see the case of s 0 = 0 in Figure 4). The results of Model 1 suggest that if Ns h , 5, there is little difference in rejection rates compared to the neutral case. Including Model 5, if the condition s h # 1/2s 0 is met, the power is not decreased. That is, the power is not highly sensitive to selection at the reference loci. Nevertheless, we recommend using synonymous sites or noncoding regions as references.
We next considered a situation in which allele frequencies at the reference loci were low at t = 0 and some allele frequencies could become 0 or 1 by t = T. When allele frequencies became 0 or 1, the statistic, t FITR , in (6) could not be defined. Therefore, it was practical to remove these loci from the calculation of t FITR . Table 4 indicates how rejection rates are changed by the data handling. Values in parentheses indicate the average number of reference loci used to calculate the FITR statistics. For L = 2 the type I error rate was inflated by a few percent (e.g., 6.69% at most, Model 1) possibly because changes in allele frequencies at reference loci are biased toward smaller values when loci are removed for which the frequencies of alleles become 0 or 1. These apparently reduced changes in allele frequencies could bring about overestimates of change at the selective locus. For L = 10 the inflation of the type I error rate becomes small. In general, to prevent inflation of the type I error rate, loci having moderate frequencies of alleles (e.g., $10%) should be used in this test.
The effects of sampling error on the type I error rate and power of the FITR are shown in Figure 5. In general, the effects of sampling error on the type I error rate were conservative. As expected, the power decreased as the number of sampled individuals increased. The degree of power reduction differed for different demographic models or values of L. This finding reflects that the power is influenced by the relative magnitudes of changes in allele frequencies at R + 1 loci and sampling errors. As L increased, the relative changes in allele frequencies to the sampling errors decreased. Thus, power was more reduced for larger L. Regarding the demographic models, the population size of Model 5 was larger than that of Model 1. Therefore, the relative changes in allele frequencies to the sampling errors were larger in Model 5, and the degree of power is large in Model 5.
In this study, we proposed a neutrality test, the FITR, to accommodate fluctuations in population size using reference loci. Our test is an extension of Feder et al.'s (2013) FIT. By computer simulation, the actual type I error rate of the FITR was nearly equal to the nominal significant level regardless of fluctuations in population size when sampling noise could be ignored. The FITR detected selection with remarkable power under conditions of rapid growth (model 4) and severe bottleneck (model 5). Even under a model of constant population size, the FITR using 10 or more reference loci had more power than the FIT.
We also discussed the performance of the FITR in practical situations. The effects of selection at the reference loci were small unless selection was strong. Our findings indicated that when R + 1 were in LE, those loci should be considered independent of each other. In addition, loci with moderate frequencies of alleles should be used as references. Our findings may facilitate the development of more sophisticated methods using independent reference loci, including a method that can quantify (estimate) the strength of selection. These methods will enable appropriate inferences about natural selection in real and dynamic populations. Figure 5.
ACKNOWLEDGMENTS I thank K. Ikeo and S. Mano for useful suggestions on the study. I also thank the associate editor and the two anonymous reviewers for their helpful comments on an earlier version of this manuscript, which improved this work substantially. This work was supported by health labor sciences research grant from The Ministry of Health Labor and Welfare (H23-jituyouka(nanbyou)-006).

LITERATURE CITED APPENDIX
Suppose that there are two segregating alleles, denoted by A h and a h , at the h-th locus ðh 2 ð0; 1; 2; . . . ; RÞÞ. The fitnesses of genotypes A h A h , A h a h , and a h a h are assumed to be 1; 1 þ 0:5s h ; and 1 þ s h , respectively. As described in the main text, the case of h = 0 corresponds to the focal locus. When Dt i ¼ ðt i 2 t i21 Þ is small compared with N i , the change in the allele frequency of a h , x h;i 2 x h;i21 , from t i21 to t i approximately follows a normal distribution, i) t FITRðiÞ is the exact LRS using data Dx i We will show that t FITRðiÞ given by (3) or (4) Under H 1 , the probability density of Dx h;i is given by where c$ ¼ ffiffiffi ffi c9 p . Therefore, t FITR(i) as given by (3) or (4) is the exact LRS using data Dx i .
ii) Ad hoc interpretation for t FITR as a test statistic using the data Dx As R grows, the t distribution for t FITR(i) approaches the normal distribution with mean s 0 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi N i Dt i x 0;i 2 1 ð1 2 x 0;i 2 1 Þ=2 p and variance 1, Unfortunately, s 0 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi N i Dt i x 0;i 2 1 ð1 2 x 0;i 2 1 Þ=2 p varies with N i , Dt i , and x 0;i 2 1 . When x 0;i 2 1 ¼ 0:5; 0:4; 0:3; 0:2; 0:1; and 0:05; ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi x 0;i21 ð1 2 x 0;i21 Þ p ¼ 0.50, 0.490, 0.458, 0.400, 0.300 and 0.218, respectively. Nevertheless, the values of ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi x 0;i21 ð1 2 x 0;i21 Þ p are roughly the same if x 0;i 2 1 are far from 0 or 1. In addition, the values of ffiffiffiffiffi N i p and ffiffiffiffiffiffi Dt i p vary with i. However, we try to consider s 0 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi N i Dt i x 0;i21 ð1 2 x 0;i21 Þ=2 p as a constant m, such that t FITRðiÞ $ Nðm; 1Þ: Consider a test H 0 :m = 0 against H 1 :m 6 ¼ 0 using L i.i.d samples from the distribution. The LRS is given by the sum of t FITRðiÞ , X L i¼1 t FITRðiÞ . That is, when R is large and the variation of the value for ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi N i Dt i x 0;i21 ð1 2 x 0;i21 Þ p is not large, t FITR given by (5) or (6) is expected to be close to the LRS.