Standard Exponential Cure Rate Model with Informative Censoring

In this article, we consider the standard cure rate model proposed by Boag (1949 Boag , J. W. ( 1949 ). Maximum likelihood estimates of the proportion of patients cured by cancer therapy . J. Roy. Statist. Soc. B 11 : 15 – 53 . [Google Scholar]) and Berkson and Gage (1952 Berkson , J. , Gage , R. ( 1952 ). Survival curve for cancer patients following treatment . J. Amer. Statist. Assoc. 47 : 501 – 515 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). We present a new definition of informative censoring similar to Lawless (1982 Lawless , J. F. ( 1982 ). Statistical Models and Methods for Lifetime Data . New York : Wiley .[Crossref] , [Google Scholar]) and the corresponding likelihood function. Under informative censoring, we obtain the Fisher information matrix of the exponential standard cure rate model. We verify, with simulated data, the impact caused by informative censoring in the coverage probabilities and in the lengths of asymptotic confidence intervals of the parameters of interest. An example with real data is analyzed.     

A Euclidean Likelihood Estimator for Bivariate Tail Dependence

The spectral measure plays a key role in the statistical modeling of multivariate extremes. Estimation of the spectral measure is a complex issue, given the need to obey a certain moment condition. We propose a Euclidean likelihood-based estimator for the spectral measure which is simple and explicitly defined, with its expression being free of Lagrange multipliers. Our estimator is shown to have the same limit distribution as the maximum empirical likelihood estimator of Einmahl and Segers (2009 Einmahl , J. H. J. , Segers , J. ( 2009 ). Maximum empirical likelihood estimation of the spectral measure of an extreme-value distribution . Ann. Statist. 37 ( 5B ): 2953 – 2989 .[Crossref], [Web of Science ®] , [Google Scholar]). Numerical experiments suggest an overall good performance and identical behavior to the maximum empirical likelihood estimator. We illustrate the method in an extreme temperature data analysis.     

Hypothesis testing in log-Birnbaum-Saunders regressions

The log-Birnbaum-Saunders regression model introduced by Rieck and Nedelman (1991 Rieck, J. R., Nedelman, J. R. (1991). A log-linear model for the Birnbaum-Saunders distribution. Technometrics 33:51–60. [Google Scholar]) is useful for modeling lifetimes of materials and equipments subject to different conditions. Our goal in this article is twofold. First, we numerically evaluate the finite sample performances of the likelihood ratio, score and Wald tests in the log-Birnbaum-Saunders regression model. Second, we introduce a RESET-like misspecification test for that model. The null hypothesis is that the model is correctly specified which is tested against the alternative hypothesis of model misspecification. The power of the test is evaluated using Monte Carlo simulations. Bootstrap-based inference is also considered. An empirical application is presented and discussed.     

Likelihood based inference for the ratio of gamma means

Inference concerning the ratio of two means based two independent two-parameter gamma models with common shape parameter was examined in Booth et al. (1999 Booth, J. G., Hobert, J. P. and Ohman, P. A. (1999). On the probable error of the ratio of two gamma means. Biometrika, 86: 439–452.  [Google Scholar]) and a computationally intensive bootstrap calibration method was developed. In this paper, a likelihood based method is proposed for small sample inference about the ratio of two means of the two-parameter gamma models when the shape parameters may or may not be equal. The proposed method is very simple to use and, as illustrated in simulation studies, gives extremely accurate results.     

A NOTE ON ‘A COMPARISON OF BAYES AND MAXIMUM LIKELIHOOD ESTIMATION OF THE INTRACLASS CORRELATION COEFFICIENT’

Palmer and Broemeling [1] Palmer, J. L. and Broemeling, L. D. 1990. A Comparison of Bayes and Maximum Likelihood Estimation of the Intraclass Correlation Coefficient. Comm. Statist.-Theory Meth, 19: 953–975. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar] compare Bayes and maximum likelihood estimates of the intraclass correlation (ICC). The prior information in their derivation of the Bayes estimator is placed on the variance components instead of the ICC itself. This paper finds a Bayes estimator of the ICC with the prior placed on the ICC. Bayes estimates based on three different priors are then compared to method of moments estimate.     

Bayesian estimation of sensitivity level and population proportion of a sensitive characteristic in a binary optional unrelated question RRT model

Sihm et al. (2016 Sihm, J. S., A. Chhabra, and S. N. Gupta. 2016. An optional unrelated question RRT model. Involve: A Journal of Mathematics 9 (2):195–209.[Crossref] , [Google Scholar]) proposed an unrelated question binary optional randomized response technique (RRT) model for estimating the proportion of population that possess a sensitive characteristic and the sensitivity level of the question. In our work, decision theoretic approach has been followed to obtain Bayes estimates of the two parameters along with their corresponding minimal Bayes posterior expected losses (BPEL) using beta prior and squared error loss function (SELF). Relative losses are also examined to compare the performances of the Bayes estimates with those of the classical estimates obtained by Sihm et al. (2016 Sihm, J. S., A. Chhabra, and S. N. Gupta. 2016. An optional unrelated question RRT model. Involve: A Journal of Mathematics 9 (2):195–209.[Crossref] , [Google Scholar]). The results obtained are illustrated with the help of real survey data using non informative prior.     

Maximum Likelihood Estimation in Compound Geometric Failure Model with Changing Parameters From Type-I Two-Stage Progressively Censored and Group Censored Samples

Boardman and Kendell (1970 Boardman , T. J. , Kendell , P. J. ( 1970 ). Estimation in compound failure models . Technometrics 12 : 891 – 908 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) considered the problem of estimation with respect to Type-I censoring when an item is subjected to only one of the two causes of failure assuming exponential model. Patel and Gajjar (1992 Patel , M. N. , Gajjar , A. V. ( 1992 ). Maximum likelihood estimation in compound exponential failure model with changing failure rates from Type-I progressively censored and group censored samples . Commun. Statist. Theor. Meth. 21 ( 10 ): 2899 – 2908 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) considered extension of the Boardman and Kendell's results in case of two-stage progressive censoring. Here we have considered geometric competing risk failure model with two independent causes of failures. Maximum likelihood estimation of the parameters is carried out using Type-I two-stage progressively censored and group censored samples. Asymptotic standard errors of the estimators are obtained for both the cases. Two illustrative examples are cited for ungroup and group competing risk models.     

Fitting spatial regressions to large datasets using unilateral approximations

Maximum likelihood estimation of a spatial model typically requires a sizeable computational capacity, even in relatively small samples, and becomes unfeasible in very large datasets. The unilateral approximation approach to spatial model estimation (suggested in Besag 1974 Besag, J. E. 1974. Spatial interaction and the statistical analysis of lattice systems. Journal of the Royal Statistical Society. Series B (Methodological) 36 (2):192–236.[Crossref], [Web of Science ®] , [Google Scholar]) provides a viable alternative to maximum likelihood estimation that reduces substantially the computing time and the storage required. In this article, we extend the method, originally proposed for conditionally specified processes, to simultaneous and to general bilateral spatial processes over rectangular lattices. We prove the estimators’ consistency and study their finite-sample properties via Monte Carlo simulations.     

Comments on the econometrics of female labor supply and children by alice and masao nokalnura

Abstract

The present paper focuses attention on the sensitivity of technical inefficiency to most commonly used one‐sided distributions of the inefficiency error term, namely the truncated normal, the half‐normal, and the exponential distributions. A generalized version of the half‐normal, which does not embody the zero‐mean restriction, is also explored. For each distribution, the likelihood function and the counterpart of the estimator of technical efficiency are explicitly stated (Jondrow, J., Lovell, C. A. K., Materov, I. S., Schmidt, P. ([1982] Jondrow, J., Lovell, C. A. K., Materov, I. S. and Schmidt, P. 1982. On estimation of technical inefficiency in the stochastic frontier production function model. J. Econometrics, 19: 233–238. [Crossref], [Web of Science ®] , [Google Scholar]), On estimation of technical inefficiency in the stochastic frontier production function model, J. Econometrics19:233–238). Based on our panel data set, related to Tunisian manufacturing firms over the period 1983–1993, formal tests lead to a strong rejection of the zero‐mean restriction embodied in the half normal distribution. Our main conclusion is that the degree of measured inefficiency is very sensitive to the postulated assumptions about the distribution of the one‐sided error term. The estimated inefficiency indices are, however, unaffected by the choice of the functional form for the production function.     

Unified multivariate survival model with a surviving fraction: an application to a Brazilian customer churn data

In this paper we propose a new lifetime model for multivariate survival data in presence of surviving fractions and examine some of its properties. Its genesis is based on situations in which there are m types of unobservable competing causes, where each cause is related to a time of occurrence of an event of interest. Our model is a multivariate extension of the univariate survival cure rate model proposed by Rodrigues et al. [37 J. Rodrigues, V.G. Cancho, M. de Castro, and F. Louzada-Neto, On the unification of long-term survival models, Statist. Probab. Lett. 79 (2009), pp. 753–759. doi: 10.1016/j.spl.2008.10.029[Crossref], [Web of Science ®] , [Google Scholar]]. The inferential approach exploits the maximum likelihood tools. We perform a simulation study in order to verify the asymptotic properties of the maximum likelihood estimators. The simulation study also focus on size and power of the likelihood ratio test. The methodology is illustrated on a real data set on customer churn data.     

Note on posterior inference for the Bingham distribution

The properties of high-dimensional Bingham distributions have been studied by Kume and Walker (2014 Kume, A., and S. G. Walker. 2014. On the Bingham distribution with large dimension. Journal of Multivariate Analysis 124:345–52.[Crossref], [Web of Science ®] , [Google Scholar]). Fallaize and Kypraios (2016 Fallaize, C. J., and T. Kypraios. 2016. Exact Bayesian inference for the Bingham distribution. Statistics and Computing 26:349–60.[Crossref], [Web of Science ®] , [Google Scholar]) propose the Bayesian inference for the Bingham distribution and they use developments in Bayesian computation for distributions with doubly intractable normalizing constants (Møller et al. 2006 Møller, J., A. N. Pettitt, R. Reeves, and K. K. Berthelsen. 2006. An efficient Markov chain Monte Carlo method for distributions with intractable normalising constants. Biometrika 93 (2):451–458.[Crossref], [Web of Science ®] , [Google Scholar]; Murray, Ghahramani, and MacKay 2006 Murray, I., Z. Ghahramani, and D. J. C. MacKay. 2006. MCMC for doubly intractable distributions. In Proceedings of the 22nd annual conference on uncertainty in artificial intelligence (UAI-06), 359–66. AUAI Press. [Google Scholar]). However, they rely heavily on two Metropolis updates that they need to tune. In this article, we propose instead a model selection with the marginal likelihood.     

The asymptotic covariance matrix of the QMLE in ARMA models

A compact analytical representation of the asymptotic covariance matrix, in terms of model parameters directly, of the quasi maximum likelihood estimator (QMLE) is derived in autoregressive moving average (ARMA) models with possible nonzero means and non-Gaussian error terms. For model parameters excluding the error variance, it is found that the Huber (1967 Huber, P. J. (1967). The behavior of maximum likelihood estimates under nonstandard conditions. Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics and Probability, Volume 1, pp. 221–233. [Google Scholar]) sandwich form for the asymptotic covariance matrix degenerates into the inverse of the associated information matrix. In comparison to the existing result that involves the second moments of some auxiliary variables for the case of zero-mean ARMA models, the analytical asymptotic covariance in this article has an advantage in that it can be conveniently estimated by plugging in the estimated model parameters directly.     

Steady-state Gibbs sampler estimation for lung cancer data

This paper is based on the application of a Bayesian model to a clinical trial study to determine a more effective treatment to lower mortality rates and consequently to increase survival times among patients with lung cancer. In this study, Qian et al. [13 J. Qian, D.K. Stangl, and S. George, A Weibull model for survival data: Using prediction to decide when to stop a clinical trial, in Bayesian Biostatistics, D. Berry and D. Stangl, eds., Marcel Dekker, New York, 1996, pp. 187–205. [Google Scholar]] strived to determine if a Weibull survival model can be used to decide whether to stop a clinical trial. The traditional Gibbs sampler was used to estimate the model parameters. This paper proposes to use the independent steady-state Gibbs sampling (ISSGS) approach, introduced by Dunbar et al. [3 M. Dunbar, H.M. Samawi, R. Vogel, and L. Yu, A more efficient Gibbs sampler estimation using steady state simulation: Application to public health studies, J. Stat. Simul. Comput. 10.1080/00949655.2013.770857.[Taylor &; Francis Online] , [Google Scholar]], to improve the original Gibbs sampler in multidimensional problems. It is demonstrated that ISSGS provides accuracy with unbiased estimation and improves the performance and convergence of the Gibbs sampler in this application.     

Frailty models power variance function with cure fraction and latent risk factors negative binomial

In this article, we propose a flexible cure rate model, which is an extension of Cancho et al. (2011 Cancho, V.G., Rodrigues, J., de Castro, M. (2011). A flexible model for survival data with a cure rate: A Bayesian approach. J. Appl. Stat. 38:57–70.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) model, by incorporating a power variance function (PVF) frailty term in latent risk. The model is more flexible in terms of dispersion and it also quantifies the unobservable heterogeneity. The parameter estimation is reached by maximum likelihood estimation procedure and Monte Carlo simulation studies are considered to evaluate the proposed model performance. The practical relevance of the model is illustrated in a real data set of preventing cancer recurrence.     

Maximum simulated likelihood estimation of the panel sample selection model

Heckman's (1976 Heckman, J. J. (1976). The common structure of statistical models of truncation, sample selection and limited dependent variables and a simple estimator for such models. Annals of Economic and Social Measurement 15:475–492. [Google Scholar], 1979 Heckman, J. J. (1979). Sample selection bias as a specification error. Econometrica 47(1):153–161.[Crossref], [Web of Science ®] , [Google Scholar]) sample selection model has been employed in many studies of linear and nonlinear regression applications. It is well known that ignoring the sample selectivity may result in inconsistency of the estimator due to the correlation between the statistical errors in the selection and main equations. In this article, we reconsider the maximum likelihood estimator for the panel sample selection model in Keane et al. (1988 Keane, M., Moffitt, R., Runkle, D. (1988). Real wages over the business cycle: Estimating the impact of heterogeneity with micro data. Journal of Political Economy 96:1232–1266.[Crossref], [Web of Science ®] , [Google Scholar]). Since the panel data model contains individual effects, such as fixed or random effects, the likelihood function is more complicated than that of the classical Heckman model. As an alternative to the existing derivation of the likelihood function in the literature, we show that the conditional distribution of the main equation follows a closed skew-normal (CSN) distribution, of which the linear transformation is still a CSN. Although the evaluation of the likelihood function involves high-dimensional integration, we show that the integration can be further simplified into a one-dimensional problem and can be evaluated by the simulated likelihood method. Moreover, we also conduct a Monte Carlo experiment to investigate the finite sample performance of the proposed estimator and find that our estimator provides reliable and quite satisfactory results.     

Comparison of Frailty Models for Acute Leukemia Data under Gompertz Baseline Distribution

In this article, we consider two different shared frailty regression models under the assumption of Gompertz as baseline distribution. Mostly assumption of gamma distribution is considered for frailty distribution. To compare the results with gamma frailty model, we consider the inverse Gaussian shared frailty model also. We compare these two models to a real life bivariate survival data set of acute leukemia remission times (Freireich et al., 1963 Freireich, E.J., Gehan, E., Frei, E., Schroeder, L.R., Wolman, I.J., Anbari, R., Burgert, E.O., Mills, S.D., Pinkel, D., Selawry, O.S., Moon, J.H., Gendel, B.R., Spurr, C.L., Storrs, R., Haurani, F., Hoogstraten, B., Lee, S. (1963). The effect of 6-mercaptopurine on the duration of steroid-induced remissions in acute leukemia: a model for evaluation of other potentially useful therapy. Blood 21:699–716.[Web of Science ®] , [Google Scholar]). Analysis is performed using Markov Chain Monte Carlo methods. Model comparison is made using Bayesian model selection criterion and a well-fitted model is suggested for the acute leukemia data.     

Model selection in regression based on pre-smoothing

In this paper, we investigate the effect of pre-smoothing on model selection. Christóbal et al 6 Christóbal Christóbal, J. A., Faraldo Roca, P. and González Manteiga, W. 1987. A class of linear regression parameter estimators constructed by nonparametric estimation. Ann. Statist.,, 15: 603–609. [Crossref], [Web of Science ®] [Google Scholar] showed the beneficial effect of pre-smoothing on estimating the parameters in a linear regression model. Here, in a regression setting, we show that smoothing the response data prior to model selection by Akaike's information criterion can lead to an improved selection procedure. The bootstrap is used to control the magnitude of the random error structure in the smoothed data. The effect of pre-smoothing on model selection is shown in simulations. The method is illustrated in a variety of settings, including the selection of the best fractional polynomial in a generalized linear model.     

Empirical likelihood ratio with doubly truncated data

Doubly truncated data appear in a number of applications, including astronomy and survival analysis. For doubly-truncated data, the lifetime T is observable only when U≤T≤V, where U and V are the left-truncated and right-truncated time, respectively. Based on the empirical likelihood approach of Zhou [21 Zhou, M. 2005. Empirical likelihood ratio with arbitrarily censored/truncated data by EM algorithm. J. Comput. Graph. Statist., 14: 643–656. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]], we propose a modified EM algorithm of Turnbull [19 Turnbull, B. W. 1976. The empirical distribution function with arbitrarily grouped censored and truncated data. J. R. Stat. Soc. Ser. B, 38: 290–295.  [Google Scholar]] to construct the interval estimator of the distribution function of T. Simulation results indicate that the empirical likelihood method can be more efficient than the bootstrap method.     

Some tests for the allometric extension model in regression

The allometric extension model is a multivariate regression model recently proposed by Tarpey and Ivey (2006 Tarpey, T., Ivey, C.T. (2006). Allometric extension for multivariate regression. J. Data Sci. 4:479–495. [Google Scholar]). This model holds when the matrix of covariances between the variables in the response vector y and the variables in the vector of regressors x has a particular structure. In this paper, we consider tests of hypotheses for this structure when (y′, x′)′ has a multivariate normal distribution. In particular, we investigate the likelihood ratio test and a Wald test.     

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]).