A generalized response surface model with varying relative potency for assessing drug interaction

When multiple drugs are administered simultaneously, investigators are often interested in assessing whether the drug combinations are synergistic, additive, or antagonistic. Based on the Loewe additivity reference model, many existing response surface models require constant relative potency and some of them use a single parameter to capture synergy, additivity, or antagonism. However, the assumption of constant relative potency is too restrictive, and these models using a single parameter to capture drug interaction are inadequate to describe the phenomenon when synergy, additivity, and antagonism are interspersed in different regions of drug combinations. We propose a generalized response surface model with a function of doses instead of one single parameter to identify and quantify departure from additivity. The proposed model can incorporate varying relative potencies among multiple drugs as well. Examples and simulations are given to demonstrate that the proposed model is effective in capturing different patterns of drug interaction.     

A semiparametric approach for the nonparametric transformation survival model with multiple covariates

The nonparametric transformation model makes no parametric assumptions on the forms of the transformation function and the error distribution. This model is appealing in its flexibility for modeling censored survival data. Current approaches for estimation of the regression parameters involve maximizing discontinuous objective functions, which are numerically infeasible to implement with multiple covariates. Based on the partial rank (PR) estimator (Khan and Tamer, 2004), we propose a smoothed PR estimator which maximizes a smooth approximation of the PR objective function. The estimator is shown to be asymptotically equivalent to the PR estimator but is much easier to compute when there are multiple covariates. We further propose using the weighted bootstrap, which is more stable than the usual sandwich technique with smoothing parameters, for estimating the standard error. The estimator is evaluated via simulation studies and illustrated with the Veterans Administration lung cancer data set.     

Lipid model membranes for drug interaction study

The present work shows a structural study on the process of incorporation of a hydrophobic drug, Ellipticine (ELPT), into lipid model membranes for drug targeting purpose. The ELPT is an alkaloid that showed an anti-proliferation activity against several types of tumor cells and against the HIV1 virus. We used the zwitterionic lipid dipalmitoyl phosphatidylcholine (DPPC) and four different anionic lipids: cardiolipin (CL), dipalmitoyl phosphatidic acid (DPPA), dipalmitoyl phosphatidylglycerol (DPPG) and dipalmitoyl phosphatidylserine (DPPS), both spread on a Langmuir monolayer and deposited on a solid substrate to mimic a model membrane and study the interaction with the drug ELPT. X-ray reflectivity results pointed toward an increase in drug loading efficiency up to 13.5% mol/mol of ELPT into mixed systems DPPC/CL. This increase in loading efficiency was also accompanied by a slight distortion in the stacking of the bilayers less evidenced after optimization of the molar ratio between the co-lipids. Grazing incidence X-ray diffraction measurements revealed an in-plane lattice distortion due to the presence of hydrocarbon chain backbone ordering in pure systems of DPPC doped with ELPT. The same was not observed in mixed membranes with DPPC/CL and DPPC/DPPA.     

A semiparametric mixture model for analyzing clustered competing risks data

A very general class of multivariate life distributions is considered for analyzing failure time clustered data that are subject to censoring and multiple modes of failure. Conditional on cluster-specific quantities, the joint distribution of the failure time and event indicator can be expressed as a mixture of the distribution of time to failure due to a certain type (or specific cause), and the failure type distribution. We assume here the marginal probabilities of various failure types are logistic functions of some covariates. The cluster-specific quantities are subject to some unknown distribution that causes frailty. The unknown frailty distribution is modeled nonparametrically using a Dirichlet process. In such a semiparametric setup, a hybrid method of estimation is proposed based on the i.i.d. Weighted Chinese Restaurant algorithm that helps us generate observations from the predictive distribution of the frailty. The Monte Carlo ECM algorithm plays a vital role for obtaining the estimates of the parameters that assess the extent of the effects of the causal factors for failures of a certain type. A simulation study is conducted to study the consistency of our methodology. The proposed methodology is used to analyze a real data set on HIV infection of a cohort of female prostitutes in Senegal.     

Accounting for multiplicities in assessing drug safety: a three-level hierarchical mixture model

Multiple comparisons and other multiplicities are among the most difficult of problems that face statisticians, frequentists, and Bayesians alike. An example is the analysis of the many types of adverse events (AEs) that are recorded in drug clinical trials. We propose a three-level hierarchical mixed model. The most basic level is type of AE. The second level is body system, each of which contains a number of types of possibly related AEs. The highest level is the collection of all body systems. Our analysis allows for borrowing across body systems, but there is greater potential-depending on the actual data-for borrowing within each body system. The probability that a drug has caused a type of AE is greater if its rate is elevated for several types of AEs within the same body system than if the AEs with elevated rates were in different body systems. We give examples to illustrate our method and we describe its application to other types of problems.     

A semiparametric cure model for interval‐censored data

There is a growing interest in the analysis of survival data with a cured proportion particularly in tumor recurrences studies. Biologically, it is reasonable to assume that the recurrence time is mainly affected by the overall health condition of the patient that depends on some covariates such as age, sex, or treatment type received. We propose a semiparametric frailty‐Cox cure model to quantify the overall health condition of the patient by a covariate‐dependent frailty that has a discrete mass at zero to characterize the cured patients, and a positive continuous part to characterize the heterogeneous health conditions among the uncured patients. A multiple imputation estimation method is proposed for the right‐censored case, which is further extended to accommodate interval‐censored data. Simulation studies show that the performance of the proposed method is highly satisfactory. For illustration, the model is fitted to a set of right‐censored melanoma incidence data and a set of interval‐censored breast cosmesis data. Our analysis suggests that patients receiving treatment of radiotherapy with adjuvant chemotherapy have a significantly higher probability of breast retraction, but also a lower hazard rate of breast retraction among those patients who will eventually experience the event with similar health conditions. The interpretation is very different to those based on models without a cure component that the treatment of radiotherapy with adjuvant chemotherapy significantly increases the risk of breast retraction.     

A Bayesian semiparametric joint hierarchical model for longitudinal and survival data

This article proposes a new semiparametric Bayesian hierarchical model for the joint modeling of longitudinal and survival data. We relax the distributional assumptions for the longitudinal model using Dirichlet process priors on the parameters defining the longitudinal model. The resulting posterior distribution of the longitudinal parameters is free of parametric constraints, resulting in more robust estimates. This type of approach is becoming increasingly essential in many applications, such as HIV and cancer vaccine trials, where patients' responses are highly diverse and may not be easily modeled with known distributions. An example will be presented from a clinical trial of a cancer vaccine where the survival outcome is time to recurrence of a tumor. Immunologic measures believed to be predictive of tumor recurrence were taken repeatedly during follow-up. We will present an analysis of this data using our new semiparametric Bayesian hierarchical joint modeling methodology to determine the association of these longitudinal immunologic measures with time to tumor recurrence.     

A semiparametric joint model for cluster size and subunit-specific interval-censored outcomes

Clustered data frequently arise in biomedical studies, where observations, or subunits, measured within a cluster are associated. The cluster size is said to be informative, if the outcome variable is associated with the number of subunits in a cluster. In most existing work, the informative cluster size issue is handled by marginal approaches based on within-cluster resampling, or cluster-weighted generalized estimating equations. Although these approaches yield consistent estimation of the marginal models, they do not allow estimation of within-cluster associations and are generally inefficient. In this paper, we propose a semiparametric joint model for clustered interval-censored event time data with informative cluster size. We use a random effect to account for the association among event times of the same cluster as well as the association between event times and the cluster size. For estimation, we propose a sieve maximum likelihood approach and devise a computationally-efficient expectation-maximization algorithm for implementation. The estimators are shown to be strongly consistent, with the Euclidean components being asymptotically normal and achieving semiparametric efficiency. Extensive simulation studies are conducted to evaluate the finite-sample performance, efficiency and robustness of the proposed method. We also illustrate our method via application to a motivating periodontal disease dataset.     

A semiparametric model for the analysis of recurrent-event panel data

In many longitudinal studies, interest focuses on the occurrence rate of some phenomenon for the subjects in the study. When the phenomenon is nonterminating and possibly recurring, the result is a recurrent-event data set. Examples include epileptic seizures and recurrent cancers. When the recurring event is detectable only by an expensive or invasive examination, only the number of events occurring between follow-up times may be available. This article presents a semiparametric model for such data, based on a multiplicative intensity model paired with a fully flexible nonparametric baseline intensity function. A random subject-specific effect is included in the intensity model to account for the overdispersion frequently displayed in count data. Estimators are determined from quasi-likelihood estimating functions. Because only first- and second-moment assumptions are required for quasi-likelihood, the method is more robust than those based on the specification of a full parametric likelihood. Consistency of the estimators depends only on the assumption of the proportional intensity model. The semiparametric estimators are shown to be highly efficient compared with the usual parametric estimators. As with semiparametric methods in survival analysis, the method provides useful diagnostics for specific parametric models, including a quasi-score statistic for testing specific baseline intensity functions. The techniques are used to analyze cancer recurrences and a pheromone-based mating disruption experiment in moths. A simulation study confirms that, for many practical situations, the estimators possess appropriate small-sample characteristics.     

A new inducible RNAi xenograft model for assessing the staged tumor response to mTOR silencing

Human xenograft tumor models are widely used for efficacy evaluation of potential cancer targets. siRNA is usually stably introduced into tumor cells prior to transplantation. However, silencing of the cancer therapeutic target usually results in reduced cell growth/survival in vitro and/or failure to establish tumors in vivo, thus hindering tumor response-based efficacy evaluation. The present study explored a new tumor response model based on regulated RNAi, which is more relevant from a clinical standpoint. As a proof of principle, an inducible lentiviral RNAi vector was used to silence the known cancer therapeutic target mTOR upon induction with Doxycycline (DOX). The responses to DOX-induced mTOR silencing were tested both in vitro and in vivo for prostate cancer PC3 models. Significant reduction in cancer cell survival was observed due to cell cycle arrest and apoptosis when mTOR silencing was induced in vitro. mTOR silencing also caused tumor regression for the early-staged PC3 tumors (100% tumor regressed and 45% became tumor-free). The advanced-staged tumors also demonstrated significant responses (100% regressed). Therefore, our results demonstrate the powerful utility of this new inducible xenograft tumor model for efficacy evaluation of cancer targets, and it provides a direct in vivo efficacy validation of mTOR as a cancer therapeutic target.     

A semiparametric transition model with latent traits for longitudinal multistate data

SUMMARY: We propose a general multistate transition model. The model is developed for the analysis of repeated episodes of multiple states representing different health status. Transitions among multiple states are modeled jointly using multivariate latent traits with factor loadings. Different types of state transition are described by flexible transition-specific nonparametric baseline intensities. A state-specific latent trait is used to capture individual tendency of the sojourn in the state that cannot be explained by covariates and to account for correlation among repeated sojourns in the same state within an individual. Correlation among sojourns across different states within an individual is accounted for by the correlation between the different latent traits. The factor loadings for a latent trait accommodate the dependence of the transitions to different competing states from a same state. We obtain the semiparametric maximum likelihood estimates through an expectation-maximization (EM) algorithm. The method is illustrated by studying repeated transitions between independence and disability states of activities of daily living (ADL) with death as an absorbing state in a longitudinal aging study. The performance of the estimation procedure is assessed by simulation studies.     

A measurement error model for time-series studies of air pollution and mortality

One barrier to interpreting the observational evidence concerning the adverse health effects of air pollution for public policy purposes is the measurement error inherent in estimates of exposure based on ambient pollutant monitors. Exposure assessment studies have shown that data from monitors at central sites may not adequately represent personal exposure. Thus, the exposure error resulting from using centrally measured data as a surrogate for personal exposure can potentially lead to a bias in estimates of the health effects of air pollution. This paper develops a multi-stage Poisson regression model for evaluating the effects of exposure measurement error on estimates of effects of particulate air pollution on mortality in time-series studies. To implement the model, we have used five validation data sets on personal exposure to PM10. Our goal is to combine data on the associations between ambient concentrations of particulate matter and mortality for a specific location, with the validation data on the association between ambient and personal concentrations of particulate matter at the locations where data have been collected. We use these data in a model to estimate the relative risk of mortality associated with estimated personal-exposure concentrations and make a comparison with the risk of mortality estimated with measurements of ambient concentration alone. We apply this method to data comprising daily mortality counts, ambient concentrations of PM10measured at a central site, and temperature for Baltimore, Maryland from 1987 to 1994. We have selected our home city of Baltimore to illustrate the method; the measurement error correction model is general and can be applied to other appropriate locations.Our approach uses a combination of: (1) a generalized additive model with log link and Poisson error for the mortality-personal-exposure association; (2) a multi-stage linear model to estimate the variability across the five validation data sets in the personal-ambient-exposure association; (3) data augmentation methods to address the uncertainty resulting from the missing personal exposure time series in Baltimore. In the Poisson regression model, we account for smooth seasonal and annual trends in mortality using smoothing splines. Taking into account the heterogeneity across locations in the personal-ambient-exposure relationship, we quantify the degree to which the exposure measurement error biases the results toward the null hypothesis of no effect, and estimate the loss of precision in the estimated health effects due to indirectly estimating personal exposures from ambient measurements.     

A microstructurally informed model for the mechanical response of three-dimensional actin networks

We propose a class of microstructurally informed models for the linear elastic mechanical behaviour of cross-linked polymer networks such as the actin cytoskeleton. Salient features of the models include the possibility to represent anisotropic mechanical behaviour resulting from anisotropic filament distributions, and a power law scaling of the mechanical properties with the filament density. Mechanical models within the class are parameterized by seven different constants. We demonstrate a procedure for determining these constants using finite element models of three-dimensional actin networks. Actin filaments and cross-links were modelled as elastic rods, and the networks were constructed at physiological volume fractions and at the scale of an image voxel. We show the performance of the model in estimating the mechanical behaviour of the networks over a wide range of filament densities and degrees of anisotropy.     

A cerebrovascular response model for functional neuroimaging including dynamic cerebral autoregulation

Functional neuroimaging techniques such as functional magnetic resonance imaging (fMRI) and near-infrared spectroscopy (NIRS) can be used to isolate an evoked response to a stimulus from significant background physiological fluctuations. Data analysis approaches typically use averaging or linear regression to remove this physiological baseline with varying degrees of success. Biophysical model-based analysis of the functional hemodynamic response has also been advanced previously with the Balloon and Windkessel models. In the present work, a biophysical model of systemic and cerebral circulation and gas exchange is applied to resting state NIRS neuroimaging data from 10 human subjects. The model further includes dynamic cerebral autoregulation, which modulates the cerebral arteriole compliance to control cerebral blood flow. This biophysical model allows for prediction, from noninvasive blood pressure measurements, of the background hemodynamic fluctuations in the systemic and cerebral circulations. Significantly higher correlations with the NIRS data were found using the biophysical model predictions compared to blood pressure regression and compared to transfer function analysis (multifactor ANOVA, p < 0.0001). This finding supports the further development and use of biophysical models for removing baseline activity in functional neuroimaging analysis. Future extensions of this work could model changes in cerebrovascular physiology that occur during development, aging, and disease.