Next Article in Journal
Effect of Atomic Size Difference on the Microstructure and Mechanical Properties of High-Entropy Alloys
Next Article in Special Issue
Composite Tests under Corrupted Data
Previous Article in Journal
First-Principles Design of Refractory High Entropy Alloy VMoNbTaW
Previous Article in Special Issue
Asymptotic Properties for Methods Combining the Minimum Hellinger Distance Estimate and the Bayesian Nonparametric Density Estimate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Likelihood Ratio Testing under Measurement Errors

1
Faculté de Mathématiques, Laboratoire de Probabilité, Statistique et Modélisation, Université Pierre et Marie Curie (Sorbonne Université), 4 place Jussieu, 75252 Paris CEDEX 05, France
2
Institute of Information Theory and Automation, The Czech Academy of Sciences, Pod Vodárenskou věží 4, 182 08 Prague 8, Czech Republic
3
Faculty of Mathematics and Physics, Charles University, Sokolovská 83, 186 75 Prague 8, Czech Republic
4
Institute of Computer Science, The Czech Academy of Sciences, Pod Vodárenskou věží 2, 182 07 Prague 8, Czech Republic
*
Author to whom correspondence should be addressed.
Submission received: 13 November 2018 / Revised: 6 December 2018 / Accepted: 7 December 2018 / Published: 13 December 2018

Abstract

:
We consider the likelihood ratio test of a simple null hypothesis (with density f 0 ) against a simple alternative hypothesis (with density g 0 ) in the situation that observations X i are mismeasured due to the presence of measurement errors. Thus instead of X i for i = 1 , , n , we observe Z i = X i + δ V i with unobservable parameter δ and unobservable random variable V i . When we ignore the presence of measurement errors and perform the original test, the probability of type I error becomes different from the nominal value, but the test is still the most powerful among all tests on the modified level. Further, we derive the minimax test of some families of misspecified hypotheses and alternatives. The test exploits the concept of pseudo-capacities elaborated by Huber and Strassen (1973) and Buja (1986). A numerical experiment illustrates the principles and performance of the novel test.

1. Introduction

Measurement technologies are often affected by random errors; if the goal of the experiment is to compare two probability distributions using data, then the conclusion can be distorted if the data are affected by some measurement errors. If the data are mismeasured due to the presence of measurement errors, the statistical inference performed with them is biased and trends or associations in the data are deformed. This is common for a broad spectrum of applications e.g., in engineering, physics, biomedicine, molecular genetics, chemometrics, econometrics etc. Some observations can be even undetected, e.g., in measurements of magnetic or luminous flux in analytical chemistry when the flux intensity falls below some flux limit. Actually, we can hardly imagine real data free of measurement errors; the question is how severe the measurement errors are and what their influence on the data analysis is [1,2,3].
A variety of functional models have been proposed for handling measurement errors in statistical inference. Technicians, geologists, and other specialists are aware of this problem, and try to reduce the effect of measurement errors with various ad hoc procedures. However, this effect cannot be completely eliminated or substantially reduced unless we have some additional knowledge on the behavior of measurement errors.
There exists a rich literature on the statistical inference in the error-in-variables (EV) models as is evidenced by the monographs of Fuller [4], Carroll et al. [5], and Cheng and van Ness [6], and the references therein. The monographs [4] and [6] deal mostly with classical Gaussian set up while [5] discusses numerous inference procedure under semi-parametric set up. Nonparametric methods in EV models are considered in [7,8] and in references therein, and in [9], among others. The regression quantile theory in the area of EV models was started by He and Liang [10]. Arias [11] used an instrumental variable estimator for quantile regression, considering biases arising from unmeasured ability and measurement errors. The papers dealing with practical aspects of measurement error models include [12,13,14,15,16], among others. Recent developments in treating the effect of measurement errors on econometric models was presented in [17] or [18] The advantage of rank and signed rank procedures in the measurement errors models was discovered recently in [19,20,21,22,23,24]. The problem of interest in the present paper is to study how the measurement errors can affect the conclusion of the likelihood ratio test.
The distribution function of measurement errors is considered unknown, up to zero expectation and unit variance. When we use the the likelihood ratio test while ignoring the possible measurement errors, we can suffer a loss in both errors of the first and second kind. However, we show that under a small variance of measurement errors, the original likelihood ratio test is still most powerful, only on a slightly changed significance level.
On the other hand, we may consider the situation that H 0 or H 1 are classes of distributions of random variables Z + δ V . Hence, both hypothesis and alternative are composite as families H 0 and H 1 ; if they are bounded by alternating Choquet capacities of order 2, then we can look for a minimax test based on the ratio of the capacities, and/over on the ratio of the pair of the least favorable distributions of H 0 and H 1 , respectively (cf. Huber and Strassen [25]).

2. Likelihood Ratio Test under Measurement Errors

Our primary goal is to test the null hypothesis H 0 that independent observations X = ( X 1 , , X n ) come from a population with a density f against the alternative H 1 that the true density is g , where f and g are fixed densities of our interest. For the identifiability, we shall assume that f and g are continuous and symmetric around 0. Although the alternative is the main concern of the experimenter, some measurement errors or just the nature may cause the situation that the true alternative should be considered as composite. Specifically, X 1 , , X n , can be affected by additive measurement errors, what appears in numerous fields, as illustrated in Section 1.
Hence the alternative is H 1 , δ under which the observations are Z i , δ = X i + δ V i , identically distributed with continuous density g δ . Here, both under the hypothesis and under the alternative, V i are independent random variables, unobservable with unknown distribution, independent of X i ; i = 1 , , n . The parameter δ > 0 is also unknown, only we assume that E V i = 0 and E V i 2 = 1 , for simplicity. The mismeasured, hence unobservable, X i are assumed to have the density g under the alternative. Quite analogously, the mismeasured observations lead to a composite hypothesis H 0 , δ under which the density of observations Z i , δ = X i + δ V i is f δ while the X i are assumed to have density f.
If we knew f δ and g δ , we would use the Neyman-Pearson critical region
W = z : i = 1 n ln g δ ( z i ) f δ ( z i ) u
with u determined so that
P f δ i = 1 n ln g δ ( z i ) f δ ( z i ) u = α ,
with a significance level α . Evidently
I i = 1 n ln g δ ( z i ) f δ ( z i ) u i = 1 n g δ ( z i ) d z i = I i = 1 n ln g ( x i ) f ( x i ) u i = 1 n g ( x i ) d x i I i = 1 n ln g δ ( z i ) f δ ( z i ) u i = 1 n f δ ( z i ) d z i = I i = 1 n ln g ( x i ) f ( x i ) u i = 1 n f ( x i ) d x i .
Indeed, notice that
E g δ I i = 1 n ln g δ ( Z i ) f δ ( Z i ) u | V 1 = v 1 , , V n = v n = E g I i = 1 n ln g ( X i ) f ( X i ) u | V 1 = v 1 , , V n = v n v i R , i = 1 , , n ,
where the expectations are considered with respect to the conditional distribution; a similar equality holds for f δ .
Combining the integration transmission in the conditional distribution, we obtain
I i = 1 n ln g δ ( x i + δ V i ) f δ ( x i + δ V i ) u i = 1 n f ( x i ) d x i I i = 1 n ln g ( x i ) f ( x i ) u i = 1 n f ( x i ) d x i = α ,
hence the size of the critical region W when used for testing H 0 against H 1 differs from α . Then we ask how the critical region W in (1) behaves when it is used as a test of H 0 . This problem we shall try to attack with an expansion of f δ , g δ in δ close to zero.

Approximations of Densities

Put f = f 0 , g = g 0 the densities of X under the hypotheses and alternative, respectively. For the identifiability, we shall assume that f 0 and g 0 are continuous and symmetric around 0. Denote f δ the density of Z δ = X + δ V . This means that X is affected by an additive measurement error δ V , where V is independent of X and E V = 0 , E V 2 = 1 , E V 4 < . Notice that if densities of X and V are strongly unimodal, then that of Z is also strongly unimodal (see [26]). Under some additional conditions on f 0 , g 0 , we shall derive approximations of f δ and g δ for small δ > 0 . More precisely, we assume that both f 0 and g 0 have differentiable and integrable derivatives up to order 5. Then we have the following expansion of f δ and a parallel result for g δ :
Theorem 1.
Assume that f 0 and g 0 are symmetric around 0, strongly unimodal with differentiable and integrable derivatives, up to the order 5. Then, as δ 0 ,
f δ ( z ) = f 0 ( x + δ V ) = f 0 ( x ) + δ 2 d 2 d z 2 f 0 ( x ) + δ 2 4 ! d 4 d z 4 f 0 ( x ) E ( V 4 ) + o ( δ 2 ) , g δ ( z ) = g 0 ( x + δ V ) = g 0 ( x ) + δ 2 d 2 d z 2 g 0 ( x ) + δ 2 4 ! d 4 d z 4 g 0 ( x ) E ( V 4 ) + o ( δ 2 )
Proof. 
Let φ ( u , δ ) = E { e i u Z } be the characteristic function of Z . Then
φ ( u , δ ) = E { e i u X } E { e i u δ V } = φ ( 0 , 0 ) φ V ( u δ ) = φ ( u , 0 ) 1 + 1 2 δ ( i u ) 2 + 1 4 ! δ 2 ( i u ) 4 E ( V 4 ) + o ( δ 2 ) = φ ( u , 0 ) 1 δ 2 u 2 + 1 4 ! δ 2 u 4 E ( V 4 ) + o ( δ 2 ) ,
where φ V denotes the characteristic function of V. Taking the inverse Fourier transform on both sides, we obtain (3), taking the above assumptions on V into account. □
Consider the problem of testing the hypothesis H 0 that the observations are distributed according to density f 0 against the alternative H 1 that they are distributed according to density g 0 . Parallelly, we consider the hypothesis H 0 , δ that observations are distributed according the g δ against the alternative H 1 , δ that the true density is g δ . Let Φ ( x ) be the likelihood ratio test with critical region W = x : i = 1 n ln g 0 ( x i ) f 0 ( x i ) > u and the significance level α , and Φ * = Φ * ( z ) be the test with critical region W * = z : i = 1 n ln g 0 ( z i ) ) f 0 ( z i ) > u based on observations z i = x i + δ V i , i = 1 , , n . We know neither δ nor V , hence the test Φ * is just an application of the critical region W for contaminated data Z 1 , , Z n . Thus, due to our lack of information, we use the test Φ even for testing H 0 , δ against H 1 , δ , and the performance of this test is of interest. This is described in the following theorem:
Theorem 2. (Assume the conditions of Theorem 1).
Then, as δ 0 , the test Φ * is the most powerful even for testing H 0 , δ against H 1 , δ , with a modified significance level satisfying
α δ α + δ 2 | f 0 ( 0 ) | + δ 2 24 E V 4 | f 0 ( 3 ) ( 0 ) | + O ( δ ) .
Proof. 
E f 0 Φ * ( X ) = I ln g 0 ( x + δ V ) f 0 ( x + δ V ) > u f 0 ( x ) d x = I ln g 0 ( x + δ V ) f 0 ( x + δ V ) > u f 0 ( x ) f 0 ( x + δ V ) f 0 ( x + δ V ) d x = I ln g 0 ( x ) f 0 ( x ) > u f 0 ( x δ V ) f 0 ( x ) f 0 ( x ) d x .
If f 0 is symmetric, then the derivative f 0 ( k ) is symmetric for k even and skew-symmetric for k odd, k = 1 , , 4 . Moreover, because | f 0 ( x ) | and | f 0 ( 3 ) ( x ) | are integrable, then lim x ± | f 0 ( x ) | = 0 and lim x ± | f 0 ( 3 ) ( x ) | = 0 . Hence, using the expansion (3), we obtain
E f 0 Φ * ( X ) = E f 0 Φ ( X ) + I ln g 0 ( x ) f 0 ( x ) > u δ 2 f 0 ( x ) + δ 2 24 E V 4 f ( 4 ) ( x ) d x + o ( δ 2 ) E f 0 Φ ( X ) + δ 2 | f 0 ( 0 ) | + δ 2 24 E V 4 | f 0 ( 3 ) ( 0 ) | + o ( δ 2 ) = α + O ( δ ) a s δ 0 .

3. Robust Testing

If the observations are missmeasured or contaminated, we observe Z δ = Z + δ V with unknown δ and unobservable V instead of Z. Hence, instead of simple f 0 and g 0 , we are led to composite hypothesis and alternative H and K . Following [25], we can try to find suitable 2-alternating capacities, dominating H and K and to construct a pertaining minimax test. As before, we assume that Z and V are independent, E V = 0 , E V 2 = 1 , and E V 4 < . Moreover, we assume that f 0 and g 0 are symmetric, strongly unimodal and differentiable up to order 5, with derivatives integrable and increasing distribution functions F 0 and G 0 , respectively. The measurement errors V are assumed to satisfy
1 E V 4 K
with a fixed K , 0 < K < . Hence the distribution of V is restricted to have the tails lighter than t-distribution with 4 degrees of freedom. We shall construct a pair of 2-alternating capacities around specific subfamilies of f 0 and g 0 .
Let us determine the capacity around g 0 ; that for f 0 is analogous. By Theorem 1 we have
g δ ( z ) = g 0 ( z ) + δ 2 d 2 d z 2 g 0 ( z ) + δ 2 4 ! d 4 d z 4 g 0 ( z ) E ( V 4 ) + o ( δ 2 ) , a s δ 0 .
We shall concentrate on the following family K * of densities (similarly for f 0 ):
K * = g δ , κ * : g δ , κ * ( z ) = g 0 ( z ) + δ 2 g 0 ( z ) + κ δ 2 24 g 0 ( 4 ) ( z ) | δ Δ , 1 κ K
with fixed suitable Δ , K > 0 .
Indeed, under our assumptions, each g δ , κ * K * is a positive and symmetric density satisfying
sup δ Δ , κ K sup z { R g δ , κ * ( z ) g 0 ( z ) C K Δ 2 + o ( Δ 2 )
for some C , 0 < C < .
Let G δ , κ * ( B ) , B B , be the probability distribution induced by density g δ , κ * K * , with B being the Borel σ -algebra. Then the set function
w ( B ) = sup G * ( B ) : G * K * if B 0 if B =
is a pseudo-capacity in the sense of Buja [27], i.e., satisfying
(a) 
w ( ) = 0 , w ( Ω ) = 1
(b) 
w ( A ) w ( B ) A B
(c) 
w ( A n ) w ( A ) A n A
(d) 
w ( A n ) w ( A ) A n A
(e) 
w ( A B ) + w ( A B ) w ( A ) + w ( B ) .
Analogously, consider a density f 0 , symmetric around 0 and satisfying the assumptions of Theorem 1 as a simple hypothesis. Construct the family H * of densities and the corresponding family of distributions F δ , κ * ( · ) , δ Δ , κ K similarly as above. Then the set function
v ( B ) = sup F * ( B ) : F * H * if B 0 if B =
is a pseudo-capacity in the sense of Buja [27].
Buja [27] showed that on any Polish space exists a (possibly different) topology which generates the same Borel algebra and on which every pseudo-capacity is a 2-alternating capacity in the sense of [25].
Let us now consider the problem of testing the hypothesis H = F * H * | F * ( · ) v ( · ) against the alternative K = G * K * | G * ( · ) w ( · ) , based on an independent random sample Z 1 , , Z n . Assume that H * and K * satisfy (5). Then, following [27] and [25], we have the main theorem providing the minimax test of H against K with significance level α ( 0 , 1 ) :
Theorem 3.
The test
ϕ ( z 1 , , z n ) = 1 if i = 1 n π ( z i ) > C γ if i = 1 n π ( z i ) = C 0 if i = 1 n π ( z i ) < C
where π ( · ) is a version of d w d v ( · ) and C and γ are chosen so that E v ϕ ( Z ) = α , is a minimax test of H against K of level α .

4. Numerical Illustration

We assume to observe independent observations Z 1 , δ , , Z n , δ for i = 1 , , n , where Z i , δ = X i + δ V i as described in Section 3, where X 1 , , X n are independent identically distributed (with a distribution function F) but unobserved. Let us further denote by Φ the distribution function of N ( 0 , 1 ) and by Φ σ * the distribution function of N ( 0 , σ 2 ) . The primary task here is to test H 0 : F Φ against
H 1 : F ( x ) = ( 1 λ ) Φ ( x ) + λ Φ σ * ( x ) , x R ,
with a fixed σ > 1 and λ ( 0 , 1 ) . We perform all the computations using the R software [28].
To describe our approach to computing the test, we will need the notation for the set of pseudo-distribution functions corresponding to the set of pseudo-densities H * denotes as
H ˜ * = F δ , κ * : F δ , κ * ( z ) = Φ ( z ) + δ 2 f 0 ( z ) + κ δ 2 24 f 0 ( 3 ) ( z ) | δ Δ , 1 κ K ,
where Φ denotes the distribution function of N ( 0 , 1 ) distribution. Under the alternative, the set analogous to K * is defined as
K ˜ * = G δ , κ * : G δ , κ * ( z ) = G 0 ( z ) + δ 2 g 0 ( z ) + κ δ 2 24 g 0 ( 3 ) ( z ) | 0 δ Δ , 1 κ K .
Our task is to approximate
v ( , z ) = sup { F δ , κ * ( z ) ; F δ , κ * H ˜ * } , z R ,
and
w ( , z ) = sup { G δ , κ * ( z ) ; G δ , κ * K ˜ * } , z R .
Here, the functions F δ , κ * ( z ) and G δ , κ * ( z ) are evaluated over a grid with step 0.05. Then, the maximization in (8) and (9) is performed for values of z over the grid and over four boundary values of ( δ , κ ) T , which are equal to ( 0 , 0 ) T , ( 0 , K ) T , ( Δ , 0 ) T , and ( Δ , K ) T . Additional computations with 10 randomly selected pairs of ( δ , κ ) T over δ [ 0 , Δ ] and κ [ 0 , K ] revealed that the optimum is attained in one of the boundary values. Further, the Radon-Nikodym derivatives of V and W are estimated by a finite difference approximation in order to compute the test statistic.
The test rejects H 0 if the test statistics i = 1 n π ( z i ) exceeds a critical value, which (as well as the p-value) can be approximated by a Monte Carlo simulation, i.e., by a repeated random generating random variables X 1 , , X n under H 0 , and we generate them 10,000 times here.
We perform the following particular numerical study. We compute the critical value of the α -test for n = 20 (or n = 40 ), λ = 0.25 , σ 2 = 3 , Δ = 0.2 , K = 1.1 , and α = 0.05 . Further, we are interested in evaluating the probability of rejecting this test for data generated from
F ( x ) = ( 1 λ ˜ ) Φ ( x ) + λ ˜ Φ σ ˜ * ( x ) , x R ,
with different values of λ ˜ and σ ˜ 2 . Its values are shown in Table 1 (for n = 20 ) and Table 2 (for n = 40 ), which are approximated using (again) 10,000 randomly generated variables from (10). The boldface numbers are equal to the power of the test (under the simple H 1 ). The proposed test seems meaningful, while its power is increased for n = 40 compared to n = 20 ; in addition, the power increases with an increasing λ ˜ if σ ˜ 2 is retained; and the power also increases with an increasing σ ˜ 2 if λ ˜ is retained.

5. Conclusions

The likelihood ratio test of f 0 against g 0 is considered in the situation that observations X i are mismeasured due to the presence of measurement errors. Thus instead of X i for i = 1 , , n , we observe Z i = X i + δ V i with unobservable parameter δ and unobservable random variable V i . When we ignore the presence of measurement errors and perform the original test, the probability of type I error becomes different from the nominal value, but the test is still the most powerful among all tests on the modified level.
Under some assumptions on f 0 and g 0 and for δ < Δ , E V 4 K , we further construct a minimax likelihood ratio test of some families of distributions of the Z i = X i + δ V i , based on the capacities of the Huber-Strassen type. The test treats the composite null and alternative hypotheses, which cover all possible measurement errors satisfying the assumptions. The advantage of the novel test is that it keeps the probability of type I error below the desired value ( α = 0.05 ) across all possible measurement errors. The test is performed in a straightforward way, while the user must specify particular (not excessively large) values of Δ and K. We do not consider this a limiting requirement, because parameters corresponding to the severity of measurement errors are commonly chosen in a similar way in numerous measurement error models [5,23] or robust optimization procedures [29]. The critical value of the test can be approximated by a simulation. The numerical experiment in Section 4 illustrates the principles and performance of the novel test.

Author Contributions

Methodology, M.B. and J.J.; Software, J.K.; Writing—Original Draft Preparation, M.B. and J.J.; Writing—Review & Editing, M.B., J.J. and J.K.; Funding Acquisition, M.B.

Funding

The research of Jana Jurečková was supported by the Grant 18-01137S of the Czech Science Foundation. The research of Jan Kalina was supported by the Grant 17-01251S of the Czech Science Foundation.

Acknowledgments

The authors would like to thank two anonymous referees for constructive advice.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boyd, A.; Lankford, H.; Loeb, S.; Wyckoff, J. Measuring test measurement error: A general approach. J. Educ. Behav. Stat. 2013, 38, 629–663. [Google Scholar] [CrossRef]
  2. Brakenhoff, T.B.; Mitroiu, M.; Keogh, R.H.; Moons, K.G.M.; Groenwold, R.H.H.; van Smeden, M. Measurement error is often neglected in medical literature: A systematic review. J. Clin. Epidemiol. 2018, 98, 89–97. [Google Scholar] [CrossRef] [PubMed]
  3. Edwards, J.K.; Cole, S.R.; Westreich, D. All your data are always missing: Incorporating bias due to measurement error into the potential outcomes framework. Int. J. Epidemiol. 2015, 44, 1452–1459. [Google Scholar] [CrossRef] [PubMed]
  4. Fuller, W.A. Measurement Error Models; John Wiley & Sons: New York, NY, USA, 1987. [Google Scholar]
  5. Carroll, R.J.; Ruppert, D.; Stefanski, L.A.; Crainiceanu, C.M. Measurement Error in Nonlinear Models: A Modern Perspective, 2nd ed.; Chapman & Hall/CRC: Boca Raton, FL, USA, 2006. [Google Scholar]
  6. Cheng, C.L.; van Ness, J.W. Statistical Regression with Measurement Error; Arnold: London, UK, 1999. [Google Scholar]
  7. Carroll, R.J.; Maca, J.D.; Ruppert, D. Nonparametric regression in the presence of measurement error. Biometrika 1999, 86, 541–554. [Google Scholar] [CrossRef] [Green Version]
  8. Carroll, R.J.; Delaigle, A.; Hall, P. Non-parametric regression estimation from data contaminated by a mixture of Berkson and classical errors. J. R. Stat. Soc. B 2007, 69, 859–878. [Google Scholar] [CrossRef]
  9. Fan, J.; Truong, Y.K. Nonparametric regression estimation involving errors-in-variables. Ann. Stat. 1993, 21, 23–37. [Google Scholar] [CrossRef]
  10. He, X.; Liang, H. Quantile regression estimate for a class of linear and partially linear errors-in-variables models. Stat. Sin. 2000, 10, 129–140. [Google Scholar]
  11. Arias, O.; Hallock, K.F.; Sosa-Escudero, W. Individual heterogeneity in the returns to schooling: Instrumental variables quantile regression using twins data. Empir. Econ. 2001, 26, 7–40. [Google Scholar] [CrossRef]
  12. Hyk, W.; Stojek, Z. Quantifying uncertainty of determination by standard additions and serial dilutions methods taking into account standard uncertainties in both axes. Anal. Chem. 2013, 85, 5933–5939. [Google Scholar] [CrossRef]
  13. Kelly, B.C. Some aspects of measurement error in linear regression of astronomical data. Astrophys. J. 2007, 665, 1489–1506. [Google Scholar] [CrossRef]
  14. Marques, T.A. Predicting and correcting bias caused by measurement error in line transect sampling using multiplicative error model. Biometrics 2004, 60, 757–763. [Google Scholar] [CrossRef] [PubMed]
  15. Rocke, D.M.; Lorenzato, S. A two-component model for measurement error in analytical chemistry. Technometrics 1995, 37, 176–184. [Google Scholar] [CrossRef]
  16. Akritas, M.G.; Bershady, M.A. Linear regression for astronomical data with measurement errors and intrinsic scatter. Astrophys. J. 1996, 470, 706–728. [Google Scholar] [CrossRef]
  17. Hausman, J. Mismeasured variables in econometric analysis: Problems from the right and problems from the left. J. Econ. Perspect. 2001, 15, 57–67. [Google Scholar] [CrossRef]
  18. Hyslop, D.R.; Imbens, Q.W. Bias from classical and other forms of measurement error. J. Bus. Econ. Stat. 2001, 19, 475–481. [Google Scholar] [CrossRef]
  19. Jurečková, J.; Picek, J.; Saleh, A.K.M.E. Rank tests and regression rank scores tests in measurement error models. Comput. Stat. Data Anal. 2010, 54, 3108–3120. [Google Scholar] [CrossRef]
  20. Jurečková, J.; Koul, H.L.; Navrátil, R.; Picek, J. Behavior of R-estimators under Measurement Errors. Bernoulli 2016, 22, 1093–1112. [Google Scholar] [CrossRef]
  21. Navrátil, R.; Saleh, A.K.M.E. Rank tests of symmetry and R-estimation of location parameter under measurement errors. Acta Univ. Palacki. Olomuc. Fac. Rerum Nat. Math. 2011, 50, 95–102. [Google Scholar]
  22. Navrátil, R. Rank tests and R-estimates in location model with measurement errors. In Proceedings of Workshop of the Jaroslav Hájek Center and Financial Mathematics in Practice I; Masaryk University: Brno, Czech Republic, 2012; pp. 37–44. [Google Scholar]
  23. Saleh, A.K.M.E.; Picek, J.; Kalina, J. R-estimation of the parameters of a multiple regression model with measurement errors. Metrika 2012, 75, 311–328. [Google Scholar] [CrossRef]
  24. Sen, P.K.; Jurečková, J.; Picek, J. Rank tests for corrupted linear models. J. Indian Stat. Assoc. 2013, 51, 201–230. [Google Scholar]
  25. Huber, P.; Strassen, V. Minimax tests and the Neyman-Pearson lemma for capacities. Ann. Stat. 1973, 2, 251–273. [Google Scholar] [CrossRef]
  26. Ibragimov, I.A. On the composition of unimodal distributions. Theor. Probab. Appl. 1956, 1, 255–260. [Google Scholar] [CrossRef]
  27. Buja, A. On the Huber-Strassen theorem. Probab. Theory Relat. Fields 1986, 73, 149–152. [Google Scholar] [CrossRef]
  28. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017; Available online: https://www.R-project.org/ (accessed on 15 September 2018).
  29. Xanthopoulos, P.; Pardalos, P.M.; Trafalis, T.B. Robust Data Mining; Springer: New York, NY, USA, 2013. [Google Scholar]
Table 1. Probability of rejecting the test in the simulation with n = 20 .
Table 1. Probability of rejecting the test in the simulation with n = 20 .
Value of λ ˜ Value of σ ˜ 2
3456
0.250.390.520.610.67
0.350.500.670.750.81
0.450.610.760.850.89
Table 2. Probability of rejecting the test in the simulation with n = 40 .
Table 2. Probability of rejecting the test in the simulation with n = 40 .
Value of λ ˜ Value of σ ˜ 2
3456
0.250.550.730.820.87
0.350.720.860.930.96
0.450.820.940.970.99

Share and Cite

MDPI and ACS Style

Broniatowski, M.; Jurečková, J.; Kalina, J. Likelihood Ratio Testing under Measurement Errors. Entropy 2018, 20, 966. https://0-doi-org.brum.beds.ac.uk/10.3390/e20120966

AMA Style

Broniatowski M, Jurečková J, Kalina J. Likelihood Ratio Testing under Measurement Errors. Entropy. 2018; 20(12):966. https://0-doi-org.brum.beds.ac.uk/10.3390/e20120966

Chicago/Turabian Style

Broniatowski, Michel, Jana Jurečková, and Jan Kalina. 2018. "Likelihood Ratio Testing under Measurement Errors" Entropy 20, no. 12: 966. https://0-doi-org.brum.beds.ac.uk/10.3390/e20120966

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop