New results on the likelihood ratio and score tests for the von Mises distribution

In this paper, we derive Bartlett and Bartlett-type corrections [G.M. Cordeiro and S.L.P. Ferrari 1991, A modified score test statistic having chi-squared distribution to order n ^?1 , Biometrika 78 (1991), pp. 573–582] to improve the likelihood ratio and Rao's score statistics for testing the mean parameter and the concentration parameter in the von Mises distribution. Simple formulae are suggested for the corrections valid for small and large values of the concentration parameter that do not depend on the modified Bessel functions and can be useful in practical applications.     

Optimal B-robust estimators for the parameters of the Burr XII distribution

The Burr XII distribution offers a flexible alternative to the distributions that play important role for modelling data in reliability, risk and process capability. However, estimating the shape parameters of the Burr XII distribution is a challenging problem. The classical estimation methods such as maximum likelihood and least squares are often used to estimate the parameters of the Burr XII distribution, but these methods are very sensitive to the outliers in the data. Thus, a robust estimation method alternative to the classical methods is needed to find robust estimators that are less sensitive to the outliers in the data. The purpose of this paper is to use the optimal B-robust estimation method [Hampel FR, Ronchetti EM, Rousseeuw PJ, Stahel WA. Robust statistics: the approach based on influence functions. New York: Wiley; 1986] to obtain robust estimators for the shape parameters of the Burr XII distribution. The simulation results show that the optimal B-robust estimators generally outperform the classical estimators in terms of the bias and root mean square errors when there are outliers in data.     

商誉和价格影响需求的易逝品动态定价与广告决策

研究需求受价格和商誉共同影响下易逝品动态定价与广告投资策略的联合决策问题。在初始库存确定的情形下,在某待定销售周期内,企业通过投入广告提升自身商誉从而拓宽市场,以企业利润最大化为目标,建立联合最优动态定价及广告投资策略模型。应用庞特里亚金极大值原理,求得不同系统参数下最优动态定价和广告投资策略及最优的销售周期。通过数值算例验证方法的有效性,结合参数灵敏度分析检验系统参数对最优策略的影响并给出相应管理意义。     

Reliability of stress–strength model for exponentiated Pareto distributions

In this paper, the reliability of a system is discussed when the strength of the system and the stress imposed on it are independent, non-identical exponentiated Pareto distributed random variables. Different point estimations and interval estimations are proposed. The point estimators obtained are maximum likelihood, uniformly minimum variance unbiased and Bayesian estimators. The interval estimations obtained are approximate, exact, bootstrap-p and bootstrap-t confidence intervals and Bayesian credible interval. Different methods and the corresponding confidence intervals are compared using Monte-carlo simulations.     

Estimation of variance in bivariate normal distribution after the preliminary test of homogeneity

A preliminary test estimator of variance in the bivariate normal distribution is proposed after the Pitman–Morgan test of homogeneity of two variances. The bias and mean square error of the estimator are derived. The relative efficiency (RE) of the preliminary test estimator is studied. Computations and 3D graphs of RE for different parameters are analyzed. In order to get the maximum RE, recommendations of the significance level for the preliminary test are given for various sample sizes by using the max–min criterion.     

Limit theorems for dependent Bernoulli variables with statistical inference

We consider a class of dependent Bernoulli variables where the conditional success probability is a linear combination of the last few trials and the original success probability. We obtain its limit theorems including the strong law of large numbers, weak invariance principle, and law of the iterated logarithm. We also derive some statistical inference results which make the model applicable. Simulation results are exhibited as well to show that with small sample size the convergence rate is satisfying and the proposed estimators behave well.     

Higher-order Bayesian Approximations for Pseudo-posterior Distributions

The theory of higher-order asymptotics provides accurate approximations to posterior distributions for a scalar parameter of interest, and to the corresponding tail area, for practical use in Bayesian analysis. The aim of this article is to extend these approximations to pseudo-posterior distributions, e.g., posterior distributions based on a pseudo-likelihood function and a suitable prior, which are proved to be particularly useful when the full likelihood is analytically or computationally infeasible. In particular, from a theoretical point of view, we derive the Laplace approximation for a pseudo-posterior distribution, and for the corresponding tail area, for a scalar parameter of interest, also in the presence of nuisance parameters. From a computational point of view, starting from these higher-order approximations, we discuss the higher-order tail area (HOTA) algorithm useful to approximate marginal posterior distributions, and related quantities. Compared to standard Markov chain Monte Carlo methods, the main advantage of the HOTA algorithm is that it gives independent samples at a negligible computational cost. The relevant computations are illustrated by two examples.     

Improving the detection of unusual observations in high‐dimensional settings

Multivariate control charts are used to monitor stochastic processes for changes and unusual observations. Hotelling's T² statistic is calculated for each new observation and an out‐of‐control signal is issued if it goes beyond the control limits. However, this classical approach becomes unreliable as the number of variables p approaches the number of observations n, and impossible when p exceeds n. In this paper, we devise an improvement to the monitoring procedure in high‐dimensional settings. We regularise the covariance matrix to estimate the baseline parameter and incorporate a leave‐one‐out re‐sampling approach to estimate the empirical distribution of future observations. An extensive simulation study demonstrates that the new method outperforms the classical Hotelling T² approach in power, and maintains appropriate false positive rates. We demonstrate the utility of the method using a set of quality control samples collected to monitor a gas chromatography–mass spectrometry apparatus over a period of 67 days.     

A method for choosing the smoothing parameter in a semi-parametric model for detecting change-points in blood flow

In a smoothing spline model with unknown change-points, the choice of the smoothing parameter strongly influences the estimation of the change-point locations and the function at the change-points. In a tumor biology example, where change-points in blood flow in response to treatment were of interest, choosing the smoothing parameter based on minimizing generalized cross-validation (GCV) gave unsatisfactory estimates of the change-points. We propose a new method, aGCV, that re-weights the residual sum of squares and generalized degrees of freedom terms from GCV. The weight is chosen to maximize the decrease in the generalized degrees of freedom as a function of the weight value, while simultaneously minimizing aGCV as a function of the smoothing parameter and the change-points. Compared with GCV, simulation studies suggest that the aGCV method yields improved estimates of the change-point and the value of the function at the change-point.     

Convergence rate of wavelet density estimator with data missing randomly when covariables are present

In this article, we study global L₂ error of non linear wavelet estimator of density in the Besov space B^s_pq for missing data model when covariables are present and prove that the estimator can achieve the optimal rate of convergence, which is similar to the result studied by Donoho et al. (1996) Donoho, D.L., Johnstone, I.M., Kerkyacharian, G., Picard, D. (1996). Density estimation by wavelet thresholding. Ann. Stat. 24:508–539.[Crossref], [Web of Science ®] , [Google Scholar] in complete independent data case with term-by-term thresholding of the empirical wavelet coefficients. Finite-sample behavior of the proposed estimator is explored via simulations.