A Bayesian approach to estimating the regression coefficients of a multinomial logit model

A Bayesian approach to estimation of the regression coefficients of a multinominal logit model with ordinal scale response categories is presented. A Monte Carlo method is used to construct the posterior distribution of the link function. The link function is treated as an arbitrary scalar function. Then the Gauss-Markov theorem is used to determine a function of the link which produces a random vector of coefficients. The posterior distribution of the random vector of coefficients is used to estimate the regression coefficients. The method described is referred to as a Bayesian generalized least square (BGLS) analysis. Two cases involving multinominal logit models are described. Case I involves a cumulative logit model and Case II involves a proportional-odds model. All inferences about the coefficients for both cases are described in terms of the posterior distribution of the regression coefficients. The results from the BGLS method are compared to maximum likelihood estimates of the regression coefficients. The BGLS method avoids the nonlinear problems encountered when estimating the regression coefficients of a generalized linear model. The method is not complex or computationally intensive. The BGLS method offers several advantages over Bayesian approaches. ^

Logistic regression models for ordinal response: A study of self-efficacy in colorectal cancer screening

The ordinal logistic regression models are used to analyze the dependant variable with multiple outcomes that can be ranked, but have been underutilized. In this study, we describe four logistic regression models for analyzing the ordinal response variable. ^ In this methodological study, the four regression models are proposed. The first model uses the multinomial logistic model. The second is adjacent-category logit model. The third is the proportional odds model and the fourth model is the continuation-ratio model. We illustrate and compare the fit of these models using data from the survey designed by the University of Texas, School of Public Health research project PCCaSO (Promoting Colon Cancer Screening in people 50 and Over), to study the patient’s confidence in the completion colorectal cancer screening (CRCS). ^ The purpose of this study is two fold: first, to provide a synthesized review of models for analyzing data with ordinal response, and second, to evaluate their usefulness in epidemiological research, with particular emphasis on model formulation, interpretation of model coefficients, and their implications. Four ordinal logistic models that are used in this study include (1) Multinomial logistic model, (2) Adjacent-category logistic model [9], (3) Continuation-ratio logistic model [10], (4) Proportional logistic model [11]. We recommend that the analyst performs (1) goodness-of-fit tests, (2) sensitivity analysis by fitting and comparing different models.^

HPV vaccine coverage in Texas counties: An application of multilevel, small area estimation

Health departments, research institutions, policy-makers, and healthcare providers are often interested in knowing the health status of their clients/constituents. Without the resources, financially or administratively, to go out into the community and conduct health assessments directly, these entities frequently rely on data from population-based surveys to supply the information they need. Unfortunately, these surveys are ill-equipped for the job due to sample size and privacy concerns. Small area estimation (SAE) techniques have excellent potential in such circumstances, but have been underutilized in public health due to lack of awareness and confidence in applying its methods. The goal of this research is to make model-based SAE accessible to a broad readership using clear, example-based learning. Specifically, we applied the principles of multilevel, unit-level SAE to describe the geographic distribution of HPV vaccine coverage among females aged 11-26 in Texas.^ Multilevel (3 level: individual, county, public health region) random-intercept logit models of HPV vaccination (receipt of ≥ 1 dose Gardasil® ) were fit to data from the 2008 Behavioral Risk Factor Surveillance System (outcome and level 1 covariates) and a number of secondary sources (group-level covariates). Sampling weights were scaled (level 1) or constructed (levels 2 & 3), and incorporated at every level. Using the regression coefficients (and standard errors) from the final models, I simulated 10,000 datasets for each regression coefficient from the normal distribution and applied them to the logit model to estimate HPV vaccine coverage in each county and respective demographic subgroup. For simplicity, I only provide coverage estimates (and 95% confidence intervals) for counties.^ County-level coverage among females aged 11-17 varied from 6.8-29.0%. For females aged 18-26, coverage varied from 1.9%-23.8%. Aggregated to the state level, these values translate to indirect state estimates of 15.5% and 11.4%, respectively; both of which fall within the confidence intervals for the direct estimates of HPV vaccine coverage in Texas (Females 11-17: 17.7%, 95% CI: 13.6, 21.9; Females 18-26: 12.0%, 95% CI: 6.2, 17.7).^ Small area estimation has great potential for informing policy, program development and evaluation, and the provision of health services. Harnessing the flexibility of multilevel, unit-level SAE to estimate HPV vaccine coverage among females aged 11-26 in Texas counties, I have provided (1) practical guidance on how to conceptualize and conduct modelbased SAE, (2) a robust framework that can be applied to other health outcomes or geographic levels of aggregation, and (3) HPV vaccine coverage data that may inform the development of health education programs, the provision of health services, the planning of additional research studies, and the creation of local health policies.^

NONLINEAR TIME SERIES/SYSTEM ANALYSIS (STATE-SPACE, DYNAMICS, INTERVENTION, RESILIENCE, KALMAN)

A discussion of nonlinear dynamics, demonstrated by the familiar automobile, is followed by the development of a systematic method of analysis of a possibly nonlinear time series using difference equations in the general state-space format. This format allows recursive state-dependent parameter estimation after each observation thereby revealing the dynamics inherent in the system in combination with random external perturbations.^ The one-step ahead prediction errors at each time period, transformed to have constant variance, and the estimated parametric sequences provide the information to (1) formally test whether time series observations y(,t) are some linear function of random errors (ELEM)(,s), for some t and s, or whether the series would more appropriately be described by a nonlinear model such as bilinear, exponential, threshold, etc., (2) formally test whether a statistically significant change has occurred in structure/level either historically or as it occurs, (3) forecast nonlinear system with a new and innovative (but very old numerical) technique utilizing rational functions to extrapolate individual parameters as smooth functions of time which are then combined to obtain the forecast of y and (4) suggest a measure of resilience, i.e. how much perturbation a structure/level can tolerate, whether internal or external to the system, and remain statistically unchanged. Although similar to one-step control, this provides a less rigid way to think about changes affecting social systems.^ Applications consisting of the analysis of some familiar and some simulated series demonstrate the procedure. Empirical results suggest that this state-space or modified augmented Kalman filter may provide interesting ways to identify particular kinds of nonlinearities as they occur in structural change via the state trajectory.^ A computational flow-chart detailing computations and software input and output is provided in the body of the text. IBM Advanced BASIC program listings to accomplish most of the analysis are provided in the appendix. ^

The prevalence of periodontal disease in Mexican American vs European American elderly in three distinct socio-economic communities in San Antonio, Texas

Purpose of the Study: This study evaluated the prevalence of periodontal disease between Mexican American elderly and European American elderly residing in three socio-economically distinct neighborhoods in San Antonio, Texas. ^ Study Group: Subjects for the original protocol were participants of the Oral Health: San Antonio Longitudinal Study of Aging (OH: SALSA), which began with National Institutes of Health (NIH) funding in 1993 (M.J. Saunders, PI). The cohort in the study was the individuals who had been enrolled in Phases I and III of the San Antonio Heart Study (SAHS). This SAHS/SALSA sample is a community-based probability sample of Mexican American and European American residents from three socio-economically distinct San Antonio neighborhoods: low-income barrio, middle-income transitional, and upper-income suburban. The OH: SALSA cohort was established between July 1993 and May 1998 by sampling two subsets of the San Antonio Heart Study (SAHS) cohort. These subsets included the San Antonio Longitudinal Study of Aging (SALSA) cohort, comprised of the oldest members of the SAHS (age 65+ yrs. old), and a younger set of controls (age 35-64 yrs. old) sampled from the remainder of the SAHS cohort. ^ Methods: The study used simple descriptive statistics to describe the sociodemographic characteristics and periodontal disease indicators of the OH: SALSA participants. Means and standard deviations were used to summarize continuous measures. Proportions were used to summarize categorical measures. Simple m x n chi square statistics was used to compare ethnic differences. A multivariable ordered logit regression was used to estimate the prevalence of periodontal disease and test ethnic group and neighborhood differences in the prevalence of periodontal disease. A multivariable model adjustment for socio-economic status (income and education), gender, and age (treated as confounders) was applied. ^ Summary: In the unadjusted and adjusted model, Mexican American elderly demonstrated the greatest prevalence for periodontitis, p < 0.05. Mexican American elderly in barrio neighborhoods demonstrated the greatest prevalence for severe periodontitis, with unadjusted prevalence rates of 31.7%, 22.3%, and 22.4% for Mexican American elderly barrio, transitional, and suburban neighborhoods, respectively. Also, Mexican American elderly had adjusted prevalence rates of 29.4%, 23.7%, and 20.4% for barrio, transitional, and suburban neighborhoods, respectively. ^ Conclusion: This study indicates that the prevalence of periodontal disease is an important oral health issue among the Mexican American elderly. The results suggest that the socioeconomic status of the residential neighborhood increased the risk for severe periodontal disease among the Mexican American elderly when compared to European American elderly. A viable approach to recognizing oral health disparities in our growing population of Mexican American elderly is imperative for the provision of special care programs that will help increase the quality of care in this minority population.^

Factors associated with HPV16 status and survival among patients with squamous cell carcinoma of the oropharynx

Background: Squamous cell carcinoma of the oropharynx (SCCOP) is characterized by local tumor aggressiveness, high recurrence rates, a high incidence of second primary tumors, and medical comorbidities. Significant trends in demographic and clinical characteristics as well as survival among SCCOP patients have been observed over time, likely owing to the changing etiology of the disease. Human papillomavirus type 16 (HPV16) infection is associated with a significant proportion of these cancers. Biomarkers that may aid in identifying patients that are at higher risk of recurrence and death are important so that these patients may be followed more closely to improve their quality of life. ^ Study population and methods: The retrospective review (Specific Aim 2) included 3891 newly diagnosed, previously untreated patients presenting to our institution between 1955 and 2004. A total of 2299 patients treated at our institution were included in survival and recursive partitioning analysis. The prospective cohort study (Specific Aim 3) included 266 patients presenting to our institution between January 2006 and September 2009. ^ Results: The results from the retrospective review showed that over time, patients presented at younger ages and were more likely to have base of tongue/tonsil tumors and to be never/former smokers and moreover survival improved significantly over time. In survival and recursive partitioning analyses, the TNM staging system was efficient in prognosticating patients prior to 1995. However, in the recent decade, the TNM staging system was completely inadequate. The factors having the greatest positive effect on overall survival since 1995 were those common to HPV-associated SCCOP. The results from the prospective cohort study indicate that patients with high nodal stage and those with late stage disease have increased levels of pretreatment serum HPV DNA. ^ Conclusions: We saw a distinct improvement in survival among SCCOP patients over the past 50 years at our institution. The main factors contributing to this were changes in clinical characteristics, in particular surrogates for HPV status. The current TNM staging system for SCCOP is inadequate and incorporation of HPV status (and perhaps smoking status) is encouraged. Furthermore, although pretreatment circulating levels of HPV DNA was associated with higher N category and overall disease stage, it has limited utility as a marker for recurrence among SCCOP patients.^

Renal Cryoablation: Investigation of Periprocedural Visualization Tools and Treatment Response Quantification

Cryoablation for small renal tumors has demonstrated sufficient clinical efficacy over the past decade as a non-surgical nephron-sparing approach for treating renal masses for patients who are not surgical candidates. Minimally invasive percutaneous cryoablations have been performed with image guidance from CT, ultrasound, and MRI. During the MRI-guided cryoablation procedure, the interventional radiologist visually compares the iceball size on monitoring images with respect to the original tumor on separate planning images. The comparisons made during the monitoring step are time consuming, inefficient and sometimes lack the precision needed for decision making, requiring the radiologist to make further changes later in the procedure. This study sought to mitigate uncertainty in these visual comparisons by quantifying tissue response to cryoablation and providing visualization of the response during the procedure. Based on retrospective analysis of MR-guided cryoablation patient data, registration and segmentation algorithms were investigated and implemented for periprocedural visualization to deliver iceball position/size with respect to planning images registered within 3.3mm with at least 70% overlap and a quantitative logit model was developed to relate perfusion deficit in renal parenchyma visualized in verification images as a result of iceball size visualized in monitoring images. Through retrospective study of 20 patient cases, the relationship between likelihood of perfusion loss in renal parenchyma and distance within iceball was quantified and iteratively fit to a logit curve. Using the parameters from the logit fit, the margin for 95% perfusion loss likelihood was found to be 4.28 mm within the iceball. The observed margin corresponds well with the clinically accepted margin of 3-5mm within the iceball. In order to display the iceball position and perfusion loss likelihood to the radiologist, algorithms were implemented to create a fast segmentation and registration module which executed in under 2 minutes, within the clinically-relevant 3 minute monitoring period. Using 16 patient cases, the average Hausdorff distance was reduced from 10.1mm to 3.21 mm with average DSC increased from 46.6% to 82.6% before and after registration.