Journal of Data Science

The Topp-Leone distribution is an attractive model for life testing and reliability studies as it acquires a bathtub shaped hazard function. In this paper, we introduce a new family of distributions, depending on Topp–Leone random variable as a generator, called the Type II generalized Topp– Leone–G (TIIGTL-G) family. Its density function can be unimodel, leftskewed, right-skewed, and reversed-J shaped, and has increasing, decreasing, upside-down, J and reversed-J hazard rates. Some special models are presented. Some of its statistical properties are studied. Explicit expressions for the ordinary and incomplete moments, quantile and generating functions, Rényi entropy and order statistics are derived. The method of maximum likelihood is used to estimate the model parameters. The importance of one special model; namely; the Type II generalized Topp–Leone exponential is illustrated through two real data sets.

Abstract: Quick identification of severe injury crashes can help Emergency Medical Services (EMS) better allocate their scarce resources to improve the survival of severely injured crash victims by providing them with a fast and timely response. Data broadcast from a vehicle’s Event Data Recorder (EDR) provide an opportunity to capture crash information and send them to EMS near real-time. A key feature of EDR data is a longitudinal measure of crash deceleration. We used functional data analysis (FDA) to ascertain key features of the deceleration trajectories (absolute integral, absolute in- tegral of its slope, and residual variance) to develop and verify a risk predic- tion model for serious (AIS 3+) injuries. We used data from the 2002-2012 EDR reports and the National Highway and National Automotive Sampling System (NASS) Crashworthiness Data System (CDS) datasets available on the National Transportation Safety Administration (NHTSA) website. We consider a variety of approaches to model deceleration data, including non- penalized and penalized splines and a variable selection method, ultimately obtaining a model with a weighted AUC of 0.93. A novel feature of our approach is the use of residual variance as a measure of predictive risk. Our model can be viewed as an important first step towards developing a real- time prediction model capable of predicting the risk of severe injury in any motor vehicle crash.

istribution of Lindley distribution constructed by combining the cumulative distribution function (cdf) of Lomax and Lindley distributions. Some mathematical properties of the new distribution are discussed including moments, quantile and moment generating function. Estimation of the model parameters is carried out using maximum likelihood method. Finally, real data examples are presented to illustrate the usefulness and applicability of this new distribution.

To date, many gene set analysis (GSA) approaches have been developed for identifying differentially expressed gene sets using microarray data. However, these methods are not directly

Inferences about the ratio of two lognormal means δ can depend on plausible values of ρ, the ratio of the normal standard deviations associated to these distributions. This aspect is not usually considered in most of the analyses carried out in some applied sciences. In this paper we propose a profile likelihood approach that allows the comparison of two independent lognormal data sets in a more exhaustive way. Inferences about δ, ρ and (δ, ρ) are jointly analyzed through a simple closed-form expression obtained for the profile likelihood function of the parameter vector (δ, ρ). A similar analysis is done for ψ and ρ, where ψ is the ratio of two lognormal medians, obtaining also a simple closed-form expression for the profile likelihood function of these parameters. These expressions allow us to construct likelihood contour plots that capture most of the information provided by the samples and become valuable to identify if a trade-off between the parameters under study occurs; in case of that, individual inferences should be analyzed carefully. A detailed series of Monte Carlo simulations are included; they illustrate the performance of profile likelihood and parametric bootstrap approaches, for different sample sizes and parameter values.

Abstract: Li and Tiwari (2008) recently developed a corrected Z-test statistic for comparing the trends in cancer age-adjusted mortality and incidence rates across overlapping geographic regions, by properly adjusting for the correlation between the slopes of the fitted simple linear regression equations. One of their key assumptions is that the error variances have unknown but common variance. However, since the age-adjusted rates are linear combinations of mortality or incidence counts, arising naturally from an underlying Poisson process, this constant variance assumption may be violated. This paper develops a weighted-least-squares based test that incorporates heteroscedastic error variances, and thus significantly extends the work of Li and Tiwari. The proposed test generally outperforms the aforementioned test through simulations and through application to the age-adjusted mortality data from the Surveillance, Epidemiology, and End Results (SEER) Program of the National Cancer Institute.

Abstract: In this paper we consider clinical trials with two treatments and a non-normally distributed response variable. In addition, we focus on ap plications which include only discrete covariates and their interactions. For such applications, the semi-parametric Area Under the ROC Curve (AUC) regression model proposed by Dodd and Pepe (2003) can be used. However, because a logistic regression procedure is used to obtain parameter estimates and a bootstrapping method is needed for computing parameter standard errors, their method may be cumbersome to implement. In this paper we propose to use a set of AUC estimates to obtain parameter estimates and combine DeLong’s method and the delta method for computing parameter standard errors. Our new method avoids heavy computation associated with the Dodd and Pepe’s method and hence is easy to implement. We conduct simulation studies to show that the two methods yield similar results. Finally, we illustrate our new method using data from urinary incontinence clinical trials.

Abstract: Of interest in this paper is the development of a model that uses inverse sampling of binary data that is subject to false-positive misclassification in an effort to estimate a proportion. From this model, both the proportion of success and false positive misclassification rate may be estimated. Also, three first-order likelihood based confidence intervals for the proportion of success are mathematically derived and studied via a Monte Carlo simulation. The simulation results indicate that the score and likelihood ratio intervals are generally preferable over the Wald interval. Lastly, the model is applied to a medical data set.

Abstract: We study the spatial distribution of clusters associated to the aftershocks of the megathrust Maule earthquake MW 8.8 of 27 February 2010. We used a recent clustering method which hinges on a nonparametric estimation of the underlying probability density function to detect subsets of points forming clusters associated with high density areas. In addition, we estimate the probability density function using a nonparametric kernel method for each of these clusters. This allows us to identify a set of regions where there is an association between frequency of events and coseismic slip. Our results suggest that high coseismic slip is spatially related to high aftershock frequency.

The power generalized Weibull distribution due to Bagdonovacius and Nikulin (2002) is an alternative,and always provides better fits than the exponentiated Weibull family for modeling lifetime data. In this paper, we consider the generalized order statistics (GOS) from this distribution. We obtain exact explicit expressions as well as recurrence relations for the single, product and conditional moments of generalized order statistics from the power generalized Weibull distribution and then we use these results to compute the means and variances of order statistics and record values for samples of different sizes for various values of the shape and scale parameters.