Estimation of a change point in the variance function based on the χ2-distribution

ABSTRACT

Let us consider that the variance function or its νth derivative in a regression model has a change/discontinuity point at an unknown location. To use the local polynomial fits, the log-variance function which break the positivity is targeted. The location and the jump size of the change point are estimated based on a one-sided kernel-weighted local-likelihood function which is provided by the χ²-distribution. The whole structure of the log-variance function is then estimated using the data sets split by the estimated location. Asymptotic results of the proposed estimators are described. Numerical works demonstrate the performances of the methods with simulated and real examples.     

A Note on Estimation of a Distribution Function in a Nonparametric Set-up Using Stein’s Shrinkage Estimation Technique

Based on Stein’s famous shrinkage estimation of a multivariate normal distribution, we propose a new type of estimators of the distribution function of a random variable in a nonparametric setup. The proposed estimators are then compared with the empirical distribution function, which is the best equivariant estimator under a well-known loss function. Our extensive simulation study shows that our proposed estimators can perform better for moderate to large sample sizes.     

Stein estimation—A review

This paper presents an expository development of Stein estimation in several distribution families. Considered are both the point estimation and confidence interval cases. Specific results for linear regression models are added. Emphasis is laid on the chronological history and on recent results.     

Minimum variance unbiased estimation of stress–strength reliability under bivariate normal and its comparisons

In many industrial and natural phenomena, we need the probability that a component is smaller than the other component. Under a stress–strength model, this is reliability of an item. Under independent setup, there are different approaches for the estimation of such reliability. Here, estimation is considered under the dependent case. Under bi-variate setup uniformly minimum variance unbiased estimator is obtained. Also comparison with available estimator based on Maximum Likelihood Estimate (MLE) is done through Mean Square Error (MSE) and bias. Also these are compared by computing L₁ distance between their distribution functions. From this idea and numerical computations, UMVUE appears to be good.     

Estimation of the parameters of a general student’s t distribution

We have observations for a t distribution with unknown mean, variance, and degrees of freedom, each of which we wish to estimate. The major problem lies in the estimate of the degrees of freedom. We show that a relatively efficient yet very simple estimator is a given function of the ratio of percentile estimates. We derive the appropriate estimator, provide equations for transformation and standard errors, contrast this with other estimators, and give examples.     

Nonparametric Estimation of Sojourn Time Distributions for Truncated Serial Event Data—a Weight-adjusted Approach

In follow-up studies, survival data often include subjects who have had a certain event at recruitment and may potentially experience a series of subsequent events during the follow-up period. This kind of survival data collected under a cross-sectional sampling criterion is called truncated serial event data. The outcome variables of interest in this paper are serial sojourn times between successive events. To analyze the sojourn times in truncated serial event data, we need to confront two potential sampling biases arising simultaneously from a sampling criterion and induced informative censoring. In this study, nonparametric estimation of the joint probability function of serial sojourn times is developed by using inverse probabilities of the truncation and censoring times as weight functions to accommodate these two sampling biases under various situations of truncation and censoring. Relevant statistical properties of the proposed estimators are also discussed. Simulation studies and two real data are presented to illustrate the proposed methods.     

Estimation of reliability of multicomponent stress–strength for a Kumaraswamy distribution

This article deals with the Bayesian and non Bayesian estimation of multicomponent stress–strength reliability by assuming the Kumaraswamy distribution. Both stress and strength are assumed to have a Kumaraswamy distribution with common and known shape parameter. The reliability of such a system is obtained by the methods of maximum likelihood and Bayesian approach and the results are compared using Markov Chain Monte Carlo (MCMC) technique for both small and large samples. Finally, two data sets are analyzed for illustrative purposes.     

An update on ‘a comparative study of tests for homogeneity of variance’

Tests for equality of variances using independent samples are widely used in data analysis. Conover et al. [A comparative study of tests for homogeneity of variance, with applications to the outer continental shelf bidding data. Technometrics. 1981;23:351–361], won the Youden Prize by comparing 56 variations of popular tests for variance on the basis of robustness and power in 60 different scenarios. None of the tests they compared were robust and powerful for the skewed distributions they considered. This study looks at 12 variations they did not consider, and shows that 10 are robust for the skewed distributions they considered plus the lognormal distribution, which they did not study. Three of these 12 have clearly superior power for skewed distributions, and are competitive in terms of robustness and power for all of the distributions considered. They are recommended for general use based on robustness, power, and ease of application.     

Local Linear Estimation of Recurrent Jump—Diffusion Models

In this article, we provide a nonparametric estimation of first and second infinitesimal moments of the underlying jump diffusion model. We show that under certain regularity conditions the nonparametric estimations of first and second infinitesimal moments based on the local linear estimator are consistent and asymptotically follow normal distributions.     

A simple generalization of the Box–Muller method for obtaining a pair of correlated standard normal variables

In the simple and widely used method of Box–Muller [G. Box and M. Muller, A note on the generation of random normal deviates, Ann. Math. Statist. 29 (1958), pp. 610–611], from a pair of uniform and independent random variables in (0,1), a pair of standard and independent normal variables is obtained. In this article, we present a very simple and elegant generalization of this method to obtain a pair of correlated standard normal variables with a given coefficient of correlation. This generalized method, which is computationally very easy, is interpreted in geometric terms, considering a translation of the uniform interval (0,1) and a rotation of a defined angle, both related to the coefficient of correlation. Some numerical results are simulated and statistically analysed, proving that the generalization is extremely simple and powerful.     

A Monte Carlo comparison of some variance estimators of the Horvitz–Thompson estimator

At the design and estimation stage of a survey, large survey organization often uses auxiliary information. This article discusses various procedures for improving variance estimation of the Horvitz–Thompson estimator of a finite population total with the aid of auxiliary information. To study the design-based properties of the proposed variance estimators relative to the standard one, a small scale Monte Carlo study is performed.     

Estimation of Sib–Sib Correlation via a Kotz-Type Density Function

We estimate sib–sib correlation by maximizing the log-likelihood of a Kotz-type distribution. Using extensive simulations we conclude that estimating sib–sib correlation using the proposed method has many advantages. Results are illustrated on a real life data set due to Galton. Testing of hypothesis about this correlation is also discussed using the three likelihood based tests and a test based on Srivastava's estimator. It is concluded that score test derived using Kotz-type density performs the best.     

Multivariate normal mean–variance mixture distribution based on Birnbaum–Saunders distribution

A multivariate normal mean–variance mixture based on a Birnbaum–Saunders (NMVMBS) distribution is introduced and several properties of this new distribution are discussed. A new robust non-Gaussian ARCH-type model is proposed in which there exists a relation between the variance of the observations, and the marginal distributions are NMVMBS. A simple EM-based maximum likelihood estimation procedure to estimate the parameters of this normal mean–variance mixture distribution is given. A simulation study and some real data are used to demonstrate the modelling strength of this new model.     

A comparison of bootstrap method and Edgeworth expansion in approximating the distribution of sample variance — one sample and two sample cases

The performance of the bootstrap method and the Edgeworth expansion in approximating the distribution of sample variance are compared when the data are from a non-normal population. Both approximations are very good. so long as the parent population is close to normal.     

Direction and collinearity factors of wilks‘s a - a review

Wilks's A was factorized by Bartlett (1951) for testing the hypothesis of goodness of fit of a hypothetical discriminant function in the case of several groups. This test has applications in various areas such as Econometrics, contingency tables, growth curves, principal components analysis, design of experiments and so on. This paper gives a consolidated account of the research done in these areas on the application of factors of Wilks's A     

Forecasting the Penetration of a New Product—A Bayesian Approach

We adopt a Bayesian approach to forecast the penetration of a new product into a market. We incorporate prior information from an existing product and/or management judgments into the data analysis. The penetration curve is assumed to be a nondecreasing function of time and may be under shape constraints. Markov-chain Monte Carlo methods are proposed and used to compute the Bayesian forecasts. An example on forecasting the penetration of color television using the information from black-and-white television is provided. The models considered can also be used to address the general bioassay and reliability stress-testing problems.     

Estimation and testing in generalized partial linear models—A comparative study

A particular semiparametric model of interest is the generalized partial linear model (GPLM) which extends the generalized linear model (GLM) by a nonparametric component.The paper reviews different estimation procedures based on kernel methods as well as test procedures on the correct specification of this model (vs. a parametric generalized linear model). Simulations and an application to a data set on East–West German migration illustrate similarities and dissimilarities of the estimators and test statistics.     

Estimation of the location of a 0-type or ∞-type singularity by Poisson observations

We consider an inhomogeneous Poisson process X on [0, T]. The intensity function of X is supposed to be strictly positive and smooth on [0, T] except at the point θ, in which it has either a 0-type singularity (tends to 0 like |x| ^p , p∈(0, 1)), or an ∞-type singularity (tends to ∞ like |x| ^p , p∈(?1, 0)). We suppose that we know the shape of the intensity function, but not the location of the singularity. We consider the problem of estimation of this location (shift) parameter θ based on n observations of the process X. We study the Bayesian estimators and, in the case p>0, the maximum-likelihood estimator. We show that these estimators are consistent, their rate of convergence is n ^1/(p+1), they have different limit distributions, and the Bayesian estimators are asymptotically efficient.     

Meta-analysis of Proportions of Rare Events–A Comparison of Exact Likelihood Methods with Robust Variance Estimation

The conventional random effects model for meta-analysis of proportions approximates within-study variation using a normal distribution. Due to potential approximation bias, particularly for the estimation of rare events such as some adverse drug reactions, the conventional method is considered inferior to the exact methods based on binomial distributions. In this article, we compare two existing exact approaches—beta binomial (B-B) and normal-binomial (N-B)—through an extensive simulation study with focus on the case of rare events that are commonly encountered in medical research. In addition, we implement the empirical (“sandwich”) estimator of variance into the two models to improve the robustness of the statistical inferences. To our knowledge, it is the first such application of sandwich estimator of variance to meta-analysis of proportions. The simulation study shows that the B-B approach tends to have substantially smaller bias and mean squared error than N-B for rare events with occurrences under 5%, while N-B outperforms B-B for relatively common events. Use of the sandwich estimator of variance improves the precision of estimation for both models. We illustrate the two approaches by applying them to two published meta-analysis from the fields of orthopedic surgery and prevention of adverse drug reactions.     

Estimation of parameters for a Birnbaum–Saunders regression model with censored data

Little work has been published on the analysis of censored data for the Birnbaum–Saunders distribution (BISA). In this article, we implement the EM algorithm to fit a regression model with censored data when the failure times follow the BISA. Three approaches to implement the E-Step of the EM algorithm are considered. In two of these implementations, the M-Step is attained by an iterative least-squares procedure. The algorithm is exemplified with a single explanatory variable in the model.