Quality Data is Key to Improving Education

The Data Quality Campaign (DQC) has been focused since 2005 on advocating for states to build robust state longitudinal data systems (SLDS). While states have made great progress in their data infrastructure, and should continue to emphasize this work, t data systems alone will not improve outcomes. It is time for both DQC and states to focus on building capacity to use the information that these systems are producing at every level – from classrooms to state houses. To impact system performance and student achievement, the ingrained culture must be replaced with one that focuses on data use for continuous improvement. The effective use of data to inform decisions, provide transparency, improve the measurement of outcomes, and fuel continuous improvement will not come to fruition unless there is a system wide focus on building capacity around the collection, analysis, dissemination, and use of this data, including through research.

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

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

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

Psychosocial and environmental predictors of calcium intake, physical activity and bone health among adolescent girls

Objectives. The purpose of this study was to identify the psychosocial and environmental predictors and the pathways they use to influence calcium intake, physical activity and bone health among adolescent girls. Methods. A secondary data analysis using a cross-sectional and longitudinal study design was implemented to examine the associations of interest. Data from the Incorporating More Physical Activity and Calcium in Teens (IMPACT) study collected in 2001-2003 were utilized for the analyses. IMPACT was a 1½ year nutrition and physical activity intervention study conducted among 718 middle-school girls in central Texas. Hierarchical regression modeling and Structural Equation Modeling (SEM) were used to determine the psychosocial predictors of calcium intake, physical activity and bone health at baseline. Hierarchical regression was used to determine if psychosocial factors at baseline were significant predictors of calcium intake and physical activity at follow-up. Data was adjusted for included BMI, lactose intolerance, ethnicity, menarchal status, intervention and participation in 7th grade PE/athletics. Results. Results of the baseline regression analysis revealed that calcium self-efficacy and milk availability at home were the strongest predictors of calcium intake. Friend engagement in physical activity, physical activity self-efficacy and participation in sports teams were the strongest predictors of physical activity. Finally, physical activity outcome expectations, social support and participation in sports teams were significant predictors of stiffness index at baseline. Results of the baseline SEM path analysis found that outcome expectations and milk availability at home directly influenced calcium intake. Knowledge and calcium self-efficacy indirectly influenced calcium intake with outcome expectations as the mediator. Physical activity self-efficacy and social support had significant direct and indirect influence on physical activity with participation in sports teams as the mediator. Participation in sports teams had a direct effect on both physical activity and stiffness index. Results of regression analysis for baseline predicting follow-up showed that participation in sports teams, self-efficacy, outcome expectations and social support at baseline were significant predictors of physical activity at follow-up. Conclusion. Results of this study reinforce the relevance of addressing both, psychosocial and environmental factors which are critical when developing interventions to improve bone health among adolescent girls. ^

Second, do less harm: The state of states' medical error reporting programs

Statement of the problem and public health significance. Hospitals were designed to be a safe haven and respite from disease and illness. However, a large body of evidence points to preventable errors in hospitals as the eighth leading cause of death among Americans. Twelve percent of Americans, or over 33.8 million people, are hospitalized each year. This population represents a significant portion of at risk citizens exposed to hospital medical errors. Since the number of annual deaths due to hospital medical errors is estimated to exceed 44,000, the magnitude of this tragedy makes it a significant public health problem. ^ Specific aims. The specific aims of this study were threefold. First, this study aimed to analyze the state of the states' mandatory hospital medical error reporting six years after the release of the influential IOM report, "To Err is Human." The second aim was to identify barriers to reporting of medical errors by hospital personnel. The third aim was to identify hospital safety measures implemented to reduce medical errors and enhance patient safety. ^ Methods. A descriptive, longitudinal, retrospective design was used to address the first stated objective. The study data came from the twenty-one states with mandatory hospital reporting programs which report aggregate hospital error data that is accessible to the public by way of states' websites. The data analysis included calculations of expected number of medical errors for each state according to IOM rates. Where possible, a comparison was made between state reported data and the calculated IOM expected number of errors. A literature review was performed to achieve the second study aim, identifying barriers to reporting medical errors. The final aim was accomplished by telephone interviews of principal patient safety/quality officers from five Texas hospitals with more than 700 beds. ^ Results. The state medical error data suggests vast underreporting of hospital medical errors to the states. The telephone interviews suggest that hospitals are working at reducing medical errors and creating safer environments for patients. The literature review suggests the underreporting of medical errors at the state level stems from underreporting of errors at the delivery level. ^

Breastfeeding the very preterm hospitalized infant: An ethnographic study

Purpose. To develop a greater understanding of the experience—including the thoughts, feelings, and actions—of mothers' initiation and maintenance of lactation within the context of the NICU following the birth of a very preterm infant. ^ Design. Mixed method with dominant focused ethnographic approach. Setting: A 76-bed neonatal intensive care unit in the largest children's hospital located in a large metropolitan city in southeast Texas. ^ Sample. Purposeful sampling resulted in 23 interviews with 14 subjects. ^ Methods. Mixed method design with a dominant qualitative approach combined with a quantitative component to further identify and expand upon the investigation of the population in question. Open-ended semi-structured interviews and fieldwork were used to explore the experience of breastfeeding in the context of the NICU for mothers of very preterm infants. Longitudinal data obtained from each subject included in-depth interviews, demographic and clinical information, milk expression patterns (including pumping frequency, duration, and milk volumes obtained), and scores obtained from the Edinburgh Postpartum Depression Scale (EPDS). ^ Findings. Thematic analysis revealed that mothers of very preterm infants experienced an interruption in the process of becoming a mother, a paradoxical experience related to aspects of their milk expression routines and patterns, and negotiating the NICU environment. Sub-themes of becoming a mother-interrupted included: attribution, separation, connection, and navigation. Additional sub-themes related to the paradoxical experience included: the pump sometimes acting as a wedge or link to the infant; diversionary thoughts/activities during pumping; and perceptions of milk flow/volume. The process of negotiation included the environment, adaptive/maladaptive strategies related to milk expression, motivating factors related to the provision of breast milk, and learning their infant's feeding cues/abilities. EPDS scores did not reveal congruent differences in those mothers scoring high compared to those scoring low. ^ Conclusions. Understanding the experiences of the mothers in this study allows for a better perspective of breastfeeding the very preterm infant in the context of the NICU. Findings from this study validate the difficult and incremental process of attaining maternal identity and the significant burden placed on these women with regards to the provision of breast milk and breastfeeding during their infant's hospitalization. ^

Racial disparities in CD4 counts at initial HIV-1 diagnosis: Analysis of the Adult Spectrum of HIV disease dataset and public health implications

Introduction. The HIV/AIDS disease burden disproportionately affects minority populations, specifically African Americans. While sexual risk behaviors play a role in the observed HIV burden, other factors including gender, age, socioeconomics, and barriers to healthcare access may also be contributory. The goal of this study was to determine how far down the HIV/AIDS disease process people of different ethnicities first present for healthcare. The study specifically analyzed the differences in CD4 cell counts at the initial HIV-1 diagnosis with respect to ethnicity. The study also analyzed racial differences in HIV/AIDS risk factors. ^ Methods. This is a retrospective study using data from the Adult Spectrum of HIV Disease (ASD), collected by the City of Houston Department of Health. The ASD database contains information on newly reported HIV cases in the Harris County District Hospitals between 1989 and 2000. Each patient had an initial and a follow-up report. The extracted variables of interest from the ASD data set were CD4 counts at the initial HIV diagnosis, race, gender, age at HIV diagnosis and behavioral risk factors. One-way ANOVA was used to examine differences in baseline CD4 counts at HIV diagnosis between racial/ethnic groups. Chi square was used to analyze racial differences in risk factors. ^ Results. The analyzed study sample was 4767. The study population was 47% Black, 37% White and 16% Hispanic [p<0.05]. The mean and median CD4 counts at diagnosis were 254 and 193 cells per ml, respectively. At the initial HIV diagnosis Blacks had the highest average CD4 counts (285), followed by Whites (233) and Hispanics (212) [p<0.001 ]. These statistical differences, however, were only observed with CD4 counts above 350 [p<0.001], even when adjusted for age at diagnosis and gender [p<0.05]. Looking at risk factors, Blacks were mostly affected by intravenous drug use (IVDU) and heterosexuality, whereas Whites and Hispanics were more affected by male homosexuality [ p<0.05]. ^ Conclusion. (1) There were statistical differences in CD4 counts with respect to ethnicity, but these differences only existed for CD4 counts above 350. These differences however do not appear to have clinical significance. Antithetically, Blacks had the highest CD4 counts followed by Whites and Hispanics. (2) 50% of this study group clinically had AIDS at their initial HIV diagnosis (median=193), irrespective of ethnicity. It was not clear from data analysis if these observations were due to failure of early HIV surveillance, HIV testing policies or healthcare access. More studies need to be done to address this question. (3) Homosexuality and bisexuality were the biggest risk factors for Whites and Hispanics, whereas for Blacks were mostly affected by heterosexuality and IVDU, implying a need for different public health intervention strategies for these racial groups. ^

The effect of genotype-by-environment interaction on longitudinal triglyceride profiles in the Bogalusa Heart Study

Triglyceride levels are a component of plasma lipids that are thought to be an important risk factor for coronary heart disease and are influenced by genetic and environmental factors, such as single nucleotide polymorphisms (SNPs), alcohol intake, and smoking. This study used longitudinal data from the Bogalusa Heart Study, a biracial community-based survey of cardiovascular disease risk factors. A sample of 1191 individuals, 4 to 38 years of age, was measured multiple times from 1973 to 2000. The study sample consisted of 730 white and 461 African American participants. Individual growth models were developed in order to assess gene-environment interactions affecting plasma triglycerides over time. After testing for inclusion of significant covariates and interactions, final models, each accounting for the effects of a different SNP, were assessed for fit and normality. After adjustment for all other covariates and interactions, LIPC -514C/T was found to interact with age3, age2, and age and a non-significant interaction of CETP -971G/A genotype with smoking status was found (p = 0.0812). Ever-smokers had higher triglyceride levels than never smokers, but persons heterozygous at this locus, about half of both races, had higher triglyceride levels after smoking cessation compared to current smokers. Since tobacco products increase free fatty acids circulating in the bloodstream, smoking cessation programs have the potential to ultimately reduce triglyceride levels for many persons. However, due to the effect of smoking cessation on the triglyceride levels of CETP -971G/A heterozygotes, the need for smoking prevention programs is also demonstrated. Both smoking cessation and prevention programs would have a great public health impact on minimizing triglyceride levels and ultimately reducing heart disease. ^

A comparative analysis of CHIP Perinatal policy in twelve states

The purpose of this comparative analysis of CHIP Perinatal policy (42 CFR § 457) was to provide a basis for understanding the variation in policy outputs across the twelve states that, as of June 2007, implemented the Unborn Child rule. This Department of Health and Human Services regulation expanded in 2002 the definition of “child” to include the period from conception to birth, allowing states to consider an unborn child a “targeted low-income child” and therefore eligible for SCHIP coverage. ^ Specific study aims were to (1) describe typologically the structural and contextual features of the twelve states that adopted a CHIP Perinatal policy; (2) describe and differentiate among the various designs of CHIP Perinatal policy implemented in the states; and (3) develop a conceptual model that links the structural and contextual features of the adopting states to differences in the forms the policy assumed, once it was implemented. ^ Secondary data were collected from publicly available information sources to describe characteristics of states’ political system, health system, economic system, sociodemographic context and implemented policy attributes. I posited that socio-demographic differences, political system differences and health system differences would directly account for the observed differences in policy output among the states. ^ Exploratory data analysis techniques, which included median polishing and multidimensional scaling, were employed to identify compelling patterns in the data. Scaled results across model components showed that economic system was most closely related to policy output, followed by health system. Political system and socio-demographic characteristics were shown to be weakly associated with policy output. Goodness-of-fit measures for MDS solutions implemented across states and model components, in one- and two-dimensions, were very good. ^ This comparative policy analysis of twelve states that adopted and implemented HHS Regulation 42 C.F.R. § 457 contributes to existing knowledge in three areas: CHIP Perinatal policy, public health policy and policy sciences. First, the framework allows for the identification of CHIP Perinatal program design possibilities and provides a basis for future studies that evaluate policy impact or performance. Second, studies of policy determinants are not well represented in the health policy literature. Thus, this study contributes to the development of the literature in public health policy. Finally, the conceptual framework for policy determinants developed in this study suggests new ways for policy makers and practitioners to frame policy arguments, encouraging policy change or reform. ^

Determinants of patients' expectations about total knee arthroplasty outcomes

Knee osteoarthritis (OA) is the most prevalent form of arthritis in the US, affecting approximately 37% of adults. Approximately 300,000 total knee arthroplasty (TKA) procedures take place in the United States each year. Total knee arthroplasty is an elective procedure available to patients as an irreversible treatment after failure of previous medical treatments. Some patients sacrifice quality of life and endure many years of pain before making the decision to undergo total knee replacement. In making their decision, it is therefore imperative for patients to understand the procedure, risks and surgical outcomes to create realistic expectations and increase outcome satisfaction. ^ From 2004-2007, 236 OA patients who underwent TKA participated in the PEAKS (Patient Expectations About Knee Surgery) study, an observational longitudinal cohort study, completed baseline and 6 month follow-up questionnaires after the surgery. We performed a secondary data analysis of the PEAKS study to: (1) determine the specific presurgical patient characteristics associated with patients’ presurgical expectations of time to functional recovery; and (2) determine the association between presurgical expectations of time to functional recovery and postsurgical patient capabilities (6 months after TKA). We utilized the WOMAC to measure knee pain and function, the SF-36 to measure health-related quality of life, and the DASS and MOS-SSS to measure psychosocial quality of life variables. Expectation and capability measures were generated from panel of experts. A list of 10 activities was used for this analysis to measure functional expectations and postoperative functional capabilities. ^ The final cohort consisted of 236 individuals, was predominately White with 154 women and 82 men. The mean age was 65 years. Patients were optimistic about their time to functional recovery. Expectation time of being able to perform the list activities per patient had a median of less than 3 months. Patients who expected to be able to perform the functional activities by 3 months had better knee function, less pain and better overall health-related quality of life. Despite expectation differences, all patients showed significant improvement 6 months after surgery. Participant expectation of time to functional recovery was not an independent predictor of capability to perform functional activities at 6 months. Better presurgical patient characteristics were, however, associated with a higher likelihood of being able to perform all activities at 6 months. ^ This study gave us initial insight on the relationship between presurgical patient characteristics and their expectations of functional recovery after total knee replacement. Future studies clarifying the relationship between patient presurgical characteristics and postsurgical functional capabilities are needed.^

The role of cardiovascular comorbidities in ovarian cancer survival

Objective. One facet of cancer care that often goes ignored is comorbidities, or diseases that exist in concert with cancer. Comorbid conditions may affect survival by influencing treatment decisions and prognosis. The purpose of this secondary data analysis was to identify whether a history of cardiovascular comorbidities among ovarian cancer patients influenced survival time at the University of Texas M. D. Anderson Cancer Center. The parent study, Project Peace, has a longitudinal design with an embedded randomized efficacy study which seeks to improve detection of depressive disorders in ovarian, peritoneal, and fallopian tube cancers. ^ Methods. Survival time was calculated for the 249 ovarian cancer patients abstracted by Project Peace staff. Cardiovascular comorbidities were documented as present, based upon information from medical records in addition to self reported comorbidities in a baseline study questionnaire. Kaplan-Meier survival curves were used to compare survival time among patients with a presence or absence of particular cardiovascular comorbidities. Cox Regression proportional models accounted for multivariable factors such as age, staging, family history of cardiovascular comorbidities, and treatment. ^ Results. Among our patient population, there was a statistically significant relationship between shorter survival time and a history of thrombosis, pericardial disease/tamponade, or COPD/pulmonary hypertension. Ovarian cancer patients with a history of thrombosis lived approximately half as long as patients without thrombosis (58.06 months vs. 121.55 months; p=.001). In addition, patients who suffered from pericardial disease/tamponade had poorer survival than those without a history of pericardial disease/tamponade (48 months vs. 80.07 months; p=.002). Ovarian cancer patients with a history of COPD or pulmonary hypertension had a median survival of 60.2 months, while the median survival for patients without these comorbidities was 80.2 months (p=.014). ^ Conclusion. Especially because of its relatively lower survival rate, greater emphasis needs to be placed on the potential influence of cardiovascular comorbid conditions in ovarian cancer.^

A path analysis of an after-school health and exercise program for elementary school children in El Paso, Texas

As schools are pressured to perform on academics and standardized examinations, schools are reluctant to dedicate increased time to physical activity. After-school exercise and health programs may provide an opportunity to engage in more physical activity without taking time away from coursework during the day. The current study is a secondary data analysis of data from a randomized trial of a 10-week after-school program (six schools, n = 903) that implemented an exercise component based on the CATCH physical activity component and health modules based on the culturally-tailored Bienestar health education program. Outcome variables included BMI and aerobic capacity, health knowledge and healthy food intentions as assessed through path analysis techniques. Both the baseline model (χ2 (df = 8) = 16.90, p = .031; RMSEA = .035 (90% CI of .010–.058), NNFI = 0.983 and the CFI = 0.995) and the model incorporating intervention participation proved to be a good fit to the data (χ2 (df = 10) = 11.59, p = .314. RMSEA = .013 (90% CI of .010–.039); NNFI = 0.996 and CFI = 0.999). Experimental group participation was not predictive of changes in health knowledge, intentions to eat healthy foods or changes in Body Mass Index, but it was associated with increased aerobic capacity, β = .067, p < .05. School characteristics including SES and Language proficiency proved to be significantly associated with changes in knowledge and physical indicators. Further effects of school level variables on intervention outcomes are recommended so that tailored interventions can be developed aimed at the specific characteristics of each participating school. ^

Secondary analysis of effectivness of tinidazole-containing-quadruple eradication therapy in asymptomatic H. pylori infected children in El Paso

Helicobacter pylori infection is frequently acquired during childhood. This microorganism is known to cause gastritis, and duodenal ulcer in pediatric patients, however most children remain completely asymptomatic to the infection. Currently there is no consensus in favor of treatment of H. pylori infection in asymptomatic children. The firstline of treatment for this population is triple medication therapy including two antibacterial agents and one proton pump inhibitor for a 2 week duration course. Decreased eradication rate of less than 75% has been documented with the use of this first-line therapy but novel tinidazole-containing quadruple sequential therapies seem worth investigating. None of the previous studies on such therapy has been done in the United States of America. As part of an iron deficiency anemia study in asymptomatic H. pylori infected children of El Paso, Texas, we conducted a secondary data analysis of study data collected in this trial to assess the effectiveness of this tinidazole-containing sequential quadruple therapy compared to placebo on clearing the infection. Subjects were selected from a group of asymptomatic children identified through household visits to 11,365 randomly selected dwelling units. After obtaining parental consent and child assent a total of 1,821 children 3-10 years of age were screened and 235 were positive to a novel urine immunoglobulin class G antibodies test for H. pylori infection and confirmed as infected using a 13C urea breath test, using a hydrolysis urea rate >10 μg/min as cut-off value. Out of those, 119 study subjects had a complete physical exam and baseline blood work and were randomly allocated to four groups, two of which received active H. pylori eradication medication alone or in combination with iron, while the other two received iron only or placebo only. Follow up visits to their houses were done to assess compliance and occurrence of adverse events and at 45+ days post-treatment, a second urea breath test was performed to assess their infection status. The effectiveness was primarily assessed on intent to treat basis (i.e., according to their treatment allocation), and the proportion of those who cleared their infection using a cut-off value >10 μg/min of for urea hydrolysis rate, was the primary outcome. Also we conducted analysis on a per-protocol basis and according to the cytotoxin associated gene A product of the H. pylori infection status. Also we compared the rate of adverse events across the two arms. On intent-to-treat and per-protocol analyses, 44.3% and 52.9%, respectively, of the children receiving the novel quadruple sequential eradication cleared their infection compared to 12.2% and 15.4% in the arms receiving iron or placebo only, respectively. Such differences were statistically significant (p<0.001). The study medications were well accepted and safe. In conclusion, we found in this study population, of mostly asymptomatically H. pylori infected children, living in the US along the border with Mexico, that the quadruple sequential eradication therapy cleared the infection in only half of the children receiving this treatment. Research is needed to assess the antimicrobial susceptibility of the strains of H. pylori infecting this population to formulate more effective therapies. ^

Demographic analysis of Certified Industrial Hygienists (CIH) with the use of survival analysis techniques

Objective. The goal of this study is to characterize the current workforce of CIHs, the lengths of professional practice careers of the past and current CIHs.^ Methods. This is a secondary data analysis of data compiled from all of the nearly 50 annual roster listings of the American Board of Industrial Hygiene (ABIH) for Certified Industrial Hygienists active in each year since 1960. Survival analysis was performed as a technique to measure the primary outcome of interest. The technique which was involved in this study was the Kaplan-Meier method for estimating the survival function.^ Study subjects: The population to be studied is all Certified Industrial Hygienists (CIHs). A CIH is defined by the ABIH as an individual who has achieved the minimum requirements for education, working experience and through examination, has demonstrated a minimum level of knowledge and competency in the prevention of occupational illnesses. ^ Results. A Cox-proportional hazards model analysis was performed by different start-time cohorts of CIHs. In this model we chose cohort 1 as the reference cohort. The estimated relative risk of the event (defined as retirement, or absent from 5 consecutive years of listing) occurred for CIHs for cohorts 2,3,4,5 relative to cohort 1 is 0.385, 0.214, 0.234, 0.299 relatively. The result show that cohort 2 (CIHs issued from 1970-1980) has the lowest hazard ratio which indicates the lowest retirement rate.^ Conclusion. The manpower of CIHs (still actively practicing up to the end of 2009) increased tremendously starting in 1980 and grew into a plateau in recent decades. This indicates that the supply and demand of the profession may have reached equilibrium. More demographic information and variables are needed to actually predict the future number of CIHs needed. ^

BEST SUBSET SELECTION FOR CATEGORICAL DATA

When choosing among models to describe categorical data, the necessity to consider interactions makes selection more difficult. With just four variables, considering all interactions, there are 166 different hierarchical models and many more non-hierarchical models. Two procedures have been developed for categorical data which will produce the "best" subset or subsets of each model size where size refers to the number of effects in the model. Both procedures are patterned after the Leaps and Bounds approach used by Furnival and Wilson for continuous data and do not generally require fitting all models. For hierarchical models, likelihood ratio statistics (G('2)) are computed using iterative proportional fitting and "best" is determined by comparing, among models with the same number of effects, the Pr((chi)(,k)('2) (GREATERTHEQ) G(,ij)('2)) where k is the degrees of freedom for ith model of size j. To fit non-hierarchical as well as hierarchical models, a weighted least squares procedure has been developed.^ The procedures are applied to published occupational data relating to the occurrence of byssinosis. These results are compared to previously published analyses of the same data. Also, the procedures are applied to published data on symptoms in psychiatric patients and again compared to previously published analyses.^ These procedures will make categorical data analysis more accessible to researchers who are not statisticians. The procedures should also encourage more complex exploratory analyses of epidemiologic data and contribute to the development of new hypotheses for study. ^