A SIMULATION STUDY OF THE SIZE AND POWER OF SOME TWO-SAMPLE TESTS FOR UNEQUALLY CENSORED SURVIVAL DATA (RANK TESTS, LOGRANK, PROPORTIONAL HAZARDS, WILCOXON TESTS, EXPONENTIAL DISTRIBUTION)

Sizes and power of selected two-sample tests of the equality of survival distributions are compared by simulation for small samples from unequally, randomly-censored exponential distributions. The tests investigated include parametric tests (F, Score, Likelihood, Asymptotic), logrank tests (Mantel, Peto-Peto), and Wilcoxon-Type tests (Gehan, Prentice). Equal sized samples, n = 18, 16, 32 with 1000 (size) and 500 (power) simulation trials, are compared for 16 combinations of the censoring proportions 0%, 20%, 40%, and 60%. For n = 8 and 16, the Asymptotic, Peto-Peto, and Wilcoxon tests perform at nominal 5% size expectations, but the F, Score and Mantel tests exceeded 5% size confidence limits for 1/3 of the censoring combinations. For n = 32, all tests showed proper size, with the Peto-Peto test most conservative in the presence of unequal censoring. Powers of all tests are compared for exponential hazard ratios of 1.4 and 2.0. There is little difference in power characteristics of the tests within the classes of tests considered. The Mantel test showed 90% to 95% power efficiency relative to parametric tests. Wilcoxon-type tests have the lowest relative power but are robust to differential censoring patterns. A modified Peto-Peto test shows power comparable to the Mantel test. For n = 32, a specific Weibull-exponential comparison of crossing survival curves suggests that the relative powers of logrank and Wilcoxon-type tests are dependent on the scale parameter of the Weibull distribution. Wilcoxon-type tests appear more powerful than logrank tests in the case of late-crossing and less powerful for early-crossing survival curves. Guidelines for the appropriate selection of two-sample tests are given. ^

Effects of informative censoring on the hazard function: Application to the proportional hazards regression model

The standard analyses of survival data involve the assumption that survival and censoring are independent. When censoring and survival are related, the phenomenon is known as informative censoring. This paper examines the effects of an informative censoring assumption on the hazard function and the estimated hazard ratio provided by the Cox model.^ The limiting factor in all analyses of informative censoring is the problem of non-identifiability. Non-identifiability implies that it is impossible to distinguish a situation in which censoring and death are independent from one in which there is dependence. However, it is possible that informative censoring occurs. Examination of the literature indicates how others have approached the problem and covers the relevant theoretical background.^ Three models are examined in detail. The first model uses conditionally independent marginal hazards to obtain the unconditional survival function and hazards. The second model is based on the Gumbel Type A method for combining independent marginal distributions into bivariate distributions using a dependency parameter. Finally, a formulation based on a compartmental model is presented and its results described. For the latter two approaches, the resulting hazard is used in the Cox model in a simulation study.^ The unconditional survival distribution formed from the first model involves dependency, but the crude hazard resulting from this unconditional distribution is identical to the marginal hazard, and inferences based on the hazard are valid. The hazard ratios formed from two distributions following the Gumbel Type A model are biased by a factor dependent on the amount of censoring in the two populations and the strength of the dependency of death and censoring in the two populations. The Cox model estimates this biased hazard ratio. In general, the hazard resulting from the compartmental model is not constant, even if the individual marginal hazards are constant, unless censoring is non-informative. The hazard ratio tends to a specific limit.^ Methods of evaluating situations in which informative censoring is present are described, and the relative utility of the three models examined is discussed. ^

The heart challenge: The effect of regular aerobic exercise on body fat distribution in third and fourth grade school children: A randomized clinical trial

The pattern of body fat distribution known as "centralized", and characterized by a predominance of subcutaneous fat on the trunk and a "pot belly", has been associated with an increased risk of chronic disease. These patterns of fat distribution, as well as the lifestyle habit variables associated with adult fatness and chronic morbidity clearly begin to develop during childhood, indicating the need for intervention and primary prevention of obesity, particularly the centralized form, during childhood or adolescence. The purpose of this study was to determine whether regular aerobic exercise could beneficially alter the distribution of body fat in 8 and 9 year old children. One hundred and eighty-eight participants were randomized into either a regular aerobic exercise treatment group or a standard physical education program control group. A variety of aerobic activities was used for intervention 5 days per week during physical education class for a period of 12 weeks. Fat distribution was measured by a number of the most commonly used indices, including ratios of body circumferences and skinfolds and indices derived from a principal components analysis. Change over time in average pulse rate was used to determine if intervention actually occurred. Approximately 10% of the students were remeasured, allowing the calculation of intra- and interexaminer measurement reliability estimates for all indices.^ This study group was comparable to the U.S. population, though the study children were slightly larger for certain measures. No effect of the exercise intervention was found. The most likely explanation for this was inadequacy of the intervention, as indicated by the lack of any change in average pulse rate with treatment. The results of the measurement reliability analysis are reported and indicate that body circumference ratios are more precise than skinfold ratios, particularly when multiple observers are used. Reliability estimates for the principal component indices were also high.^ It remains unclear whether the distribution of body fat can be altered with exercise. It is likely that this issue will remain undecided until one highly reliable, valid, and sensitive measure of fat distribution can be found. ^

Cell surface galactosyltransferase structure and function during mesenchymal cell migration

In order to more fully understand the function of surface GalTase on mesenchymal cells, anti-GalTase IgG was used to (a) examine the role of surface GalTase during mouse mesenchymal cell migration on laminin and fibronectin; (b) define the plasma membrane distribution of GalTase by indirect immunofluorescence on migrating cells; (c) quantitate the level of surface GalTase on migrating cells; and (d) determine whether GalTase is associated with the cytoskeleton.^ Results show that anti-GalTase IgG was able to inhibit migration (48-80% as compared to basal rate) when cells were migrating on laminin-containing matrices. Monovalent Fab fragments inhibited migration on laminin by 90% after 4 hours. On the other hand, anti-GalTase IgG had no effect on cells migrating on fibronectin. This illustrates the substrate specificity of GalTase mediated-migration. When anti-GalTase IgG was used to localize surface GalTase on cells migratory on laminin, the enzyme was restricted to the leading and trailing edges of the cell. Assays indicate that GalTase is elevated approximately 3-fold when cells are migrating on laminin-containing matrices as compared to migratory cells on plastic or fibronectin, or as compared to stationary cells on any substrate. Laminin appears to recruit GalTase from preexisting intracellular pools to the growing lamellipodia.^ Double-label indirect immunofluorescence studies indicate that there is an apparent co-localization between some of the surface GalTase and some actin filaments. This relationship was explored by extracting cells prelabeled with anti-GalTase IgG and quantitated by radiolabeled second antibodies. Results show that 79% of the surface GalTase is associated with the cytoskeleton (as judged by detergent insolubility) when monovalent antibodies (Fab) are used. However virtually all (80-100%) of the surface GalTase can be induced to associate with the cytoskeleton when cross-linked with bivalent antibodies. Furthermore, when cells in suspension are incubated with divalent antibodies, an additional 66% of the surface GalTase can be induced to associate with the cytoskeleton. The elevated levels of surface GalTase detectable on cells migrating on laminin also appear to be associated with the cytoskeleton.^ Several lines of evidence suggest that GalTase is associated with F-actin. Data suggest that laminin induces the expression of surface GalTase to the growing lamellipodia where it becomes associated with the cytoskeleton leading to cell spreading and migration. (Abstract shortened with permission of author.) ^

Cytochrome P450 in mammalian brain: Purification, anatomical distribution and regulation of the cytochrome P450 monooxygenase system in rat brain

The cytochrome P450 (P450) monooxygenase system plays a major role in metabolizing a wide variety of xenobiotic as well as endogenous compounds. In performing this function, it serves to protect the body from foreign substances. However, in a number of cases, P450 activates procarcinogens to cause harm. In most animals, the highest level of activity is found in the liver. Virtually all tissues demonstrate P450 activity, though, and the role of the P450 monooxygenase system in these other organs is not well understood. In this project I have studied the P450 system in rat brain; purifying NADPH-cytochrome P450 reductase (reductase) from that tissue. In addition, I have examined the distribution and regulation of expression of reductase and P450 in various anatomical regions of the rat brain.^ NADPH-cytochrome P450 reductase was purified to apparent homogeneity and cytochrome P450 partially purified from whole rat brain. Purified reductase from brain was identical to liver P450 reductase by SDS-PAGE and Western blot techniques. Kinetic studies utilizing cerebral P450 reductase reveal Km values in close agreement with those determined with enzyme purified from rat liver. Moreover, the brain P450 reductase was able to function successfully in a reconstituted microsomal system with partially purified brain cytochrome P450 and with purified hepatic P4501A1 as measured by 7-ethoxycoumarin and 7-ethoxyresorufin O-deethylation. These results indicate that the reductase and P450 components may interact to form a competent drug metabolism system in brain tissue.^ Since the brain is not a homogeneous organ, dependent upon the well orchestrated interaction of numerous parts, pathology in one nucleus may have a large impact upon its overall function. Hence, the anatomical distribution of the P450 monooxygenase system in brain is important in elucidating its function in that organ. Related to this is the regulation of P450 expression in brain. In order to study these issues female rats--both ovariectomized and not--were treated with a number of xenobiotic compounds and sex steroids. The brains from these animals were dissected into 8 discrete regions and the presence and relative level of message for P4502D and reductase determined using polymerase chain reaction. Results of this study indicate the presence of mRNA for reductase and P4502D isoforms throughout the rat brain. In addition, quantitative PCR has allowed the determination of factors affecting the expression of message for these enzymes. ^

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

Bayesian dual threshold design with Dirichlet distribution: An alternative way to frequentist multi-stage phase II designs in oncology

Many phase II clinical studies in oncology use two-stage frequentist design such as Simon's optimal design. However, they have a common logistical problem regarding the patient accrual at the interim. Strictly speaking, patient accrual at the end of the first stage may have to be suspended until all patients have events, success or failure. For example, when the study endpoint is six-month progression free survival, patient accrual has to be stopped until all outcomes from stage I is observed. However, study investigators may have concern when accrual is suspended after the first stage due to the loss of accrual momentum during this hiatus. We propose a two-stage phase II design that resolves the patient accrual problem due to an interim analysis, and it can be used as an alternative way to frequentist two-stage phase II studies in oncology. ^

Changing epidemiology of trauma deaths leads to a bimodal distribution

Introduction. Injury mortality was classically described with a tri-modal distribution, with immediate deaths at the scene, early deaths due to hemorrhage, and late deaths from organ failure. We hypothesized that trauma systems development have improved pre-hospital care, early resuscitation, and critical care, and altered this pattern. ^ Methods. This is a population-based study of all trauma deaths in an urban county with a mature trauma system (n=678, median age 33 years, 81% male, 43% gunshot, 20% motor vehicle crashes). Deaths were classified as immediate (scene), early (in hospital, ≤ 4 hours from injury), or late (>4 hours post injury). Multinomial regression was used to identify independent predictors of immediate and early vs. late deaths, adjusted for age, gender, race, intention, mechanism, toxicology and cause of death. ^ Results. There were 416 (61%) immediate, 199 (29%) early, and 63 (10%) late deaths. Immediate deaths remained unchanged and early deaths occurred much earlier (median 52 minutes vs. 120). However, unlike the classic trimodal distribution, there was no late peak. Intentional injuries, alcohol intoxication, asphyxia, and injuries to the head and chest were independent predictors of immediate deaths. Alcohol intoxication and injuries to the chest were predictors of early deaths, while pelvic fractures and blunt assaults were associated with late deaths. ^ Conclusion. Trauma deaths now have a bimodal distribution. Elimination of the late peak likely represents advancements in resuscitation and critical care that have reduced organ failure. Further reductions in mortality will likely come from prevention of intentional injuries, and injuries associated with alcohol intoxication. ^

A community-based assessment of the spatiotemporal distribution of influenza in a U.S.-Mexico border county during the spring 2009 novel H1N1 swine-origin influenza A pandemic

This study retrospectively evaluated the spatial and temporal disease patterns associated with influenza-like illness (ILI), positive rapid influenza antigen detection tests (RIDT), and confirmed H1N1 S-OIV cases reported to the Cameron County Department of Health and Human Services between April 26 and May 13, 2009 using the space-time permutation scan statistic software SaTScan in conjunction with geographical information system (GIS) software ArcGIS 9.3. The rate and age-adjusted relative risk of each influenza measure was calculated and a cluster analysis was conducted to determine the geographic regions with statistically higher incidence of disease. A Poisson distribution model was developed to identify the effect that socioeconomic status, population density, and certain population attributes of a census block-group had on that area's frequency of S-OIV confirmed cases over the entire outbreak. Predominant among the spatiotemporal analyses of ILI, RIDT and S-OIV cases in Cameron County is the consistent pattern of a high concentration of cases along the southern border with Mexico. These findings in conjunction with the slight northward space-time shifts of ILI and RIDT cluster centers highlight the southern border as the primary site for public health interventions. Finally, the community-based multiple regression model revealed that three factors—percentage of the population under age 15, average household size, and the number of high school graduates over age 25—were significantly associated with laboratory-confirmed S-OIV in the Lower Rio Grande Valley. Together, these findings underscore the need for community-based surveillance, improve our understanding of the distribution of the burden of influenza within the community, and have implications for vaccination and community outreach initiatives.^

SOCIO-ENVIRONMENTAL FACTORS ASSOCIATED WITH THE DISTRIBUTION OF DRACUNCULUS MEDINENSIS (GUINEA WORM DISEASE) AND SUGGESTED METHODS FOR CONTROL, IBARAPA, NIGERIA, 1982 (AFRICA)

A population-based cross-sectional survey of socio-environmental factors associated with the prevalence of Dracunculus medinensis (guinea worm disease) was conducted in Idere, a rural agricultural community in Ibarapa, Oyo state, Nigeria, during 1982.^ The epidemiologic data were collected by household interview of all 501 households. The environmental data were collected by analysis of water samples collected from all domestic water sources and rainfall records.^ The specific objectives of this research were to: (a) Describe the prevalence of guinea worm disease in Idere during 1982 by age, sex, area of residence, drinking water source, religion and weekly amount of money spent by the household to collect potable drinking water. (b) Compare the characteristics of cases with non-cases of guinea worm in order to identify factors associated with high risk of infection. (c) Investigate domestic water sources for the distribution of Cyclops. (d) Determine the extent of potable water shortage with a view to identifying factors responsible for such shortage in the community. (e) Describe the effects of guinea worm on school attendance during 1980/1982 school years by class and location of school from piped water supply.^ The findings of this research indicate that during 1982, 31.8 percent of Idere's 6,527 residents experienced guinea worm infection, with higher prevalence of infection recorded in males in their most productive years and females in their teenage years. The role of sex and age to risk of higher infection rate was explained in the context of water related exposure and water intake due to dehydration from physical occupational actitives of subgroups.^ Potable water available to residents was considerably below the minimum recommended by WHO for tropical climates, with sixty-eight percent of water needs of the residents coming from unprotected surface water which harbour Cyclops, the obligatory intermediate host of Dracunculus medinensis. An association was found between periods of relative high density of Cyclops in domestic water and rainfall.^ Impact of guinea worm infection on educational activities was considerable and its implications were discussed, including the implications of the research findings in relation to control of guinea worm disease in Ibarapa. ^

Comparison of Tumor Shrinkage and Cumulative Dose Distribution for Lung Cancers

Background: The physical characteristic of protons is that they deliver most of their radiation dose to the target volume and deliver no dose to the normal tissue distal to the tumor. Previously, numerous studies have shown unique advantages of proton therapy over intensity-modulated radiation therapy (IMRT) in conforming dose to the tumor and sparing dose to the surrounding normal tissues and the critical structures in many clinical sites. However, proton therapy is known to be more sensitive to treatment uncertainties such as inter- and intra-fractional variations in patient anatomy. To date, no study has clearly demonstrated the effectiveness of proton therapy compared with the conventional IMRT under the consideration of both respiratory motion and tumor shrinkage in non-small cell lung cancer (NSCLC) patients. Purpose: This thesis investigated two questions for establishing a clinically relevant comparison of the two different modalities (IMRT and proton therapy). The first question was whether or not there are any differences in tumor shrinkage between patients randomized to IMRT versus passively scattered proton therapy (PSPT). Tumor shrinkage is considered a standard measure of radiation therapy response that has been widely used to gauge a short-term progression of radiation therapy. The second question was whether or not there are any differences between the planned dose and 5D dose under the influence of inter- and intra-fractional variations in the patient anatomy for both modalities. Methods: A total of 45 patients (25 IMRT patients and 20 PSPT patients) were used to quantify the tumor shrinkage in terms of the change of the primary gross tumor volume (GTVp). All patients were randomized to receive either IMRT or PSPT for NSCLC. Treatment planning goals were identical for both groups. All patients received 5 to 8 weekly repeated 4-dimensional computed tomography (4DCT) scans during the course of radiation treatments. The original GTVp contours were propagated to T50 of weekly 4DCT images using deformable image registration and their absolute volumes were measured. Statistical analysis was performed to compare the distribution of tumor shrinkage between the two population groups. In order to investigate the difference between the planned dose and the 5D dose with consideration of both breathing motion and anatomical change, we re-calculated new dose distributions at every phase of the breathing cycle for all available weekly 4DCT data sets which resulted 50 to 80 individual dose calculations for each of the 7 patients presented in this thesis. The newly calculated dose distributions were then deformed and accumulated to T50 of the planning 4DCT for comparison with the planned dose distribution. Results: At the end of the treatment, both IMRT and PSPT groups showed mean tumor volume reductions of 23.6% ( 19.2%) and 20.9% ( 17.0 %) respectively. Moreover, the mean difference in tumor shrinkage between two groups is 3% along with the corresponding 95% confidence interval, [-8%, 14%]. The rate of tumor shrinkage was highly correlated with the initial tumor volume size. For the planning dose and 5D dose comparison study, all 7 patients showed a mean difference of 1 % in terms of target coverage for both IMRT and PSPT treatment plans. Conclusions: The results of the tumor shrinkage investigation showed no statistically significant difference in tumor shrinkage between the IMRT and PSPT patients, and the tumor shrinkage between the two modalities is similar based on the 95% confidence interval. From the pilot study of comparing the planned dose with the 5D dose, we found the difference to be only 1%. Overall impression of the two modalities in terms of treatment response as measured by the tumor shrinkage and 5D dose under the influence of anatomical change that were designed under the same protocol (i.e. randomized trial) showed similar result.