A Bootstrap-Regression Procedure to Capture Unit Specific Effects in Data Envelopment Analysis

The Data Envelopment Analysis (DEA) efficiency score obtained for an individual firm is a point estimate without any confidence interval around it. In recent years, researchers have resorted to bootstrapping in order to generate empirical distributions of efficiency scores. This procedure assumes that all firms have the same probability of getting an efficiency score from any specified interval within the [0,1] range. We propose a bootstrap procedure that empirically generates the conditional distribution of efficiency for each individual firm given systematic factors that influence its efficiency. Instead of resampling directly from the pooled DEA scores, we first regress these scores on a set of explanatory variables not included at the DEA stage and bootstrap the residuals from this regression. These pseudo-efficiency scores incorporate the systematic effects of unit-specific factors along with the contribution of the randomly drawn residual. Data from the U.S. airline industry are utilized in an empirical application.

A coalescent analysis for modeling the mutation process in colorectal cancer

Colorectal cancer is the forth most common diagnosed cancer in the United States. Every year about a hundred forty-seven thousand people will be diagnosed with colorectal cancer and fifty-six thousand people lose their lives due to this disease. Most of the hereditary nonpolyposis colorectal cancer (HNPCC) and 12% of the sporadic colorectal cancer show microsatellite instability. Colorectal cancer is a multistep progressive disease. It starts from a mutation in a normal colorectal cell and grows into a clone of cells that further accumulates mutations and finally develops into a malignant tumor. In terms of molecular evolution, the process of colorectal tumor progression represents the acquisition of sequential mutations. ^ Clinical studies use biomarkers such as microsatellite or single nucleotide polymorphisms (SNPs) to study mutation frequencies in colorectal cancer. Microsatellite data obtained from single genome equivalent PCR or small pool PCR can be used to infer tumor progression. Since tumor progression is similar to population evolution, we used an approach known as coalescent, which is well established in population genetics, to analyze this type of data. Coalescent theory has been known to infer the sample's evolutionary path through the analysis of microsatellite data. ^ The simulation results indicate that the constant population size pattern and the rapid tumor growth pattern have different genetic polymorphic patterns. The simulation results were compared with experimental data collected from HNPCC patients. The preliminary result shows the mutation rate in 6 HNPCC patients range from 0.001 to 0.01. The patients' polymorphic patterns are similar to the constant population size pattern which implies the tumor progression is through multilineage persistence instead of clonal sequential evolution. The results should be further verified using a larger dataset. ^

A cross-sectional analysis of physical activity and obesity in inner city Houston Texas Mexican American youth

Obesity prevalence among children and adolescents is rising. It is one of the most attributable causes of hospitalization and death. Overweight and obese children are more likely to suffer from associated conditions such as hypertension, dyslipidemia, chronic inflammation, increased blood clotting tendency, endothelial dysfunction, hyperinsulinemia, and asthma. These children and adolescents are also more likely to be overweight and obese in adulthood. Interestingly, rates of obesity and overweight are not evenly distributed across racial and ethnic groups. Mexican American youth have higher rates of obesity and are at higher risk of becoming obese than non-Hispanic black and non-Hispanic white children. ^ Methods. This cross-sectional study describes the association between rates of obesity and physical activity in a sample of 1313 inner-city Mexican American children and adolescents (5-19 years of age) in Houston, Texas. This study is important because it will contribute to our understanding of childhood and adolescent obesity in this at-risk population. ^ Data from the Mexican American Feasibility Cohort using the Mano a Mano questionnaire are used to describe this population's status of obesity and physical activity. An initial sample taken from 5000 households in inner city Houston Texas was used as the baseline for this prospective cohort. The questionnaire was given in person to the participants to complete (or to parents for younger children) at a home visit by two specially trained bilingual interviewers. Analysis comprised prevalence estimates of obesity represented as percentile rank (<85%= normal weight, >85%= at risk, >95%= obese) by age and gender. The association between light, moderate, strenuous activity, and obesity was also examined using linear regression. ^ Results. Overall, 46% of this Mexican American Feasibility cohort is overweight or obese. The prevalence for children in the 6-11 age range (53.2%) was significantly greater than that reported from NHANES, 1999–2002 data (39.4%). Although the percentage of overweight and obese among the 12-19 year olds was greater than that reported in NHANES (38.5% versus 38.6%) this difference was not statistically significant. ^ A significant association between BMI and sit time and moderate physical activity (both p < 0.05) found in this sample. For males, this association was significant for moderate physical activity (p < 0.01). For the females, this association was significant for BMI and sit time (p < 0.05). These results need to be interpreted in the light of design and measurement limitations. ^ Conclusion. This study supports observations that the inner city Houston Texas Mexican American child and adolescent population is more overweight and obese than nationally reported figures, and that there are positive relationships between BMI, activity levels, and sit time in this population. This study supports the need for public health initiatives within the Houston Hispanic community. ^

Functional analysis of BLAP75, a BLM-associated protein

Bloom syndrome (BS) is an autosomal recessive disorder characterized by dwarfism, immunodeficiency, impaired fertility, and most importantly, early development of a broad range of cancers. The hallmark of BS cells is hyper-recombination, characterized by a drastically elevated frequency of sister chromatid exchange (SCE). BLM, the gene mutated in BS, encodes a DNA helicase of the RecQ protein family. BLM is thought to participate in several DNA transactions and to interact with many proteins involved in DNA replication, recombination, and repair. However, the precise function of BLM and the BLM-dependent anti-tumor mechanism remain obscure. ^ A novel protein, BLAP75 (BLM-associated polypeptide, 75KD), was identified to form an evolutionarily conserved complex with BLM and DNA topoisomerase IIIα (Topo IIIα). Our work demonstrates that loss of BLAP75 destabilized BLM and Topo IIIα proteins. BLAP75 colocalized with BLM in subnuclear foci in response to DNA damage and the recruitment of BLM to these foci was BLAP75-dependent. Moreover, depletion of BLAP75 by siRNA resulted in an elevated SCE rate similar to cells depleted of BLM by siRNA. In addition, RNAi-mediated silencing of BLAP75 greatly diminished cell viability. This cellular deficiency was rescued by expression of wild type BLAP75 but not BLAP75 with mutated conserved domain III, which abrogated the interaction between BLAP75, BLM and Topo IIIα, suggesting that the integrity of BLM-Topo IIIα-BLAP75 complex might be critical for cell survival. Finally, I found that BLAP75 was phosphorylated during mitosis and upon various DNA-damaging agents, implying that BLAP75 might also function in mitosis and DNA damage response. ^ Taken together, this study has defined BLAP75 as an integral component of the BLM complex to maintain genome stability. Our findings provide insights into the molecular mechanisms of the BLM helicase pathway and tumorigenesis process associated with these mechanisms. ^

Health care policy analysis of Medicaid's early and periodic screening, diagnosis, and treatment (EPSDT) benefit in Colorado

The purpose of this thesis is to identify "best practice" recommendations for successful implementation of the EPSDT outreach program at Memorial Health System's Hospital for Children in Colorado Springs through a policy analysis of Medicaid EPSDT services in Colorado. A successful program at Memorial will increase education and awareness of EPSDT services, enrollment, and access to and utilization of health care services for eligible children. Methodology utilized in this study included questionnaires designed for the EPSDT contract administrator and outreach coordinators/workers; analysis of current federal and state policies; and studies conducted at the federal and state level, and by various advocacy groups. The need for this analysis of EPSDT came about in part through an awareness of increasingly high numbers of children in poverty and who are uninsured. Though the percentage of children living in poverty in Colorado is slightly below the national average (see Table 2), according to data analyzed by The Annie E. Casey Foundation, the percentage of children (0-18) living in poverty in Colorado increased from 10% in 2000 to 16% in 2006, a dramatic increase of 60% surpassed by only one other state in the nation (The Annie E. Casey Foundation, 2008). By comparison, the U.S. percentage of children in poverty during the same time frame rose from 17% to 18% (The Annie E. Casey Foundation, 2008). What kind of health care services are available to this vulnerable and growing group of Coloradans, and what are the barriers that affect their enrollment in, access to and utilization of these health care services? Barriers identified included difficulty with the application process; system and process issues; a lack of providers; and a lack of awareness and knowledge of EPSDT. Fiscal restraints and legislation at the federal and state level are also barriers to increasing enrollment and access to services. Outreach services are a critical component of providing EPSDT services, and there were several recommendations regarding outreach and case management that will benefit the program in the future. Through this analysis and identification of a broad range of barriers, a clearer picture emerged of current challenges within the EPSDT program as well as a broad range of strategies and recommendations to address these challenges. Through increased education and advocacy for EPSDT and the services it encompasses; stronger collaboration and cooperation between all groups involved, including providing a Medical Home for all eligible children; and new legislation putting more money and focus on comprehensive health care for low-income uninsured children; enrollment, access to and utilization of developmentally appropriate and quality health care services can be achieved. ^

Interim analysis of clinical trials: Simulation studies of fractional Brownian motion

Interim clinical trial monitoring procedures were motivated by ethical and economic considerations. Classical Brownian motion (Bm) techniques for statistical monitoring of clinical trials were widely used. Conditional power argument and α-spending function based boundary crossing probabilities are popular statistical hypothesis testing procedures under the assumption of Brownian motion. However, it is not rare that the assumptions of Brownian motion are only partially met for trial data. Therefore, I used a more generalized form of stochastic process, called fractional Brownian motion (fBm), to model the test statistics. Fractional Brownian motion does not hold Markov property and future observations depend not only on the present observations but also on the past ones. In this dissertation, we simulated a wide range of fBm data, e.g., H = 0.5 (that is, classical Bm) vs. 0.5< H <1, with treatment effects vs. without treatment effects. Then the performance of conditional power and boundary-crossing based interim analyses were compared by assuming that the data follow Bm or fBm. Our simulation study suggested that the conditional power or boundaries under fBm assumptions are generally higher than those under Bm assumptions when H > 0.5 and also matches better with the empirical results. ^

ANALYSIS OF ESTROGEN-INDUCED ANEUPLOIDY AND CHROMOSOME DISPLACEMENT

Diethylstilbestrol (DES) is a known human carcinogen and teratogen whose mechanism of action remains undetermined. As essentially diploid Chinese hamster cell line (Don) was used to test diethylstilbestrol (DES), dienestrol, hexestrol and the naturally occurring estrogens, estradiol and estriol for their ability to cause metaphase arrest and to induce aneuploidy. These compounds arrest mitosis within a narrow range of high concentrations and induce aneuploidy in recovering cell populations. DES was the most effective arrestant on a comparative molar basis. Estradiol and estriol were less potent as arrestants but were effective inducers of aneuploidy. Aneuploidy was induced in a non-random manner. The smallest chromosomes were most frequently recorded in aneuploid cells. Using anti-tubulin antibody and indirect immunofluorescence, it was found that DES inhibits bi-polar spindle assembly and disrupts the cytoplasmic microtubule complex (CMTC). Estradiol arrests mitosis in a manner that allows spindle assembly. Estradiol has no apparent effect on the CMTC. The naturally occurring estrogens caused chromosome displacement during mitotic arrest. Electron microscopy confirmed that the displaced chromosomes appeared at the polar regions of arrested cells. The arresting effect of estradiol, and to some extent DES, was reduced by the addition of dibutyryl cyclic adenosine monophosphate (db-cAMP). Aneuploidy induction by DES and similar compounds may be related to their carcinogenic and/or teratogenic potential. ^

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

DEVELOPMENT OF A BEAM-SPECIFIC PLANNING TARGET VOLUME AND A ROBUST PLAN ANALYSIS TOOL FOR PROTON THERAPY

Proton therapy is growing increasingly popular due to its superior dose characteristics compared to conventional photon therapy. Protons travel a finite range in the patient body and stop, thereby delivering no dose beyond their range. However, because the range of a proton beam is heavily dependent on the tissue density along its beam path, uncertainties in patient setup position and inherent range calculation can degrade thedose distribution significantly. Despite these challenges that are unique to proton therapy, current management of the uncertainties during treatment planning of proton therapy has been similar to that of conventional photon therapy. The goal of this dissertation research was to develop a treatment planning method and a planevaluation method that address proton-specific issues regarding setup and range uncertainties. Treatment plan designing method adapted to proton therapy: Currently, for proton therapy using a scanning beam delivery system, setup uncertainties are largely accounted for by geometrically expanding a clinical target volume (CTV) to a planning target volume (PTV). However, a PTV alone cannot adequately account for range uncertainties coupled to misaligned patient anatomy in the beam path since it does not account for the change in tissue density. In order to remedy this problem, we proposed a beam-specific PTV (bsPTV) that accounts for the change in tissue density along the beam path due to the uncertainties. Our proposed method was successfully implemented, and its superiority over the conventional PTV was shown through a controlled experiment.. Furthermore, we have shown that the bsPTV concept can be incorporated into beam angle optimization for better target coverage and normal tissue sparing for a selected lung cancer patient. Treatment plan evaluation method adapted to proton therapy: The dose-volume histogram of the clinical target volume (CTV) or any other volumes of interest at the time of planning does not represent the most probable dosimetric outcome of a given plan as it does not include the uncertainties mentioned earlier. Currently, the PTV is used as a surrogate of the CTV’s worst case scenario for target dose estimation. However, because proton dose distributions are subject to change under these uncertainties, the validity of the PTV analysis method is questionable. In order to remedy this problem, we proposed the use of statistical parameters to quantify uncertainties on both the dose-volume histogram and dose distribution directly. The robust plan analysis tool was successfully implemented to compute both the expectation value and its standard deviation of dosimetric parameters of a treatment plan under the uncertainties. For 15 lung cancer patients, the proposed method was used to quantify the dosimetric difference between the nominal situation and its expected value under the uncertainties.

Grain size analysis on CRP sediment cores from the Ross Sea, Antarctica

The upper 1200 m of pre-Pliocene sediment recovered by Cape Roberts Project (CRP) drilling off the Victoria Land coast of Antarctica between 1997-1999 has been subdivided into 54 unconformity-bound stratigraphic sequences, spanning the period c. 32 to 17 Ma. The sequences are recognised on the basis of the cyclical vertical stacking of their constituent lithofacies, which are enclosed by erosion surfaces produced during the grounding of the advancing ice margin onto the sea floor. Each sequence represents deposition in a range of offshore shelf to coastal glacimarine sedimentary environments during oscillations in the ice margin across the Western Ross Sea shelf, and coeval fluctuations in water depth. This paper applies spectral analysis techniques to depth- and time-series of sediment grain size (500 samples) for intervals of the core with adequate chronological data. Time series analysis of 0.5-l.0m-spaced grainsize data spanning sequences 9-11 (CRP-2/2A) and sequences 1-7 (CRP-3) suggests that the length of individual sequences correspond to Milankovitch frequencies, probably 41 k.y., but possibly as low as 100 k.y. Higher frequency periodic components at 23 k.y. (orbital precession) and 15-10 k.y. (sub-orbital) are recognised at the intrasequence-scale, and may represent climatic cycles akin to the ice rafting episodes described in the North Atlantic Ocean during the Quaternary. The cyclicity recorded by glacimarine sequences in CRP core provides direct evidence from the periphery of Antarctica for orbital oscillations in the size of the Oligocene-Early Miocene East Antarctic Ice Sheet.

Planktonic and benthic stables Isotopes, age model and carbon analysis of HUD91/039 from the northernmost Baffin Bay

A multiproxy study of palaeoceanographic and climatic changes in northernmost Baffin Bay shows that major environmental changes have occurred since the deglaciation of the area at about 12 500 cal. yr BP. The interpretation is based on sedimentology, benthic and planktonic foraminifera and their isotopic composition, as well as diatom assemblages in the sedimentary records at two core sites, one located in the deeper central part of northernmost Baffin Bay and one in a separate trough closer to the Greenland coast. A revised chronology for the two records is established on the basis of 15 previously published AMS 14C age determinations. A basal diamicton is overlain by laminated, fossil-free sediments. Our data from the early part of the fossiliferous record (12 300 - 11 300 cal. yr BP), which is also initially laminated, indicate extensive seasonal sea-ice cover and brine release. There is indication of a cooling event between 11 300 and 10 900 cal. yr BP, and maximum Atlantic Water influence occurred between 10 900 and 8200 cal. yr BP (no sediment recovery between 8200 and 7300 cal. yr BP). A gradual, but fluctuating, increase in sea-ice cover is seen after 7300 cal. yr BP. Sea-ice diatoms were particularly abundant in the central part of northernmost Baffin Bay, presumably due to the inflow of Polar waters from the Arctic Ocean, and less sea ice occurred at the near-coastal site, which was under continuous influence of the West Greenland Current. Our data from the deep, central part show a fluctuating degree of upwelling after c. 7300 cal. yr BP, culminating between 4000 and 3050 cal. yr BP. There was a gradual increase in the influence of cold bottom waters from the Arctic Ocean after about 3050 cal. yr BP, when agglutinated foraminifera became abundant. A superimposed short-term change in the sea-surface proxies is correlated with the Little Ice Age cooling.

Grain size analysis of sediment coring site CRP-2 from the Ross Sea, Antarctica

The purpose of this note is to present results of grain size analyses from 118 samples of the CRP-2/2A core using sieve and Sedigraph techniques. The samples were selected to represent the range of facies encountered, and tend to become more widely spaced with depth. Fifteen came from the upper 27 m of Quaternary and Pliocene sediments, 62 from the early Miocene-late Oligocene strata (27 to 307 mbsf), and 41 from the early Oligocene strata beneath (307 to 624 mbsf). The results are intended to provide reference data for lithological descriptions in the core logs (Cape Roberts Science Team, 1999), and to help with facies interpretation. The analytical technique used for determining size frequency of the sand fraction in our samples (sieving) is simple, physical and widely practised for over a century. Thus it provides a useful reference point for analyses produced by other faster and more sophisticated techniques, such as the Malvern laser particle size analysis system (Woolfe et al., 2000), and estimates derived from measurements taken with down-hole logging tools (Bücker, pers. com., 1999).

Sedimentation rate calculations and isotopic analysis on sediment cores along the IODP Expedition 311 transect

The oceanographic and tectonic conditions of accretionary margins are well-suited for several potential processes governing methane generation, storage and release. To identify the relevant methane evolution pathways in the northern Cascadia accretionary margin, a four-site transect was drilled during Integrated Ocean Drilling Program Expedition 311. The d13C values of methane range from a minimum value of -82.2 per mil on an uplifted ridge of accreted sediment near the deformation front (Site U1326, 1829 mbsl, meters below sea level) to a maximum value of -39.5 per mil at the most landward location within an area of steep canyons near the shelf edge (Site U1329, 946 mbsl). An interpretation based solely on methane isotope values might conclude the 13C-enrichment of methane indicates a transition from microbially- to thermogenically-sourced methane. However, the co-existing CO2 exhibits a similar trend of 13C-enrichment along the transect with values ranging from -22.5 per mil to +25.7 per mil. The magnitude of the carbon isotope separation between methane and CO2 (Ec = 63.8 ± 5.8) is consistent with isotope fractionation during microbially mediated carbonate reduction. These results, in conjunction with a transect-wide gaseous hydrocarbon content composed of > 99.8% (by volume) methane and uniform dDCH4 values (-172 per mil ± 8) that are distinct from thermogenic methane at a seep located 60 km from the Expedition 311 transect, suggest microbial CO2 reduction is the predominant methane source at all investigated sites. The magnitude of the intra-site downhole 13C-enrichment of CO2 within the accreted ridge (Site U1326) and a slope basin nearest the deformation front (Site U1325, 2195 mbsl) is ~ 5 per mil. At the mid-slope site (Site U1327, 1304 mbsl) the downhole 13C-enrichment of the CO2 is ~ 25 per mil and increases to ~ 40 per mil at the near-shelf edge Site U1329. This isotope fractionation pattern is indicative of more extensive diagenetic alteration at sites with greater 13C-enrichment. The magnitude of the 13C-enrichment of CO2 correlates with decreasing sedimentation rates and a diminishing occurrence of stratigraphic gas hydrate. We suggest the decreasing sedimentation rates increase the exposure time of sedimentary organic matter to aerobic and anaerobic degradation, during burial, thereby reducing the availability of metabolizable organic matter available for methane production. This process is reflected in the occurrence and distribution of gas hydrate within the northern Cascadia margin accretionary prism. Our observations are relevant for evaluating methane production and the occurrence of stratigraphic gas hydrate within other convergent margins.

AMS radiocarbon dates and pollen analysis of ODP Hole 175-1078C

ODP Site 1078 situated under the coast of Angola provides the first record of the vegetation history for Angola. The upper 11 m of the core covers the past 30 thousand years, which has been analysed palynologically in decadal to centennial resolution. Alkenone sea surface temperature estimates were analysed in centennial resolution. We studied sea surface temperatures and vegetation development during full glacial, deglacial, and interglacial conditions. During the glacial the vegetation in Angola was very open consisting of grass and heath lands, deserts and semi-deserts, which suggests a cool and dry climate. A change to warmer and more humid conditions is indicated by forest expansion starting in step with the earliest temperature rise in Antarctica, 22 thousand years ago. We infer that around the period of Heinrich Event 1, a northward excursion of the Angola Benguela Front and the Congo Air Boundary resulted in cool sea surface temperatures but rain forest remained present in the northern lowlands of Angola. Rain forest and dry forest area increase 15 thousand years ago. During the Holocene, dry forests and Miombo woodlands expanded. Also in Angola globally recognised climate changes at 8 thousand and 4 thousand years ago had an impact on the vegetation. During the past 2 thousand years, savannah vegetation became dominant.

Seagrass time-series analysis products, in Moreton Bay, Australia, with links to GeoTIFFs

The spatial and temporal dynamics of seagrasses have been well studied at the leaf to patch scales, however, the link to large spatial extent landscape and population dynamics is still unresolved in seagrass ecology. Traditional remote sensing approaches have lacked the temporal resolution and consistency to appropriately address this issue. This study uses two high temporal resolution time-series of thematic seagrass cover maps to examine the spatial and temporal dynamics of seagrass at both an inter- and intra-annual time scales, one of the first globally to do so at this scale. Previous work by the authors developed an object-based approach to map seagrass cover level distribution from a long term archive of Landsat TM and ETM+ images on the Eastern Banks (~200 km**2), Moreton Bay, Australia. In this work a range of trend and time-series analysis methods are demonstrated for a time-series of 23 annual maps from 1988 to 2010 and a time-series of 16 monthly maps during 2008-2010. Significant new insight was presented regarding the inter- and intra-annual dynamics of seagrass persistence over time, seagrass cover level variability, seagrass cover level trajectory, and change in area of seagrass and cover levels over time. Overall we found that there was no significant decline in total seagrass area on the Eastern Banks, but there was a significant decline in seagrass cover level condition. A case study of two smaller communities within the Eastern Banks that experienced a decline in both overall seagrass area and condition are examined in detail, highlighting possible differences in environmental and process drivers. We demonstrate how trend and time-series analysis enabled seagrass distribution to be appropriately assessed in context of its spatial and temporal history and provides the ability to not only quantify change, but also describe the type of change. We also demonstrate the potential use of time-series analysis products to investigate seagrass growth and decline as well as the processes that drive it. This study demonstrates clear benefits over traditional seagrass mapping and monitoring approaches, and provides a proof of concept for the use of trend and time-series analysis of remotely sensed seagrass products to benefit current endeavours in seagrass ecology.