Some properties of modified integrated squared error estimators for the stable laws

The estimation procedure of Paulson, Holcomb and Leitch (1975) for the parameters of the stable laws is shown to be similar in spirit to the modified X²minimum procedure. This observation suggests that a class of modified integrated squared error procedures may be developed for the stable laws as well as much more generally. For the stable case, some influence curves, asymptotic covariances, and efficiencies are given, and the robustness of maximum likelihood estimators is discussed.     

On MISE of a Non linear Wavelet Estimator of the Regression Function Based on Biased Data under Strong Mixing

In this paper, we consider the adaptation of the non linear wavelet-based estimator of the regression function for the biased data setup under strong mixing. We provide an asymptotic sharp bound for the mean integrated squared error (MISE) of the estimator, that is nearly optimal in the minimax sense over a large range of Besov function classes.     

A fourier integral density estimate: a Monte Carlo study

Davis (1977) proposed the use of a kernel density estimate which is the sample characteristic function integrated over (-A(n) , A(n)), where A(n) is chosen to minimize the mean integrated square error of the estimate. The scalar, A(n), is determined by the sample size and the population characteristic function. This paper investigates, in a Monte Carlo study, the mean integrated square error obtained under a procedure suggested by Davis (1977) for estimating A(n) when the population characteristic function is unknown.     

Identifiability of two‐component skew normal mixtures with one known component

We give sufficient identifiability conditions for estimating mixing proportions in two‐component mixtures of skew normal distributions with one known component. We consider the univariate case and two multivariate extensions: a multivariate skew normal distribution (MSN) and the canonical fundamental skew normal distribution (CFUSN). The characteristic function of the CFUSN distribution is additionally derived.     

A Bezier curve method in interval-censored data

In this paper we propose a Bezier curve method to estimate the survival function and the median survival time in interval-censored data. We compare the proposed estimator with other existing methods such as the parametric method, the single point imputation method, and the nonparametric maximum likelihood estimator through extensive numerical studies, and it is shown that the proposed estimator performs better than others in the sense of mean squared error and mean integrated squared error. An illustrative example based on a real data set is given.     

Omnibus tests for the error distribution in the linear regression model

Test procedures are constructed for testing the goodness-of-fit of the error distribution in the regression context. The test statistic is based on an L ²-type distance between the characteristic function of the (assumed) error distribution and the empirical characteristic function of the residuals. The asymptotic null distribution as well as the behavior of the test statistic under contiguous alternatives is investigated, while the issue of the choice of suitable estimators has been particularly emphasized. Theoretical results are accompanied by a simulation study.     

Testing for normality with panel data

Classical omnibus and more recent methods are adapted to panel data situations in order to jointly test for normality of the error components. The test statistics incorporate either the empirical distribution function or the empirical characteristic function, these functions resulting from estimation of the fixed and random components. Monte Carlo results show that the new procedure based on the empirical characteristic function compares favorably with classical methods.     

Goodness-of-fit test based on empirical characterisitc function

This paper discusses a goodness-of-fit test that uses the integral of the squared modulus of the difference between the empirical characteristic function of the sample data and the characteristic function of the hypothesized distribution. Monte Carlo procedures are employed to obtain the empirical percentage points for testing the fit of normal, logistic and exponential distributions with unknown location and scale parameters. Results of Monte Carlo power comparisons with other well-developed goodness-of-fit tests are summarized. Tne proposed test is shown to have superior power for testing the fit of the logistic distibotion (for moderate sample sizes) against a wide range of alternative distributions.     

Optimal Sequential Design in a Controlled Non-parametric Regression

Abstract. In a non‐parametric regression, the heteroscedasticity (dependence of the variance of the regression error on the predictor) can be a serious complication in estimation or visualization of an underlying regression function. If a controlled sampling is permitted, then the statistician can choose the design of predictors which attenuates the effect of heteroscedasticity. It is proposed to use a design which minimizes the mean integrated squared error of the regression function estimation. Then the corresponding optimal design density is proportional to the standard deviation of the regression error (the so‐called scale function). Because in general the statistician does not know an underlying scale function, the natural question is as follows: is it possible to suggest a sequential design which performs as well as an oracle that knows the underlying scale function? The answer is ‘yes’, and a corresponding sequential procedure is developed. It is proved, for the first time in the literature, that a data‐driven sequential design, together with an adaptive regression estimator, can mimic the oracle and be sharp minimax. Further, it is shown that the suggested method is feasible for small samples.     

Comparison of linex and quadratic bayes estimators foe the rayleigh distribution

In this paper, the Bayes estimators for the parameter, the reliability function, and failure rate function of the Rayleigh distribution are obtained when based on complete or type II censored samples. Some types of the linex loss function are used. Comparieons in terms of risks of those under linex loss and squared error loss function with Bayes estimators relative to squared error loss function are made, Numerical example and simulation example are included.     

Subsampling-extrapolation bandwidth selection in bivariate kernel density estimation

This paper focuses on bivariate kernel density estimation that bridges the gap between univariate and multivariate applications. We propose a subsampling-extrapolation bandwidth matrix selector that improves the reliability of the conventional cross-validation method. The proposed procedure combines a U-statistic expression of the mean integrated squared error and asymptotic theory, and can be used in both cases of diagonal bandwidth matrix and unconstrained bandwidth matrix. In the subsampling stage, one takes advantage of the reduced variability of estimating the bandwidth matrix at a smaller subsample size m (m < n); in the extrapolation stage, a simple linear extrapolation is used to remove the incurred bias. Simulation studies reveal that the proposed method reduces the variability of the cross-validation method by about 50% and achieves an expected integrated squared error that is up to 30% smaller than that of the benchmark cross-validation. It shows comparable or improved performance compared to other competitors across six distributions in terms of the expected integrated squared error. We prove that the components of the selected bivariate bandwidth matrix have an asymptotic multivariate normal distribution, and also present the relative rate of convergence of the proposed bandwidth selector.     

A generalization of the compound rayleigh distribution: using a bayesian method on cancer survival times

In this paper the generalized compound Rayleigh model, exhibiting flexible hazard rate, is high¬lighted. This makes it attractive for modelling survival times of patients showing characteristics of a random hazard rate. The Bayes estimators are derived for the parameters of this model and some survival time parameters from a right censored sample. This is done with respect to conjugate and discrete priors on the parameters of this model, under the squared error loss function, Varian's asymmetric linear-exponential (linex) loss function and a weighted linex loss function. The future survival time of a patient is estimated under these loss functions. A Monte Carlo simu¬lation procedure is used where closed form expressions of the estimators cannot be obtained. An example illustrates the proposed estimators for this model.     

Stein rule prediction of the composite target function in a general linear regression model

This paper considers the problem of simultaneous prediction of the actual and average values of the dependent variable in a general linear regression model. Utilizing the philosophy of Stein rule procedure, a family of improved predictors for a linear function of the actual and expected value of the dependent variable for the forecast period has been proposed. An unbiased estimator for the mean squared error (MSE) matrix of the proposed family of predictors has been obtained and dominance of the family of Stein rule predictors over the best linear unbiased predictor (BLUP) has been established under a quadratic loss function.     

Estimation of the Stochastic Volatility Model by the Empirical Characteristic Function Method

The stochastic volatility model has no closed form for its likelihood and hence the maximum likelihood estimation method is difficult to implement. However, it can be shown that the model has a known characteristic function. As a consequence, the model is estimable via the empirical characteristic function. In this paper, the characteristic function of the model is derived and the estimation procedure is discussed. An application is considered for daily returns of Australian/New Zealand dollar exchange rate. Model checking suggests that the stochastic volatility model together with the empirical characteristic function estimates fit the data well.     

Nonparametric density estimation based on beta prime kernel

Abstract

In this work, we propose beta prime kernel estimator for estimation of a probability density functions defined with nonnegative support. For the proposed estimator, beta prime probability density function used as a kernel. It is free of boundary bias and nonnegative with a natural varying shape. We obtained the optimal rate of convergence for the mean squared error (MSE) and the mean integrated squared error (MISE). Also, we use adaptive Bayesian bandwidth selection method with Lindley approximation for heavy tailed distributions and compare its performance with the global least squares cross-validation bandwidth selection method. Simulation studies are performed to evaluate the average integrated squared error (ISE) of the proposed kernel estimator against some asymmetric competitors using Monte Carlo simulations. Moreover, real data sets are presented to illustrate the findings.     

Regression with an infinite number of observations applied to estimating the parameters of the stable distribution using the empirical characteristic function

ABSTRACT

The parameters of stable law parameters can be estimated using a regression based approach involving the empirical characteristic function. One approach is to use a fixed number of points for all parameters of the distribution to estimate the characteristic function. In this work the results are derived where all points in an interval is used to estimate the empirical characteristic function, thus least squares estimators of a linear function of the parameters, using an infinite number of observations. It was found that the procedure performs very good in small samples.     

Nonparametric smooth estimation of the expected inactivity time function

The expected inactivity time (EIT) function (also known as the mean past lifetime function) is a well known reliability function which has application in many disciplines such as survival analysis, actuarial studies and forensic science, to name but a few. In this paper, we use a fixed design local polynomial fitting technique to obtain estimators for the EIT function when the lifetime random variable has an unknown distribution. It will be shown that the proposed estimators are asymptotically unbiased, consistent and also, when standardized, has an asymptotic normal distribution. An optimal bandwidth, which minimizes the AMISE (asymptotic mean integrated squared error) of the estimator, is derived. Numerical examples based on simulated samples from various lifetime distributions common in reliability studies will be presented to evaluate the performances of these estimators. Finally, three real life applications will also be presented to further illustrate the wide applicability of these estimators.     

Estimation of the Spectral Density with Assigned Risk

It is well known that adaptive sequential nonparametric estimation of differentiable functions with assigned mean integrated squared error and minimax expected stopping time is impossible. In other words, no sequential estimator can compete with an oracle estimator that knows how many derivatives an estimated curve has. Differentiable functions are typical in probability density and regression models but not in spectral density models, where considered functions are typically smoother. This paper shows that for a large class of spectral densities, which includes spectral densities of classical autoregressive moving average processes, an adaptive minimax sequential estimation with assigned mean integrated squared error is possible. Furthermore, a two‐stage sequential procedure is proposed, which is minimax and adaptive to smoothness of an underlying spectral density.     

A note on Bayesian nonparametric regression function estimation

In this note the problem of nonparametric regression function estimation in a random design regression model with Gaussian errors is considered from the Bayesian perspective. It is assumed that the regression function belongs to a class of functions with a known degree of smoothness. A prior distribution on the given class can be induced by a prior on the coefficients in a series expansion of the regression function through an orthonormal system. The rate of convergence of the resulting posterior distribution is employed to provide a measure of the accuracy of the Bayesian estimation procedure defined by the posterior expected regression function. We show that the Bayes’ estimator achieves the optimal minimax rate of convergence under mean integrated squared error over the involved class of regression functions, thus being comparable to other popular frequentist regression estimators.     

On the effect of misspecifying the error density in a deconvolution problem

The author studies the effect of a misspecification of the error density on the mean integrated squared error (MISE) of the deconvolution estimator. He shows that the MISE converges to a certain functional which he defines. He also illustrates the fact that the limit can sometimes be infinite. Finally, he derives some guidelines for selecting the error density in order to ensure robustness properties of the procedure.