Bayesian Estimation for Linear Functions of Poisson Rates

The probability matching prior for linear functions of Poisson parameters is derived. A comparison is made between the confidence intervals obtained by Stamey and Hamilton (2006 Stamey, J., Hamilton, C. (2006). A note on confidence intervals for a linear function of Poisson rates. Commun. Statist. Simul. &; Computat. 35(4):849–856.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), and the intervals derived by us when using the Jeffreys’ and probability matching priors. The intervals obtained from the Jeffreys’ prior are in some cases fiducial intervals (Krishnamoorthy and Lee, 2010 Krishnamoorthy, K., Lee, M. (2010). Inference for functions of parameters in discrete distributions based on fiducial approach: Binomial and Poisson cases. J. Statist. Plann. Infere. 140(5):1182–1192.[Crossref], [Web of Science ®] , [Google Scholar]). A weighted Monte Carlo method is used for the probability matching prior. The power and size of the test, using Bayesian methods, is compared to tests used by Krishnamoorthy and Thomson (2004 Krishnamoorthy, K., Thomson, J. (2004). A more powerful test for comparing two Poisson means. J. Statist. Plann. Infere. 119(1):23–35.[Crossref], [Web of Science ®] , [Google Scholar]). The Jeffreys’, probability matching and two other priors are used.     

Credibility premium for rate-making systems

Credibility formula has been developed in many fields of actuarial sciences. Based upon Payandeh (2010 Payandeh, A.T. (2010). A new approach to the credibility formula. Insur.: Math. Econ. 46(2):334–338.[Crossref], [Web of Science ®] , [Google Scholar]), this article extends concept of credibility formula to relatively premium of a given rate-making system. More precisely, it calculates Payandeh’s (2010 Payandeh, A.T. (2010). A new approach to the credibility formula. Insur.: Math. Econ. 46(2):334–338.[Crossref], [Web of Science ®] , [Google Scholar]) credibility factor for zero-inflated Poisson gamma distributions with respect to several loss functions. A comparison study has been given.     

An objective Bayesian estimation in a two-period crossover design

Two-period crossover design is one of the commonly used designs in clinical trials. But, the estimation of treatment effect is complicated by the possible presence of carryover effect. It is known that ignoring the carryover effect when it exists can lead to poor estimates of the treatment effect. The classical approach by Grizzle (1965 Grizzle, J.E. (1965). The two-period change-over design and its use in clinical trials. Biometrics 21:467–480. See Grizzle (1974) for corrections.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) consists of two stages. First, a preliminary test is conducted on carryover effect. If the carryover effect is significant, analysis is based only on data from period one; otherwise, analysis is based on data from both periods. A Bayesian approach with improper priors was proposed by Grieve (1985 Grieve, A.P. (1985). A Bayesian analysis of the two-period crossover design for clinical trials. Biometrics 41:979–990.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) which uses a mixture of two models: a model with carryover effect and another without. The indeterminacy of the Bayes factor due to the arbitrary constant in the improper prior was addressed by assigning a minimally discriminatory value to the constant. In this article, we present an objective Bayesian estimation approach to the two-period crossover design which is also based on a mixture model, but using the commonly recommended Zellner–Siow g-prior. We provide simulation studies and a real data example and compare the numerical results with Grizzle (1965 Grizzle, J.E. (1965). The two-period change-over design and its use in clinical trials. Biometrics 21:467–480. See Grizzle (1974) for corrections.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])’s and Grieve (1985 Grieve, A.P. (1985). A Bayesian analysis of the two-period crossover design for clinical trials. Biometrics 41:979–990.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])’s approaches.     

Buffered Autoregressive Models With Conditional Heteroscedasticity: An Application to Exchange Rates

This article introduces a new model called the buffered autoregressive model with generalized autoregressive conditional heteroscedasticity (BAR-GARCH). The proposed model, as an extension of the BAR model in Li et al. (2015 Li, G.D., Guan, B., Li, W.K., and Yu, P. L.H. (2015), “Hysteretic Autoregressive Time Series Models,” Biometrika, 102, 717–723.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), can capture the buffering phenomena of time series in both the conditional mean and variance. Thus, it provides us a new way to study the nonlinearity of time series. Compared with the existing AR-GARCH and threshold AR-GARCH models, an application to several exchange rates highlights the importance of the BAR-GARCH model.     

A Bayesian control chart for a common coefficient of variation

By using the medical data analyzed by Kang et al. (2007 Kang, C.W., Lee, M.S., Seong, Y.J., Hawkins, D.M. (2007). A control chart for the coefficient of variation. J. Qual. Technol. 39(2):151–158.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), a Bayesian procedure is applied to obtain control limits for the coefficient of variation. Reference and probability matching priors are derived for a common coefficient of variation across the range of sample values. By simulating the posterior predictive density function of a future coefficient of variation, it is shown that the control limits are effectively identical to those obtained by Kang et al. (2007 Kang, C.W., Lee, M.S., Seong, Y.J., Hawkins, D.M. (2007). A control chart for the coefficient of variation. J. Qual. Technol. 39(2):151–158.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) for the specific dataset they used. This article illustrates the flexibility and unique features of the Bayesian simulation method for obtaining posterior distributions, predictive intervals, and run-lengths in the case of the coefficient of variation. A simulation study shows that the 95% Bayesian confidence intervals for the coefficient of variation have the correct frequentist coverage.     

Weighted geometric distribution with new characterizations of geometric distribution

In this article, we derive a new generalized geometric distribution through a weight function, which can also be viewed as a discrete analog of weighted exponential distribution introduced by Gupta and Kundu (2009 Gupta, R. D., and D. Kundu. 2009. A new class of weighted exponential distributions. Statistics 43:621–34.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). We derive some distributional properties like moments, generating functions, hazard function, and infinite divisibility followed by different estimation methods to estimate the parameters. New characterizations of the geometric distribution are presented using the proposed generalized geometric distribution. The superiority of the proposed distribution to other competing models is demonstrated with the help of two real count datasets.     

Linear transformation models for survival analysis with tumor growth information in cancer screening study

The complication in analyzing tumor data is that the tumors detected in a screening program tend to be slowly progressive tumors, which is the so-called left-truncated sampling that is inherent in screening studies. Under the assumption that all subjects have the same tumor growth function, Ghosh (2008 Ghosh, D. (2008). Proportional hazards regression for cancer studies. Biometrics 64:141–148.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) developed estimation procedures for the Cox proportional hazards model. Shen (2011a Shen, P.-S. (2011a). Proportional hazards regression for cancer screening data. J. Stat. Comput. Simul. 18:367–377.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) demonstrated that Ghosh (2008 Ghosh, D. (2008). Proportional hazards regression for cancer studies. Biometrics 64:141–148.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])'s approach can be extended to the case when each subject has a specific growth function. In this article, under linear transformation model, we present a general framework to the analysis of data from cancer screening studies. We developed estimation procedures under linear transformation model, which includes Cox's model as a special case. A simulation study is conducted to demonstrate the potential usefulness of the proposed estimators.     

Extrapolation techniques in U-statistic variance estimation

This article considers the problem of variance estimation of a U-statistic. Following the proposal of a linearly extrapolated variance estimator in Wang and Chen (2015 Wang, Q., Chen, S. (2015). A general class of linearly extrapolated variance estimators. Stat. Probab. Lett. 98:29–38.[Crossref], [Web of Science ®] , [Google Scholar]), we consider a second-order extrapolation technique and devise a variance estimator that is nearly second-order unbiased. Simulation studies confirm that the second-order extrapolated variance estimator has smaller bias than the linearly extrapolated variance estimator and the jackknife variance estimator across a wide selection of distributions. In addition, the proposal also yields a smaller mean squared error than its counterparts. In the end, we discuss the advantages of the proposed variance estimator in regression analysis and model selection.     

Non parametric confidence intervals for the extinction probability of a Galton–Watson branching process

It is demonstrated that the confidence intervals (CIs) for the probability of eventual extinction and other parameters of a Galton–Watson branching process based upon the maximum likelihood estimators can often have substantially lower coverage when compared to the desired nominal confidence coefficient, especially in small, more realistic sample sizes. The same conclusion holds for the traditional bootstrap CIs. We propose several adjustments to these CIs, which greatly improves coverage in most cases. We also make a correction in an asymptotic variance formula given in Stigler (1971 Stigler, S.M. (1971). The estimation of the probability of extinction and other parameters associated with branching processes. Biometrika 58(3):499–508.[Crossref], [Web of Science ®] , [Google Scholar]). The focus here is on implementation of the CIs which have good coverage, in a wide variety of cases. We also consider expected CI lengths. Some recommendations are made.     

Bimodal symmetric-asymmetric power-normal families

This article proposes new symmetric and asymmetric distributions applying methods analogous as the ones in Kim (2005 Kim, H.J. (2005). On a class of two-piece skew-normal distributions. Statist.: J. Theoret. Appl. Statist. 39:537–553.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and Arnold et al. (2009 Arnold, B.C., H.W. Gómez, and H.S. Salinas. (2009). On multiple constraint skewed models. Statist. J. Theoret. Appl. Statist. 43: 279–293.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) to the exponentiated normal distribution studied in Durrans (1992 Durrans, S.R. (1992). Distributions of fractional order statistics in hydrology. Water Resour. Res. 28:1649–1655.[Crossref], [Web of Science ®] , [Google Scholar]), that we call the power-normal (PN) distribution. The proposed bimodal extension, the main focus of the paper, is called the bimodal power-normal model and is denoted by BPN(α) model, where α is the asymmetry parameter. The authors give some properties including moments and maximum likelihood estimation. Two important features of the model proposed is that its normalizing constant has closed and simple form and that the Fisher information matrix is nonsingular, guaranteeing large sample properties of the maximum likelihood estimators. Finally, simulation studies and real applications reveal that the proposed model can perform well in both situations.     

Bahadur Representation of Linear Kernel Quantile Estimator for Stationary Processes

In this article, we discuss the method of linear kernel quantile estimator proposed by Parzen (1979 Parzen, E. (1979). Nonparametric statistical data modeling. J. Amer. Statist. Assoc. 74:105121.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). We establish a Bahadur representation in sense of almost surely convergence with the rate log^{? α}n under the case of S-mixing random variable sequence which was proposed by Berkes (2009 Berkes, I., Hörmann, S., (2009). Asymptotic results for the itpirical process of stationary sequences. Stoch. Process. Their Applic. 119:12981324.[Crossref], [Web of Science ®] , [Google Scholar]). We also obtain the strong consistence of this estimator and its convergence rate.     

On the optimal allocation of active redundancies in series system

This paper studies the allocations of two non identical active redundancies in series systems in terms of the reversed hazard rate order and hazard rate order, which generalizes some results built in Valdés and Zequeira (2003 Valdés, J. E., and R. I. Zequeira 2003. On the optimal allocation of an active redundancy in a two-component series system. Stat. Probab. Lett. 63:325–32.[Crossref], [Web of Science ®] , [Google Scholar], 2006 Valdés, J. E., and R. I. Zequeira 2006. On the optimal allocation of two active redundancies in a two-component series system. Oper. Res. Lett. 34:49–52.[Crossref], [Web of Science ®] , [Google Scholar]).     

A new efficient unrelated randomized response model

The present paper suggests an interesting and useful ramification of the unrelated randomized response model due to Pal and Singh (2012 Pal, S., and S. Singh. 2012. A new unrelated question randomized response model. Statistics 46 (1):99–109.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) [A new unrelated question randomized response model. Statistics 46 (1), 99–109] that can be used for any sampling scheme. We have shown theoretically and numerically that the proposed model is more efficient than Pal and Singh (2012 Pal, S., and S. Singh. 2012. A new unrelated question randomized response model. Statistics 46 (1):99–109.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) model.     

On weighted cumulative residual entropy

In analogy with the weighted Shannon entropy proposed by Belis and Guiasu (1968 Belis, M., Guiasu, S. (1968). A quantitative-qualitative measure of information in cybernetic systems. IEEE Trans. Inf. Th. IT-4:593–594.[Crossref], [Web of Science ®] , [Google Scholar]) and Guiasu (1986 Guiasu, S. (1986). Grouping data by using the weighted entropy. J. Stat. Plann. Inference 15:63–69.[Crossref], [Web of Science ®] , [Google Scholar]), we introduce a new information measure called weighted cumulative residual entropy (WCRE). This is based on the cumulative residual entropy (CRE), which is introduced by Rao et al. (2004 Rao, M., Chen, Y., Vemuri, B.C., Wang, F. (2004). Cumulative residual entropy: a new measure of information. IEEE Trans. Info. Theory 50(6):1220–1228.[Crossref], [Web of Science ®] , [Google Scholar]). This new information measure is “length-biased” shift dependent that assigns larger weights to larger values of random variable. The properties of WCRE and a formula relating WCRE and weighted Shannon entropy are given. Related studies of reliability theory is covered. Our results include inequalities and various bounds to the WCRE. Conditional WCRE and some of its properties are discussed. The empirical WCRE is proposed to estimate this new information measure. Finally, strong consistency and central limit theorem are provided.     

Optimal allocations of coverage limits for two independent random losses of insurance policy

This article recasts the optimal allocations of coverage limits for two independent random losses. Under some regularity conditions on the two concerned probability density functions, we build the sufficient and necessary condition for the existence of the optimal allocation of coverage limits, and derive the optimal allocation whenever they do exist. The results supplement Lu and Meng (2011 Lu, Z.Y., Meng, L.L. (2011). Stochastic comparisons for allocations of upper limits and deductibles with applications. Insur.: Math. Econ. 48:338–343.[Crossref], [Web of Science ®] , [Google Scholar], Proposition 5.2) and Hu and Wang (2014 Hu, S., Wang, R. (2014). Stochastic comparisons and optimal allocation for policy limits and deductibles. Commun. Stat. – Theory Methods 43:151–164.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar], Theorem 5.1).     

An extension of almost sure central limit theorem for the maximum of stationary Gaussian random fields

Let X₁, X₂, … be a sequence of stationary standardized Gaussian random fields. The almost sure limit theorem for the maxima of stationary Gaussian random fields is established. Our results extend and improve the results in Csáki and Gonchigdanzan (2002 Csáki, E., Gonchigdanzan, K. (2002). Almost sure limit theorems for the maximum of stationary Gaussian sequences. Stat. Probab. Lett. 58:195–203.[Crossref], [Web of Science ®] , [Google Scholar]) and Choi (2010 Choi, H. (2010). Almost sure limit theorem for stationary Gaussian random fields. J. Korean Stat. Soc. 39:449–454.[Crossref], [Web of Science ®] , [Google Scholar]).     

Improved Exponential Type Estimators of the Mean of a Sensitive Variable in the Presence of Nonsensitive Auxiliary Information

Recently, Koyuncu et al. (2013 Koyuncu, N., Gupta, S., Sousa, R. (2014). Exponential type estimators of the mean of a sensitive variable in the presence of non-sensitive auxiliary information. Communications in Statistics- Simulation and Computation[PubMed], [Web of Science ®] , [Google Scholar]) proposed an exponential type estimator to improve the efficiency of mean estimator based on randomized response technique. In this article, we propose an improved exponential type estimator which is more efficient than the Koyuncu et al. (2013 Koyuncu, N., Gupta, S., Sousa, R. (2014). Exponential type estimators of the mean of a sensitive variable in the presence of non-sensitive auxiliary information. Communications in Statistics- Simulation and Computation[PubMed], [Web of Science ®] , [Google Scholar]) estimator, which in turn was shown to be more efficient than the usual mean estimator, ratio estimator, regression estimator, and the Gupta et al. (2012 Gupta, S., Shabbir, J., Sousa, R., Corte-Real, P. (2012). Regression estimation of the mean of a sensitive variable in the presence of auxiliary information. Communications in Statistics – Theory and Methods 41:2394–2404.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) estimator. Under simple random sampling without replacement (SRSWOR) scheme, bias and mean square error expressions for the proposed estimator are obtained up to first order of approximation and comparisons are made with the Koyuncu et al. (2013 Koyuncu, N., Gupta, S., Sousa, R. (2014). Exponential type estimators of the mean of a sensitive variable in the presence of non-sensitive auxiliary information. Communications in Statistics- Simulation and Computation[PubMed], [Web of Science ®] , [Google Scholar]) estimator. A simulation study is used to observe the performances of these two estimators. Theoretical findings are also supported by a numerical example with real data. We also show how to, extend the proposed estimator to the case when more than one auxiliary variable is available.     

Matrix mean squared error comparisons of some biased estimators with two biasing parameters

To deal with multicollinearity problem, the biased estimators with two biasing parameters have recently attracted much research interest. The aim of this article is to compare one of the last proposals given by Yang and Chang (2010 Yang, H., and X. Chang. 2010. A new two-parameter estimator in linear regression. Communications in Statistics: Theory and Methods 39 (6):923–34.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) with Liu-type estimator (Liu 2003 Liu, K. 2003. Using Liu-type estimator to combat collinearity. Communications in Statistics: Theory and Methods 32 (5):1009–20.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and k ? d class estimator (Sakallioglu and Kaciranlar 2008 Sakallioglu, S., and S. Kaciranlar. 2008. A new biased estimator based on ridge estimation. Statistical Papers 49:669–89.[Crossref], [Web of Science ®] , [Google Scholar]) under the matrix mean squared error criterion. As well as giving these comparisons theoretically, we support the results with the extended simulation studies and real data example, which show the advantages of the proposal given by Yang and Chang (2010 Yang, H., and X. Chang. 2010. A new two-parameter estimator in linear regression. Communications in Statistics: Theory and Methods 39 (6):923–34.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) over the other proposals with increasing multicollinearity level.     

Successive comparisons between ordered normal means using isotonic regression estimators

In this study we discuss multiple comparison procedures for checking differences among a sequence of normal means with ordered restriction. Lee and Spurrier (1995 Lee, R.E., Spurrier, J.D. (1995). Successive comparisons between ordered treatments. J. Stat. Plann. Inference 43:323–330.[Crossref], [Web of Science ®] , [Google Scholar]) proposed a multiple comparison procedure which tests the difference between two adjacent means using the difference of sample means. In this study we propose a multiple comparison procedure modifying Lee and Spurrier's (1995 Lee, R.E., Spurrier, J.D. (1995). Successive comparisons between ordered treatments. J. Stat. Plann. Inference 43:323–330.[Crossref], [Web of Science ®] , [Google Scholar]) procedure using isotonic regression estimators instead of sample means. We determine the critical value for pairwise comparisons for a specified significance level. Furthermore, we formulate the power of the test. Finally, we give some numerical examples regarding critical values and power of the test intended to compare our procedure with Lee and Spurrier's (1995 Lee, R.E., Spurrier, J.D. (1995). Successive comparisons between ordered treatments. J. Stat. Plann. Inference 43:323–330.[Crossref], [Web of Science ®] , [Google Scholar]) procedure.     

Shared inverse Gaussian frailty models based on additive hazards

Frailty models are used in the survival analysis to account for the unobserved heterogeneity in individual risks to disease and death. To analyze the bivariate data on related survival times (e.g., matched pairs experiments, twin, or family data), the shared frailty models were suggested. These models are based on the assumption that frailty acts multiplicatively to hazard rate. In this article, we assume that frailty acts additively to hazard rate. We introduce the shared inverse Gaussian frailty models with three different baseline distributions, namely the generalized log-logistic, the generalized Weibull, and exponential power distribution. We introduce the Bayesian estimation procedure using Markov chain Monte Carlo technique to estimate the parameters involved in these models. We apply these models to a real-life bivariate survival dataset of McGilchrist and Aisbett (1991 McGilchrist, C.A., Aisbett, C.W. (1991). Regression with frailty in survival analysis. Biometrics 47:461–466.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) related to the kidney infection data, and a better model is suggested for the data.