Age -period -cohort analysis: An extension of Cox's lifetable regression with time -dependent covariates

It is well known that an identification problem exists in the analysis of age-period-cohort data because of the relationship among the three factors (date of birth + age at death = date of death). There are numerous suggestions about how to analyze the data. No one solution has been satisfactory. The purpose of this study is to provide another analytic method by extending the Cox's lifetable regression model with time-dependent covariates. The new approach contains the following features: (1) It is based on the conditional maximum likelihood procedure using a proportional hazard function described by Cox (1972), treating the age factor as the underlying hazard to estimate the parameters for the cohort and period factors. (2) The model is flexible so that both the cohort and period factors can be treated as dummy or continuous variables, and the parameter estimations can be obtained for numerous combinations of variables as in a regression analysis. (3) The model is applicable even when the time period is unequally spaced.^ Two specific models are considered to illustrate the new approach and applied to the U.S. prostate cancer data. We find that there are significant differences between all cohorts and there is a significant period effect for both whites and nonwhites. The underlying hazard increases exponentially with age indicating that old people have much higher risk than young people. A log transformation of relative risk shows that the prostate cancer risk declined in recent cohorts for both models. However, prostate cancer risk declined 5 cohorts (25 years) earlier for whites than for nonwhites under the period factor model (0 0 0 1 1 1 1). These latter results are similar to the previous study by Holford (1983).^ The new approach offers a general method to analyze the age-period-cohort data without using any arbitrary constraint in the model. ^

TIME DEPENDENT ALTERATIONS IN CENTRAL NERVOUS SYSTEM DOPAMINE RECEPTORS INDUCED BY DOPAMINE AGONIST TREATMENT

Levodopa, the precursor of dopamine, is currently the drug of choice in the treatment of Parkinson's disease. Recently, two direct dopamine agonists, bromocriptine and pergolide, have been tested for the treatment of Parkinson's disease because of reduced side effects compared to levodopa. Few studies have evaluated the effects of long-term treatment of dopamine agonists on dopamine receptor regulation in the central nervous system. Thus, the purpose of this study was to determine whether chronic dopamine agonist treatment produces a down-regulation of striatal dopamine receptor function and to compare the results of the two classes of dopaminergic drugs.^ Levodopa with carbidopa, a peripheral decarboxylase inhibitor, was administered orally to rats whereas bromocriptine and pergolide were injected intraperitoneally once daily. Several neurochemical parameters were examined from 1 to 28 days.^ Levodopa minimally decreased striatal D-1 receptor activity but increased the number of striatal D-2 binding sites. Levodopa increased the V(,max) of tyrosine hydroxylase (TH) in all brain regions tested. Protein blot analysis of striatal TH indicated a significant increase in the amount of TH present. Dopamine-beta-hydroxylase (DBH) activity was markedly decreased in all brain regions studied and mixing experiments of control and drug-treated cortices did not show the presence of an increased level of endogenous inhibitors.^ Bromocriptine treatment decreased the number of D-2 binding sites. Striatal TH activity was decreased and protein blot analysis indicated no change in TH quantity. The specificity of bromocriptine for striatal TH suggested that bromocriptine preferentially interacts with dopamine autoreceptors.^ Combination levodopa-bromocriptine was administered for 12 days. There was a decrease in both D-1 receptor activity and D-2 binding sites, and a decrease in brain HVA levels suggesting a postsynaptic receptor action. Pergolide produced identical results to the combination levodopa-bromocriptine studies.^ In conclusion, combination levodopa-bromocriptine and pergolide treatments exhibited the expected down-regulation of dopamine receptor activity. In contrast, levodopa appeared to up-regulate dopamine receptor activity. Thus, these data may help to explain, on a biochemical basis, the decrease in the levodopa-induced side effects noted with combination levodopa-bromocriptine or pergolide therapies in the treatment of Parkinson's disease. ^

Analysis of recurrent failure times: A time-dependent Yule process approach

Analysis of recurrent events has been widely discussed in medical, health services, insurance, and engineering areas in recent years. This research proposes to use a nonhomogeneous Yule process with the proportional intensity assumption to model the hazard function on recurrent events data and the associated risk factors. This method assumes that repeated events occur for each individual, with given covariates, according to a nonhomogeneous Yule process with intensity function λx(t) = λ 0(t) · exp( x′β). One of the advantages of using a non-homogeneous Yule process for recurrent events is that it assumes that the recurrent rate is proportional to the number of events that occur up to time t. Maximum likelihood estimation is used to provide estimates of the parameters in the model, and a generalized scoring iterative procedure is applied in numerical computation. ^ Model comparisons between the proposed method and other existing recurrent models are addressed by simulation. One example concerning recurrent myocardial infarction events compared between two distinct populations, Mexican-American and Non-Hispanic Whites in the Corpus Christi Heart Project is examined. ^

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.) ^

AN ILLNESS-DEATH PROCESS WITH TIME-DEPENDENT COVARIATES

A general model for the illness-death stochastic process with covariates has been developed for the analysis of survival data. This model incorporates important baseline and time-dependent covariates to make proper adjustment for the transition probabilities and survival probabilities. The follow-up period is subdivided into small intervals and a constant hazard is assumed for each interval. An approximation formula is derived to estimate the transition parameters when the exact transition time is unknown.^ The method developed is illustrated by using data from a study on the prevention of the recurrence of a myocardial infarction and subsequent mortality, the Beta-Blocker Heart Attack Trial (BHAT). This method provides an analytical approach which simultaneously includes provision for both fatal and nonfatal events in the model. According to this analysis, the effectiveness of the treatment can be compared between the Placebo and Propranolol treatment groups with respect to fatal and nonfatal events. ^

The number and timing of time -dependent covariate measurements for the Cox regression model

The problem of analyzing data with updated measurements in the time-dependent proportional hazards model arises frequently in practice. One available option is to reduce the number of intervals (or updated measurements) to be included in the Cox regression model. We empirically investigated the bias of the estimator of the time-dependent covariate while varying the effect of failure rate, sample size, true values of the parameters and the number of intervals. We also evaluated how often a time-dependent covariate needs to be collected and assessed the effect of sample size and failure rate on the power of testing a time-dependent effect.^ A time-dependent proportional hazards model with two binary covariates was considered. The time axis was partitioned into k intervals. The baseline hazard was assumed to be 1 so that the failure times were exponentially distributed in the ith interval. A type II censoring model was adopted to characterize the failure rate. The factors of interest were sample size (500, 1000), type II censoring with failure rates of 0.05, 0.10, and 0.20, and three values for each of the non-time-dependent and time-dependent covariates (1/4,1/2,3/4).^ The mean of the bias of the estimator of the coefficient of the time-dependent covariate decreased as sample size and number of intervals increased whereas the mean of the bias increased as failure rate and true values of the covariates increased. The mean of the bias of the estimator of the coefficient was smallest when all of the updated measurements were used in the model compared with two models that used selected measurements of the time-dependent covariate. For the model that included all the measurements, the coverage rates of the estimator of the coefficient of the time-dependent covariate was in most cases 90% or more except when the failure rate was high (0.20). The power associated with testing a time-dependent effect was highest when all of the measurements of the time-dependent covariate were used. An example from the Systolic Hypertension in the Elderly Program Cooperative Research Group is presented. ^

Time to renal disease and end-stage renal disease in PROFILE: a multiethnic lupus cohort.

BACKGROUND: Renal involvement is a serious manifestation of systemic lupus erythematosus (SLE); it may portend a poor prognosis as it may lead to end-stage renal disease (ESRD). The purpose of this study was to determine the factors predicting the development of renal involvement and its progression to ESRD in a multi-ethnic SLE cohort (PROFILE). METHODS AND FINDINGS: PROFILE includes SLE patients from five different United States institutions. We examined at baseline the socioeconomic-demographic, clinical, and genetic variables associated with the development of renal involvement and its progression to ESRD by univariable and multivariable Cox proportional hazards regression analyses. Analyses of onset of renal involvement included only patients with renal involvement after SLE diagnosis (n = 229). Analyses of ESRD included all patients, regardless of whether renal involvement occurred before, at, or after SLE diagnosis (34 of 438 patients). In addition, we performed a multivariable logistic regression analysis of the variables associated with the development of renal involvement at any time during the course of SLE.In the time-dependent multivariable analysis, patients developing renal involvement were more likely to have more American College of Rheumatology criteria for SLE, and to be younger, hypertensive, and of African-American or Hispanic (from Texas) ethnicity. Alternative regression models were consistent with these results. In addition to greater accrued disease damage (renal damage excluded), younger age, and Hispanic ethnicity (from Texas), homozygosity for the valine allele of FcgammaRIIIa (FCGR3A*GG) was a significant predictor of ESRD. Results from the multivariable logistic regression model that included all cases of renal involvement were consistent with those from the Cox model. CONCLUSIONS: Fcgamma receptor genotype is a risk factor for progression of renal disease to ESRD. Since the frequency distribution of FCGR3A alleles does not vary significantly among the ethnic groups studied, the additional factors underlying the ethnic disparities in renal disease progression remain to be elucidated.

Clinical, microbiological, and immunological markers of gingival and periodontal disease induction and progression

The potential impact of periodontal disease, a suspected risk factor for systemic diseases, presents challenges for health promotion and disease prevention strategies. This study examined clinical, microbiological, and immunological factors in a disease model to identify potential biomarkers that may be useful in predicting the onset and severity of both inflammatory and destructive periodontal disease. This project used an historical cohort design based on data obtained from 47 adult, female nonhuman primates followed over a 6-year period for 5 unique projects where the ligature-induced model of periodontitis was utilized. Standardization of protocols for sample collection allowed for comparison over time. Bleeding and pocket depth measures were selected as the dependent variables of relevance to humans based upon the literature and historical observations. Exposure variables included supragingival plaque, attachment level, total bacteria, black-pigmented bacteria, Gram-negative and Gram-positive bacteria, total IgG and IgA in crevicular fluid, specific IgG antibody in both crevicular fluid and serum, and IgG antibody to four select pathogenic microorganisms. Three approaches were used to analyze the data from this study. The first approach tested for differences in the means of the response variables within the group and among longitudinal observations within the group at each time point. The second approach examined the relationship among the clinical, microbiological, and immunological variables using correlation coefficients and stratified analyses. Multivariable models using GEE for repeated measures were produced as a predictive description of the induction and progression of gingivitis and periodontal disease. The multivariable models for bleeding (gingivitis) include supragingival plaque, total bacteria and total IgG while the second also contains supragingival plaque, Gram-positive bacteria, and total IgG. Two multivariable models emerged for periodontal disease. One multivariable model contains plaque, total bacteria, total IgG and attachment level. The second model includes black-pigmented bacteria, total bacteria, antibody to Campylobacter rectus, and attachment level. Utilization of the nonhuman primate model to prospectively examine causal hypotheses can provide a focus for human research on the mechanisms of progression from health to gingivitis to periodontitis. Ultimately, causal theories can guide strategies to prevent disease initiation and reduce disease severity. ^

The characterization of residential fungal spores and the relationship with housing characteristics in the Houston area

A number of indoor environmental factors, including bioaerosol or aeroallergen concentrations have been identified as exacerbators for asthma and allergenic conditions of the respiratory system. People generally spend 90% to 95% of their time indoors. Therefore, understanding the environmental factors that affect the presence of aeroallergens indoors as well as outdoors is important in determining their health impact, and in identifying potential intervention methods. This study aimed to assess the relationship between indoor airborne fungal spore concentrations and indoor surface mold levels, indoor versus outdoor airborne fungal spore concentrations and the effect of previous as well as current water intrusion. Also, the association between airborne concentration of indoor fungal spores and surface mold levels and the age of the housing structure were examined. Further, the correlation between indoor concentrations of certain species was determined as well. ^ Air and surface fungal measurements and related information were obtained from a Houston-area data set compiled from visits to homes filing insurance claims. During the sampling visit these complaint homes exhibited either visible mold or a combination of visible mold and water intrusion problems. These data were examined to assess the relationships between the independent and dependent variables using simple linear regression analysis, and independent t-tests. To examine the correlation between indoor concentrations of certain species, Spearman correlation coefficients were used. ^ There were 126 houses sampled, with spring, n=43 (34.1%), and winter, n=42 (33.3%), representing the seasons with the most samples. The summer sample illustrated the highest geometric mean concentration of fungal spores, GM=5,816.5 relative to winter, fall and spring (GM=1,743.4, GM=3,683.5 and GM=2,507.4, respectively). In all seasons, greater concentrations of fungal spores were observed during the cloudy weather conditions. ^ The results indicated no statistically significant association between outdoor total airborne fungal spore concentration and total living room airborne fungal spore concentration (β = 0.095, p = 0.491). Second, living room surface mold levels were not associated with living room airborne fungal spore concentration, (β= 0.011, p = 0.669). Third, houses with and without previous water intrusion did not differ significantly with respect to either living room (t(111) = 0.710, p = 0.528) or bedroom (t(111) =1.673, p = 0.162) airborne fungal spore concentrations. Likewise houses with and without current water intrusion did not differ significantly with respect to living room (t(109)=0.716, p = 0.476) or bedroom (t(109) = 1.035, p = 0.304) airborne fungal spore concentration. Fourth, houses with and without current water intrusion did not differ significantly with respect to living room (χ 2 (5) = 5.61, p = 0.346), or bedroom (χ 2 (5) = 1.80, p = 0.875) surface mold levels. Fifth, the age of the house structure did not predict living room (β = 0.023, p = 0.102) and bedroom (β = 0.023, p = 0.065) surface mold levels nor living room (β = 0.002, p = 0.131) and bedroom (β = 0.001, p = 0.650) fungal spore airborne concentration. Sixth, in houses with visually observed mold growth there was statistically significant differences between the mean living room concentrations and mean outdoor concentrations for Cladosporium (t (107) = 11.73, p < 0.0001), Stachybotrys (t (106)=2.288, p = 0.024, and Nigrosporia (t (102) = 2.267, p = 0.025). Finally, there was a significant correlation between several living room fungal species pairs, namely, Cladosporium and Stachybotrys (r = 0.373, p <0.01, n=65), Curvularia and Aspergillus/Penicillium (r = 0.205, p < 0.05, n= 111)), Curvularia and Stachybotrys (r = 0.205, p < 0.05, n=111), Nigrospora and Chaetomium (r = 0.254, p < 0.01, n=105) and Stachybotrys and Nigrospora (r = 0.269, p < 0.01, n=105). ^ This study has demonstrated several positive findings, i.e., significant pairwise correlations of concentrations of several fungal species in living room air, and significant differences between indoor and outdoor concentrations of three fungal species in homes with visible mold. No association was observed between indoor and outdoor fungal spore concentrations. Neither living room nor bedroom airborne spore concentrations and surface mold levels were related to the age of the house or to water intrusion, either previous or current. Therefore, these findings suggest the need for evaluating additional parameters, as well as combinations of factors such as humidity, temperature, age of structure, ventilation, and room size to better understand the determinants of airborne fungal spore concentrations and surface mold levels in homes. ^

Factors affecting prenatal care use and pregnancy outcome in Texas: A comparative study among and within Black, Anglo, and Hispanic populations

This research project is a study in the field of public health to test the relationships of demographic, socioeconomic, behavioral, and biological factors with (1) prenatal care use and (2) pregnancy outcome, measured by birth weight. It has been postulated that demographic, socioeconomic, and behavioral factors are associated with differences in the use of prenatal care services. It has also been postulated that differences in demographic, socioeconomic, behavioral, and biological factors result in differences in birth weight. This research attempts to test these two basic conceptual frameworks. At the same time, an attempt is made to determine the population groups and subgroups that are at increased risk (1) of using fewer prenatal care visits, and (2) of displaying a higher incidence of low birth weight babies. An understanding of these relationships of the demographic, socioeconomic, behavioral, and biological factors in the use of prenatal care visits and pregnancy outcome, measured by birth weight, will potentially offer guidance in the planning and policy development of maternal and child health services. The research considers four major components of maternal characteristics: (1) Demographic factors. Ethnicity, household size, maternal parity, and maternal age; (2) Socioeconomic factors. Maternal education, family income, maternal employment, health insurance coverage, and household dwelling; (3) Behavioral factors. Maternal smoking, attendance at child development classes, mother's first prenatal care visit, total number of prenatal care visits, and adequacy of care; and, (4) Biological factors. Maternal weight gain during pregnancy.^ The research considers 16 independent variables and two dependent variables.^ It was concluded that: (1) Generally, differences in demographic, socioeconomic, and behavioral factors were associated with differences in the average number of prenatal care visits between and within population groups and subgroups. The Hispanic mothers were the lowest users of prenatal care services. (2) In some cases, differences in demographic, socioeconomic, behavioral, and biological factors demonstrated differences in the average birth weight of infants between and within population groups and subgroups. (3) Differences in demographic, socioeconomic, behavioral, and biological factors resulted in differences in the rates of low birth weight babies between and within population groups and subgroups. The Black mothers delivered the highest incidence of low birth weight infants.^ These findings could provide guidance in the formulation of public health policies such as MCH services, an increase in the use of prenatal care services by prospective mothers, resulting in reduction of the incidence of low birth weight babies, and consequently aid in reducing the rates of infant mortality. ^

THE RELATIONSHIP BETWEEN DEGREE OF BLOOD PRESSURE REDUCTION AND MORTALITY AMONG HYPERTENSIVES IN THE HYPERTENSION DETECTION AND FOLLOW-UP PROGRAM

The relationship between degree of diastolic blood pressure (DBP) reduction and mortality was examined among hypertensives, ages 30-69, in the Hypertension Detection and Follow-up Program (HDFP). The HDFP was a multi-center community-based trial, which followed 10,940 hypertensive participants for five years. One-year survival was required for inclusion in this investigation since the one-year annual visit was the first occasion where change in blood pressure could be measured on all participants. During the subsequent four years of follow-up on 10,052 participants, 568 deaths occurred. For levels of change in DBP and for categories of variables related to mortality, the crude mortality rate was calculated. Time-dependent life tables were also calculated so as to utilize available blood pressure data over time. In addition, the Cox life table regression model, extended to take into account both time-constant and time-dependent covariates, was used to examine the relationship change in blood pressure over time and mortality.^ The results of the time-dependent life table and time-dependent Cox life table regression analyses supported the existence of a quadratic function which modeled the relationship between DBP reduction and mortality, even after adjusting for other risk factors. The minimum mortality hazard ratio, based on a particular model, occurred at a DBP reduction of 22.6 mm Hg (standard error = 10.6) in the whole population and 8.5 mm Hg (standard error = 4.6) in the baseline DBP stratum 90-104. After this reduction, there was a small increase in the risk of death. There was not evidence of the quadratic function after fitting the same model using systolic blood pressure. Methodologic issues involved in studying a particular degree of blood pressure reduction were considered. The confidence interval around the change corresponding to the minimum hazard ratio was wide and the obtained blood pressure level should not be interpreted as a goal for treatment. Blood pressure reduction was attributed, not only to pharmacologic therapy, but also to regression to the mean, and to other unknown factors unrelated to treatment. Therefore, the surprising results of this study do not provide direct implications for treatment, but strongly suggest replication in other populations. ^

Effect of patient accrual patterns on power and size of log-rank test for time-to-event outcome in clinical trials

The determination of size as well as power of a test is a vital part of a Clinical Trial Design. This research focuses on the simulation of clinical trial data with time-to-event as the primary outcome. It investigates the impact of different recruitment patterns, and time dependent hazard structures on size and power of the log-rank test. A non-homogeneous Poisson process is used to simulate entry times according to the different accrual patterns. A Weibull distribution is employed to simulate survival times according to the different hazard structures. The current study utilizes simulation methods to evaluate the effect of different recruitment patterns on size and power estimates of the log-rank test. The size of the log-rank test is estimated by simulating survival times with identical hazard rates between the treatment and the control arm of the study resulting in a hazard ratio of one. Powers of the log-rank test at specific values of hazard ratio (≠1) are estimated by simulating survival times with different, but proportional hazard rates for the two arms of the study. Different shapes (constant, decreasing, or increasing) of the hazard function of the Weibull distribution are also considered to assess the effect of hazard structure on the size and power of the log-rank test. ^

TRIM24-REGULATED ESTROGEN RESPONSE IS DEPENDENT ON SPECIFIC HISTONE MODIFICATIONS IN BREAST CANCER CELLS

In this dissertation, I discovered that function of TRIM24 as a co-activator of ERα-mediated transcriptional activation is dependent on specific histone modifications in tumorigenic human breast cancer-derived MCF7 cells. In the first part, I proved that TRIM24-PHD finger domain, which recognizes unmethylated histone H3 lysine K4 (H3K4me0), is critical for ERα-regulated transcription. Therefore, when LSD1-mediated demethylation of H3K4 is inhibited, activation of TRIM24-regulated ERα target genes is greatly impaired. Importantly, I demonstrated that TRIM24 and LSD1 are cyclically recruited to estrogen responsive elements (EREs) in a time-dependent manner upon estrogen induction, and depletion of their expression exert corresponding time-dependent effect on target gene activation. I also identified that phosphorylation of histone H3 threonine T6 disrupts TRIM24 from binding to the chromatin and from activating ERα-regulated targets. In the second part, I revealed that TRIM24 depletion has additive effect to LSD1 inhibitor- and Tamoxifen-mediated reduction in survival and proliferation in breast cancer cells.