beran.est: The ocation estimator of Beran (1974)

Description Usage Arguments Details Value Author(s) References Examples

View source: R/Beran_estimator.R

Description

Suppose the observations are sampled from a density f(x-m) where m is the location parameter and f is an unknown density. This function, based on Beran (1974), gives a nonparametric estimator of m. This estimator uses Fourier basis expansion to estimate the score function corresponding to the location model, which has the form

φ_F(t)=f'(F^{-1}(t))/f(F^{-1}(t)).

The Fourier coefficients of the score function are computed by score.coeff.

Usage

1
beran.est(x, theta, M, init, alpha)

Arguments

x

An array of length n; represents the data.

theta

A small number, should be of order O_p(n^{-1/2}). The default is 4n^{-1/2}. This is the tuning parameter for Fourier coefficient estimation.

M

The number of Fourier basis to be used to approximate φ_F.

init

Optional. A vector of initial estimators of m. The default value is the sample median.

alpha

The confidence level for the confidence bands. An (1-alpha) percent confidence interval is constructed. Alpha should lie in the interval (0, 0.50). The default value is 0.05.

Details

theta: theta do not need to depend on the range of the data because the estimators depend only on the rank of the dataponts. If theta=z_n which equals z n^{-1/2}, then the estimated coefficients are root-n consistent by Theorem 2.1 of Beran (1974).

M: A higher value of M decreases bias, but increases the variance. Theorem 4.1 of Beran (1974) states that if M\to∞ as the sample size n grows and \lim_{n\to ∞} M^6/n=0, then this estimator is asymptotically efficient. A larger M always gives a more conservative confidence interval.

init: Beran (1974) recommends using sample median and warns that the method will be sensitive to the choice of the preliminary estimator. This should be a root-n consistent estimator of m. Some other choices are the mean, or the trimmed mean.

Value

A list of length two.

Author(s)

Nilanjana Laha (maintainer), nlaha@hsph.harvard.edu.

References

Beran, R. (1974). Asymptotically efficient adaptive rank estimates in location models. Ann. Statist., 2, 63-74.

Examples

1
2
x <- rlogis(100); beran.est(x, M=10)
beran.est(x, M=20)

nilanjanalaha/log.location documentation built on Dec. 31, 2020, 12:07 a.m.