Bayesian joint modeling of longitudinal and survival data

The joint modeling of longitudinal and survival data is a new approach to many applications such as HIV, cancer vaccine trials and quality of life studies. There are recent developments of the methodologies with respect to each of the components of the joint model as well as statistical processes that link them together. Among these, second order polynomial random effect models and linear mixed effects models are the most commonly used for the longitudinal trajectory function. In this study, we first relax the parametric constraints for polynomial random effect models by using Dirichlet process priors, then three longitudinal markers rather than only one marker are considered in one joint model. Second, we use a linear mixed effect model for the longitudinal process in a joint model analyzing the three markers. In this research these methods were applied to the Primary Biliary Cirrhosis sequential data, which were collected from a clinical trial of primary biliary cirrhosis (PBC) of the liver. This trial was conducted between 1974 and 1984 at the Mayo Clinic. The effects of three longitudinal markers (1) Total Serum Bilirubin, (2) Serum Albumin and (3) Serum Glutamic-Oxaloacetic transaminase (SGOT) on patients' survival were investigated. Proportion of treatment effect will also be studied using the proposed joint modeling approaches. ^ Based on the results, we conclude that the proposed modeling approaches yield better fit to the data and give less biased parameter estimates for these trajectory functions than previous methods. Model fit is also improved after considering three longitudinal markers instead of one marker only. The results from analysis of proportion of treatment effects from these joint models indicate same conclusion as that from the final model of Fleming and Harrington (1991), which is Bilirubin and Albumin together has stronger impact in predicting patients' survival and as a surrogate endpoints for treatment. ^

Analysis of life expectancy using linear and nonlinear models

Life expectancy has consistently increased over the last 150 years due to improvements in nutrition, medicine, and public health. Several studies found that in many developed countries, life expectancy continued to rise following a nearly linear trend, which was contrary to a common belief that the rate of improvement in life expectancy would decelerate and was fit with an S-shaped curve. Using samples of countries that exhibited a wide range of economic development levels, we explored the change in life expectancy over time by employing both nonlinear and linear models. We then observed if there were any significant differences in estimates between linear models, assuming an auto-correlated error structure. When data did not have a sigmoidal shape, nonlinear growth models sometimes failed to provide meaningful parameter estimates. The existence of an inflection point and asymptotes in the growth models made them inflexible with life expectancy data. In linear models, there was no significant difference in the life expectancy growth rate and future estimates between ordinary least squares (OLS) and generalized least squares (GLS). However, the generalized least squares model was more robust because the data involved time-series variables and residuals were positively correlated. ^

Heart rate variability in response to pain stimulus in VLBW infants followed longitudinally during NICU stay.

The objective of this longitudinal study, conducted in a neonatal intensive care unit, was to characterize the response to pain of high-risk very low birth weight infants (<1,500 g) from 23 to 38 weeks post-menstrual age (PMA) by measuring heart rate variability (HRV). Heart period data were recorded before, during, and after a heel lanced or wrist venipunctured blood draw for routine clinical evaluation. Pain response to the blood draw procedure and age-related changes of HRV in low-frequency and high-frequency bands were modeled with linear mixed-effects models. HRV in both bands decreased during pain, followed by a recovery to near-baseline levels. Venipuncture and mechanical ventilation were factors that attenuated the HRV response to pain. HRV at the baseline increased with post-menstrual age but the growth rate of high-frequency power was reduced in mechanically ventilated infants. There was some evidence that low-frequency HRV response to pain improved with advancing PMA.

Assessment of the effect on statistical power of regression model misspecification by using techniques of mathematical statistics and simulation study

Objectives. This paper seeks to assess the effect on statistical power of regression model misspecification in a variety of situations. ^ Methods and results. The effect of misspecification in regression can be approximated by evaluating the correlation between the correct specification and the misspecification of the outcome variable (Harris 2010).In this paper, three misspecified models (linear, categorical and fractional polynomial) were considered. In the first section, the mathematical method of calculating the correlation between correct and misspecified models with simple mathematical forms was derived and demonstrated. In the second section, data from the National Health and Nutrition Examination Survey (NHANES 2007-2008) were used to examine such correlations. Our study shows that comparing to linear or categorical models, the fractional polynomial models, with the higher correlations, provided a better approximation of the true relationship, which was illustrated by LOESS regression. In the third section, we present the results of simulation studies that demonstrate overall misspecification in regression can produce marked decreases in power with small sample sizes. However, the categorical model had greatest power, ranging from 0.877 to 0.936 depending on sample size and outcome variable used. The power of fractional polynomial model was close to that of linear model, which ranged from 0.69 to 0.83, and appeared to be affected by the increased degrees of freedom of this model.^ Conclusion. Correlations between alternative model specifications can be used to provide a good approximation of the effect on statistical power of misspecification when the sample size is large. When model specifications have known simple mathematical forms, such correlations can be calculated mathematically. Actual public health data from NHANES 2007-2008 were used as examples to demonstrate the situations with unknown or complex correct model specification. Simulation of power for misspecified models confirmed the results based on correlation methods but also illustrated the effect of model degrees of freedom on power.^

Improving adherence in African American women with uncontrolled primary hypertension

Background and Objectives: African American (AA) women are disproportionately affected with hypertension (HTN). The aim of this randomized controlled trial was to evaluate the effectiveness of a 6-week culturally-tailored educational intervention for AA women with primary HTN who lived in rural Northeast Texas. ^ Methods: Sixty AA women, 29 to 86 years (M 57.98 ±12.37) with primary HTN were recruited from four rural locations and randomized to intervention (n =30) and wait-list control groups ( n =30) to determine the effectiveness of the intervention on knowledge, attitudes, beliefs, social support, adherence to a hypertension regimen, and blood pressure (BP) control. Survey and BP measurements were collected at baseline, 3 weeks, 6 weeks (post intervention) and 6 months post intervention. Culturally-tailored educational classes were provided for 90 minutes once a week for 6 weeks in two local churches and a community center. The wait-list control group received usual care and were offered education at the conclusion of the data collection six months post-intervention. Linear mixed models were used to test for differences between the groups. ^ Results: A significant overall main effect (Time) was found for systolic blood pressure, F(3, 174) =11.104, p=.000, and diastolic blood pressure. F(3, 174) =4.781, p=.003 for both groups. Age was a significant covariate for diastolic blood pressure. F(1, 56) =6.798 p=.012. Participants 57 years or older (n=30) had lower diastolic BPS than participants younger than 57 (n=30). No significant differences were found between groups on knowledge, adherence, or attitudes. Participants with lower incomes had significantly less knowledge about HBP Prevention (r=.036, p=.006). ^ Conclusion: AA women who participated in a 6 week intervention program demonstrated a significant decrease in BP over a 6 month period regardless of whether they were in the intervention or control group. These rural AA women had a relatively good knowledge of HTN and reported an average level of compliance, compared to other populations. Satisfaction with the program was high and there was no attrition, suggesting that AA women will participate in research studies that are culturally tailored to them, held in familiar community locations, and conducted by a trusted person with whom they can identify. Future studies using a different program with larger sample sizes are warranted to try to decrease the high level of HTN-related complications in AA women. ^

Statistical issues on studies of mixed longitudinal design

Mixed longitudinal designs are important study designs for many areas of medical research. Mixed longitudinal studies have several advantages over cross-sectional or pure longitudinal studies, including shorter study completion time and ability to separate time and age effects, thus are an attractive choice. Statistical methodology used in general longitudinal studies has been rapidly developing within the last few decades. Common approaches for statistical modeling in studies with mixed longitudinal designs have been the linear mixed-effects model incorporating an age or time effect. The general linear mixed-effects model is considered an appropriate choice to analyze repeated measurements data in longitudinal studies. However, common use of linear mixed-effects model on mixed longitudinal studies often incorporates age as the only random-effect but fails to take into consideration the cohort effect in conducting statistical inferences on age-related trajectories of outcome measurements. We believe special attention should be paid to cohort effects when analyzing data in mixed longitudinal designs with multiple overlapping cohorts. Thus, this has become an important statistical issue to address. ^ This research aims to address statistical issues related to mixed longitudinal studies. The proposed study examined the existing statistical analysis methods for the mixed longitudinal designs and developed an alternative analytic method to incorporate effects from multiple overlapping cohorts as well as from different aged subjects. The proposed study used simulation to evaluate the performance of the proposed analytic method by comparing it with the commonly-used model. Finally, the study applied the proposed analytic method to the data collected by an existing study Project HeartBeat!, which had been evaluated using traditional analytic techniques. Project HeartBeat! is a longitudinal study of cardiovascular disease (CVD) risk factors in childhood and adolescence using a mixed longitudinal design. The proposed model was used to evaluate four blood lipids adjusting for age, gender, race/ethnicity, and endocrine hormones. The result of this dissertation suggest the proposed analytic model could be a more flexible and reliable choice than the traditional model in terms of fitting data to provide more accurate estimates in mixed longitudinal studies. Conceptually, the proposed model described in this study has useful features, including consideration of effects from multiple overlapping cohorts, and is an attractive approach for analyzing data in mixed longitudinal design studies.^

Health outcomes, quality improvement and cost-effectiveness of the Community Diabetes Education (CoDE) program: Results of a randomized controlled trial

The objectives of this dissertation were to evaluate health outcomes, quality improvement measures, and the long-term cost-effectiveness and impact on diabetes-related microvascular and macrovascular complications of a community health worker-led culturally tailored diabetes education and management intervention provided to uninsured Mexican Americans in an urban faith-based clinic. A prospective, randomized controlled repeated measures design was employed to compare the intervention effects between: (1) an intervention group (n=90) that participated in the Community Diabetes Education (CoDE) program along with usual medical care; and (2) a wait-listed comparison group (n=90) that received only usual medical care. Changes in hemoglobin A1c (HbA1c) and secondary outcomes (lipid status, blood pressure and body mass index) were assessed using linear mixed-models and an intention-to-treat approach. The CoDE group experienced greater reduction in HbA1c (-1.6%, p<.001) than the control group (-.9%, p<.001) over the 12 month study period. After adjusting for group-by-time interaction, antidiabetic medication use at baseline, changes made to the antidiabetic regime over the study period, duration of diabetes and baseline HbA1c, a statistically significant intervention effect on HbA1c (-.7%, p=.02) was observed for CoDE participants. Process and outcome quality measures were evaluated using multiple mixed-effects logistic regression models. Assessment of quality indicators revealed that the CoDE intervention group was significantly more likely to have received a dilated retinal examination than the control group, and 53% achieved a HbA1c below 7% compared with 38% of control group subjects. Long-term cost-effectiveness and impact on diabetes-related health outcomes were estimated through simulation modeling using the rigorously validated Archimedes Model. Over a 20 year time horizon, CoDE participants were forecasted to have less proliferative diabetic retinopathy, fewer foot ulcers, and reduced numbers of foot amputations than control group subjects who received usual medical care. An incremental cost-effectiveness ratio of $355 per quality-adjusted life-year gained was estimated for CoDE intervention participants over the same time period. The results from the three areas of program evaluation: impact on short-term health outcomes, quantification of improvement in quality of diabetes care, and projection of long-term cost-effectiveness and impact on diabetes-related health outcomes provide evidence that a community health worker can be a valuable resource to reduce diabetes disparities for uninsured Mexican Americans. This evidence supports formal integration of community health workers as members of the diabetes care team.^

Changes in the PQRST intervals and heart rate variability associated with rewarming in two newborns undergoing hypothermia therapy.

BACKGROUND: Little is known about the effects of hypothermia therapy and subsequent rewarming on the PQRST intervals and heart rate variability (HRV) in term newborns with hypoxic-ischemic encephalopathy (HIE). OBJECTIVES: This study describes the changes in the PQRST intervals and HRV during rewarming to normal core body temperature of 2 newborns with HIE after hypothermia therapy. METHODS: Within 6 h after birth, 2 newborns with HIE were cooled to a core body temperature of 33.5 degrees C for 72 h using a cooling blanket, followed by gradual rewarming (0.5 degrees C per hour) until the body temperature reached 36.5 degrees C. Custom instrumentation recorded the electrocardiogram from the leads used for clinical monitoring of vital signs. Generalized linear mixed models were calculated to estimate temperature-related changes in PQRST intervals and HRV. Results: For every 1 degrees C increase in body temperature, the heart rate increased by 9.2 bpm (95% CI 6.8-11.6), the QTc interval decreased by 21.6 ms (95% CI 17.3-25.9), and low and high frequency HRV decreased by 0.480 dB (95% CI 0.052-0.907) and 0.938 dB (95% CI 0.460-1.416), respectively. CONCLUSIONS: Hypothermia-induced changes in the electrocardiogram should be monitored carefully in future studies.

The volume effect in irradiated mouse colorectum

Damage of the colorectum is the dose-limiting normal tissue complication following radiotherapy of prostate and cervical cancers. One approach for decreasing complications is to physically reduce the treatment volume. Mathematical models have been previously developed to describe the change in associated toxicity with a change in irradiated volume, i.e. the "volume effect", for serial-type normal tissues including the colorectum. The first goal of this thesis was to test the hypothesis that there would not be a threshold length in the development of obstruction after irradiation of mouse colorectum, as predicted by the Probability model of the volume effect. The second goal was to examine if there were differences in the threshold and in the incidence of colorectal obstruction after irradiation of two mouse strains, C57B1/6 (C57) and C3Hf/Kam (C3H), previously found to be fibrosis-prone and-resistant, respectively, after lung irradiation due, in part, to genetic differences. The hypothesis examined was that differences in incidence between strains were due to the differential expression of the fibrogenic cytokines $\rm TGF\beta$ and $\rm TNF\alpha.$ Various lengths of C57 and C3H mouse colorectum were irradiated and the incidence of colorectal obstruction was followed up to 15 months. A threshold length was observed for both mouse strains, in contradiction of model predictions. The mechanism of the threshold was epithelial regeneration after irradiation. C57 mice had significantly higher incidence of colorectal obstruction compared to C3H mice, especially at smaller irradiated lengths. Colorectal tissue was obtained at various times after irradiation and prepared for histology, immunohistochemistry and RNase protection assay for measurement of $\rm TGF\beta 1,$ 2, 3 and $\rm TNF\alpha$ mRNA. Distinct strain differences in the histological time of appearance and spatial locations of fibrosis were observed. However, there were no consistent strain difference in mRNA levels or immunolocalization for any of the cytokines examined. The data indicate the need for volume effect models that account for biologically important processes, such as the effect of epithelial regeneration after irradiation. As well, changes in fibrogenic cytokines at the mRNA level do not contribute to the strain difference in radiation-induced colorectal obstruction. ^

New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

Genes to blood pressure: The effects of variation in genes of the renin-angiotensin system on the hormonal and renal control of blood pressure

Objective. Essential hypertension affects 25% of the US adult population and is a leading contributor to morbidity and mortality. Because BP is a multifactorial phenotype that resists simple genetic analysis, intermediate phenotypes within the complex network of BP regulatory systems may be more accessible to genetic dissection. The Renin-Angiotensin System (RAS) is known to influence intermediate and long-term blood pressure regulation through alterations in vascular tone and renal sodium and fluid resorption. This dissertation examines associations between renin (REN), angiotensinogen (AGT), angiotensin-converting enzyme (ACE) and angiotensin II type 1 receptor (AT1) gene variation and interindividual differences in plasma hormone levels, renal hemodynamics, and BP homeostasis.^ Methods. A total of 150 unrelated men and 150 unrelated women, between 20.0 and 49.9 years of age and free of acute or chronic illness except for a history of hypertension (11 men and 7 women, all off medications), were studied after one week on a controlled sodium diet. RAS plasma hormone levels, renal hemodynamics and BP were determined prior to and during angiotensin II (Ang II) infusion. Individuals were genotyped by PCR for a variable number tandem repeat (VNTR) polymorphism in REN, and for the following restriction fragment length polymorphisms (RFLP): AGT M235T, ACE I/D, and AT1 A1166C. Associations between clinical measurements and allelic variation were examined using multiple linear regression statistical models.^ Results. Women homozygous for the AT1 1166C allele demonstrated higher intracellular levels of sodium (p = 0.044). Men homozygous for the AGT T235 allele demonstrated a blunted decrement in renal plasma flow in response to Ang II infusion (p = 0.0002). There were no significant associations between RAS gene variation and interindividual variation in RAS plasma hormone levels or BP.^ Conclusions. Rather than identifying new BP controlling genes or alleles, the study paradigm employed in this thesis (i.e., measured genes, controlled environments and interventions) may provide mechanistic insight into how candidate genes affect BP homeostasis. ^

Assessing homogeneity in workplace exposure under stationary and non-stationary scenarios

In most epidemiological studies, historical monitoring data are scant and must be pooled to identify occupational groups with homogeneous exposures. Homogeneity of exposure is generally assessed in a group of workers who share a common job title or work in a common area. While published results suggest that the degree of homogeneity varies widely across job groups, less is known whether such variation differs across industrial sectors, classes of contaminants, or in the methods used to group workers. Relying upon a compilation of results presented in the literature, patterns of homogeneity among nearly 500 occupational groups of workers were evaluated on the basis of type of industry and agent. Additionally, effects of the characteristics of the sampling strategy on estimated indicators of homogeneity of exposure were assessed. ^ Exposure profiles for occupational groups of workers have typically been assessed under the assumption of stationarity, i.e., the mean exposure level and variance of the distribution that describes the underlying population of exposures are constant over time. Yet, the literature has shown that occupational exposures have declined in the last decades. This renders traditional methods for the description of exposure profiles inadequate. Thus, work was needed to develop appropriate methods to assess homogeneity for groups of workers whose exposures have changed over time. A study was carried out applying mixed effects models with a term for temporal trend to appropriately describe exposure profiles of groups of workers in the nickel-producing industry over a 20-year period. Using a sub-set of groups of nickel-exposed workers, another study was conducted to develop and apply a framework to evaluate the assumption of stationarity of the variances in the presence of systematic changes in exposure levels over time. ^

The impact of a child's unexpected hospitalization upon Taiwanese parents' situational anxiety

Background. The parents of a sick child likely experience situational anxiety due to their young child being unexpectedly hospitalized. The emotional upheaval may be great enough that their anxiety inhibits them in providing positive support to their hospitalized child. Because anxiety affects psychological distress as well as behavioral distress, identifying parental distress helps parents improving their coping mechanisms. ^ Purpose. The study compared situational anxiety levels between Taiwanese fathers and mothers and focused on differences between parental anxiety levels at the beginning of the child's unplanned hospitalization and at time of discharge. The study also identified factors related to the parents' distress and use of coping mechanisms. ^ Methods. A descriptive, comparative research design was used to determine the difference between the anxiety levels of 62 Taiwanese father-mother dyads during the situational crisis of their child's unexpected hospitalization. The Mandarin version (M) of Visual Analog Scale (VAS-M), State-Trait Anxiety Inventory (STAI-M), and the Index of Parent Participation/Hospitalized Child (IPP/HC-M) were used to differentiate maternal and paternal anxiety levels and identify factors related to the parents' distress. Questionnaires were completed by parents within 24-36 hours of the child's hospital admission and within 24 hours prior to discharge. A paired t-test, two sample t-test, and linear mixed regression model were used to test and support the study hypothesis. ^ Results. The findings reveal that the mothers' anxiety levels did not significantly differ from the fathers' anxiety level when their child had a sudden admission to the hospital. In particular, parental state anxiety levels did not decrease during the child's hospital stay and subsequent discharge. Moreover, anxiety levels did not differ between parents regardless of whether the child's disease was acute or chronic. The most effective factor related to parental situational anxiety was parental perception of the severity of the child's illness. ^ Conclusions. Parental anxiety was found to be significantly related to changes in their perception of the severity of their child's illness. However, the study was not able to illustrate how parental involvement in the child's hospital care was related to parental perception of the severity of their child's illness. Future studies, using a qualitative approach to gamer more information as to what variables influence parental anxiety during a situational crisis, may provide a richer database from which to modify key variables as well as the instruments used to improve the quality of the data obtained. ^

Understanding disparities in emergency department visits for asthma

Asthma is the most common chronic disorder in childhood, affecting an estimated 6.2 million children under 18 years (1). The purpose of this study was to look at individual- and community-level characteristics simultaneously to examine and explain the factors that contribute to the use of emergency department services by children 18 years old or less and to determine if there was an association between air quality and ED visits in the same population, from 2005-2007 in Houston/Harris County. Data were collected from the Houston Safety Net Hospital Emergency Department Use Study and the 2000 US Census. Bivariate and multivariate logistic regression models and mixed effects models were used to analyze data that was collected during the study period.^ There were 704,902 ED visits made by children 18 and younger, who were living in Houston from January 1, 2005 to December 31, 2007. Of those, 19,098 had a primary discharge diagnosis of asthma. Asthma ED visits varied by season, with proportions of ED visits for asthma highest from September-December. African-American children were 2.6 (95% CI, 2.43-2.66) times more likely to have an ED visit for asthma compared to White children. Poverty, single parent headed households, and younger age all a greater likelihood of having gone to the ED for asthma treatment. Compared to Whites living in lightly-monitored pollution areas, African-Americans and Hispanics living in heavily monitored areas were 1.15 (95% CI, 1.04-1.28) times more likely to have an ED visit for asthma.^ Race and poverty seem to account for a large portion of the disparities in ED use found among children. This was true even after accounting for multiple individual- and community-level variables. These results suggest that racial disparities in asthma continue to pose risks for African American children, and they point to the need for additional research into potential explanations and remedies. Programs to reduce inappropriate ED use must be sensitive to an array of complex socioeconomic issues within minority and income populations. ^

Determinants of CD4 cell count decline after treatment interruption

Patients who had started HAART (Highly Active Anti-Retroviral Treatment) under previous aggressive DHHS guidelines (1997) underwent a life-long continuous HAART that was associated with many short term as well as long term complications. Many interventions attempted to reduce those complications including intermittent treatment also called pulse therapy. Many studies were done to study the determinants of rate of fall in CD4 count after interruption as this data would help guide treatment interruptions. The data set used here was a part of a cohort study taking place at the Johns Hopkins AIDS service since January 1984, in which the data were collected both prospectively and retrospectively. The patients in this data set consisted of 47 patients receiving via pulse therapy with the aim of reducing the long-term complications. ^ The aim of this project was to study the impact of virologic and immunologic factors on the rate of CD4 loss after treatment interruption. The exposure variables under investigation included CD4 cell count and viral load at treatment initiation. The rates of change of CD4 cell count after treatment interruption was estimated from observed data using advanced longitudinal data analysis methods (i.e., linear mixed model). Using random effects accounted for repeated measures of CD4 per person after treatment interruption. The regression coefficient estimates from the model was then used to produce subject specific rates of CD4 change accounting for group trends in change. The exposure variables of interest were age, race, and gender, CD4 cell counts and HIV RNA levels at HAART initiation. ^ The rate of fall of CD4 count did not depend on CD4 cell count or viral load at initiation of treatment. Thus these factors may not be used to determine who can have a chance of successful treatment interruption. CD4 and viral load were again studied by t-tests and ANOVA test after grouping based on medians and quartiles to see any difference in means of rate of CD4 fall after interruption. There was no significant difference between the groups suggesting that there was no association between rate of fall of CD4 after treatment interruption and above mentioned exposure variables. ^