Joint modeling of correlated binary outcomes: HIV-1 and HSV-2 co-infection

Herpes Simplex Virus Type 2 (HSV-2) facilitates the sexual acquisition and transmission of HIV-1 infection and is highly prevalent in most regions experiencing severe HIV epidemics. In sub-Saharan Africa, where HIV infection is a public health burden, the prevalence of HSV-2 is substantially high. The high prevalence of HSV-2 and the association between HSV-2 infection and HIV-1 acquisition could play a significant role in the spread of HIV-1 in the region. The objective of our study was to identify risk factors for HSV-2 and HIV-1 infections among men in sub-Saharan Africa. We used a joint response model that accommodates the interdependence between the two infections in assessing their risk factors. Simulation studies show superiority of the joint response model compared to the traditional models which ignore the dependence between the two infections. We found higher odds of having HSV-2/HIV-1 among older men, in men who had multiple sexual partners, abused alcohol, or reported symptoms of sexually transmitted infections. These findings suggest that interventions that identify and control the risk factors of the two infections should be part of HIV-1 prevention programs in sub-Saharan Africa where antiretroviral therapy is not readily available.     

On the probability of ruin in a continuous risk model with two types of delayed claims

ABSTRACT

This article studies a risk model involving one type of main claims and two types of by-claims, which is an extension of the general risk model with delayed claims. We suppose that every main claim may not induce any by-claims or may induce one by-claim belonging to one of the two types of by-claims with a certain probability. In addition, assume that the by-claim and its associated main claim may occur at the same time and that the occurrence of the by-claim may be delayed. An integro-differential equation system for survival probabilities is derived by using two auxiliary risk models. The expression of the survival probability is obtained by applying Laplace transforms and Rouché theorem. Furthermore, we provide a method for solving the survival probability when the two by-claim amounts satisfy different exponential distributions. As a special case, an explicit expression of survival probability is given when all the claim amounts obey the same exponential distribution. Finally, numerical results are provided to examine the proposed method.     

A litter-based approach to risk assessment in developmental toxicity studies via a power family of completely monotone functions

Summary.  A new class of distributions for exchangeable binary data is proposed that originates from modelling the joint success probabilities of all orders by a power family of completely monotone functions. The distribution proposed allows flexible modelling of the dose–response relationship for both the marginal response probability and the pairwise odds ratio and is especially well suited for a litter-based approach to risk assessment. Specifically, the risk of at least one adverse response within a litter takes on a simple form under the distribution proposed and can be reduced further to a generalized linear model if a complementary log–log-link function is used. Existing distributions such as the beta–binomial or folded logistic functions have a tendency to assign too much probability to zero, leading to an underestimation of the risk that at least one foetus is affected and an overestimation of the safe dose. The distribution proposed does not suffer from this problem. With the aid of symbolic differentiation, the distribution proposed can be fitted easily and quickly via the method of scoring. The usefulness of the class of distributions proposed and its superiority over existing distributions are demonstrated in a series of examples involving developmental toxicology and teratology data.     

Variable-Selection-Based Epidemic Disease Diagnosis

Epidemic surveillance in a community involves monitoring infection trend, triggering alarms before outbreaks, and identifying sources and paths of disease transmission. Algorithms for outbreak detection that are derived from industrial statistical process control (SPC) and scan statistics have been reported in the literature, but there are relatively few methods reported for identifying transmission paths. In this work, we propose an expanded spatial-temporal (EST) model for identifying infection sources. Three dimensional information, subject, location, and time, are expanded into a two-dimensional space by dividing the time horizon into segments and multiplying each segment by the locations. Based on the EST model, we further propose a variable-selection algorithm to identify potential location/time combinations as sources of infection, and thus achieve diagnosis. Numerical simulations show that the proposed scheme is effective in locating infection sources.     

Estimating the transmission probability of human immunodeficiency virus in injecting drug users in Thailand

We estimate the transmission probability for the human immunodeficiency virus from seroconversion data of a cohort of injecting drug users (IDUs) in Thailand. The transmission probability model developed accounts for interval censoring and incorporates each IDU's reported frequency of needle sharing and injecting acts. Using maximum likelihood methods, the per needle sharing act transmission probability estimate between infectious and susceptible IDUs is 0.008. The effects of covariates, disease dynamics, mismeasured exposure information and the uncertainty of the disease prevalence on the transmission probability estimate are considered.     

Semiparametric Proportional Mean Residual Life Model With Censoring Indicators Missing at Random

For right-censored survival data, the information that whether the observed time is survival or censoring time is frequently lost. This is the case for the competing risk data. In this article, we consider statistical inference for the right-censored survival data with censoring indicators missing at random under the proportional mean residual life model. Simple and augmented inverse probability weighted estimating equation approaches are developed, in which the nonmissingness probability and some unknown conditional expectations are estimated by the kernel smoothing technique. The asymptotic properties of all the proposed estimators are established, while extensive simulation studies demonstrate that our proposed methods perform well under the moderate sample size. At last, the proposed method is applied to a data set from a stage II breast cancer trial.     

ESTIMATING WITHIN-SCHOOL CONTACT NETWORKS TO UNDERSTAND INFLUENZA TRANSMISSION

Many epidemic models approximate social contact behavior by assuming random mixing within mixing groups (e.g., homes, schools, and workplaces). The effect of more realistic social network structure on estimates of epidemic parameters is an open area of exploration. We develop a detailed statistical model to estimate the social contact network within a high school using friendship network data and a survey of contact behavior. Our contact network model includes classroom structure, longer durations of contacts to friends than non-friends and more frequent contacts with friends, based on reports in the contact survey. We performed simulation studies to explore which network structures are relevant to influenza transmission. These studies yield two key findings. First, we found that the friendship network structure important to the transmission process can be adequately represented by a dyad-independent exponential random graph model (ERGM). This means that individual-level sampled data is sufficient to characterize the entire friendship network. Second, we found that contact behavior was adequately represented by a static rather than dynamic contact network. We then compare a targeted antiviral prophylaxis intervention strategy and a grade closure intervention strategy under random mixing and network-based mixing. We find that random mixing overestimates the effect of targeted antiviral prophylaxis on the probability of an epidemic when the probability of transmission in 10 minutes of contact is less than 0.004 and underestimates it when this transmission probability is greater than 0.004. We found the same pattern for the final size of an epidemic, with a threshold transmission probability of 0.005. We also find random mixing overestimates the effect of a grade closure intervention on the probability of an epidemic and final size for all transmission probabilities. Our findings have implications for policy recommendations based on models assuming random mixing, and can inform further development of network-based models.     

A Discrete-Time Multilevel Mixture Model for Event History Data with Long-Term Survivors,with an Application to an Analysis of Contraceptive Sterilization in Bangladesh

Event history models typically assume that the entire population is at risk of experiencing the event of interest throughout the observation period. However, there will often be individuals, referred to as long-term survivors, who may be considered a priori to have a zero hazard throughout the study period. In this paper, a discrete-time mixture model is proposed in which the probability of long-term survivorship and the timing of event occurrence are modelled jointly. Another feature of event history data that often needs to be considered is that they may come from a population with a hierarchical structure. For example, individuals may be nested within geographical regions and individuals in the same region may have similar risks of experiencing the event of interest due to unobserved regional characteristics. Thus, the discrete-time mixture model is extended to allow for clustering in the likelihood and timing of an event within regions. The model is further extended to allow for unobserved individual heterogeneity in the hazard of event occurrence. The proposed model is applied in an analysis of contraceptive sterilization in Bangladesh. The results show that a woman's religion and education level affect her probability of choosing sterilization, but not when she gets sterilized. There is also evidence of community-level variation in sterilization timing, but not in the probability of sterilization.     

Selection of the splined variables and convergence rates in a partial spline model

A method based on the principle of unbiased risk estimation is used to select the splined variables in an exploratory partial spline model proposed by Wahba (1985). The probability of correct selection based on the proposed procedure is discussed under regularity conditions. Furthermore, the resulting estimate of the regression function achieves the optimal rates of convergence over a general class of smooth regression functions (Stone 1982) when its underlying smoothness condition is not known.     

Nonparametric estimation of the sojourn time distributions for a multipath model

Summary. We use a multipath (multistate) model to describe data with multiple end points. Statistical inference based on the intermediate end point is challenging because of the problems of nonidentifiability and dependent censoring. We study nonparametric estimation for the path probability and the sojourn time distributions between the states. The methodology proposed can be applied to analyse cure models which account for the competing risk of death. Asymptotic properties of the estimators proposed are derived. Simulation shows that the methods proposed have good finite sample performance. The methodology is applied to two data sets.     

信用风险、融资担保和动态保费估值

基于市场化信息和风险中性的概念,度量了有担保贷款的边际违约概率和累积违约概率,确定了贷款担保风险的精算现值,根据市场信用价差的变化给出了动态保费每期的调整幅度,并利用数值模拟进行了担保费率的比较静态分析,最后根据实际的担保数据给出了动态保费的实证检验.结果显示,实际的违约支付非常接近于动态保费的估值,证明动态保费估值模型是一个简单、可行和实用的定价模型.     

Modeling potential time to event data with competing risks

Patients receiving radical prostatectomy are at risk of metastasis or prostate cancer related death, and often need repeated clinical evaluations to determine whether additional adjuvant or salvage therapies are needed. Since the prostate cancer is a slowly progressing disease, and these additional therapies come with significant side effects, it is important for clinical decision making purposes to estimate a patient’s risk of cancer metastasis, in the presence of a competing risk by death, under the hypothetical condition that the patient does not receive any additional therapy. In observational studies, patients may receive additional therapy by choice; the time to metastasis without any therapy is often a potential outcome and not always observed. We study the competing risks model of Fine and Gray (J Am Stat Assoc, 94:496–509, 1999) with adjustment for treatment choice by inverse probability censoring weighting (IPCW). The model can be fit using standard software for partial likelihood with double IPCW weights. The proposed methodology is used in a prostate cancer study to predict the post-prostatectomy cumulative incidence probability of cancer metastasis without additional adjuvant or salvage therapies.     

Proportional hazards model for competing risks data with missing cause of failure

We consider the semiparametric proportional hazards model for the cause-specific hazard function in analysis of competing risks data with missing cause of failure. The inverse probability weighted equation and augmented inverse probability weighted equation are proposed for estimating the regression parameters in the model, and their theoretical properties are established for inference. Simulation studies demonstrate that the augmented inverse probability weighted estimator is doubly robust and the proposed method is appropriate for practical use. The simulations also compare the proposed estimators with the multiple imputation estimator of Lu and Tsiatis (2001). The application of the proposed method is illustrated using data from a bone marrow transplant study.     

A new policy for designing acceptance sampling plan based on Bayesian inference in the presence of inspection errors

The purpose of this article is to present a new policy for designing an acceptance sampling plan based on the minimum proportion of the lot that should be inspected in the presence of inspection errors. It is assumed that inspection is not perfect and every defective item cannot be detected with complete certainty. The Bayesian method is used for obtaining the probability distribution function of the number of defective items in the lot. To design this model, two constraints of producer risk and consumer risk are considered during the inspection process by using two specified points on operating characteristic curve. In order to illustrate the application of the proposed model, an example is presented. In addition, a sensitivity analysis is performed to analyze the model performance under different scenarios of process parameters and the results are elaborated. Finally, the efficiency of the proposed model is compared with the sampling method of Spencer and Kevan de Lopez (2017) at the same conditions.     

Efficient heterogeneous sampling for stochastic simulation with an illustration in health care applications

In modeling disease transmission, contacts are assumed to have different infection rates. A proper simulation must model the heterogeneity in the transmission rates. In this article, we present a computationally efficient algorithm that can be applied to a population with heterogeneous transmission rates. We conducted a simulation study to show that the algorithm is more efficient than other algorithms for sampling the disease transmission in a subset of the heterogeneous population. We use a valid stochastic model of pandemic influenza to illustrate the algorithm and to estimate the overall infection attack rates of influenza A (H1N1) in a Canadian city.     

Individual-level modeling of the spread of influenza within households

A class of individual-level models (ILMs) outlined by R. Deardon et al., [Inference for individual level models of infectious diseases in large populations, Statist. Sin. 20 (2010), pp. 239–261] can be used to model the spread of infectious diseases in discrete time. The key feature of these ILMs is that they take into account covariate information on susceptible and infectious individuals as well as shared covariate information such as geography or contact measures. Here, such ILMs are fitted in a Bayesian framework using Markov chain Monte Carlo techniques to data sets from two studies on influenza transmission within households in Hong Kong during 2008 to 2009 and 2009 to 2010. The focus of this paper is to estimate the effect of vaccination on infection risk and choose a model that best fits the infection data.     

Bayeslan Credit Ratings

In this contribution we aim at improving ordinal variable selection in the context of causal models for credit risk estimation. In this regard, we propose an approach that provides a formal inferential tool to compare the explanatory power of each covariate and, therefore, to select an effective model for classification purposes. Our proposed model is Bayesian nonparametric thus keeps the amount of model specification to a minimum. We consider the case in which information from the covariates is at the ordinal level. A noticeable instance of this regards the situation in which ordinal variables result from rankings of companies that are to be evaluated according to different macro and micro economic aspects, leading to ordinal covariates that correspond to various ratings, that entail different magnitudes of the probability of default. For each given covariate, we suggest to partition the statistical units in as many groups as the number of observed levels of the covariate. We then assume individual defaults to be homogeneous within each group and heterogeneous across groups. Our aim is to compare and, therefore select, the partition structures resulting from the consideration of different explanatory covariates. The metric we choose for variable comparison is the calculation of the posterior probability of each partition. The application of our proposal to a European credit risk database shows that it performs well, leading to a coherent and clear method for variable averaging of the estimated default probabilities.     

Time-Varying Systemic Risk: Evidence From a Dynamic Copula Model of CDS Spreads

This article proposes a new class of copula-based dynamic models for high-dimensional conditional distributions, facilitating the estimation of a wide variety of measures of systemic risk. Our proposed models draw on successful ideas from the literature on modeling high-dimensional covariance matrices and on recent work on models for general time-varying distributions. Our use of copula-based models enables the estimation of the joint model in stages, greatly reducing the computational burden. We use the proposed new models to study a collection of daily credit default swap (CDS) spreads on 100 U.S. firms over the period 2006 to 2012. We find that while the probability of distress for individual firms has greatly reduced since the financial crisis of 2008–2009, the joint probability of distress (a measure of systemic risk) is substantially higher now than in the precrisis period. Supplementary materials for this article are available online.     

基于高斯谱聚类的风险商户聚类分析

随着电子支付的普及,市场涌现出越来越多的第三方支付平台,而当前关于第三方支付平台商户风险方面的研究相对较少。故本文提出基于高斯谱聚类的风险商户聚类方法,首先使用高斯混合模型构建交易-交易群体的双模网络;其次借助网络中信息传递的思想构建“商户-交易群体网络”的双模网络;再次使用双模网络聚类方法中的谱聚类方法同时对网络中的两类节点聚类,对商户节点聚类的结果可区分出不同风险级别的商户,对交易群体节点聚类的结果可以进一步描述风险商户的交易特征;最后本文分别在模拟数据和某第方支付平台的实际数据中验证了模型的有效性。实验结果表明,本文提出的方法不仅可以准确地区分出不同风险级别的商户群体,而且能总结归纳风险商户的交易特征,为风险商户的监管提供参考。     

A discrete-time and finite-state Markov Chain model for association football matches

A birth process model proposed by Dixon and Robinson has been widely used in football spread betting market. However, multiple goals in a minute are permitted in the model, which does not conform to historical record. Moreover, it is difficult to calculate the outcome probability of the process accurately. The article presents a discrete-time and finite-state Markov chain model for real-time forecast of football matches and a recursive algorithm is derived to calculate the outcome probability accurately. The empirical study shows that the proposed model outperforms the models of Dixon and Robinson and Dixon and Coles.