Search

In this paper, a new multivariate zero-inflated binomial (MZIB) distribution is proposed to analyse the correlated proportional data with excessive zeros. The distributional properties of purposed model are studied. The Fisher scoring algorithm and EM algorithm are given for the computation of estimates of parameters in the proposed MZIB model with/without covariates. The score tests and the likelihood ratio tests are derived for assessing both the zero-inflation and the equality of multiple binomial probabilities in correlated proportional data. A limited simulation study is performed to evaluate the performance of derived EM algorithms for the estimation of parameters in the model with/without covariates and to compare the nominal levels and powers of both score tests and likelihood ratio tests. The whitefly data is used to illustrate the proposed methodologies. [ABSTRACT FROM AUTHOR]

The widely used Fellegi–Sunter model for probabilistic record linkage does not leverage information contained in field values and consequently leads to identical classification of match status regardless of whether records agree on rare or common values. Since agreement on rare values is less likely to occur by chance than agreement on common values, records agreeing on rare values are more likely to be matches. Existing frequency-based methods typically rely on knowledge of error probabilities associated with field values and frequencies of agreed field values among matches, often derived using prior studies or training data. When such information is unavailable, applications of these methods are challenging. In this paper, we propose a simple two-step procedure for frequency-based matching using the Fellegi–Sunter framework to overcome these challenges. Matching weights are adjusted based on frequency distributions of the agreed field values among matches and non-matches, estimated by the Fellegi–Sunter model without relying on prior studies or training data. Through a real-world application and simulation, our method is found to produce comparable or better performance than the unadjusted method. Furthermore, frequency-based matching provides greater improvement in matching accuracy when using poorly discriminating fields with diminished benefit as the discriminating power of matching fields increases. [ABSTRACT FROM AUTHOR]

In this paper, we introduce a new first-order generalized Poisson integer-valued autoregressive process, for modeling integer-valued time series exhibiting a piecewise structure and overdispersion. Basic probabilistic and statistical properties of this model are discussed. Conditional least squares and conditional maximum likelihood estimators are derived. The asymptotic properties of the estimators are established. Moreover, two special cases of the process are discussed. Finally, some numerical results of the estimates and a real data example are presented. [ABSTRACT FROM AUTHOR]

In this paper, we propose a new kind of multivariate t distribution by allowing different degrees of freedom for each univariate component. Compared with the classical multivariate t distribution, it is more flexible in the model specification that can be used to deal with the variant amounts of tail weights on marginals in multivariate data modeling. In particular, it could include components following the multivariate normal distribution, and it contains the product of independent t-distributions as a special case. Subsequently, it is extended to the regression model as the joint distribution of the error terms. Important distributional properties are explored and useful statistical methods are developed. The flexibility of the specified structure in better capturing the characteristic of data is exemplified by both simulation studies and real data analyses. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the mean change-point models based on associated sequences. Under some weak conditions, we obtain a limit distribution of CUSUM statistic which can be used to judge the mean change-mount δ n is satisfied or dissatisfied n 1 / 2 δ n = o (1). We also study the consistency of sample covariances and change-point location statistics. Based on Normality and Lognormality data, some simulations such as empirical sizes, empirical powers and convergence are presented to test our results. As an important application, we use CUSUM statistics to do the mean change-point analysis for a financial series. [ABSTRACT FROM AUTHOR]

The models used to describe the kinetics of ruminal degradation are usually nonlinear models where the dependent variable is the proportion of degraded food. The method of least squares is the standard approach used to estimate the unknown parameters but this method can lead to unacceptable predictions. To solve this issue, a beta nonlinear model and the Bayesian perspective is proposed in this article. The application of standard methodologies to obtain prior distributions, such as the Jeffreys prior or the reference priors, involves serious difficulties here because this model is a nonlinear non-normal regression model, and the constrained parameters appear in the log-likelihood function through the Gamma function. This paper proposes an objective method to obtain the prior distribution, which can be applied to other models with similar complexity, can be easily implemented in OpenBUGS, and solves the problem of unacceptable predictions. The model is generalized to a larger class of models. The methodology was applied to real data with three models that were compared using the Deviance Information Criterion and the root mean square prediction error. A simulation study was performed to evaluate the coverage of the credible intervals. [ABSTRACT FROM AUTHOR]

Geyser eruption is one of the most popular signature attractions at the Yellowstone National Park. The interdependence of geyser eruptions and impacts of covariates are of interest to researchers in geyser studies. In this paper, we propose a parametric covariate-adjusted recurrent event model for estimating the eruption gap time. We describe a general bivariate recurrent event process, where a bivariate lognormal distribution and a Gumbel copula with different marginal distributions are used to model an interdependent dual-type event system. The maximum likelihood approach is used to estimate model parameters. The proposed method is applied to analyzing the Yellowstone geyser eruption data for a bivariate geyser system and offers a deeper understanding of the event occurrence mechanism of individual events as well as the system as a whole. A comprehensive simulation study is conducted to evaluate the performance of the proposed method. [ABSTRACT FROM AUTHOR]

In this paper, we estimate the dynamic parameters of a time-varying coefficient model through radial kernel functions in the context of a longitudinal study. Our proposal is based on a linear combination of weighted kernel functions involving a bandwidth, centered around a given set of time points. In addition, we study different alternatives of estimation and inference including a Frequentist approach using weighted least squares along with bootstrap methods, and a Bayesian approach through both Markov chain Monte Carlo and variational methods. We compare the estimation strategies mention above with each other, and our radial kernel functions proposal with an expansion based on regression spline, by means of an extensive simulation study considering multiples scenarios in terms of sample size, number of repeated measurements, and subject-specific correlation. Our experiments show that the capabilities of our proposal based on radial kernel functions are indeed comparable with or even better than those obtained from regression splines. We illustrate our methodology by analyzing data from two AIDS clinical studies. [ABSTRACT FROM AUTHOR]

For banks using the Advanced Internal Ratings-Based Approach in accordance with Basel III requirements, the amount of required regulatory capital relies on the banks' estimates of the probability of default, the loss given default and the conversion factor for their credit risk portfolio. Therefore, for both model development and validation, assessing the models' predictive and discriminatory abilities is of key importance in order to ensure an adequate quantification of risk. This paper compares different measures of discriminatory power suitable for multi-class target variables such as in loss given default (LGD) models, which are currently used among banks and supervisory authorities. This analysis highlights the disadvantages of using measures that solely rely on pairwise comparisons when applied in a multi-class setting. Thus, for multi-class classification problems, we suggest using a generalisation of the well-known area under the receiver operating characteristic (ROC) curve known as the volume under the ROC surface (VUS). Furthermore, we present the R-package VUROCS, which allows for a time-efficient computation of the VUS as well as associated (co)variance estimates and illustrate its usage based on real-world loss data and validation principles. [ABSTRACT FROM AUTHOR]

The length of the gap is the key factor affecting its reliability. Based on the mechanism of the gap null gate, this paper regards the two endpoint thresholds of the gap length as bivariate random variables and establishes successful response models. Score test statistic is presented to test the correlation coefficient. The DIC criterion is also provided to compare the models. With the experimental data of the gap null gate, we build Probit model and Logit model as the successful response models, and prove that the correlation coefficients in the both models can be regarded as 0. By comparing the DIC value, we find that the Probit model is more suitable to describe the distribution of the endpoint thresholds of the reliability window. Finally, both the point estimation and interval estimation results of the reliability window are given to illustrate the feasibility of the method shown in the paper. [ABSTRACT FROM AUTHOR]

In clinical researches, various clinical scales with ordered categories are used to evaluate the efficacy and safety/toxicity of treatments. Such the clinical scales are sometimes summarized on the transition between the baseline and the study end point as a square contingency table. Also, clinical scales may be reclassified into three groups. However, the cutpoints can be varied depending on clinical researches or clinicians. Hence, this paper proposes a measure which is expressed by using same weights for collapsed tables and which can see the directionality for two kinds of asymmetries. Also, this paper shows an application of the proposed measure to clinical data, and that the proposed measure is a useful statistical method for analyzing ordered categorical data. [ABSTRACT FROM AUTHOR]

This paper is motivated by the analytical challenges we encounter when analyzing the ChAMP (Chondral Lesions And Meniscus Procedures) study, a randomized controlled trial to compare debridement to observation of chondral lesions in arthroscopic knee surgery. The main outcome, WOMAC (Western Ontario and McMaster Universities Osteoarthritis Index) pain score, is derived from the patient's responses to the questionnaire collected in the study. The major goal is to identify potentially important variables that contribute to this outcome. In this paper, the model of interest is a semiparametric model for the pain score. To address the missing data issue, we adopt a flexible missingness mechanism that is much more versatile in practice than a single parametric model. Then we propose a pairwise conditional likelihood approach to estimate the unknown parameter in the semiparametric model without the need of modeling its nonparametric counterpart nor the missingness mechanism. For variable selection, we apply a regularization approach with a variety of stability enhanced tuning parameter selection methods. We conduct comprehensive simulation studies to evaluate the performance of the proposed method. We also apply the proposed method to the ChAMP study to demonstrate its usefulness. [ABSTRACT FROM AUTHOR]

The paper explores the relationships between subjective and contextual covariates as determinants of discrimination perception in a large sample survey. A modelling approach is implemented to detect the perception of discrimination on the ground of sexual orientation through expressed ratings, in a cross-country perspective. Exploiting the capabilities of a mixture model, which interprets ordinal evaluations as the combined results of two latent components, the cognitive process of selection among discrete ordered alternatives is discussed in terms of attractiveness towards the item and uncertainty in the response pattern. Gender, age and political orientation turn out to be significant variables at the individual level. Furthermore, the ILGA country score and education level, as contextual variables, exert prominent effects at country level. Empirical evidence from the Special Eurobarometer 2015 on Opinions as well as a simulation experiment support the usefulness of the approach. [ABSTRACT FROM AUTHOR]

In this paper, we consider a regression model under the generalized Waring distribution for modeling count data. We develop and implement local influence diagnostic techniques based on likelihood displacement. Also we develop case-deletion methods. The generalized Waring regression model is presented as a mixture of the Negative Binomial and the Beta II distributions, and it is compared to the Negative Binomial and Waring regression models. Estimation is performed by maximum likelihood function. The influence measures developed in this paper are applied to a Spanish football league data set. Empirical results show that the generalized Waring regression model performs better when compared to the Negative Binomial and Waring regression models. Technical details are presented in the Appendix. [ABSTRACT FROM AUTHOR]

Investigating the main determinants of the mechanical performance of metals is not a simple task. Already known physically inspired qualitative relations between 2D microstructure characteristics and 3D mechanical properties can act as the starting point of the investigation. Isotonic regression allows to take into account ordering relations and leads to more efficient and accurate results when the underlying assumptions actually hold. The main goal in this paper is to test order relations in a model inspired by a materials science application. The statistical estimation procedure is described considering three different scenarios according to the knowledge of the variances: known variance ratio, completely unknown variances, and variances under order restrictions. New likelihood ratio tests are developed in the last two cases. Both parametric and non-parametric bootstrap approaches are developed for finding the distribution of the test statistics under the null hypothesis. Finally an application on the relation between geometrically necessary dislocations and number of observed microstructure precipitations is shown. [ABSTRACT FROM AUTHOR]

Health care audits are crucial in managing the government insurance programs that are estimated to have losses amounting to billions of dollars every year. Statistical methods such as sampling have long been used to handle their size and complexity. Sampling from health care claims data can benefit from multi-stage approaches, especially when the evaluation of the tradeoffs between precision and cost is important. The use of decision models could facilitate health care auditors and policy makers make the best use of these sampling outputs. This paper proposes an integrated multi-stage sampling and decision-making framework that enables auditors address the tradeoffs between audit costs and expected overpayment recovery. We illustrate the framework and discuss insights utilizing a variety of overpayment scenarios for payment populations including U.S. Medicare Part B claims payment data. [ABSTRACT FROM AUTHOR]

In this paper, we introduced a new probability distribution named as inverse power Maxwell distribution. The proposal distribution can be seen as an extension of the Maxwell distribution with more flexibility in modeling upside-down lifetime data. Some statistical properties of this distribution are derived. In estimation viewpoint, five methods are used for estimating the unknown parameters of the distribution and these methods are performed through the simulation study. Finally, two real data sets were analyzed to illustrate the applicability of the proposed distribution, proving that it fits each real data set much better than some other existing distributions. [ABSTRACT FROM AUTHOR]

Linear regression models are useful statistical tools to analyze data sets in different fields. There are several methods to estimate the parameters of a linear regression model. These methods usually perform under normally distributed and uncorrelated errors. If error terms are correlated the Conditional Maximum Likelihood (CML) estimation method under normality assumption is often used to estimate the parameters of interest. The CML estimation method is required a distributional assumption on error terms. However, in practice, such distributional assumptions on error terms may not be plausible. In this paper, we propose to estimate the parameters of a linear regression model with autoregressive error term using Empirical Likelihood (EL) method, which is a distribution free estimation method. A small simulation study is provided to evaluate the performance of the proposed estimation method over the CML method. The results of the simulation study show that the proposed estimators based on EL method are remarkably better than the estimators obtained from CML method in terms of mean squared errors (MSE) and bias in almost all the simulation configurations. These findings are also confirmed by the results of the numerical and real data examples. [ABSTRACT FROM AUTHOR]

In obstetrics and gynecology, knowledge about how women's features are associated with childbirth is important. This leads to establishing guidelines and can help managers to describe the dynamics of pregnant women's hospital stays. Then, time is a variable of great importance and can be described by survival models. An issue that should be considered in the modeling is the inclusion of women for whom the duration of labor cannot be observed due to fetal death, generating a proportion of times equal to zero. Additionally, another proportion of women's time may be censored due to some intervention. The aim of this paper was to present the Log-Normal zero-inflated cure regression model and to evaluate likelihood-based parameter estimation by a simulation study. In general, the inference procedures showed a better performance for larger samples and low proportions of zero inflation and cure. To exemplify how this model can be an important tool for investigating the course of the childbirth process, we considered the Better Outcomes in Labor Difficulty project dataset and showed that parity and educational level are associated with the main outcomes. We acknowledge the World Health Organization for granting us permission to use the dataset. [ABSTRACT FROM AUTHOR]

In this paper, we develop a variable selection framework with the spike-and-slab prior distribution via the hazard function of the Cox model. Specifically, we consider the transformation of the score and information functions for the partial likelihood function evaluated at the given data from the parameter space into the space generated by the logarithm of the hazard ratio. Thereby, we reduce the nonlinear complexity of the estimation equation for the Cox model and allow the utilization of a wider variety of stable variable selection methods. Then, we use a stochastic variable search Gibbs sampling approach via the spike-and-slab prior distribution to obtain the sparsity structure of the covariates associated with the survival outcome. Additionally, we conduct numerical simulations to evaluate the finite-sample performance of our proposed method. Finally, we apply this novel framework on lung adenocarcinoma data to find important genes associated with decreased survival in subjects with the disease. [ABSTRACT FROM AUTHOR]

Health data are often not symmetric to be adequately modeled through the usual normal distributions; most of them exhibit skewed patterns. They can indeed be modeled better through the larger family of skew-normal distributions covering both skewed and symmetric cases. Since outliers are not uncommon in complex real-life experimental datasets, a robust methodology automatically taking care of the noises in the data would be of great practical value to produce stable and more precise research insights leading to better policy formulation. In this paper, we develop a class of robust estimators and testing procedures for the family of skew-normal distributions using the minimum density power divergence approach with application to health data. In particular, a robust procedure for testing of symmetry is discussed in the presence of outliers. Two efficient computational algorithms are discussed. Besides deriving the asymptotic and robustness theory for the proposed methods, their advantages and utilities are illustrated through simulations and a couple of real-life applications for health data of athletes from Australian Institute of Sports and AIDS clinical trial data. [ABSTRACT FROM AUTHOR]

This paper proposes a differing methodology from the Brazilian Electricity Regulatory Agency on the efficiency estimation for the Brazilian electricity distribution sector. Our proposal combines robust state-space models and stochastic frontier analysis to measure the operational cost efficiency in a panel data set from 60 Brazilian electricity distribution utilities. The modeling joins the main literature in energy economics with advanced econometric and statistic techniques in order to estimate the efficiencies. Moreover, the suggested model is able to deal with changes in the inefficiencies across time whilst the Bayesian paradigm – through Markov chain Monte Carlo techniques – facilitates the inference on all unknowns. The method enables a significant degree of flexibility in the resultant efficiencies and a complete photography about the distribution sector. [ABSTRACT FROM AUTHOR]

In this paper, inference for a competing risks model is studied when latent failure times follow Kumaraswamy distribution and causes of failure are partially observed. Under generalized progressive hybrid censoring, existence and uniqueness of maximum likelihood estimators of model parameters are established. The confidence intervals are obtained by using asymptotic distribution theory. We further compute Bayes estimators along with credible intervals. In addition, inference is also discussed when there is order restricted shape parameters. The performance of all estimates is investigated using Monte-Carlo simulations. Finally, analysis of a real data set is presented for illustration purposes. [ABSTRACT FROM AUTHOR]

Sparse functional data are commonly observed in real-data analyzes. For such data, we propose a new classification method based on functional principal component analysis (FPCA) and bootstrap aggregating. Bootstrap aggregating is believed to improve the single classifier. In this paper, we apply this belief to an FPCA based classification, and compare the classification performance with that of the single classifiers. The simulation results show that the proposed method performs better than the conventional single classifiers. We then conduct two real-data analyzes. [ABSTRACT FROM AUTHOR]

Investigations of multivariate population are pretty common in applied researches, and the two-way crossed factorial design is a common design used at the exploratory phase in industrial applications. When assumptions such as multivariate normality and covariance homogeneity are violated, the conventional wisdom is to resort to nonparametric tests for hypotheses testing. In this paper we compare the performances, and in particular the power, of some nonparametric and semi-parametric methods that have been developed in recent years. Specifically, we examined resampling methods and robust versions of classical multivariate analysis of variance (MANOVA) tests. In a simulation study, we generate data sets with different configurations of factor's effect, number of replicates, number of response variables under null hypothesis, and number of response variables under alternative hypothesis. The objective is to elicit practical advice and guides to practitioners regarding the sensitivity of the tests in the various configurations, the tradeoff between power and type I error, the strategic impact of increasing number of response variables, and the favourable performance of one test when the alternative is sparse. A real case study from an industrial engineering experiment in thermoformed packaging production is used to compare and illustrate the application of the various methods. [ABSTRACT FROM AUTHOR]

In some applications, the clustered survival data are arranged spatially such as clinical centers or geographical regions. Incorporating spatial variation in these data not only can improve the accuracy and efficiency of the parameter estimation, but it also investigates the spatial patterns of survivorship for identifying high-risk areas. Competing risks in survival data concern a situation where there is more than one cause of failure, but only the occurrence of the first one is observable. In this paper, we considered Bayesian subdistribution hazard regression models with spatial random effects for the clustered HIV/AIDS data. An intrinsic conditional autoregressive (ICAR) distribution was employed to model the areal spatial random effects. Comparison among competing models was performed by the deviance information criterion. We illustrated the gains of our model through application to the HIV/AIDS data and the simulation studies. [ABSTRACT FROM AUTHOR]

In many real-world applications of monitoring multivariate spatio-temporal data that are non-stationary over time, one is often interested in detecting hot-spots with spatial sparsity and temporal consistency, instead of detecting system-wise changes as in traditional statistical process control (SPC) literature. In this paper, we propose an efficient method to detect hot-spots through tensor decomposition, and our method has three steps. First, we fit the observed data into a Smooth Sparse Decomposition Tensor (SSD-Tensor) model that serves as a dimension reduction and de-noising technique: it is an additive model decomposing the original data into: smooth but non-stationary global mean, sparse local anomalies, and random noises. Next, we estimate model parameters by the penalized framework that includes Least Absolute Shrinkage and Selection Operator (LASSO) and fused LASSO penalty. An efficient recursive optimization algorithm is developed based on Fast Iterative Shrinkage Thresholding Algorithm (FISTA). Finally, we apply a Cumulative Sum (CUSUM) Control Chart to monitor model residuals after removing global means, which helps to detect when and where hot-spots occur. To demonstrate the usefulness of our proposed SSD-Tensor method, we compare it with several other methods including scan statistics, LASSO-based, PCA-based, T2-based control chart in extensive numerical simulation studies and a real crime rate dataset. [ABSTRACT FROM AUTHOR]

Bayesian bandwidth selections in multivariate associated kernel estimation of probability density functions are known to improve classical methods such as cross-validation techniques in terms of execution time and smoothing quality. The paper focuses on a basic multivariate gamma kernel which is appropriated to estimate densities with support [ 0 , ∞) d . For this purpose, we consider a Bayesian adaptive estimation of the bandwidths vector under the usual quadratic loss function. The exact expression of the posterior distribution and the vector of bandwidths are obtained. Simulations studies highlight the excellent performance of the proposed approach, comparing to the global cross-validation bandwidth selection, and under integrated squared errors. Two bivariate and trivariate applications to the Old Faithful geyser data and new ones on drinking water pumps in the Sahel, respectively, are made. [ABSTRACT FROM AUTHOR]

The recent controversy about the size of crowds at candlelight protests in Korea raises an interesting question regarding the methods used to estimate crowd size. Protest organizers tend to count all participants in the event from its start to finish, while the police usually report the crowd size at its peak. While several counting methods are available to estimate the size of a crowd at a given time, counting the total number of the participants at a protest is not straightforward. In this paper, we propose a new estimator to count the total number of participants that we call the size of a dynamic crowd. We assume that the arrival and departure times of the crowd are randomly observed and that the number of the attendees in the crowd at a specific time is estimable. We estimate the number of total attendees during the entire gathering based on the capture-recapture model. We also propose a bootstrap procedure to construct a confidence interval for the crowd size. We demonstrate the performance of the proposed method with simulation studies and the data from Korea's March for Science, a global event across the world on Earth Day, April 22, 2017. [ABSTRACT FROM AUTHOR]

Partially paired data, either with incompleteness in one or both arms, are common in practice. For testing equality of means of two arms, practitioners often use only the portion of data with complete pairs and perform paired tests. Although such tests (referred as 'naive paired tests') are legitimate, their powers might be low as only partial data are utilized. The recently proposed 'P-value pooling methods', based on combining P-values from two tests, use all data, have reasonable type-I error control and good power property. While it is generally believed that 'P-value pooling methods' are superior to 'naive paired tests' in terms of power as the former use more data than the latter, no detailed power comparison has been done. This paper aims to compare powers of 'naive paired tests' and 'P-value pooling methods' analytically and our findings are counterintuitive, i.e. the 'P-value pooling methods' do not always outperform the naive paired tests in terms of power. Based on these results, we present guidance on how to select the best test for testing equality of means with partially paired data. [ABSTRACT FROM AUTHOR]

Rotatable designs that are available for process/ product optimization trials are mostly symmetric in nature. In many practical situations, response surface designs (RSDs) with mixed factor (unequal) levels are more suitable as these designs explore more regions in the design space but it is hard to get rotatable designs with a given level of asymmetry. When experimenting with unequal factor levels via asymmetric second order rotatable design (ASORDs), the lack of fit of the model may become significant which ultimately leads to the estimation of parameters based on a higher (or third) order model. Experimenting with a new third order rotatable design (TORD) in such a situation would be expensive as the responses observed from the first stage runs would be kept underutilized. In this paper, we propose a method of constructing asymmetric TORD by sequentially augmenting some additional points to the ASORDs without discarding the runs in the first stage. The proposed designs will be more economical to obtain the optimum response as the design in the first stage can be used to fit the second order model and with some additional runs, third order model can be fitted without discarding the initial design. [ABSTRACT FROM AUTHOR]

Statistical surveillance for covariance matrices has attracted increasing attention recently. Many approaches have been developed for monitoring general shifts that are arbitrary deviations, as well as sparse shifts occurring in only a few elements. This paper considers directional shifts that occur in only one independent parameter, which is common if the process is relatively stable. A directional covariance matrix control chart is proposed, which fully exploits directional shift information and borrows the strong power of likelihood ratio test. Therefore, this chart provides a powerful tool for monitoring covariance matrices. In addition, the proposed chart does not require specifying the regularisation parameter, and it enjoys a concise quadratic form, thereby easy to implement. Furthermore, this chart naturally leads to a diagnostic prescription for identifying the shifting element in the covariance matrix. Simulation results have demonstrated the efficiency of the suggested control chart and its accompanying diagnostic scheme. [ABSTRACT FROM AUTHOR]

Classical continuous goodness-of-fit (GOF) testing is employed for examining whether the data come from an assumed parametric model. In many cases, GOF tests assume a uniform null distribution and examine extreme values of the order statistics of the samples. Many of these statistics can be expressed by a function of the order statistics and the p-values amount to a joint probability statement based on the uniform order statistics. In this paper, we utilize Steck's recursion method and propose two high precision computing algorithms to compute the p-values for these GOF statistics. The numerical difficulties in implementing Steck's method are discussed and compared with solutions provided in high precision libraries. [ABSTRACT FROM AUTHOR]

We consider the problem of estimating the maximum posterior probability (MAP) state sequence for a finite state and finite emission alphabet hidden Markov model (HMM) in the Bayesian setup, where both emission and transition matrices have Dirichlet priors. We study a training set consisting of thousands of protein alignment pairs. The training data is used to set the prior hyperparameters for Bayesian MAP segmentation. Since the Viterbi algorithm is not applicable any more, there is no simple procedure to find the MAP path, and several iterative algorithms are considered and compared. The main goal of the paper is to test the Bayesian setup against the frequentist one, where the parameters of HMM are estimated using the training data. [ABSTRACT FROM AUTHOR]

In this application note paper, we propose and examine the performance of a Bayesian approach for a homoscedastic nonlinear regression (NLR) model assuming errors with two-piece scale mixtures of normal (TP-SMN) distributions. The TP-SMN is a large family of distributions, covering both symmetrical/ asymmetrical distributions as well as light/heavy tailed distributions, and provides an alternative to another well-known family of distributions, called scale mixtures of skew-normal distributions. The proposed family and Bayesian approach provides considerable flexibility and advantages for NLR modelling in different practical settings. We examine the performance of the approach using simulated and real data. [ABSTRACT FROM AUTHOR]

In this paper, we discuss the bivariate Birnbaum-Saunders accelerated lifetime model, in which we have modeled the dependence structure of bivariate survival data through the use of frailty models. Specifically, we propose the bivariate model Birnbaum-Saunders with the following frailty distributions: gamma, positive stable and logarithmic series. We present a study of inference and diagnostic analysis for the proposed model, more concisely, are proposed a diagnostic analysis based in local influence and residual analysis to assess the fit model, as well as, to detect influential observations. In this regard, we derived the normal curvatures of local influence under different perturbation schemes and we performed some simulation studies for assessing the potential of residuals to detect misspecification in the systematic component, the presence in the stochastic component of the model and to detect outliers. Finally, we apply the methodology studied to real data set from recurrence in times of infections of 38 kidney patients using a portable dialysis machine, we analyzed these data considering independence within the pairs and using the bivariate Birnbaum-Saunders accelerated lifetime model, so that we could make a comparison and verify the importance of modeling dependence within the times of infection associated with the same patient. [ABSTRACT FROM AUTHOR]

In this paper, we consider the ultrahigh-dimensional sufficient dimension reduction (SDR) for censored data and measurement error in covariates. We first propose the feature screening procedure based on censored data and the covariates subject to measurement error. With the suitable correction of mismeasurement, the error-contaminated variables detected by the proposed feature screening procedure are the same as the truly important variables. Based on the selected active variables, we develop the SDR method to estimate the central subspace and the structural dimension with both censored data and measurement error incorporated. The theoretical results of the proposed method are established. Simulation studies are reported to assess the performance of the proposed method. The proposed method is implemented to NKI breast cancer data. [ABSTRACT FROM AUTHOR]

This paper deals with the problem of maximum likelihood and Bayesian estimation of stress–strength reliability involving paired observation with ties using bivariate exponentiated half-logistic distribution. This problem is of importance because in some real applications the strength of the component is highly dependent on the stress experienced by it. A bivariate extension of exponentiated half-logistic is discussed and an expression for the stress–strength reliability is obtained. This model is also useful to analyse data having the unusual feature of having a number of pairs with tied scores, even when the scores are continuous. The maximum likelihood estimate and interval estimate of the stress–strength reliability has been developed. The Bayes estimates of the stress–strength reliability under squared error loss function are obtained using importance sampling technique. Simulation studies are conducted to evaluate the performance of maximum likelihood and Bayes estimates. Two real-life data sets are analysed to numerically illustrate the usefulness of the developed method. [ABSTRACT FROM AUTHOR]

In this paper, 91 different tests for exponentiality are reviewed. Some of the tests are universally consistent while others are against some special classes of life distributions. Power performances of 40 of these different tests for exponentiality of datasets are compared through extensive Monte Carlo simulations. The comparisons are conducted for different sample sizes of 10, 25, 50 and 100 for different groups of distributions according to the shape of their hazard functions at 5 percent level of significance. Also, the techniques are applied to two real-world datasets and a measure of power is employed for the comparison of the tests. The results show that some tests which are very good under one group of alternative distributions are not so under another group. Also, some tests maintained relatively high power over all the groups of alternative distributions studied while some others maintained poor power performances over all the groups of alternative distributions. Again, the result obtained from real-world datasets agree completely with those of the simulation studies. [ABSTRACT FROM AUTHOR]

In this paper, we suggest a technique to quantify model risk, particularly model misspecification for binary response regression problems found in financial risk management, such as in credit risk modelling. We choose the probability of default model as one instance of many other credit risk models that may be misspecified in a financial institution. By way of illustrating the model misspecification for probability of default, we carry out quantification of two specific statistical predictive response techniques, namely the binary logistic regression and complementary log–log. The maximum likelihood estimation technique is employed for parameter estimation. The statistical inference, precisely the goodness of fit and model performance measurements, are assessed. Using the simulation dataset and Taiwan credit card default dataset, our finding reveals that with the same sample size and very small simulation iterations, the two techniques produce similar goodness-of-fit results but completely different performance measures. However, when the iterations increase, the binary logistic regression technique for balanced dataset reveals prominent goodness of fit and performance measures as opposed to the complementary log–log technique for both simulated and real datasets. [ABSTRACT FROM AUTHOR]

This paper is concerned with testing and dating structural breaks in the dependence structure of multivariate time series. We consider a cumulative sum (CUSUM) type test for constant copula-based dependence measures, such as Spearman's rank correlation and quantile dependencies. The asymptotic null distribution is not known in closed form and critical values are estimated by an i.i.d. bootstrap procedure. We analyze size and power properties in a simulation study under different dependence measure settings, such as skewed and fat-tailed distributions. To date breakpoints and to decide whether two estimated break locations belong to the same break event, we propose a pivot confidence interval procedure. Finally, we apply the test to the historical data of 10 large financial firms during the last financial crisis from 2002 to mid-2013. [ABSTRACT FROM AUTHOR]

This paper discusses methods for clustering a continuous covariate in a survival analysis model. The advantages of using a categorical covariate defined from discretizing a continuous covariate (via clustering) is (i) enhanced interpretability of the covariate's impact on survival and (ii) relaxing model assumptions that are usually required for survival models, such as the proportional hazards model. Simulations and an example are provided to illustrate the methods. [ABSTRACT FROM AUTHOR]

Students' migration mobility is the new form of migration: students migrate to improve their skills and become more valued for the job market. The data regard the migration of Italian Bachelors who enrolled at Master Degree level, moving typically from poor to rich areas. This paper investigates the migration and other possible determinants on the Master Degree students' performance. The Clustering of Effects approach for Quantile Regression Coefficients Modelling has been used to cluster the effects of some variables on the students' performance for three Italian macro-areas. Results show evidence of similarity between Southern and Centre students, with respect to the Northern ones. [ABSTRACT FROM AUTHOR]

Modeling cyber threats, such as the computer malicious software (malware) propagation dynamics in cyberspace, is an important research problem because models can deepen our understanding of dynamical cyber threats. In this paper, we study the statistical modeling of the macro-level evolution of dynamical cyber attacks. Specifically, we propose a Bayesian structural time series approach for modeling the computer malware propagation dynamics in cyberspace. Our model not only possesses the parsimony property (i.e. using few model parameters) but also can provide the predictive distribution of the dynamics by accommodating uncertainty. Our simulation study shows that the proposed model can fit and predict the computer malware propagation dynamics accurately, without requiring to know the information about the underlying attack-defense interaction mechanism and the underlying network topology. We use the model to study the propagation of two particular kinds of computer malware, namely the Conficker and Code Red worms, and show that our model has very satisfactory fitting and prediction accuracies. [ABSTRACT FROM AUTHOR]

Ranked set sampling (RSS) design as a cost-effective sampling is a powerful tool in situations where measuring the variable of interest is costly and time-consuming; however, ranking information about sampling units can be obtained easily through inexpensive and easy to measure characteristics at little or no cost. In this paper, we study RSS data for analysis of an ordinal population. First, we compare the problem of non-representative extreme samples under RSS and commonly-used simple random sampling. Using RSS data with tie information, we propose non-parametric and maximum likelihood estimators for population parameters. Through extensive numerical studies, we investigate the effect of various factors including ranking ability, tie generating mechanisms, the number of categories and population setting on the performance of the estimators. Finally, we apply the proposed methods to the bone disorder data to estimate the proportions of patients with osteopenia and osteoporosis status. [ABSTRACT FROM AUTHOR]

In this paper, we first present a characterization of exponential distribution based on the extropy of record values and next introduce a goodness-of-fit test for exponentiality. Monte Carlo simulation is used to compute the critical values of our proposed test for different sample sizes and significance levels. To show the advantage of the proposed test, we adopt 58 competitor tests and compute the adjusted power against different alternatives with distinct types of hazard function. The power results show that our proposed test has superior adjusted power if the alternatives have increasing failure rates or bathtub decreasing-increasing failure rates, especially when the sample size is small. Finally, three real examples are used to illustrate the applicability and robustness of our proposed test by monitoring the p-values of the tests. [ABSTRACT FROM AUTHOR]

In this paper, we consider the classification problem of functional data including the sleeping energy expenditure (SEE) data, focusing on functional classification. Many existing classification rules are not effective in distinguishing the two classes of SEE data, because the trajectories of each observation have very different patterns for each class. It is often observed that some aspect of data such as the variability of paths is helpful in classification of functional data. To reflect this issue, we introduce a variable measuring the length of path in functional data and then propose a logistic model with fused lasso that considers the behavior of fluctuation of path as well as local correlations within each path. Our proposed model shows a significant improvement over some models used in the existing literature on the classification accuracy rate of functional data such as SEE data. We carry out simulation studies to show the finite sample performance and the gain that it makes in comparison with fused lasso without considering path length. With two more real datasets studied in some existing literature, we demonstrate that the new model achieves better or similar accuracy rate than the best accuracy rates reported in those studies. [ABSTRACT FROM AUTHOR]

This paper describes a comprehensive survival analysis for the inverse Gaussian distribution employing Bayesian and Fiducial approaches. It focuses on making inferences on the inverse Gaussian (IG) parameters μ and λ and the average remaining time of censored units. A flexible Gibbs sampling approach applicable in the presence of censoring is discussed and illustrations with Type II, progressive Type II, and random rightly censored observations are included. The analyses are performed using both simulated IG data and empirical data examples. Further, the bootstrap comparisons are made between the Bayesian and Fiducial estimates. It is concluded that the shape parameter ( ϕ = λ / μ) of the inverse Gaussian distribution has the most impact on the two analyses, Bayesian vs. Fiducial, and so does the size of censoring in data to a lesser extent. Overall, both these approaches are effective in estimating IG parameters and the average remaining lifetime. The suggested Gibbs sampler allowed a great deal of flexibility in implementation for all types of censoring considered. [ABSTRACT FROM AUTHOR]

Seoul, the capital city of Korea with over 10 million residents, has been experiencing serious air pollution problems. Previous studies on source apportionment of PM2.5 in Seoul are based on measurements of chemical compositions of PM2.5 from a single monitoring site. In this paper, we analyse PM2.5 concentration data collected from multiple sites in 24 districts of Seoul and estimate regional source profiles using Bayesian multivariate receptor model. The regional source profiles provide information for the identification of major PM2.5 sources as well as the regions relatively more seriously affected by each source than other regions. These regional characteristics relevant to PM2.5 can help establish effective, customised, region-specific PM2.5 control strategies for each region rather than general strategies that apply to every region of Seoul. [ABSTRACT FROM AUTHOR]

In this paper, we investigate robust parameter estimation and variable selection for binary regression models with grouped data. We investigate estimation procedures based on the minimum-distance approach. In particular, we employ minimum Hellinger and minimum symmetric chi-squared distances criteria and propose regularized minimum-distance estimators. These estimators appear to possess a certain degree of automatic robustness against model misspecification and/or for potential outliers. We show that the proposed non-penalized and penalized minimum-distance estimators are efficient under the model and simultaneously have excellent robustness properties. We study their asymptotic properties such as consistency, asymptotic normality and oracle properties. Using Monte Carlo studies, we examine the small-sample and robustness properties of the proposed estimators and compare them with traditional likelihood estimators. We also study two real-data applications to illustrate our methods. The numerical studies indicate the satisfactory finite-sample performance of our procedures. [ABSTRACT FROM AUTHOR]