Search

The formation history of binary black hole systems is imprinted on the distribution of their masses, spins, and eccentricity. While much has been learned studying these parameters in turn, recent studies have explored the joint distribution of binary black hole parameters in two or more dimensions. Most notably, it has recently been argued that binary black hole mass ratio and effective inspiral spin χeff are anticorrelated. We point out a previously overlooked subtlety in such 2D population studies: in order to conduct a controlled test for correlation, one ought to fix the two marginal distributions – lest the purported correlation be driven by improved fit in just one dimension. We address this subtlety using a tool from applied statistics: the copula density function. We use the previous work correlating mass ratio and χeff as a case study to demonstrate the power of copulas in gravitational-wave astronomy while scrutinizing their astrophysical inferences. Our findings, however, affirm their conclusions that binary black holes with unequal component masses exhibit larger χeff (98.7 per cent credibility). We conclude by discussing potential astrophysical implications of these findings as well as prospects for future studies using copulas. [ABSTRACT FROM AUTHOR]

In dose–response analysis, it is a challenge to choose appropriate linear or curvilinear shapes when considering multiple, differently scaled endpoints. It has been proposed to fit several marginal regression models that try sets of different transformations of the dose levels as explanatory variables for each endpoint. However, the multiple testing problem underlying this approach, involving correlated parameter estimates for the dose effect between and within endpoints, could only be adjusted heuristically. An asymptotic correction for multiple testing can be derived from the score functions of the marginal regression models. Based on a multivariate t‐distribution, the correction provides a one‐step adjustment of p‐values that accounts for the correlation between estimates from different marginal models. The advantages of the proposed methodology are demonstrated through three example datasets, involving generalized linear models with differently scaled endpoints, differing covariates, and a mixed effect model and through simulation results. The methodology is implemented in an R package. [ABSTRACT FROM AUTHOR]

Suppose HT ht . One may try to remove such ambiguity by restricting the minimum support of HT ht to be greater than I m i , but this essentially excludes the possibility that the primary event censors the mediating event. Therefore, the usual notion of having well-defined interventions for I A i and I M i for defining HT ht does not go through because HT ht may happen before I m i , and in that case, the role of I m i is ambiguous. Under a composition assumption, HT ht . [Extracted from the article]

We propose a new highly flexible and tractable Bayesian approach to undertake variable selection in non‐Gaussian regression models. It uses a copula decomposition for the joint distribution of observations on the dependent variable. This allows the marginal distribution of the dependent variable to be calibrated accurately using a nonparametric or other estimator. The family of copulas employed are "implicit copulas" that are constructed from existing hierarchical Bayesian models widely used for variable selection, and we establish some of their properties. Even though the copulas are high dimensional, they can be estimated efficiently and quickly using Markov chain Monte Carlo. A simulation study shows that when the responses are non‐Gaussian, the approach selects variables more accurately than contemporary benchmarks. A real data example in the Web Appendix illustrates that accounting for even mild deviations from normality can lead to a substantial increase in accuracy. To illustrate the full potential of our approach, we extend it to spatial variable selection for fMRI. Using real data, we show our method allows for voxel‐specific marginal calibration of the magnetic resonance signal at over 6000 voxels, leading to an increase in the quality of the activation maps. [ABSTRACT FROM AUTHOR]

Dynamic prediction uses patient information collected during follow-up to produce individualized survival predictions at given time points beyond treatment or diagnosis. This allows clinicians to obtain updated predictions of a patient's prognosis that can be used in making personalized treatment decisions. Two commonly used approaches for dynamic prediction are landmarking and joint modeling. Landmarking does not constitute a comprehensive probability model, and joint modeling often requires strong distributional assumptions and computationally intensive methods for estimation. We introduce an alternative approximate approach for dynamic prediction that aims to overcome the limitations of both methods while achieving good predictive performance. We separately specify the marker and failure time distributions conditional on surviving up to a prediction time of interest and use standard variable selection and goodness-of-fit techniques to identify the best-fitting models. Taking advantage of its analytic tractability and easy two-stage estimation, we use a Gaussian copula to link the marginal distributions smoothly at each prediction time with an association function. With simulation studies, we examine the proposed method's performance. We illustrate its use for dynamic prediction in an application to predicting death for heart valve transplant patients using longitudinal left ventricular mass index information. [ABSTRACT FROM AUTHOR]

This research is motivated by discovering and underpinning genetic causes for the progression of a bilateral eye disease, age-related macular degeneration (AMD), of which the primary outcomes, progression times to late-AMD, are bivariate and interval-censored due to intermittent assessment times. We propose a novel class of copula-based semiparametric transformation models for bivariate data under general interval censoring, which includes the case 1 interval censoring (current status data) and case 2 interval censoring. Specifically, the joint likelihood is modeled through a two-parameter Archimedean copula, which can flexibly characterize the dependence between the two margins in both tails. The marginal distributions are modeled through semiparametric transformation models using sieves, with the proportional hazards or odds model being a special case. We develop a computationally efficient sieve maximum likelihood estimation procedure for the unknown parameters, together with a generalized score test for the regression parameter(s). For the proposed sieve estimators of finite-dimensional parameters, we establish their asymptotic normality and efficiency. Extensive simulations are conducted to evaluate the performance of the proposed method in finite samples. Finally, we apply our method to a genome-wide analysis of AMD progression using the Age-Related Eye Disease Study data, to successfully identify novel risk variants associated with the disease progression. We also produce predicted joint and conditional progression-free probabilities, for patients with different genetic characteristics. [ABSTRACT FROM AUTHOR]

The need for a method to construct multidimensional distribution function is increasing recently, in the era of huge multiwavelength surveys. We have proposed a systematic method to build a bivariate luminosity or mass function of galaxies by using a copula. It allows us to construct a distribution function when only its marginal distributions are known, and we have to estimate the dependence structure from data. A typical example is the situation that we have univariate luminosity functions at some wavelengths for a survey, but the joint distribution is unknown. Main limitation of the copula method is that it is not easy to extend a joint function to higher dimensions (d > 2), except some special cases like multidimensional Gaussian. Even if we find such a multivariate analytic function in some fortunate case, it would often be inflexible and impractical. In this work, we show a systematic method to extend the copula method to unlimitedly higher dimensions by a vine copula. This is based on the pair-copula decomposition of a general multivariate distribution. We show how the vine copula construction is flexible and extendable. We also present an example of the construction of a stellar mass–atomic gas–molecular gas three-dimensional mass function. We demonstrate the maximum likelihood estimation of the best functional form for this function, as well as a proper model selection via vine copula. [ABSTRACT FROM AUTHOR]

Several gene‐based association tests for time‐to‐event traits have been proposed recently to detect whether a gene region (containing multiple variants), as a set, is associated with the survival outcome. However, for bivariate survival outcomes, to the best of our knowledge, there is no statistical method that can be directly applied for gene‐based association analysis. Motivated by a genetic study to discover the gene regions associated with the progression of a bilateral eye disease, age‐related macular degeneration (AMD), we implement a novel functional regression (FR) method under the copula framework. Specifically, the effects of variants within a gene region are modeled through a functional linear model, which then contributes to the marginal survival functions within the copula. Generalized score test statistics are derived to test for the association between bivariate survival traits and the genetic region. Extensive simulation studies are conducted to evaluate the type I error control and power performance of the proposed approach, with comparisons to several existing methods for a single survival trait, as well as the marginal Cox FR model using the robust sandwich estimator for bivariate survival traits. Finally, we apply our method to a large AMD study, the Age‐related Eye Disease Study, and to identify the gene regions that are associated with AMD progression. [ABSTRACT FROM AUTHOR]

Summary: The outcome of high‐throughput biological experiments is affected by many operational factors in the experimental and data‐analytical procedures. Understanding how these factors affect the reproducibility of the outcome is critical for establishing workflows that produce replicable discoveries. In this article, we propose a regression framework, based on a novel cumulative link model, to assess the covariate effects of operational factors on the reproducibility of findings from high‐throughput experiments. In contrast to existing graphical approaches, our method allows one to succinctly characterize the simultaneous and independent effects of covariates on reproducibility and to compare reproducibility while controlling for potential confounding variables. We also establish a connection between our model and certain Archimedean copula models. This connection not only offers our regression framework an interpretation in copula models, but also provides guidance on choosing the functional forms of the regression. Furthermore, it also opens a new way to interpret and utilize these copulas in the context of reproducibility. Using simulations, we show that our method produces calibrated type I error and is more powerful in detecting difference in reproducibility than existing measures of agreement. We illustrate the usefulness of our method using a ChIP‐seq study and a microarray study. [ABSTRACT FROM AUTHOR]

Summary: It is desirable to adjust Spearman's rank correlation for covariates, yet existing approaches have limitations. For example, the traditionally defined partial Spearman's correlation does not have a sensible population parameter, and the conditional Spearman's correlation defined with copulas cannot be easily generalized to discrete variables. We define population parameters for both partial and conditional Spearman's correlation through concordance–discordance probabilities. The definitions are natural extensions of Spearman's rank correlation in the presence of covariates and are general for any orderable random variables. We show that they can be neatly expressed using probability‐scale residuals (PSRs). This connection allows us to derive simple estimators. Our partial estimator for Spearman's correlation between X and Y adjusted for Z is the correlation of PSRs from models of X on Z and of Y on Z, which is analogous to the partial Pearson's correlation derived as the correlation of observed‐minus‐expected residuals. Our conditional estimator is the conditional correlation of PSRs. We describe estimation and inference, and highlight the use of semiparametric cumulative probability models, which allow preservation of the rank‐based nature of Spearman's correlation. We conduct simulations to evaluate the performance of our estimators and compare them with other popular measures of association, demonstrating their robustness and efficiency. We illustrate our method in two applications, a biomarker study and a large survey. [ABSTRACT FROM AUTHOR]

This article proposes a modeling strategy to infer the impact of a covariate on the dependence structure of right-censored clustered event time data. The joint survival function of the event times is modeled using a conditional copula whose parameter depends on a cluster-level covariate in a functional way. We use a local likelihood approach to estimate the form of the copula parameter and outline a generalized likelihood ratio-type test strategy to formally test its constancy. A bootstrap procedure is employed to obtain an approximate $p$-value for the test. The performance of the proposed estimation and testing methods is evaluated in simulations under different rates of right-censoring and for various parametric copula families, considering both parametrically and nonparametrically estimated margins. We apply the methods to data from the Diabetic Retinopathy Study to assess the impact of age at diabetes onset on the time to loss of visual acuity. [ABSTRACT FROM AUTHOR]

A common approach towards modelling bivariate scatterplots is decomposition into Gaussian components, i.e. Gaussian mixture modelling. This implicitly assumes linear relationships between the variables within each of the components in the mixture. An alternative, namely dependence modelling by mixtures of copulas, is advocated in this paper. This approach allows separate modelling of the univariate marginal distributions and the dependence which can possibly be non-linear and/or asymmetric. It also accommodates the use of a variety of parametric families for modelling each component and for each variable. The variety of dependence structures can be extended by introducing rotated versions of the copulas. Gaussian mixture modelling on the one hand, and separate modelling of univariate marginal distributions and dependence on the other hand, are illustrated by application to pulsar period – period-derivative observations. Parameter estimation for mixtures of copulas is performed using the method of maximum likelihood and selected copula models are subjected to non-parametric goodness-of-fit testing. [ABSTRACT FROM AUTHOR]

In many biomedical studies that involve correlated data, an outcome is often repeatedly measured for each individual subject along with the number of these measurements, which is also treated as an observed outcome. This type of data has been referred as multivariate random length data by Barnhart and Sampson (1995). A common approach to handling such type of data is to jointly model the multiple measurements and the random length. In previous literature, a key assumption is the multivariate normality for the multiple measurements. Motivated by a reproductive study, we propose a new copula-based joint model which relaxes the normality assumption. Specifically, we adopt the Clayton-Oakes model for multiple measurements with flexible marginal distributions specified as semi-parametric transformation models. The random length is modeled via a generalized linear model. We develop an approximate EM algorithm to derive parameter estimators and standard errors of the estimators are obtained through bootstrapping procedures and the finite-sample performance of the proposed method is investigated using simulation studies. We apply our method to the Mount Sinai Study of Women Office Workers (MSSWOW), where women were prospectively followed for 1 year for studying fertility. [ABSTRACT FROM AUTHOR]

In this article, we propose an association model to estimate the penetrance (risk) of successive cancers in the presence of competing risks. The association between the successive events is modeled via a copula and a proportional hazards model is specified for each competing event. This work is motivated by the analysis of successive cancers for people with Lynch Syndrome in the presence of competing risks. The proposed inference procedure is adapted to handle missing genetic covariates and selection bias, induced by the data collection protocol of the data at hand. The performance of the proposed estimation procedure is evaluated by simulations and its use is illustrated with data from the Colon Cancer Family Registry (Colon CFR). [ABSTRACT FROM AUTHOR]

The heritability of chronic diseases can be effectively studied by examining the nature and extent of within-family associations in disease onset times. Families are typically accrued through a biased sampling scheme in which affected individuals are identified and sampled along with their relatives who may provide right-censored or current status data on their disease onset times. We develop likelihood and composite likelihood methods for modeling the within-family association in these times through copula models in which dependencies are characterized by Kendall's [Formula: see text] Auxiliary data from independent individuals are exploited by augmentating composite likelihoods to increase precision of marginal parameter estimates and consequently increase efficiency in dependence parameter estimation. An application to a motivating family study in psoriatic arthritis illustrates the method and provides some evidence of excessive paternal transmission of risk. [ABSTRACT FROM AUTHOR]

Current status data occur in contexts including demographic studies and tumorigenicity experiments. In such cases, each subject is observed only once and the failure time of interest is either left- or right-censored (Kalbfleisch & Prentice, 2002). Many methods have been developed for the analysis of such data (Huang, 1996; Sun, 2006), most of which assume that the failure time and the observation time are independent completely or given covariates. In this paper, we present a sieve maximum likelihood approach for current status data when independence does not hold. A copula model and monotone I-splines are used and the asymptotic properties of the resulting estimators are established. In particular, the estimated regression parameters are shown to be semiparametrically efficient. An illustrative example is provided. [ABSTRACT FROM AUTHOR]

In capture-recapture studies, the use of individual covariates has been recommended to get stable population estimates. However, some residual heterogeneity might still exist and ignoring such heterogeneity could lead to underestimating the population size ( N). In this work, we explore two new models with capture probabilities depending on both covariates and unobserved random effects, to estimate the size of a population. Inference techniques including Horvitz-Thompson estimate and confidence intervals for the population size, are derived. The selection of a particular model is carried out using the Akaike information criterion (AIC). First, we extend the random effect model of Darroch et al. (1993, Journal of American Statistical Association 88, 1137-1148) to handle unit level covariates and discuss its limitations. The second approach is a generalization of the traditional zero-truncated binomial model that includes a random effect to account for an unobserved heterogeneity. This approach provides useful tools for inference about N, since key quantities such as moments, likelihood functions and estimates of N and their standard errors have closed form expressions. Several models for the unobserved heterogeneity are available and the marginal capture probability is expressed using the Logit and the complementary Log-Log link functions. The sensitivity of the inference to the specification of a model is also investigated through simulations. A numerical example is presented. We compare the performance of the proposed estimator with that obtained under model [ABSTRACT FROM AUTHOR]

This article develops a Bayesian semiparametric approach to the extended hazard model, with generalization to high-dimensional spatially grouped data. County-level spatial correlation is accommodated marginally through the normal transformation model of Li and Lin (2006, Journal of the American Statistical Association 101, 591-603), using a correlation structure implied by an intrinsic conditionally autoregressive prior. Efficient Markov chain Monte Carlo algorithms are developed, especially applicable to fitting very large, highly censored areal survival data sets. Per-variable tests for proportional hazards, accelerated failure time, and accelerated hazards are efficiently carried out with and without spatial correlation through Bayes factors. The resulting reduced, interpretable spatial models can fit significantly better than a standard additive Cox model with spatial frailties. [ABSTRACT FROM AUTHOR]

A copula-based model is described which enables joint analysis of multiple progressive multistate processes. Unlike intensity-based or frailty-based approaches to joint modeling, the copula formulation proposed herein ensures that a wide range of marginal multistate processes can be specified and the joint model will retain these marginal features. The copula formulation also facilitates a variety of approaches to estimation and inference including composite likelihood and two-stage estimation procedures. We consider processes with Markov margins in detail, which are often suitable when chronic diseases are progressive in nature. We give special attention to the setting in which individuals are examined intermittently and transition times are consequently interval-censored. Simulation studies give empirical insight into the different methods of analysis and an application involving progression in joint damage in psoriatic arthritis provides further illustration. [ABSTRACT FROM AUTHOR]

Random sums play an important role in insurance claim theory and risk models. This paper aims to explore the numerical characteristics of dependent random sums using copulas. We obtain theoretical formulae for the expectation and variance of random sums in the case of a general copula. The expectation is an extension of the Wald equality and can be used to estimate prospective value, while variance provides a theoretical guarantee for measuring the risk premium. As corollaries, the corresponding results for the compound Poisson process and random sums using Farlie-Gumbel-Morgenstern copulas have been obtained. We also conduct stochastic simulation to verify the accuracy of our theoretical results. [ABSTRACT FROM AUTHOR]

We discuss sample size determination for clinical trials evaluating the joint effects of an intervention on two potentially correlated co-primary time-to-event endpoints. For illustration, we consider the most common case, a comparison of two randomized groups, and use typical copula families to model the bivariate endpoints. A correlation structure of the bivariate logrank statistic is specified to account for the correlation among the endpoints, although the between-group comparison is performed using the univariate logrank statistic. We propose methods to calculate the required sample size to compare the two groups and evaluate the performance of the methods and the behavior of required sample sizes via simulation. [ABSTRACT FROM AUTHOR]

We propose a novel class of models for functional data exhibiting skewness or other shape characteristics that vary with spatial or temporal location. We use copulas so that the marginal distributions and the dependence structure can be modeled independently. Dependence is modeled with a Gaussian or t-copula, so that there is an underlying latent Gaussian process. We model the marginal distributions using the skew t family. The mean, variance, and shape parameters are modeled nonparametrically as functions of location. A computationally tractable inferential framework for estimating heterogeneous asymmetric or heavy-tailed marginal distributions is introduced. This framework provides a new set of tools for increasingly complex data collected in medical and public health studies. Our methods were motivated by and are illustrated with a state-of-the-art study of neuronal tracts in multiple sclerosis patients and healthy controls. Using the tools we have developed, we were able to find those locations along the tract most affected by the disease. However, our methods are general and highly relevant to many functional data sets. In addition to the application to one-dimensional tract profiles illustrated here, higher-dimensional extensions of the methodology could have direct applications to other biological data including functional and structural magnetic resonance imaging (MRI). [ABSTRACT FROM AUTHOR]

There are copula-based statistical models in the literature for regression with dependent data such as clustered and longitudinal overdispersed counts, for which parameter estimation and inference are straightforward. For situations where the main interest is in the regression and other univariate parameters and not the dependence, we propose a “weighted scores method”, which is based on weighting score functions of the univariate margins. The weight matrices are obtained initially fitting a discretized multivariate normal distribution, which admits a wide range of dependence. The general methodology is applied to negative binomial regression models. Asymptotic and small-sample efficiency calculations show that our method is robust and nearly as efficient as maximum likelihood for fully specified copula models. An illustrative example is given to show the use of our weighted scores method to analyze utilization of health care based on family characteristics. [ABSTRACT FROM PUBLISHER]

Summary The goal of mediation analysis is to assess direct and indirect effects of a treatment or exposure on an outcome. More generally, we may be interested in the context of a causal model as characterized by a directed acyclic graph (DAG), where mediation via a specific path from exposure to outcome may involve an arbitrary number of links (or 'stages'). Methods for estimating mediation (or pathway) effects are available for a continuous outcome and a continuous mediator related via a linear model, while for a categorical outcome or categorical mediator, methods are usually limited to two-stage mediation. We present a method applicable to multiple stages of mediation and mixed variable types using generalized linear models. We define pathway effects using a potential outcomes framework and present a general formula that provides the effect of exposure through any specified pathway. Some pathway effects are nonidentifiable and their estimation requires an assumption regarding the correlation between counterfactuals. We provide a sensitivity analysis to assess the impact of this assumption. Confidence intervals for pathway effect estimates are obtained via a bootstrap method. The method is applied to a cohort study of dental caries in very low birth weight adolescents. A simulation study demonstrates low bias of pathway effect estimators and close-to-nominal coverage rates of confidence intervals. We also find low sensitivity to the counterfactual correlation in most scenarios. [ABSTRACT FROM AUTHOR]

The study of dependence between random variables is a mainstay in statistics. In many cases, the strength of dependence between two or more random variables varies according to the values of a measured covariate. We propose inference for this type of variation using a conditional copula model where the copula function belongs to a parametric copula family and the copula parameter varies with the covariate. In order to estimate the functional relationship between the copula parameter and the covariate, we propose a nonparametric approach based on local likelihood. Of importance is also the choice of the copula family that best represents a given set of data. The proposed framework naturally leads to a novel copula selection method based on cross-validated prediction errors. We derive the asymptotic bias and variance of the resulting local polynomial estimator, and outline how to construct pointwise confidence intervals. The finite-sample performance of our method is investigated using simulation studies and is illustrated using a subset of the Matched Multiple Birth data. [ABSTRACT FROM AUTHOR]

In life history studies, interest often lies in the analysis of the interevent, or gap times and the association between event times. Gap time analyses are challenging however, even when the length of follow-up is determined independently of the event process, because associations between gap times induce dependent censoring for second and subsequent gap times. This article discusses nonparametric estimation of the association between consecutive gap times based on Kendall's τ in the presence of this type of dependent censoring. A nonparametric estimator that uses inverse probability of censoring weights is provided. Estimates of conditional gap time distributions can be obtained following specification of a particular copula function. Simulation studies show the estimator performs well and compares favorably with an alternative estimator. Generalizations to a piecewise constant Clayton copula are given. Several simulation studies and illustrations with real data sets are also provided. [ABSTRACT FROM AUTHOR]

This paper discusses copula model selection procedures and goodness-of-fit tests under censoring. The proposed methodology is based on a comparison of nonparametric and model-based estimators of the probability integral transformation, K. New weighted estimators for K are introduced. The resulting tests are compared to an existing approach by simulation and illustrated with an example involving bleeding changes in a woman’s reproductive history. [ABSTRACT FROM PUBLISHER]

There has been a recent emphasis on the identification of biomarkers and other biologic measures that may be potentially used as surrogate endpoints in clinical trials. We focus on the setting of data from a single clinical trial. In this article, we consider a framework in which the surrogate must occur before the true endpoint. This suggests viewing the surrogate and true endpoints as semicompeting risks data; this approach is new to the literature on surrogate endpoints and leads to an asymmetrical treatment of the surrogate and true endpoints. However, such a data structure also conceptually complicates many of the previously considered measures of surrogacy in the literature. We propose novel estimation and inferential procedures for the relative effect and adjusted association quantities proposed by Buyse and Molenberghs (1998, Biometrics 54, 1014–1029). The proposed methodology is illustrated with application to simulated data, as well as to data from a leukemia study. [ABSTRACT FROM AUTHOR]

Many longitudinal studies generate both the time to some event of interest and repeated measures data. This article is motivated by a study on patients with a renal allograft, in which interest lies in the association between longitudinal proteinuria (a dichotomous variable) measurements and the time to renal graft failure. An interesting feature of the sample at hand is that nearly half of the patients were never tested positive for proteinuria (≥1g/day) during follow-up, which introduces a degenerate part in the random-effects density for the longitudinal process. In this article we propose a two-part shared parameter model framework that effectively takes this feature into account, and we investigate sensitivity to the various dependence structures used to describe the association between the longitudinal measurements of proteinuria and the time to renal graft failure. [ABSTRACT FROM AUTHOR]

In many follow-up studies, patients are subject to concurrent events. In this article, we consider semicompeting risks data as defined by Fine, Jiang, and Chappell (2001, Biometrika 88, 907–919) where one event is censored by the other but not vice versa. The proposed model involves marginal survival functions for the two events and a parametric family of copulas for their dependency. This article suggests a general method for estimating the dependence parameter when the dependency is modeled with an Archimedean copula. It uses the copula-graphic estimator of Zheng and Klein (1995, Biometrika 82, 127–138) for estimating the survival function of the nonterminal event, subject to dependent censoring. Asymptotic properties of these estimators are derived. Simulations show that the new methods work well with finite samples. The copula-graphic estimator is shown to be more accurate than the estimator proposed by Fine et al. (2001) ; its performances are similar to those of the self-consistent estimator of Jiang, Fine, Kosorok, and Chappell (2005, Scandinavian Journal of Statistics 33, 1–20). The analysis of a data set, emphasizing the estimation of characteristics of the observable region, is presented as an illustration. [ABSTRACT FROM AUTHOR]

The product-limit estimator calculated from data subject to random left-truncation relies on the testable assumption of quasi-independence between the failure time and the truncation time. In this paper, we propose a model for a truncated sample of pairs (Xi,Yi) satisfying Yi > Xi. A possible dependency between the truncation time and the variable of interest is modelled with a parametric family of copulas. The model also features a distribution function FX(.) and a survival distribution SY(.) associated with the marginal behaviours of X and Y in the observable region Y > X. Semiparametric estimators for these two functions are proposed; they do not make any parametric assumption about either FX(.) or SY(.). We derive an estimator for the copula parameter α based on the conditional Kendall's tau. We generalise the copula-graphic estimators of Zheng & Klein (1995) to truncated variables. The asymptotic distributions of all these estimators are then investigated. The methods are illustrated with a real dataset on HIV infection by transfusion and by simulations. [ABSTRACT FROM AUTHOR]

A Gaussian copula regression model gives a tractable way of handling a multivariate regression when some of the marginal distributions are non-Gaussian. Our paper presents a general Bayesian approach for estimating a Gaussian copula model that can handle any combination of discrete and continuous marginals, and generalises Gaussian graphical models to the Gaussian copula framework. Posterior inference is carried out using a novel and efficient simulation method. The methods in the paper are applied to simulated and real data. [ABSTRACT FROM AUTHOR]

The objective of this paper is to estimate a bivariate density nonparametrically from a dataset from the joint distribution and datasets from one or both marginal distributions. We develop a copula-based solution, which has potential benefits even when the marginal datasets are empty. For example, if the copula density is sufficiently smooth in the region where we wish to estimate it, the joint density can be estimated with a high degree of accuracy. Similar improvements in performance are available if the marginals are close to being independent. We use wavelet estimators to approximate the copula density, which in cases of statistical interest can be unbounded along boundaries. Our techniques are also useful for solving recently-considered related problems, for example where the marginal distributions are determined by parametric models. The methodology is also readily extended to more general multivariate settings. [ABSTRACT FROM AUTHOR]

This article presents methodology for multivariate proportional hazards (PH)regression models. The methods employ flexible piecewise constant or spline specifications for baseline hazard functions in either marginal or conditional PH models, along with assumptions about the association among lifetimes. Because the models are parametric, ordinary maximum likelihood can be applied; it is able to deal easily with such data features as interval censoring or sequentially observed lifetimes, unlike existing semiparametric methods. A bivariate Clayton model (1978, Biometrika 65,141-151)is used to illustrate the approach taken. Because a parametric assumption about association is made, efficiency and robustness comparisons are made between estimation based on the bivariate Clayton model and "working independence" methods that specify only marginal distributions for each lifetime variable. [ABSTRACT FROM AUTHOR]

We propose a new archimedean copula model for bivariate survival data that is motivated by Dabrowska's (1988) measure of association. The model can represent negatively correlated or moderately positively correlated data but not highly positively correlated data. Local and global measures of association are calculated. A generalisation is presented. [ABSTRACT FROM AUTHOR]

Assuming that the two failure times of interest with bivariate current status data follow a bivariate copula model, we propose a two‐stage estimation procedure to estimate the association parameter which is related to Kendall's tau. Asymptotic properties of the proposed semiparametric estimator show that, although the first‐stage marginal estimators have a convergence rate of only n1/3, the resulting parameter estimator still converges to a normal random variable with the usual n1/2 rate. The variance of the proposed estimator can be consistently estimated. Simulation results are presented, and a community‐based study of cardiovascular diseases in Taiwan provides an illustrative example. [ABSTRACT FROM PUBLISHER]

We consider estimation of the cross ratio in Clayton's (1978) copula in which covariates are incorporated into the marginal distributions via semiparametric accelerated life regression models. Generalisations of Oakes' (1982, 1986) concordance estimating equations yield a closed form for the association parameter under right censoring. Joint inferences about covariate effects in the marginal models and the cross ratio are obtained with U-statistics, martingales and a resampling technique for nonsmooth estimating equations. Simulating under an exponential model, we find that our procedure may provide a more precise estimate of association than other methods for proportional hazards models. A goodness-of-fit test for the assumed copula is presented and used in an analysis of a diabetic retinopathy dataset. [ABSTRACT FROM PUBLISHER]

We consider a bivariate regression model in which the error structure is arbitrary, except for equality of its marginal distributions. Two simple but inefficient approaches to analysis are (i) to work with the differences between the two components of each response and (ii) to derive working independence estimators as if the two component responses were independent. The former strategy does well if the explanatory variables are balanced over pairs and the correlation between the two components of the response is high; the latter strategy does well when the response correlation is low or the covariates are highly unbalanced. We show that combining the two estimators can give a procedure with high efficiency in all situations. We show that an adjusted marginal score can be calculated that is uncorrelated with the within-pair score, so that separate intrapair and interpair estimators can be calculated. For bivariate normal regression models our approach is identical to the usual procedure for combining intra- and interblock information. In general the approaches are not the same. Detailed calculations and a numerical example are given for a bivariate extreme value model due to Gumbel that is also important in survival analysis. Extensions from bivariate to general multivariate responses are indicated. [ABSTRACT FROM PUBLISHER]

SUMMARY Using a model in which treatment affects only the failure type of primary clinical interest, the use of a cause-specific analysis is examined when the latent failure times are derived from a frailty model. The relationship between the cause-specific hazard ratio and the marginal hazard ratio is expressed in terms of a local dependence measure of the latent failure times. [ABSTRACT FROM AUTHOR]

Focuses on bivariate or multivariate survival modeling in presence of cure fraction. Methodology in the very general setting of copula model; Analysis of clustered data when one wants to fit and compare different parametric correlation models; Application of methodology using breast cancer data from Swedish twin registry data of twins.

Frailty and copula models specify a parametric dependence structure for multivariate failure-time data. Estimation of some joint quantities can be highly sensitive to the assumed parametric form, and hence model fit is an important issue. This paper lays out a general diagnostic framework for evaluating and selecting frailty and copula models. The approach is based on the cumulative sum of residuals that are calculated in bivariate time. The residuals reflect the difference between the observed and expected bivariate association structures. The proposed model-checking process is interpretable with a limiting distribution which can be approximated using the bootstrap. Simulations and a data example illustrate the practical application of the method. [ABSTRACT FROM AUTHOR]