A note on the equality of the OLSE and the BLUE of the parametric function in the general Gauss–Markov model

In this note we consider the equality of the ordinary least squares estimator (OLSE) and the best linear unbiased estimator (BLUE) of the estimable parametric function in the general Gauss–Markov model. Especially we consider the structures of the covariance matrix V for which the OLSE equals the BLUE. Our results are based on the properties of a particular reparametrized version of the original Gauss–Markov model.      

Response to the “letter to the editor”

First of all, the parameters involved in Chen and Hsu (1995) are unknown. For convenience, we discussed the two possibilities of the unknown parameter μ namely μ =M and μ M, separately. In other words, we do not known wheter μ = M or μ M. Therefore, we must estimate it before drawing any statistical conclusion about the index C_pmk, The conclusion should be applied only after th ehypothesis if μ = M is tested.     

Estimation of the Parameters of the Birnbaum–Saunders Distribution

Some alternative estimators to the maximum likelihood estimators of the two parameters of the Birnbaum–Saunders distribution are proposed. Most have high efficiencies as measured by root mean square error and are robust to departure from the model as well as to outliers. In addition, the proposed estimators are easy to compute. Both complete and right-censored data are discussed. Simulation studies are provided to compare the performance of the estimators.     

On the Restricted Liu Estimator in the Gauss–Markov Model

In this paper, we compare two estimators, the RLE (restricted Liu estimator) and the RLSE (restricted least squares estimator) of parameters in linear models under Gauss–Markov models. Using generalized inverse of matrices, we found some equivalency conditions for the superiority of the RLE with respect to the MSE criterion.     

Increasing the Power of the Test Through Pre-Test – A Robust Method

Robust test procedures are developed for testing the intercept of a simple regression model when the slope is (i) completely unspecified, (ii) specified to a fixed value, or (iii) suspected to be a fixed value. Defining (i) unrestricted (UT), (ii) restricted (RT), and (iii) pre-test test (PTT) functions for the intercept parameter under the three choices of the slope, tests are formulated using the M-estimation methodology. The asymptotic distributions of the test statistics and their asymptotic power functions are derived. The analytical and graphical comparisons of the tests reveal that the PTT achieves a reasonable dominance over the other tests.     

Distribution of the mean reversion estimator in the Ornstein–Uhlenbeck process

ABSTRACT

We derive the exact distribution of the maximum likelihood estimator of the mean reversion parameter (κ) in the Ornstein–Uhlenbeck process using numerical integration through analytical evaluation of a joint characteristic function. Different scenarios are considered: known or unknown drift term, fixed or random start-up value, and zero or positive κ. Monte Carlo results demonstrate the remarkably reliable performance of our exact approach across all the scenarios. In comparison, misleading results may arise under the asymptotic distributions, including the advocated infill asymptotic distribution, which performs poorly in the tails when there is no intercept in the regression and the starting value of the process is nonzero.     

How risky is the optimal portfolio which maximizes the Sharpe ratio?

In this paper, we investigate the properties of the optimal portfolio in the sense of maximizing the Sharpe ratio (SR) and develop a procedure for the calculation of the risk of this portfolio. This is achieved by constructing an optimal portfolio which minimizes the Value-at-Risk (VaR) and at the same time coincides with the tangent (market) portfolio on the efficient frontier which is related to the SR portfolio. The resulting significance level of the minimum VaR portfolio is then used to determine the risk of both the market portfolio and the corresponding SR portfolio. However, the expression of this significance level depends on the unknown parameters which have to be estimated in practice. It leads to an estimator of the significance level whose distributional properties are investigated in detail. Based on these results, a confidence interval for the suggested risk measure of the SR portfolio is constructed and applied to real data. Both theoretical and empirical findings document that the SR portfolio is very risky since the corresponding significance level is smaller than 90 % in most of the considered cases.     

Linear sufficiency and completeness in the context of estimating the parametric function in the general Gauss–Markov model

In this paper we consider linear sufficiency and linear completeness in the context of estimating the estimable parametric function K^′β${\mathbf{K}}^{\prime }\mathit{\beta }$ under the general Gauss–Markov model {y,Xβ,σ²V}$\{\mathbf{y},\mathbf{X}\mathit{\beta },{\sigma }^2\mathbf{V}\}$. We give new characterizations for linear sufficiency, and define and characterize linear completeness in a case of estimation of K^′β${\mathbf{K}}^{\prime }\mathit{\beta }$. Also, we consider a predictive approach for obtaining the best linear unbiased estimator of K^′β${\mathbf{K}}^{\prime }\mathit{\beta }$, and subsequently, we give the linear analogues of the Rao–Blackwell and Lehmann–Scheffé Theorems in the context of estimating K^′β${\mathbf{K}}^{\prime }\mathit{\beta }$.     

A note on the posterior inference for the Yule–Simon distribution

The Yule–Simon distribution has been out of the radar of the Bayesian community, so far. In this note, we propose an explicit Gibbs sampling scheme when a Gamma prior is chosen for the shape parameter. The performance of the algorithm is illustrated with simulation studies, including count data regression, and a real data application to text analysis. We compare our proposal to the frequentist counterparts showing better performance of our algorithm when a small sample size is considered.     

Bootstrapping the augmented Dickey–Fuller test for unit root using the MDIC

In this paper, we consider the bootstrap procedure for the augmented Dickey–Fuller (ADF) unit root test by implementing the modified divergence information criterion (MDIC, Mantalos et al. [An improved divergence information criterion for the determination of the order of an AR process, Commun. Statist. Comput. Simul. 39(5) (2010a), pp. 865–879; Forecasting ARMA models: A comparative study of information criteria focusing on MDIC, J. Statist. Comput. Simul. 80(1) (2010b), pp. 61–73]) for the selection of the optimum number of lags in the estimated model. The asymptotic distribution of the resulting bootstrap ADF/MDIC test is established and its finite sample performance is investigated through Monte-Carlo simulations. The proposed bootstrap tests are found to have finite sample sizes that are generally much closer to their nominal values, than those tests that rely on other information criteria, like the Akaike information criterion [H. Akaike, Information theory and an extension of the maximum likelihood principle, in Proceedings of the 2nd International Symposium on Information Theory, B.N. Petrov and F. Csáki, eds., Akademiai Kaido, Budapest, 1973, pp. 267–281]. The simulations reveal that the proposed procedure is quite satisfactory even for models with large negative moving average coefficients.     

A Cautionary Note on the Use of the Grønnesby and Borgan Goodness-of-Fit Test for the Cox Proportional Hazards Model

Grønnesby and Borgan (1996, Lifetime Data Analysis 2, 315–328) propose an omnibus goodness-of-fit test for the Cox proportional hazards model. The test is based on grouping the subjects by their estimated risk score and comparing the number of observed and a model based estimated number of expected events within each group. We show, using extensive simulations, that even for moderate sample sizes the choice of number of groups is critical for the test to attain the specified size. In light of these results we suggest a grouping strategy under which the test attains the correct size even for small samples. The power of the test statistic seems to be acceptable when compared to other goodness-of-fit tests.     

On selecting from k finite populations the population with the largest α -quantile

Several procedures for ranking populations according to the quantile of a given order have been discussed in the literature. These procedures deal with continuous distributions. This paper deals with the problem of selecting a population with the largest α-quantile from k ≥ 2 finite populatins, where the size of each population is known. A selection rule is given based on the sample quantiles, where he samples are drawn without replacement. A formula for the minimum probability of a correct selection for the given rule, for a certain configuration of the population α-quantiles, is given in terms of the sample numbers.     

Bayesian Model Selection: Measuring the χ2 Discrepancy with the Uniform Distribution

In the last years, many articles have been written about Bayesian model selection. In this article, a different and easier method is proposed and analyzed. The key idea of this method is based on the well-known property that, under the true model, the cumulative distribution function is distributed as a uniform distribution over the interval (0, 1). The method is first introduced for the continuous case and then for the discrete case by smoothing the cumulative distribution function. Some asymptotical properties of the method are obtained by developing an alternative to Helly's theorems. Finally, the performance of the method is evaluated by simulation, showing a good behavior.     

Computing the distribution of sequential Hölder norms of the Brownian motion

ABSTRACT

In this article, we give explicit formulas and study practical computations for the distribution function of sequential Hölder norms of a Brownian motion and of a Brownian bridge. We also discuss some statistical applications in the detection of some short “epidemic” changes in a sample.     

On the Behrens–Fisher problem from the spatial median point of view

We present a non-parametric affine-invariant test for the multivariate Behrens–Fisher problem. The proposed method based on the spatial medians is asymptotic and does not require normality of the data. To improve its finite sample performance, we apply a correction of the type which was already used in a similar test based on trimmed means, however, our simulations show that in the case of heavy-tailed distributions our method performs better. Also in a simulation comparison with a recently published rank-based test our test yields satisfactory results.     

Is the Selected Population the Best?—Location and Scale Parameter Cases

We observe X ₁,…,X _k , where X _i has density f(x,θ _i ) possessing monotone likelihood ratio. The best population corresponds to the largest θ _i . We select the population corresponding to the largest X _i . The goal is to attach the best possible p-value to the inference: the selected population has the uniquely largest θ _i . Gutmann and Maymin (1987 Gutmann , S. , Maymin , Z. ( 1987 ). Is the selected population the best? Ann. Statist . 15 : 456 – 461 .[Crossref], [Web of Science ®] , [Google Scholar]) considered the location parameter case and derived the supremum of the error probability by conditioning on S, the index of the largest X _i . Using this conditioning approach, Kannan and Panchapakesan (2009 Kannan , N. , Panchapakesan , S. ( 2009 ). Does the selected normal population have the smallest variance? Amer. J. Math. Management Sci . 29 : To appear . [Google Scholar]) considered the problem for the gamma family. We consider here a unified approach to both the location and scale parameter cases, and obtain the supremum of the error probability without using conditioning.