Assessment of modeling longitudinal binary data based on graphical methods

Longitudinal categorical data are commonly applied in a variety of fields and are frequently analyzed by generalized estimating equation (GEE) method. Prior to making further inference based on the GEE model, the assessment of model fit is crucial. Graphical techniques have long been in widespread use for assessing the model adequacy. We develop alternative graphical approaches utilizing plots of marginal model-checking condition and local mean deviance to assess the GEE model with logit link for longitudinal binary responses. The applications of the proposed procedures are illustrated through two longitudinal binary datasets.     

Asymptotic relative efficiency of the linear discriminant function under partial nonrandom classification of the training data

This paper considers the problem where the linear discriminant rule is formed from training data that are only partially classified with respect to the two groups of origin. A further complication is that the data of unknown origin do not constitute an observed random sample from a mixture of the two under- lying groups. Under the assumption of a homoscedastic normal model, the overall error rate of the sample linear discriminant rule formed by maximum likelihood from the partially classified training data is derived up to and including terms of the first order in the case of univariate feature data. This first- order expansion of the sample rule so formed is used to define its asymptotic efficiency relative to the rule formed from a completely classified random training set and also to the rule formed from a completely unclassified random set.     

Identifying Skeleton Curves in Noisy Data

This article presents a Bayesian method to reconstruct the centerline in noisy data using B-spline curves. The method is illustrated on simulated two- and three-dimensional data and is applied to recover the centerline of the colon in single photon emission computed tomography images.     

On a system of nonhomogeneous components sharing a common frailty

The components of a reliability system subjected to a common random environment usually have dependent lifetimes. This paper studies the stochastic properties of such a system with lifetimes of the components following multivariate frailty models and multivariate mixed proportional reversed hazard rate (PRHR) models, respectively. Through doing stochastic comparison, we devote to throwing a new light on how the random environment affects the number of working components of a reliability system and on assessing the performance of a k-out-of-n system.     

New calibration estimator based on two auxiliary variables in stratified sampling

In recent years, calibration estimation has become an important field of research in survey sampling. This paper proposes a new calibration estimator for the population mean in the presence of two auxiliary variables in stratified sampling. The theory of new calibration estimator is given and optimum calibration weights are derived. A simulation study is carried out to performance of the proposed calibration estimator over other existing calibration estimators. The results reveal that the proposed calibration estimators are more efficient than other existing calibration estimators in stratified sampling.     

Robust block bootstrap panel predictability tests*

This article develops two block bootstrap-based panel predictability test procedures that are valid under very general conditions. Some of the allowable features include cross-sectional dependence, heterogeneous predictive slopes, persistent predictors, and complex error dynamics, including cross-unit endogeneity. While the first test procedure tests if there is any predictability at all, the second procedure determines the units for which predictability holds in case of a rejection by the first. A weak unit root framework is adopted to allow persistent predictors, and a novel theory is developed to establish asymptotic validity of the proposed bootstrap. Simulations are used to evaluate the performance of our tests in small samples, and their implementation is illustrated through an empirical application to stock returns.     

An Extended Single‐index Model with Missing Response at Random

An extended single‐index model is considered when responses are missing at random. A three‐step estimation procedure is developed to define an estimator for the single‐index parameter vector by a joint estimating equation. The proposed estimator is shown to be asymptotically normal. An algorithm for computing this estimator is proposed. This algorithm only involves one‐dimensional nonparametric smoothers, thereby avoiding the data sparsity problem caused by high model dimensionality. Some simulation studies are conducted to investigate the finite sample performances of the proposed estimators.     

Multivariate log-skewed distributions with normal kernel and their applications

We introduce two classes of multivariate log-skewed distributions with normal kernel: the log canonical fundamental skew-normal (log-CFUSN) and the log unified skew-normal. We also discuss some properties of the log-CFUSN family of distributions. These new classes of log-skewed distributions include the log-normal and multivariate log-skew normal families as particular cases. We discuss some issues related to Bayesian inference in the log-CFUSN family of distributions, mainly we focus on how to model the prior uncertainty about the skewing parameter. Based on the stochastic representation of the log-CFUSN family, we propose a data augmentation strategy for sampling from the posterior distributions. This proposed family is used to analyse the US national monthly precipitation data. We conclude that a high-dimensional skewing function lead to a better model fit.     

Conditional and marginal estimates in case-control family data – extensions and sensitivity analyses

This work considers two specific estimation techniques for the family-specific proportional hazards model and for the population-averaged proportional hazards model. So far, these two estimation procedures were presented and studied under the gamma frailty distribution mainly because of its simple interpretation and mathematical tractability. Modifications of both procedures for other frailty distributions, such as the inverse Gaussian, positive stable and a specific case of discrete distribution, are presented. By extensive simulations, it is shown that under the family-specific proportional hazards model, the gamma frailty model appears to be robust to frailty distribution mis-specification in both bias and efficiency loss in the marginal parameters. The population-averaged proportional hazards model, is found to be robust under the gamma frailty model mis-specification only under moderate or weak dependency within cluster members.     

A goodness of fit test for the multilevel proportional odds model

It is essential to test the goodness of fit of the model before making inferences based on it. Multilevel modeling of ordinal categorical responses is not as developed as for continuous responses. Assessing model adequacy in terms of the goodness of fit with ordinal categorical responses is still being developed and no satisfactory tests are available so far. As a consequence of that, this study concentrates on developing such a goodness of fit test for Multilevel Proportional Odds models and to study the properties of the test.