Barriers to Managing Coronary Artery Disease Among Older Women

Numerous barriers to managing coronary artery disease (CAD) among older women are reported in the literature; however, few studies adjust for demographic and health status differences. A survey assessing barriers and other factors was distributed to a stratified random sampling of older women with CAD. Factor analysis and multiple logistic regression procedures were used to estimate the impact of these issues on receiving a CAD-related office visit. The most problematic barriers included denial and low health literacy. Efforts to promote patient awareness of heart health and better communication between patients and clinicians may alleviate these barriers.     

Cities evolution tree and applications to predicting urban growth

China’s unprecedented urbanization and economic modernization have brought about significant growth and land use changes in cities and their surroundings. Traditional analyses of the relationship between urbanization and land use change have not fully incorporated important factors like city function or developmental stage. We define and apply a new visualization and analysis tool—an evolution tree—showing how an understanding of evolutionary pathways taken by growing cities contributes to visualization and forecasting of urban development and spatial expansion. We clustered 253 Chinese cities by type and developmental stage, arranging them in a “tree” structure, with types grouped as branches and leaf position along branches reflecting developmental stage. Urban growth rates are higher for business-oriented, manufacturing-oriented, and functionally comprehensive cities than for those oriented toward tourism. Extent of spatial growth is dominated by city type and developmental stage. Cities evolve along typical pathways, punctuated by occasional leap developments across functional types.     

Adaptive modeling of the human-environment relationship applied to estimation of the population carrying capacity in an earthquake zone

Natural catastrophes such as earthquakes can, in addition to causing loss of life, disrupt the urbanization process through the need for forced population redistribution and the modification of resource and environmental carrying capacity. The population carrying capacity (PCC) of an altered environment following an earthquake is a crucial determinant in the relocation of displaced persons. We use data adaptive methods to model the correlation between the physical environment and human population density in estimating PCC in areas affected by the 2008 Wenchuan earthquake. Comparing actual population distributions with ideal population distributions allows for the identification of villages where population exceeds PCC, or conversely, areas where the environment can support a higher population. Such a comparison can provide the basis for a relocation plan, a critical element of post-catastrophe policy-making.     

Feature screening in ultrahigh-dimensional additive Cox model

The additive Cox model is flexible and powerful for modelling the dynamic changes of regression coefficients in the survival analysis. This paper is concerned with feature screening for the additive Cox model with ultrahigh-dimensional covariates. The proposed screening procedure can effectively identify active predictors. That is, with probability tending to one, the selected variable set includes the actual active predictors. In order to carry out the proposed procedure, we propose an effective algorithm and establish the ascent property of the proposed algorithm. We further prove that the proposed procedure possesses the sure screening property. Furthermore, we examine the finite sample performance of the proposed procedure via Monte Carlo simulations, and illustrate the proposed procedure by a real data example.     

Likelihood ratio order of parallel systems under multiple-outlier models

This paper studies the likelihood ratio ordering of parallel systems under multiple-outlier models. We introduce a partial order, the so-called θ-order, and show that the θ-order between the parameter vectors of the parallel systems implies the likelihood ratio order between the systems.     

A new algorithm to estimate monotone nonparametric link functions and a comparison with parametric approach

The generalized linear model (GLM) is a class of regression models where the means of the response variables and the linear predictors are joined through a link function. Standard GLM assumes the link function is fixed, and one can form more flexible GLM by either estimating the flexible link function from a parametric family of link functions or estimating it nonparametically. In this paper, we propose a new algorithm that uses P-spline for nonparametrically estimating the link function which is guaranteed to be monotone. It is equivalent to fit the generalized single index model with monotonicity constraint. We also conduct extensive simulation studies to compare our nonparametric approach for estimating link function with various parametric approaches, including traditional logit, probit and robit link functions, and two recently developed link functions, the generalized extreme value link and the symmetric power logit link. The simulation study shows that the link function estimated nonparametrically by our proposed algorithm performs well under a wide range of different true link functions and outperforms parametric approaches when they are misspecified. A real data example is used to illustrate the results.     

Estimation of Sobol's sensitivity indices under generalized linear models

We derive explicit formulas for Sobol's sensitivity indices (SSIs) under the generalized linear models (GLMs) with independent or multivariate normal inputs. We argue that the main-effect SSIs provide a powerful tool for variable selection under GLMs with identity links under polynomial regressions. We also show via examples that the SSI-based variable selection results are similar to the ones obtained by the random forest algorithm but without the computational burden of data permutation. Finally, applying our results to the problem of gene network discovery, we identify through the SSI analysis of a public microarray dataset several novel higher-order gene–gene interactions missed out by the more standard inference methods. The relevant functions for SSI analysis derived here under GLMs with identity, log, and logit links are implemented and made available in the R package Sobol sensitivity.     

HDDA: DataSifter: statistical obfuscation of electronic health records and other sensitive datasets

There are no practical and effective mechanisms to share high-dimensional data including sensitive information in various fields like health financial intelligence or socioeconomics without compromising either the utility of the data or exposing private personal or secure organizational information. Excessive scrambling or encoding of the information makes it less useful for modelling or analytical processing. Insufficient preprocessing may compromise sensitive information and introduce a substantial risk for re-identification of individuals by various stratification techniques. To address this problem, we developed a novel statistical obfuscation method (DataSifter) for on-the-fly de-identification of structured and unstructured sensitive high-dimensional data such as clinical data from electronic health records (EHR). DataSifter provides complete administrative control over the balance between risk of data re-identification and preservation of the data information. Simulation results suggest that DataSifter can provide privacy protection while maintaining data utility for different types of outcomes of interest. The application of DataSifter on a large autism dataset provides a realistic demonstration of its promise practical applications.     

Bayesian bridge-randomized penalized quantile regression estimation for linear regression model with AP(q) perturbation

Bridge penalized regression has many desirable statistical properties such as unbiasedness, sparseness as well as ‘oracle’. In Bayesian framework, bridge regularized penalty can be implemented based on generalized Gaussian distribution (GGD) prior. In this paper, we incorporate Bayesian bridge-randomized penalty and its adaptive version into the quantile regression (QR) models with autoregressive perturbations to conduct Bayesian penalization estimation. Employing the working likelihood of the asymmetric Laplace distribution (ALD) perturbations, the Bayesian joint hierarchical models are established. Based on the mixture representations of the ALD and generalized Gaussian distribution (GGD) priors of coefficients, the hybrid algorithms based on Gibbs sampler and Metropolis-Hasting sampler are provided to conduct fully Bayesian posterior estimation. Finally, the proposed Bayesian procedures are illustrated by some simulation examples and applied to a real data application of the electricity consumption.     

Tail asymptotic of discounted aggregate claims with compound dependence under risky investment

This paper considers the tail asymptotic of discounted aggregate claims with compound dependence under risky investment. The price of risky investment is modeled by a geometric Lévy process, while claims are modeled by a one-sided linear process whose innovations further obeying a so-called upper tail asymptotic independence. When the innovations are heavy tailed, we derive some uniform asymptotic formulas. The results show that the linear dependence has significant impact on the tail asymptotic of discounted aggregate claims but the upper tail asymptotic independence is negligible.