Development of automated image analysis tools for verification of radiotherapy field accuracy with an electronic portal imaging device

The successful management of cancer with radiation relies on the accurate deposition of a prescribed dose to a prescribed anatomical volume within the patient. Treatment set-up errors are inevitable because the alignment of field shaping devices with the patient must be repeated daily up to eighty times during the course of a fractionated radiotherapy treatment. With the invention of electronic portal imaging devices (EPIDs), patient's portal images can be visualized daily in real-time after only a small fraction of the radiation dose has been delivered to each treatment field. However, the accuracy of human visual evaluation of low-contrast portal images has been found to be inadequate. The goal of this research is to develop automated image analysis tools to detect both treatment field shape errors and patient anatomy placement errors with an EPID. A moments method has been developed to align treatment field images to compensate for lack of repositioning precision of the image detector. A figure of merit has also been established to verify the shape and rotation of the treatment fields. Following proper alignment of treatment field boundaries, a cross-correlation method has been developed to detect shifts of the patient's anatomy relative to the treatment field boundary. Phantom studies showed that the moments method aligned the radiation fields to within 0.5mm of translation and 0.5$\sp\circ$ of rotation and that the cross-correlation method aligned anatomical structures inside the radiation field to within 1 mm of translation and 1$\sp\circ$ of rotation. A new procedure of generating and using digitally reconstructed radiographs (DRRs) at megavoltage energies as reference images was also investigated. The procedure allowed a direct comparison between a designed treatment portal and the actual patient setup positions detected by an EPID. Phantom studies confirmed the feasibility of the methodology. Both the moments method and the cross-correlation technique were implemented within an experimental radiotherapy picture archival and communication system (RT-PACS) and were used clinically to evaluate the setup variability of two groups of cancer patients treated with and without an alpha-cradle immobilization aid. The tools developed in this project have proven to be very effective and have played an important role in detecting patient alignment errors and field-shape errors in treatment fields formed by a multileaf collimator (MLC). ^

HUMAN PROPHASE CHROMOSOME UNIQUE BAND SEQUENCES; DEFINITION, CHARACTERIZATION, AND APPLICATION (CYTOGENETICS, IMAGE PROCESSING)

Extensive experience with the analysis of human prophase chromosomes and studies into the complexity of prophase GTG-banding patterns have suggested that at least some prophase chromosomal segments can be accurately identified and characterized independently of the morphology of the chromosome as a whole. In this dissertation the feasibility of identifying and analyzing specified prophase chromosome segments was thus investigated as an alternative approach to prophase chromosome analysis based on whole chromosome recognition. Through the use of prophase idiograms at the 850-band-stage (FRANCKE, 1981) and a comparison system based on the calculation of cross-correlation coefficients between idiogram profiles, we have demonstrated that it is possible to divide the 24 human prophase idiograms into a set of 94 unique band sequences. Each unique band sequence has a banding pattern that is recognizable and distinct from any other non-homologous chromosome portion.^ Using chromosomes 11p and 16 thru 22 to demonstrate unique band sequence integrity at the chromosome level, we found that prophase chromosome banding pattern variation can be compensated for and that a set of unique band sequences very similar to those at the idiogram level can be identified on actual chromosomes.^ The use of a unique band sequence approach in prophase chromosome analysis is expected to increase efficiency and sensitivity through more effective use of available banding information. The use of a unique band sequence approach to prophase chromosome analysis is discussed both at the routine level by cytogeneticists and at an image processing level with a semi-automated approach to prophase chromosome analysis. ^

Intensive care noise and mean arterial blood pressure in ELBW neonates

Objective. Loud noises in neonatal intensive care units (NICUs) may impede growth and development for extremely low birthweight (ELBW, < 1000 grams) newborns. The objective of this study was to measure the association between NICU sound levels and ELBW neonates' arterial blood pressure to determine whether these newborns experience noise-induced stress. ^ Methods. Noise and arterial blood pressure recordings were collected for 9 ELBW neonates during the first week of life. Sound levels were measured inside the incubator, and each subject's arterial blood pressures were simultaneously recorded for 15 minutes (at 1 sec intervals). Time series cross-correlation functions were calculated for NICU noise and mean arterial blood pressure (MABP) recordings for each subject. The grand mean noise-MABP cross-correlation was calculated for all subjects and for lower and higher birthweight groups for comparison. ^ Results. The grand mean noise-MABP cross-correlation for all subjects was mostly negative (through 300 sec lag time) and nearly reached significance at the 95% level at 111 sec lag (mean r = -0.062). Lower birthweight newborns (454-709 g) experienced significant decreases in blood pressure with increasing NICU noise after 145 sec lag (peak r = -0.074). Higher birthweight newborns had an immediate negative correlation with NICU sound levels (at 3 sec lag, r = -0.071), but arterial blood pressures increased to a positive correlation with noise levels at 197 sec lag (r = 0.075). ^ Conclusions. ELBW newborns' arterial blood pressure was influenced by NICU noise levels during the first week of life. Lower birthweight newborns may have experienced an orienting reflex to NICU sounds. Higher birthweight newborns experienced an immediate orienting reflex to increasing sound levels, but arterial blood pressure increased approximately 3 minutes after increases in noise levels. Increases in arterial blood pressure following increased NICU sound levels may result from a stress response to noise. ^

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

Correlates of tobacco cessation counseling among Hispanic physicians in the US: a cross-sectional survey study.

BACKGROUND: Physician advice is an important motivator for attempting to stop smoking. However, physicians' lack of intervention with smokers has only modestly improved in the last decade. Although the literature includes extensive research in the area of the smoking intervention practices of clinicians, few studies have focused on Hispanic physicians. The purpose of this study was to explore the correlates of tobacco cessation counseling practices among Hispanic physicians in the US. METHODS: Data were collected through a validated survey instrument among a cross-sectional sample of self-reported Hispanic physicians practicing in New Mexico, and who were members of the New Mexico Hispanic Medical Society in the year 2001. Domains of interest included counseling practices, self-efficacy, attitudes/responsibility, and knowledge/skills. Returned surveys were analyzed to obtain frequencies and descriptive statistics for each survey item. Other analyses included: bivariate Pearson's correlation, factorial ANOVAs, and multiple linear regressions. RESULTS: Respondents (n = 45) reported a low level of compliance with tobacco control guidelines and recommendations. Results indicate that physicians' familiarity with standard cessation protocols has a significant effect on their tobacco-related practices (r = .35, variance shared = 12%). Self-efficacy and gender were both significantly correlated to tobacco related practices (r = .42, variance shared = 17%). A significant correlation was also found between self-efficacy and knowledge/skills (r = .60, variance shared = 36%). Attitudes/responsibility was not significantly correlated with any of the other measures. CONCLUSION: More resources should be dedicated to training Hispanic physicians in tobacco intervention. Training may facilitate practice by increasing knowledge, developing skills and, ultimately, enhancing feelings of self-efficacy.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

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

A novel method for evaluating an interaction effect of correlated continuous covariates in a linear model

Interaction effect is an important scientific interest for many areas of research. Common approach for investigating the interaction effect of two continuous covariates on a response variable is through a cross-product term in multiple linear regression. In epidemiological studies, the two-way analysis of variance (ANOVA) type of method has also been utilized to examine the interaction effect by replacing the continuous covariates with their discretized levels. However, the implications of model assumptions of either approach have not been examined and the statistical validation has only focused on the general method, not specifically for the interaction effect.^ In this dissertation, we investigated the validity of both approaches based on the mathematical assumptions for non-skewed data. We showed that linear regression may not be an appropriate model when the interaction effect exists because it implies a highly skewed distribution for the response variable. We also showed that the normality and constant variance assumptions required by ANOVA are not satisfied in the model where the continuous covariates are replaced with their discretized levels. Therefore, naïve application of ANOVA method may lead to an incorrect conclusion. ^ Given the problems identified above, we proposed a novel method modifying from the traditional ANOVA approach to rigorously evaluate the interaction effect. The analytical expression of the interaction effect was derived based on the conditional distribution of the response variable given the discretized continuous covariates. A testing procedure that combines the p-values from each level of the discretized covariates was developed to test the overall significance of the interaction effect. According to the simulation study, the proposed method is more powerful then the least squares regression and the ANOVA method in detecting the interaction effect when data comes from a trivariate normal distribution. The proposed method was applied to a dataset from the National Institute of Neurological Disorders and Stroke (NINDS) tissue plasminogen activator (t-PA) stroke trial, and baseline age-by-weight interaction effect was found significant in predicting the change from baseline in NIHSS at Month-3 among patients received t-PA therapy.^

Selection bias and the cross-validation of regression models for prediction

Strategies are compared for the development of a linear regression model with stochastic (multivariate normal) regressor variables and the subsequent assessment of its predictive ability. Bias and mean squared error of four estimators of predictive performance are evaluated in simulated samples of 32 population correlation matrices. Models including all of the available predictors are compared with those obtained using selected subsets. The subset selection procedures investigated include two stopping rules, C$\sb{\rm p}$ and S$\sb{\rm p}$, each combined with an 'all possible subsets' or 'forward selection' of variables. The estimators of performance utilized include parametric (MSEP$\sb{\rm m}$) and non-parametric (PRESS) assessments in the entire sample, and two data splitting estimates restricted to a random or balanced (Snee's DUPLEX) 'validation' half sample. The simulations were performed as a designed experiment, with population correlation matrices representing a broad range of data structures.^ The techniques examined for subset selection do not generally result in improved predictions relative to the full model. Approaches using 'forward selection' result in slightly smaller prediction errors and less biased estimators of predictive accuracy than 'all possible subsets' approaches but no differences are detected between the performances of C$\sb{\rm p}$ and S$\sb{\rm p}$. In every case, prediction errors of models obtained by subset selection in either of the half splits exceed those obtained using all predictors and the entire sample.^ Only the random split estimator is conditionally (on $\\beta$) unbiased, however MSEP$\sb{\rm m}$ is unbiased on average and PRESS is nearly so in unselected (fixed form) models. When subset selection techniques are used, MSEP$\sb{\rm m}$ and PRESS always underestimate prediction errors, by as much as 27 percent (on average) in small samples. Despite their bias, the mean squared errors (MSE) of these estimators are at least 30 percent less than that of the unbiased random split estimator. The DUPLEX split estimator suffers from large MSE as well as bias, and seems of little value within the context of stochastic regressor variables.^ To maximize predictive accuracy while retaining a reliable estimate of that accuracy, it is recommended that the entire sample be used for model development, and a leave-one-out statistic (e.g. PRESS) be used for assessment. ^

Motivational interviewing and feedback for college drinkers: Handling missing values with complete case analysis and multiple imputation

Objective: In this secondary data analysis, three statistical methodologies were implemented to handle cases with missing data in a motivational interviewing and feedback study. The aim was to evaluate the impact that these methodologies have on the data analysis. ^ Methods: We first evaluated whether the assumption of missing completely at random held for this study. We then proceeded to conduct a secondary data analysis using a mixed linear model to handle missing data with three methodologies (a) complete case analysis, (b) multiple imputation with explicit model containing outcome variables, time, and the interaction of time and treatment, and (c) multiple imputation with explicit model containing outcome variables, time, the interaction of time and treatment, and additional covariates (e.g., age, gender, smoke, years in school, marital status, housing, race/ethnicity, and if participants play on athletic team). Several comparisons were conducted including the following ones: 1) the motivation interviewing with feedback group (MIF) vs. the assessment only group (AO), the motivation interviewing group (MIO) vs. AO, and the intervention of the feedback only group (FBO) vs. AO, 2) MIF vs. FBO, and 3) MIF vs. MIO.^ Results: We first evaluated the patterns of missingness in this study, which indicated that about 13% of participants showed monotone missing patterns, and about 3.5% showed non-monotone missing patterns. Then we evaluated the assumption of missing completely at random by Little's missing completely at random (MCAR) test, in which the Chi-Square test statistic was 167.8 with 125 degrees of freedom, and its associated p-value was p=0.006, which indicated that the data could not be assumed to be missing completely at random. After that, we compared if the three different strategies reached the same results. For the comparison between MIF and AO as well as the comparison between MIF and FBO, only the multiple imputation with additional covariates by uncongenial and congenial models reached different results. For the comparison between MIF and MIO, all the methodologies for handling missing values obtained different results. ^ Discussions: The study indicated that, first, missingness was crucial in this study. Second, to understand the assumptions of the model was important since we could not identify if the data were missing at random or missing not at random. Therefore, future researches should focus on exploring more sensitivity analyses under missing not at random assumption.^

Rabies control in the developing world: Results of a cross-sectional study in Bohol, Philippines

Rabies remains a significant problem in much of the developed world, where canine rabies is not well controlled, and the bite of an infected dog is the most common means of transmission. The Philippines continues to report several hundred cases of human rabies every year, and many more cases go undetected. In recent years, the province of Bohol has been targeted by the Philippine government and the World Health Organization for a rabies eradication program. ^ The primary objective of this dissertation research was to describe factors associated with dog vaccination coverage and knowledge, attitudes, and practices regarding rabies among households in Bohol, Philippines. Utilizing a cross-sectional cluster survey design, we sampled 460 households and 541 dogs residing within dog-owning households. ^ Multivariate linear regression was used to examine potential associations between knowledge, attitudes, and practices (KAPs) and variables of interest. Forty-six percent of households knew that rabies was spread through the bite of an infected dog. The mean knowledge score was 8.36 (SD: ± 3.4; range: 1–24). We found that having known someone with rabies was significantly associated with an almost one point increase in the knowledge score (β = 0.88; p = 0.02). The mean attitudes score was 5.65 (SD: ± 0.63; range: 2–6), and the mean practices score was 7.07 (SD: ± 1.7; range: 2–9). Both the attitudes score and the practices score were positively and significantly associated with only the knowledge score and no other covariates. ^ Multivariate logistic regression was used to examine associations between dog vaccination coverage and variables of interest. Approximately 71% of owned dogs in Bohol were reported as vaccinated at some time during their lives. We found that a dog's age was significantly associated with vaccination, and the odds of vaccination increased in a linear fashion with age. We also found that dogs had approximately twice the odds of being vaccinated if they were confined both day and night to the household premises or if the owner was employed; however, these results were only marginally significant (p = 0.07) in the multivariate model. ^ Finally, a systematic review was conducted on canine rabies vaccination and dog population demographics in the developing world. We found few studies on this topic, especially in countries where the burden of rabies is greatest. Overall, dog ownership is high. Dogs are quite young and do not live very long due to disease and accidents. The biggest deterrent to vaccination is the rapid dog population turnover. ^ It is our hope that this work will be used to improve dog rabies vaccination programs around the world and save lives, both human and canine.^

Robust effect sizes and their confidence intervals for group difference between trajectories in hierarchical linear growth model

Hierarchical linear growth model (HLGM), as a flexible and powerful analytic method, has played an increased important role in psychology, public health and medical sciences in recent decades. Mostly, researchers who conduct HLGM are interested in the treatment effect on individual trajectories, which can be indicated by the cross-level interaction effects. However, the statistical hypothesis test for the effect of cross-level interaction in HLGM only show us whether there is a significant group difference in the average rate of change, rate of acceleration or higher polynomial effect; it fails to convey information about the magnitude of the difference between the group trajectories at specific time point. Thus, reporting and interpreting effect sizes have been increased emphases in HLGM in recent years, due to the limitations and increased criticisms for statistical hypothesis testing. However, most researchers fail to report these model-implied effect sizes for group trajectories comparison and their corresponding confidence intervals in HLGM analysis, since lack of appropriate and standard functions to estimate effect sizes associated with the model-implied difference between grouping trajectories in HLGM, and also lack of computing packages in the popular statistical software to automatically calculate them. ^ The present project is the first to establish the appropriate computing functions to assess the standard difference between grouping trajectories in HLGM. We proposed the two functions to estimate effect sizes on model-based grouping trajectories difference at specific time, we also suggested the robust effect sizes to reduce the bias of estimated effect sizes. Then, we applied the proposed functions to estimate the population effect sizes (d ) and robust effect sizes (du) on the cross-level interaction in HLGM by using the three simulated datasets, and also we compared the three methods of constructing confidence intervals around d and du recommended the best one for application. At the end, we constructed 95% confidence intervals with the suitable method for the effect sizes what we obtained with the three simulated datasets. ^ The effect sizes between grouping trajectories for the three simulated longitudinal datasets indicated that even though the statistical hypothesis test shows no significant difference between grouping trajectories, effect sizes between these grouping trajectories can still be large at some time points. Therefore, effect sizes between grouping trajectories in HLGM analysis provide us additional and meaningful information to assess group effect on individual trajectories. In addition, we also compared the three methods to construct 95% confident intervals around corresponding effect sizes in this project, which handled with the uncertainty of effect sizes to population parameter. We suggested the noncentral t-distribution based method when the assumptions held, and the bootstrap bias-corrected and accelerated method when the assumptions are not met.^