Optimal, minimax, and Bayesian two-stage designs in two-dose phase II clinical trials

Background: For most cytotoxic and biologic anti-cancer agents, the response rate of the drug is commonly assumed to be non-decreasing with an increasing dose. However, an increasing dose does not always result in an appreciable increase in the response rate. This may especially be true at high doses for a biologic agent. Therefore, in a phase II trial the investigators may be interested in testing the anti-tumor activity of a drug at more than one (often two) doses, instead of only at the maximum tolerated dose (MTD). This way, when the lower dose appears equally effective, this dose can be recommended for further confirmatory testing in a phase III trial under potential long-term toxicity and cost considerations. A common approach to designing such a phase II trial has been to use an independent (e.g., Simon's two-stage) design at each dose ignoring the prior knowledge about the ordering of the response probabilities at the different doses. However, failure to account for this ordering constraint in estimating the response probabilities may result in an inefficient design. In this dissertation, we developed extensions of Simon's optimal and minimax two-stage designs, including both frequentist and Bayesian methods, for two doses that assume ordered response rates between doses. ^ Methods: Optimal and minimax two-stage designs are proposed for phase II clinical trials in settings where the true response rates at two dose levels are ordered. We borrow strength between doses using isotonic regression and control the joint and/or marginal error probabilities. Bayesian two-stage designs are also proposed under a stochastic ordering constraint. ^ Results: Compared to Simon's designs, when controlling the power and type I error at the same levels, the proposed frequentist and Bayesian designs reduce the maximum and expected sample sizes. Most of the proposed designs also increase the probability of early termination when the true response rates are poor. ^ Conclusion: Proposed frequentist and Bayesian designs are superior to Simon's designs in terms of operating characteristics (expected sample size and probability of early termination, when the response rates are poor) Thus, the proposed designs lead to more cost-efficient and ethical trials, and may consequently improve and expedite the drug discovery process. The proposed designs may be extended to designs of multiple group trials and drug combination trials.^

Neural analogs of associative learning in {\it Aplysia\/}: Experimental and theoretical approaches

A change in synaptic strength arising from the activation of two neuronal pathways at approximately the same time is a form of associative plasticity and may underlie classical conditioning. Previously, a cellular analog of a classical conditioning protocol has been demonstrated to produce short-term associative plasticity at the connections between sensory and motor neurons in Aplysia. A similar training protocol produced long-term (24 hour) enhancement of excitatory postsynaptic potentials (EPSPs). EPSPs produced by sensory neurons in which activity was paired with a reinforcing stimulus were significantly larger than unpaired controls 24 hours after training. To examined whether the associative plasticity observed at these synapses may be involved in higher-order forms of classical conditioning, a neural analog of contingency was developed. In addition, computer simulations were used to analyze whether the associative plasticity observed in Aplysia could, in theory, account for second-order conditioning and blocking. ^

Is sensitization to stimulants time dependent and mediated via NMDA receptors?

Chronic administration of psychomotor stimulants has been reported to produce behavioral sensitization to its effects on motor activity. This adaptation may be related to the pathophysiology of recurrent psychiatric disorders. Since disturbances in circadian rhythms are also found in many of these disorders, the relationship between sensitization and chronobiological factors became of interest. Therefore, a computerized monitoring system investigated the following: whether repeated exposure to methylphenidate (MPD) and amphetamine (AMP) could produce sensitization to its locomotor effects in the rat; whether sensitization to MPD and AMP was dependent on the circadian time of drug administration; whether the baseline levels of locomotor activity would be effected by repeated exposure to MPD and AMP; whether the expression of a sensitized response could be affected by the photoperiod; and whether MK-801, a non-competitive NMDA antagonist, could disrupt the development of sensitization to MPD. Dawley rats were housed in test cages and motor activity was recorded continuously for 16 days. The first 2 days served as baseline for each rat, and on day 3 each rat received a saline injection. The locomotor response to 0.6, 2.5, or 10 mg/kg of MPD was tested on day 4, followed by five days of single injections of 2.5 mg/kg MPD (days 5–9). After five days without injection (days 10–14) rats were re-challenged (day 15) with the same doses they received on day 4. There were three separate dose groups ran at four different times of administration, 08:00, 14:00, 20:00, or 02:00 (i.e. 12 groups). The same protocol was conducted with AMP with the doses of 0.3, 0.6, and 1.2 mg/kg given on day 4 and 15, and 0.6 mg/kg AMP as the repeated dose on days 5 to 9. In the second set of experiments only sensitization to MPD was investigated. The expression of the sensitized response was dose-dependent and mainly observed with challenge of the lower dose groups. The development of sensitization to MPD and ANT was differentially time-dependent. For MPD, the most robust sensitization occurred during the light phase, with no sensitization during the middle of the dark phase. (Abstract shortened by UMI.) ^

Follow-up care after a health fair screening in three Houston neighborhoods in 2008

Background. Each year thousands of people participate in mass health screenings for diabetes and hypertension, but little is known about whether or not those who receive higher than normal screening results obtain the recommended follow-up medical care, or what barriers they perceive to doing so. ^ Methods. Study participants were recruited from attendees at three health fairs in low-income neighborhoods in Houston, Texas Potential participants had higher than normal blood pressure (> 90/140 mgHg) or blood glucose readings (100 mm/dL fasting or 140 mm/dL random). Study participants were called at one, two, and three months and asked if they had obtained follow-up medical care; those who had not yet obtained follow-up care were asked to identify barriers. Using a modified Aday-Andersen model of health service access, the independent variables were individual and community characteristics and self-perceived need. The dependent variable was obtaining follow-up care, with barriers to care a secondary outcome. ^ Results. Eighty-two study participants completed the initial questionnaire and 59 participants completed the study protocol. Forty-eight participants (59% under an intent to treat analysis, 81% of those completing the study protocol) obtained follow-up care. Those who completed the initial questionnaire and who reported a regular source of care were significantly more likely to obtain follow-up care. For those who completed the study protocol the relationship between having a regular source of care and obtaining follow-up care approached but did not reach significance. For those who completed the initial questionnaire, self-described health status, when examined as a binary variable (good, very good, excellent, or poor, fair, not sure) was associated with obtaining follow-up care for those who rated their health as poor, fair, or not sure. While the group who completed the study protocol did not reach statistical significance, the same relationship between self-described health status of poor, fair, or not sure and obtaining follow-up care was present. The participants who completed the study protocol and described their blood pressure as OK or a little high were statistically more likely to get follow-up care than those who described it as high or very high. All those on oral medications for hypertension (12/12) and diabetes (4/4) who were told to obtain follow-up care did so; however, the small sample size allows this correlation to be of statistical significance only for those treating hypertension. ^ The variables significantly associated with obtaining follow-up care were having a regular source of care, self-described health status of poor, fair, or not sure, self-described blood pressure of OK or a little high, and taking medication for blood pressure. ^ At the follow-up telephone calls, 34 participants identified barriers to care; cost was a significant barrier reported by 16 participants, and 10 reported that they didn’t have time because they were working long hours after Hurricane Ike. ^ The study included the offer of access assistance: information about nearby safety-net providers, a visit to or information from the Health Information Center at their Neighborhood Center location, or information from Project Safety Net (a searchable web site for safety net providers). Access assistance was offered at the health fairs and then again at follow-up telephone calls to those who had not yet obtained follow-up care. Of the 48 participants who reported obtaining follow-up care, 26 said they had made use of the access assistance to do so. The use of access assistance was associated with being Hispanic, not having health insurance or a regular source of care, and speaking Spanish. It was also associated with being worried about blood glucose. ^ Conclusion. Access assistance, as a community enabling characteristic, may be useful in aiding low-income people in obtaining medical care. ^

A promotora de salud CVD intervention: Participant responses after Su Corazon, Su Vida education

In this study we sought to identify and understand feelings, benefits and barriers to making heart healthy behavioral changes by reviewing and analyzing participant responses to a follow-up telephone survey conducted as part of the HEART project (Health Education Awareness Research Team). Individuals who participated in HEART attended classes and received “Su Corazon, Su Vida” education. The HEART follow-up telephone survey was conducted only on those participants who were part of the experimental group. A total of 93 individuals from this group were successfully contacted for participation in the telephone survey after the classes ended. Quantitative data regarding ‘feelings’ and ‘difficulty making heart healthy behavioral changes’ were analyzed by calculating frequencies of each category of response for post-intervention weeks 9, 13, and 15. In addition, Wilcoxon rank-sum tests were conducted for post-intervention at weeks 9, 13, and 15 to measure associations between feelings and difficulties making heart healthy behavioral changes. Changes in responses over time for feelings and difficulties making heart healthy behavioral changes were looked at by counting differences in responses between pairs of follow up weeks. Qualitative responses to the survey were analyzed by categorizing content of responses under themes in order to identify factors related to feelings and difficulties making heart healthy behavioral changes. Telephone survey participants showed positive attitudes towards making nutritional and physical activity changes. Out of the 93 telephone survey respondents, 53 (57%) reported some type of physical activity change during the follow-up period while 46 (49%) reported specific changes in nutrition. Data from the “difficulty to making changes” responses were categorized under constructs from the Health Belief Model, perceived benefits and barriers. Overall, the barriers for physical activity were health issues, individual habits and time. Barriers to eating healthy were family support, individual habits, and knowledge. This study suggests that with respect to nutritional knowledge barriers, educational programs should explore other ways of teaching and familiarizing individuals with information sources that may be more appropriate for those populations not accustomed to them. For example, nutrition labels, portions, recipes, and use of photonovelas. Our findings of the barriers to changes in food preparation due to lack of family support may also suggest the need for the development of programs where influential partners or relatives are involved in order to create a more supportive environment which may provide more opportunity for change toward healthier lifestyle behaviors. Finally, the physical activity barriers found suggest that it may be beneficial to recommend appropriate exercises for those with specific health problems or those with time restrictions due to work or travel so that physical activity is not completely avoided.^

NONLINEAR TIME SERIES/SYSTEM ANALYSIS (STATE-SPACE, DYNAMICS, INTERVENTION, RESILIENCE, KALMAN)

A discussion of nonlinear dynamics, demonstrated by the familiar automobile, is followed by the development of a systematic method of analysis of a possibly nonlinear time series using difference equations in the general state-space format. This format allows recursive state-dependent parameter estimation after each observation thereby revealing the dynamics inherent in the system in combination with random external perturbations.^ The one-step ahead prediction errors at each time period, transformed to have constant variance, and the estimated parametric sequences provide the information to (1) formally test whether time series observations y(,t) are some linear function of random errors (ELEM)(,s), for some t and s, or whether the series would more appropriately be described by a nonlinear model such as bilinear, exponential, threshold, etc., (2) formally test whether a statistically significant change has occurred in structure/level either historically or as it occurs, (3) forecast nonlinear system with a new and innovative (but very old numerical) technique utilizing rational functions to extrapolate individual parameters as smooth functions of time which are then combined to obtain the forecast of y and (4) suggest a measure of resilience, i.e. how much perturbation a structure/level can tolerate, whether internal or external to the system, and remain statistically unchanged. Although similar to one-step control, this provides a less rigid way to think about changes affecting social systems.^ Applications consisting of the analysis of some familiar and some simulated series demonstrate the procedure. Empirical results suggest that this state-space or modified augmented Kalman filter may provide interesting ways to identify particular kinds of nonlinearities as they occur in structural change via the state trajectory.^ A computational flow-chart detailing computations and software input and output is provided in the body of the text. IBM Advanced BASIC program listings to accomplish most of the analysis are provided in the appendix. ^

Survival analysis of circulating tumor cells in triple-negative breast cancer

Background. Breast cancer is the most frequently diagnosed cancer and the leading cause of cancer death among females, accounting for 23% (1.38 million) of the total new cancer cases and 14% (458,400) of the total cancer deaths in 2008. [1] Triple-negative breast cancer (TNBC) is an aggressive phenotype comprising 10–20% of all breast cancers (BCs). [2-4] TNBCs show absence of estrogen, progesterone and HER2/neu receptors on the tumor cells. Because of the absence of these receptors, TNBCs are not candidates for targeted therapies. Circulating tumor cells (CTCs) are observed in blood of breast cancer patients even at early stages (Stage I & II) of the disease. Immunological and molecular analysis can be used to detect the presence of tumor cells in the blood (Circulating tumor cells; CTCs) of many breast cancer patients. These cells may explain relapses in early stage breast cancer patients even after adequate local control. CTC detection may be useful in identifying patients at risk for disease progression, and therapies targeting CTCs may improve outcome in patients harboring them. Methods . In this study we evaluated 80 patients with TNBC who are enrolled in a larger prospective study conducted at M D Anderson Cancer Center in order to determine whether the presence of circulating tumor cells is a significant prognostic factor in relapse free and overall survival . Patients with metastatic disease at the time of presentation were excluded from the study. CTCs were assessed using CellSearch System™ (Veridex, Raritan, NJ). CTCs were defined as nucleated cells lacking the presence of CD45 but expressing cytokeratins 8, 18 or 19. The distribution of patient and tumor characteristics was analyzed using chi square test and Fisher's exact test. Log rank test and Cox regression analysis was applied to establish the association of circulating tumor cells with relapse free and overall survival. Results. The median age of the study participants was 53years. The median duration of follow-up was 40 months. Eighty-eight percent (88%) of patients were newly diagnosed (without a previous history of breast cancer), and (60%) of patients were chemo naïve (had not received chemotherapy at the time of their blood draw for CTC analysis). Tumor characteristics such as stage (P=0.40), tumor size (P=69), sentinel nodal involvement (P=0.87), axillary lymph node involvement (P=0.13), adjuvant therapy (P=0.83), and high histological grade of tumor (P=0.26) did not predict the presence of CTCs. However, CTCs predicted worse relapse free survival (1 or more CTCs log rank P value = 0.04, at 2 or more CTCs P = 0.02 and at 3 or more CTCs P < 0.0001) and overall survival (at 1 or more CTCs log rank P value = 0.08, at 2 or more CTCs P = 0.01 and at 3 or more CTCs P = 0.0001. Conclusions. The number of circulating tumor cells predicted worse relapse free survival and overall survival in TNBC patients.^

TIME SERIES ANALYSIS AS INPUT FOR PREDICTIVE MODELING: PREDICTING CARDIAC ARREST IN A PEDIATRIC INTENSIVE CARE UNIT

The first manuscript, entitled "Time-Series Analysis as Input for Clinical Predictive Modeling: Modeling Cardiac Arrest in a Pediatric ICU" lays out the theoretical background for the project. There are several core concepts presented in this paper. First, traditional multivariate models (where each variable is represented by only one value) provide single point-in-time snapshots of patient status: they are incapable of characterizing deterioration. Since deterioration is consistently identified as a precursor to cardiac arrests, we maintain that the traditional multivariate paradigm is insufficient for predicting arrests. We identify time series analysis as a method capable of characterizing deterioration in an objective, mathematical fashion, and describe how to build a general foundation for predictive modeling using time series analysis results as latent variables. Building a solid foundation for any given modeling task involves addressing a number of issues during the design phase. These include selecting the proper candidate features on which to base the model, and selecting the most appropriate tool to measure them. We also identified several unique design issues that are introduced when time series data elements are added to the set of candidate features. One such issue is in defining the duration and resolution of time series elements required to sufficiently characterize the time series phenomena being considered as candidate features for the predictive model. Once the duration and resolution are established, there must also be explicit mathematical or statistical operations that produce the time series analysis result to be used as a latent candidate feature. In synthesizing the comprehensive framework for building a predictive model based on time series data elements, we identified at least four classes of data that can be used in the model design. The first two classes are shared with traditional multivariate models: multivariate data and clinical latent features. Multivariate data is represented by the standard one value per variable paradigm and is widely employed in a host of clinical models and tools. These are often represented by a number present in a given cell of a table. Clinical latent features derived, rather than directly measured, data elements that more accurately represent a particular clinical phenomenon than any of the directly measured data elements in isolation. The second two classes are unique to the time series data elements. The first of these is the raw data elements. These are represented by multiple values per variable, and constitute the measured observations that are typically available to end users when they review time series data. These are often represented as dots on a graph. The final class of data results from performing time series analysis. This class of data represents the fundamental concept on which our hypothesis is based. The specific statistical or mathematical operations are up to the modeler to determine, but we generally recommend that a variety of analyses be performed in order to maximize the likelihood that a representation of the time series data elements is produced that is able to distinguish between two or more classes of outcomes. The second manuscript, entitled "Building Clinical Prediction Models Using Time Series Data: Modeling Cardiac Arrest in a Pediatric ICU" provides a detailed description, start to finish, of the methods required to prepare the data, build, and validate a predictive model that uses the time series data elements determined in the first paper. One of the fundamental tenets of the second paper is that manual implementations of time series based models are unfeasible due to the relatively large number of data elements and the complexity of preprocessing that must occur before data can be presented to the model. Each of the seventeen steps is analyzed from the perspective of how it may be automated, when necessary. We identify the general objectives and available strategies of each of the steps, and we present our rationale for choosing a specific strategy for each step in the case of predicting cardiac arrest in a pediatric intensive care unit. Another issue brought to light by the second paper is that the individual steps required to use time series data for predictive modeling are more numerous and more complex than those used for modeling with traditional multivariate data. Even after complexities attributable to the design phase (addressed in our first paper) have been accounted for, the management and manipulation of the time series elements (the preprocessing steps in particular) are issues that are not present in a traditional multivariate modeling paradigm. In our methods, we present the issues that arise from the time series data elements: defining a reference time; imputing and reducing time series data in order to conform to a predefined structure that was specified during the design phase; and normalizing variable families rather than individual variable instances. The final manuscript, entitled: "Using Time-Series Analysis to Predict Cardiac Arrest in a Pediatric Intensive Care Unit" presents the results that were obtained by applying the theoretical construct and its associated methods (detailed in the first two papers) to the case of cardiac arrest prediction in a pediatric intensive care unit. Our results showed that utilizing the trend analysis from the time series data elements reduced the number of classification errors by 73%. The area under the Receiver Operating Characteristic curve increased from a baseline of 87% to 98% by including the trend analysis. In addition to the performance measures, we were also able to demonstrate that adding raw time series data elements without their associated trend analyses improved classification accuracy as compared to the baseline multivariate model, but diminished classification accuracy as compared to when just the trend analysis features were added (ie, without adding the raw time series data elements). We believe this phenomenon was largely attributable to overfitting, which is known to increase as the ratio of candidate features to class examples rises. Furthermore, although we employed several feature reduction strategies to counteract the overfitting problem, they failed to improve the performance beyond that which was achieved by exclusion of the raw time series elements. Finally, our data demonstrated that pulse oximetry and systolic blood pressure readings tend to start diminishing about 10-20 minutes before an arrest, whereas heart rates tend to diminish rapidly less than 5 minutes before an arrest.