Free Knot Splines with RJMCMC in Survival Data Analysis

The B-spline representation is a common tool to improve the fitting of smooth nonlinear functions, it offers a fitting as a piecewise polynomial. The regions that define the pieces are separated by a sequence of knots. The main difficulty in this type of modeling is the choice of the number and the locations of these knots. The Reversible Jump Markov Chain Monte Carlo (RJMCMC) algorithm provides a solution to simultaneously select these two parameters by considering the knots as free parameters. This algorithm belongs to the MCMC techniques that allow simulations from target distributions on spaces of varying dimension. The aim of the present investigation is to use this algorithm in the framework of the analysis of survival time, for the Cox model in particular. In fact, the relation between the hazard ratio function and the covariates being assumed to be log-linear, this assumption is too restrictive. Thus, we propose to use the RJMCMC algorithm to model the log hazard ratio function by a B-spline representation with an unknown number of knots at unknown locations. This method is illustrated with two real data sets: the Stanford heart transplant data and lung cancer survival data. Another application of the RJMCMC is selecting the significant covariates, and a simulation study is performed.     

Application of a non-homogeneous Markov chain with seasonal transition probabilities to ozone data

In this work, we assume that the sequence recording whether or not an ozone exceedance of an environmental threshold has occurred in a given day is ruled by a non-homogeneous Markov chain of order one. In order to account for the possible presence of cycles in the empirical transition probabilities, a parametric form incorporating seasonal components is considered. Results show that even though some covariates (namely, relative humidity and temperature) are not included explicitly in the model, their influence is captured in the behavior of the transition probabilities. Parameters are estimated using the Bayesian point of view via Markov chain Monte Carlo algorithms. The model is applied to ozone data obtained from the monitoring network of Mexico City, Mexico. An analysis of how the methodology could be used as an aid in the decision-making is also given.     

Financial Well-Being Index in México

Social Indicators Research - This paper aims to construct a financial well-being index and analyze its determinants among adults in Mexico. It employs multiple correspondence analysis, multiple...     

HOW DOES PRODUCT MARKET COMPETITION SHAPE INCENTIVE CONTRACTS?

This paper studies the effect of product market competition on the explicit compensation packages that firms offer to their CEOs, executives and workers. We use a large sample of both traded and nontraded UK firms and exploit a quasi‐natural experiment associated to an increase in competition. The sudden appreciation of the pound in 1996 implied different changes in competition for sectors with different degrees of openness. Our difference in differences estimates show that a higher level of product market competition increases the performance pay sensitivity of compensation schemes, in particular for executives. (JEL: J33, L20, G34)     

Simultaneous marginal survival estimators when doubly censored data is present

A doubly censoring scheme occurs when the lifetimes T being measured, from a well-known time origin, are exactly observed within a window [L, R] of observational time and are otherwise censored either from above (right-censored observations) or below (left-censored observations). Sample data consists on the pairs (U, δ) where U = min{R, max{T, L}} and δ indicates whether T is exactly observed (δ = 0), right-censored (δ = 1) or left-censored (δ = −1). We are interested in the estimation of the marginal behaviour of the three random variables T, L and R based on the observed pairs (U, δ). We propose new nonparametric simultaneous marginal estimators [^(S)]_T, [^(S)]_L{\hat S_{T}, \hat S_{L}} and [^(S)]_R{\hat S_{R}} for the survival functions of T, L and R, respectively, by means of an inverse-probability-of-censoring approach. The proposed estimators [^(S)]_T, [^(S)]_L{\hat S_{T}, \hat S_{L}} and [^(S)]_R{\hat S_{R}} are not computationally intensive, generalize the empirical survival estimator and reduce to the Kaplan-Meier estimator in the absence of left-censored data. Furthermore, [^(S)]_T{\hat S_{T}} is equivalent to a self-consistent estimator, is uniformly strongly consistent and asymptotically normal. The method is illustrated with data from a cohort of drug users recruited in a detoxification program in Badalona (Spain). For these data we estimate the survival function for the elapsed time from starting IV-drugs to AIDS diagnosis, as well as the potential follow-up time. A simulation study is discussed to assess the performance of the three survival estimators for moderate sample sizes and different censoring levels.     

Asymptotic Properties for a Class of Partially Identified Models

We propose inference procedures for partially identified population features for which the population identification region can be written as a transformation of the Aumann expectation of a properly defined set valued random variable (SVRV). An SVRV is a mapping that associates a set (rather than a real number) with each element of the sample space. Examples of population features in this class include interval‐identified scalar parameters, best linear predictors with interval outcome data, and parameters of semiparametric binary models with interval regressor data. We extend the analogy principle to SVRVs and show that the sample analog estimator of the population identification region is given by a transformation of a Minkowski average of SVRVs. Using the results of the mathematics literature on SVRVs, we show that this estimator converges in probability to the population identification region with respect to the Hausdorff distance. We then show that the Hausdorff distance and the directed Hausdorff distance between the population identification region and the estimator, when properly normalized by , converge in distribution to functions of a Gaussian process whose covariance kernel depends on parameters of the population identification region. We provide consistent bootstrap procedures to approximate these limiting distributions. Using similar arguments as those applied for vector valued random variables, we develop a methodology to test assumptions about the true identification region and its subsets. We show that these results can be used to construct a confidence collection and a directed confidence collection. Those are (respectively) collection of sets that, when specified as a null hypothesis for the true value (a subset of values) of the population identification region, cannot be rejected by our tests.     

A generalized Fleming and Harrington's class of tests for interval‐censored data

The class $G^{\rho,\lambda }$ of weighted log‐rank tests proposed by Fleming & Harrington [Fleming & Harrington (1991) Counting Processes and Survival Analysis, Wiley, New York] has been widely used in survival analysis and is nowadays, unquestionably, the established method to compare, nonparametrically, k different survival functions based on right‐censored survival data. This paper extends the $G^{\rho,\lambda }$ class to interval‐censored data. First we introduce a new general class of rank based tests, then we show the analogy to the above proposal of Fleming & Harrington. The asymptotic behaviour of the proposed tests is derived using an observed Fisher information approach and a permutation approach. Aiming to make this family of tests interpretable and useful for practitioners, we explain how to interpret different choices of weights and we apply it to data from a cohort of intravenous drug users at risk for HIV infection. The Canadian Journal of Statistics 40: 501–516; 2012 © 2012 Statistical Society of Canada     

Freedom,Identity, and the Social History of Empire in Atlantic Cuba, 1795–1817

This essay examines a maritime case involving a Charleston-based French privateer in 1795 as well as members, subjects, and citizens of multiple societies from Western Africa, Cuba, Spain, and Great Britain in order to reveal the fluid legal foundations of trans-Atlantic regimes with respect to slavery and freedom. The French maritime seizure of a British slaving vessel in the Caribbean and the subsequent transfer of enslaved African peoples from British to French to Spanish authority in Cuba set the stage for nearly 40 years of litigation in both hemispheres. The case led to legal disputes in Philadelphia, Madrid, and Havana as members of this Atlantic society sought to confront the most prominent political and social issues of this era. While officials in Philadelphia sought to avoid entanglement in the imperial wars between France and Great Britain, officials in Madrid and Havana competed with each other in order to define the legal identities of a group of Carabalí-Oru peoples onboard the British vessel. For their part, this group of Carabalí-Oru peoples exploited the ambiguities created by a system of Atlantic jurisprudence while relying on their memory of the Middle Passage in order to challenge their unjust enslavement in Cuba.     

Italian Adult Gambling Behavior: At Risk and Problem Gambler Profiles

The present study examined adult gambling behaviours from a local perspective in order to assess the adult at risk and problem gambler’s profile stratified by genre and by different forms of game. 4773 Italian adults from 18 to 94 years old were administered a survey to assess socio-cultural information related to gambling behaviour and the SOGS to evaluate gambling behaviour severity. Logistic regression evidenced that both at risk and problem gamblers are associated with male gender, players that use to play to more than one game, gambling with strategy-based games. People with a gambler father or both parents who used to gamble were significantly more associated with problem gambling behaviour than participants with non-gambler parents. These results present adult profiles of at risk and problem providing a more clear understanding about the relationships between gambling behavior severity and type of gambling.     

Interval censoring: Model characterizations for the validity of the simplified likelihood

In survival data analysis, the interval censoring problem has generally been treated via likelihood methods. Because this likelihood is complex, it is often assumed that the censoring mechanisms do not affect the mortality process. The authors specify conditions that ensure the validity of such a simplified likelihood. They prove the equivalence between different characterizations of noninformative censoring and define a constant‐sum condition analogous to the one derived in the context of right censoring. They also prove that when the noninformative or constant‐sum condition holds, the simplified likelihood can be used to obtain the nonparametric maximum likelihood estimator of the death time distribution function.