Inverse Regression in Binary Response LDV Model

Sliced inverse regression, a link-free and distribution-free method, is applied to binary response limited dependent variable models. An inverse regression property of binary response LDV model is found. Based on this property, if the distributions of X _j (j = 1, 2,…, p) satisfy the linearity condition, then β can be estimated up to a positive multiplicative scalar without any assumptions on the distribution of error ε. Moreover, the estimator can be proved to be asymptotically normal based on which testing hypotheses are considered. Simulations results are reported.     

Generalized Linear Failure Rate Distribution

The exponential and Rayleigh are the two most commonly used distributions for analyzing lifetime data. These distributions have several desirable properties and nice physical interpretations. Unfortunately, the exponential distribution only has constant failure rate and the Rayleigh distribution has increasing failure rate. The linear failure rate distribution generalizes both these distributions which may have non increasing hazard function also. This article introduces a new distribution, which generalizes linear failure rate distribution. This distribution generalizes the well-known (1) exponential distribution, (2) linear failure rate distribution, (3) generalized exponential distribution, and (4) generalized Rayleigh distribution. The properties of this distribution are discussed in this article. The maximum likelihood estimates of the unknown parameters are obtained. A real data set is analyzed and it is observed that the present distribution can provide a better fit than some other very well-known distributions.     

Estimation of Survivorship Function Based on Doubly Censored Data with Application to HIV/AIDS

Left censoring concept has been defined in different ways in statistical applications. Turnbull (1974 Turnbull , B. W. ( 1974 ). Nonparametric estimation of a survivorship function with doubly censored data . J. Amer. Statist. Assoc. 69 : 169 – 173 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) defines it in a particular way. Whereas in recent literature, especially in epidemiological studies, it has been defined in another way. This difference between the two approaches is the main reason that despite simplicity, Turnbull method cannot be applicable in all cases of doubly censored data. In this article we present a modified Turnbull method for analysis of doubly censored data adequate with recent definition. Comparison has been done with other statistical methods, including imputation estimator, full likelihood-based and conditional likelihood-based approach using Iranian HIV data.     

Inference for a Simple Step-Stress Model with Type-I Censoring and Lognormally Distributed Lifetimes

Accelerated life-testing (ALT) is a very useful technique for examining the reliability of highly reliable products. It allows the experimenter to obtain failure data more quickly at increased stress levels than under normal operating conditions. A step-stress model is one special class of ALT, and in this article we consider a simple step-stress model under the cumulative exposure model with lognormally distributed lifetimes in the presence of Type-I censoring. We then discuss inferential methods for the unknown parameters of the model by the maximum likelihood estimation method. Some numerical methods, such as the Newton–Raphson and quasi-Newton methods, are discussed for solving the corresponding non-linear likelihood equations. Next, we discuss the construction of confidence intervals for the unknown parameters based on (i) the asymptotic normality of the maximum likelihood estimators (MLEs), and (ii) parametric bootstrap resampling technique. A Monte Carlo simulation study is carried out to examine the performance of these methods of inference. Finally, a numerical example is presented in order to illustrate all the methods of inference developed here.     

Free Knot Splines with RJMCMC in Survival Data Analysis

The B-spline representation is a common tool to improve the fitting of smooth nonlinear functions, it offers a fitting as a piecewise polynomial. The regions that define the pieces are separated by a sequence of knots. The main difficulty in this type of modeling is the choice of the number and the locations of these knots. The Reversible Jump Markov Chain Monte Carlo (RJMCMC) algorithm provides a solution to simultaneously select these two parameters by considering the knots as free parameters. This algorithm belongs to the MCMC techniques that allow simulations from target distributions on spaces of varying dimension. The aim of the present investigation is to use this algorithm in the framework of the analysis of survival time, for the Cox model in particular. In fact, the relation between the hazard ratio function and the covariates being assumed to be log-linear, this assumption is too restrictive. Thus, we propose to use the RJMCMC algorithm to model the log hazard ratio function by a B-spline representation with an unknown number of knots at unknown locations. This method is illustrated with two real data sets: the Stanford heart transplant data and lung cancer survival data. Another application of the RJMCMC is selecting the significant covariates, and a simulation study is performed.     

Adaptive Lasso Variable Selection for the Accelerated Failure Models

This article considers the adaptive lasso procedure for the accelerated failure time model with multiple covariates based on weighted least squares method, which uses Kaplan-Meier weights to account for censoring. The adaptive lasso method can complete the variable selection and model estimation simultaneously. Under some mild conditions, the estimator is shown to have sparse and oracle properties. We use Bayesian Information Criterion (BIC) for tuning parameter selection, and a bootstrap variance approach for standard error. Simulation studies and two real data examples are carried out to investigate the performance of the proposed method.     

Delta Method,Moment Convergence,and Inference

Statistics, as functions of the observations, are usually given by well-behaved functions. This fact is used to obtain limit distributions for statistics whose components are given by asymptotically linear functions. These results are then extended to the moments of distributions, covariance matrices and confidence regions for parameters of interest. These regions may be used to test, through duality, hypothesis on these parameters. A theoretical application is presented.     

Kernel Density-Based Linear Regression Estimate

For linear regression models with non normally distributed errors, the least squares estimate (LSE) will lose some efficiency compared to the maximum likelihood estimate (MLE). In this article, we propose a kernel density-based regression estimate (KDRE) that is adaptive to the unknown error distribution. The key idea is to approximate the likelihood function by using a nonparametric kernel density estimate of the error density based on some initial parameter estimate. The proposed estimate is shown to be asymptotically as efficient as the oracle MLE which assumes the error density were known. In addition, we propose an EM type algorithm to maximize the estimated likelihood function and show that the KDRE can be considered as an iterated weighted least squares estimate, which provides us some insights on the adaptiveness of KDRE to the unknown error distribution. Our Monte Carlo simulation studies show that, while comparable to the traditional LSE for normal errors, the proposed estimation procedure can have substantial efficiency gain for non normal errors. Moreover, the efficiency gain can be achieved even for a small sample size.     

A Note on Wald's Type Chi-squared Tests of Fit

The Wald's method for constructing chi-squared tests of fit has been formulated more accurately. It is shown that Wald's type statistics will follow the central chi-squared distribution if and only if the limit covariance matrix of standardized frequencies will not depend on unknown parameters. Several examples that illustrate this important fact are presented. In particular, it is shown that the goodness-of-fit statistic developed by Moore and Stubblebine does not follow the chi-squared limit distribution, and, hence, cannot be used for testing multivariate normality.     

An Investigation of Agreement Between the Item Difficulty Coefficient Calculated in Accordance With Classical Test Theory and Item Response Theory With Bland–Altman Method

The purpose of this study is to investigate agreement between item difficulty coefficients calculated relying on classical test theory and item response theory with Bland–Altman method. According to results, although there is a high correlation between Pj and b coefficient estimated with HGLM (hierarchical generalized linear model), 1P, and 3P models, it can be said that there is no agreement between two methods and cannot be used interchangeably. It is observed that the confidence limit is wide according to Bland–Altman graphics. Therefore, it can be said that there is no agreement between item difficulty values obtained from two methods. Bland–Altman method which is used in clinical studies mostly is suggested to be used in the comparison of methods used especially in the evaluation of student performance in education, in agreement studies among specialist considerations especially in terms of providing additional information to the studies in which correlation coefficient is calculated.