Statistics and Probability Commons^™

A Bayesian Hierarchical Approach To Multirater Correlated Roc Analysis, Tim Johnson, Valen Johnson

The University of Michigan Department of Biostatistics Working Paper Series

In a common ROC study design, several readers are asked to rate diagnostics of the same cases processed under different modalities. We describe a Bayesian hierarchical model that facilitates the analysis of this study design by explicitly modeling the three sources of variation inherent to it. In so doing, we achieve substantial reductions in the posterior uncertainty associated with estimates of the differences in areas under the estimated ROC curves and corresponding reductions in the mean squared error (MSE) of these estimates. Based on simulation studies, both the widths of confidence intervals and MSE of estimates of differences in the …

A Bayesian Chi-Squared Test For Goodness Of Fit, Valen Johnson

The University of Michigan Department of Biostatistics Working Paper Series

This article describes an extension of classical x 2 goodness-of-fit tests to Bayesian model assessment. The extension, which essentially involvesevaluating Pearson's goodness-of-fit statistic at a parameter value drawn from its posterior distribution, has the important property that it is asymptoti-cally distributed as a x2 random variable on K-1 degrees of freedom, indepen-dently of the dimension of the underlying parameter vector. By averaging over the posterior distribution of this statistic, a global goodness-of-fit diagnostic is obtained. Advantages of this diagnostic{which may be interpreted as the area under an ROC curve{include ease of interpretation, computational conve-nience, and favorable power properties. The proposed …

Individualized Predictions Of Disease Progression Following Radiation Therapy For Prostate Cancer., Jeremy Taylor, Menggang Yu, Howard M. Sandler

The University of Michigan Department of Biostatistics Working Paper Series

Background: Following treatment for localized prostate cancer, men are monitored with serial PSA measurements. Refining the predictive value of post-treatment PSA determinations may add to clinical management and we have developed a model that predicts for an individual patient future PSA values and estimates the time to future clinical recurrence.

Methods: Data from 934 patients treated for prostate cancer between 1987 and 2000 were used to develop a comprehensive statistical model to fit the clinical recurrence events and pattern of PSA data. A logistic regression model was used for the probability of cure, non-linear hierarchical mixed models were used for …

Piecewise Constant Cross-Ratio Estimation For Association In Bivariate Survival Data With Application To Studying Markers Of Menopausal Transition, Bin Nan, Xihong Lin, Lynda D. Lisabet, Sioban Harlow

The University of Michigan Department of Biostatistics Working Paper Series

A question of significant interest in female reproductive aging is to identify bleeding criteria for the menopausal transition. Although various bleeding criteria, or markers, have been proposed for the menopausal transition, their validity has not been adequately examined. The Tremin Trust data are collected from a long-term cohort study that followed a group of women throughout their whole reproductive life, and provide a unique opportunity for assessing the association between age at onset of a bleeding marker and age onset of menopause. Formal statistical analysis of this dependence is challenging give the fact that both the marker event and menopause …

Individual Prediction In Prostate Cancer Studies Using A Joint Longitudinal-Survival-Cure Model, Menggang Yu, Jeremy Taylor, Howard M. Sandler

The University of Michigan Department of Biostatistics Working Paper Series

For monitoring patients treated for prostate cancer, Prostate Specific Antigen (PSA) is measured periodically after they receive treatment. Increases in PSA are suggestive of recurrence of the cancer and are used in making decisions about possible new treatments. The data from studies of such patients typically consist of longitudinal PSA measurements, censored event times and baseline covariates. Methods for the combined analysis of both longitudinal and survival data have been developed in recent years, with the main emphasis being on modeling and estimation. We analyze data from a prostate cancer study that has been extended by adding a mixture structure …

Incorporating Death Into Health-Related Variables In Longitudinal Studies, Paula Diehr, Laura Lee Johnson, Donald L. Patrick, Bruce Psaty

UW Biostatistics Working Paper Series

Background: The aging process can be described as the change in health-related variables over time. Unfortunately, simple graphs of available data may be misleading if some people die, since they may confuse patterns of mortality with patterns of change in health. Methods have been proposed to incorporate death into self-rated health (excellent to poor) and the SF-36 profile scores, but not for other variables.

Objectives: (1) To incorporate death into the following variables: ADLs, IADLs, mini-mental state examination, depressive symptoms, body mass index (BMI), blocks walked per week, bed days, hospitalization, systolic blood pressure, and the timed walk. (2) To …

Overlap Bias In The Case-Crossover Design, With Application To Air Pollution Exposures, Holly Janes, Lianne Sheppard, Thomas Lumley

UW Biostatistics Working Paper Series

The case-crossover design uses cases only, and compares exposures just prior to the event times to exposures at comparable control, or “referent” times, in order to assess the effect of short-term exposure on the risk of a rare event. It has commonly been used to study the effect of air pollution on the risk of various adverse health events. Proper selection of referents is crucial, especially with air pollution exposures, which are shared, highly seasonal, and often have a long term time trend. Hence, careful referent selection is important to control for time-varying confounders, and in order to ensure that …

Uncertainty And The Value Of Diagnostic Information With Application To Axillary Lymph Node Dissection In Breast Cancer, Giovanni Parmigiani

Johns Hopkins University, Dept. of Biostatistics Working Papers

In clinical decision making, it is common to ask whether, and how much, a diagnostic procedure is contributing to subsequent treatment decisions. Statistically, quantification of the value of the information provided by a diagnostic procedure can be carried out using decision trees with multiple decision points, representing both the diagnostic test and the subsequent treatments that may depend on the test's results. This article investigates probabilistic sensitivity analysis approaches for exploring and communicating parameter uncertainty in such decision trees. Complexities arise because uncertainty about a model's inputs determines uncertainty about optimal decisions at all decision nodes of a tree. We …

Survival Model Predictive Accuracy And Roc Curves, Patrick Heagerty, Yingye Zheng

UW Biostatistics Working Paper Series

The predictive accuracy of a survival model can be summarized using extensions of the proportion of variation explained by the model, or R^2, commonly used for continuous response models, or using extensions of sensitivity and specificity which are commonly used for binary response models.

In this manuscript we propose new time-dependent accuracy summaries based on time-specific versions of sensitivity and specificity calculated over risk sets. We connect the accuracy summaries to a previously proposed global concordance measure which is a variant of Kendall's tau. In addition, we show how standard Cox regression output can be used to obtain estimates of …

Semiparametric Estimation Of Time-Dependent: Roc Curves For Longitudinal Marker Data, Yingye Zheng, Patrick Heagerty

UW Biostatistics Working Paper Series

One approach to evaluating the strength of association between a longitudinal marker process and a key clinical event time is through predictive regression methods such as a time-dependent covariate hazard model. For example, a time-varying covariate Cox model specifies the instantaneous risk of the event as a function of the time-varying marker and additional covariates. In this manuscript we explore a second complementary approach which characterizes the distribution of the marker as a function of both the measurement time and the ultimate event time. Our goal is to flexibly extend the standard diagnostic accuracy concepts of sensitivity and specificity to …

Optimization Of Breast Cancer Screening Modalities, Yu Shen, Giovanni Parmigiani

Johns Hopkins University, Dept. of Biostatistics Working Papers

Mathematical models and decision analyses based on microsimulations have been shown to be useful in evaluating relative merits of various screening strategies in terms of cost and mortality reduction. Most investigations regarding the balance between mortality reduction and costs have focused on a single modality, mammography. A systematic evaluation of the relative expenses and projected benefit of combining clinical breast examination and mammography is not at present available. The purpose of this report is to provide methodologic details including assumptions and data used in the process of modeling for complex decision analyses, when searching for optimal breast cancer screening strategies …

Modeling The Incubation Period Of Anthrax, Ron Brookmeyer, Elizabeth Johnson, Sarah Barry

Johns Hopkins University, Dept. of Biostatistics Working Papers

Models of the incubation period of anthrax are important to public health planners because they can be used to predict the delay before outbreaks are detected, the size of an outbreak and the duration of time that persons should remain on antibiotics to prevent disease. The difficulty is that there is little direct data about the incubation period in humans. The objective of this paper is to develop and apply models for the incubation period of anthrax. Mechanistic models that account for the biology of spore clearance and germination are developed based on a competing risks formulation. The models predict …

Clinical Effectiveness Of A Subperiosteal Anchorage Device, Monica Anne Witte

Loma Linda University Electronic Theses, Dissertations & Projects

The purpose of this pilot study was to evaluate the clinical efficacy of a subperiosteal anchorage device, the palatal OnPlant™, during orthodontic retraction of protruding anterior teeth in cases requiring maxillary premolar extraction. Seven subjects (5 female, 2 male), ages 13 to 55, were selected for the study. The OnPlant was surgically placed in the mid-palatal region through a well-defined subperiosteal tunnel. Following the manufacturer recommended osseointegration period of four months, the OnPlant was uncovered and attached to the first molars by means of a transpalatal bar. Standard orthodontic treatment then commenced to retract the anterior teeth after the first …

Comparison Of Bracket Bond Strength By Total, Self-Etch And Laser Treatment, Kyo Sung Shawn Kim

Loma Linda University Electronic Theses, Dissertations & Projects

Laser etching of enamel surfaces alters the physical and chemical characteristics of the enamel. These changes in characteristics enhance the bonding to enamel. The purpose of this study was to compare the shear and tensile bond strength of Er,Cr:YSGG hydrokinetic laser system (Biolase) with 37% phosphoric acid and self etching primer 30 minutes and 72 hours after bonding. Four different laser power output setting was observed: 1.5W, 2.0W, 2.5W, and 3.0W. Two hundred forty bovine teeth free of defect, caries, and dentin exposure were mounted in acrylic resin and divided into 24 groups of 10 teeth. Sixteen groups of 10 …

A Bond Strength Comparison Of Led And Halogen Light Curing Units, John Richard Kavanagh

Loma Linda University Electronic Theses, Dissertations & Projects

The purpose of this research was to compare the tensile bond strength of orthodontic flat based buttons bonded to bovine teeth with four commercial light-emitting diode (LED) curing lights and one conventional quartz tungsten halogen (QTH) curing light.

The dental market has recently been introduced to a number of commercially available light-emitting diode (LED) curing lights. Tensile bond strength was evaluated for LED curing lights (Rembrandt® Allegro^TM, Den-Mat Corp, Santa Maria, CA), (LEDemetron, Kerr/Demetron Corp, Danbury, Conn), (Ortholux LED, 3M^TMESPE^TM, St. Paul, MN),and (FLASH-lite 1001, Discus Dental, Culver City, CA) are compared with one …

Time-Series Studies Of Particulate Matter, Michelle L. Bell, Jonathan M. Samet, Francesca Dominici

Johns Hopkins University, Dept. of Biostatistics Working Papers

Studies of air pollution and human health have evolved from descriptive studies of the early phenomena of large increases in adverse health effects following extreme air pollution episodes, to time-series analyses and the development of sophisticated regression models. In fact, advanced statistical methods are necessary to address the many challenges inherent in the detection of a small pollution risk in the presence of many confounders. This paper reviews the history, methods, and findings of the time-series studies estimating health risks associated with short-term exposure to particulate matter, though much of the discussion is applicable to epidemiological studies of air pollution …

Smooth Quantile Ratio Estimation With Regression: Estimating Medical Expenditures For Smoking Attributable Diseases, Francesca Dominici, Scott L. Zeger

Johns Hopkins University, Dept. of Biostatistics Working Papers

In this paper we introduce a semi-parametric regression model for estimating the difference in the expected value of two positive and highly skewed random variables as a function of covariates. Our method extends Smooth Quantile Ratio Estimation (SQUARE), a novel estimator of the mean difference of two positive random variables, to a regression model.

The methodological development of this paper is motivated by a common problem in econometrics where we are interested in estimating the difference in the average expenditures between two populations, say with and without a disease, taking covariates into account. Let Y1 and Y2 be two positive …

A Corrected Pseudo-Score Approach For Additive Hazards Model With Longitudinal Covariates Measured With Error, Xiao Song, Yijian Huang

UW Biostatistics Working Paper Series

In medical studies, it is often of interest to characterize the relationship between a time-to-event and covariates, not only time-independent but also time-dependent. Time-dependent covariates are generally measured intermittently and with error. Recent interests focus on the proportional hazards framework, with longitudinal data jointly modeled through a mixed effects model. However, approaches under this framework depend on the normality assumption of the error, and might encounter intractable numerical difficulties in practice. This motivates us to consider an alternative framework, that is, the additive hazards model, under which little has been done when time-dependent covariates are measured with error. We propose …

Smooth Quantile Ratio Estimation, Francesca Dominici, Leslie Cope, Daniel Q. Naiman, Scott L. Zeger

Johns Hopkins University, Dept. of Biostatistics Working Papers

In a study of health care expenditures attributable to smoking, we seek to compare the distribution of medical costs for persons with lung cancer or chronic obstructive pulmonary disease (cases) to those without (controls) using a national survey which includes hundreds of cases and thousands of controls. The distribution of costs is highly skewed toward larger values, making estimates of the mean from the smaller sample dependent on a small fraction of the biggest values. One approach to deal with the smaller sample is to rely on a simple parametric model such as the log-normal, but this makes the undesirable …

A Varying-Coefficient Cox Model For The Effect Of Age At A Marker Event On Age At Menopause, Bin Nan, Xihong Lin, Lynda D. Lisabeth, Sioban D. Harlow

The University of Michigan Department of Biostatistics Working Paper Series

. It is of recent interest in reproductive health research to investigate the validity of a marker event for the onset of menopausal transition and to estimate age at menopause using age at the marker event. We propose a varying coefficient Cox model to investigate the association between age at a marker event, denned as a specific bleeding pattern change, and age at menopause, where both events are subject to censoring and their association varies with age at the marker event. Estimation proceeds using the regression spline method. The proposed method is applied to the Tremin Trust Data to evaluate …

The Psychological Characteristics Of Obese Children, Thomas Taylor Mitchell

Loma Linda University Electronic Theses, Dissertations & Projects

Several studies along with the surgeon general of the United States have identified obesity as a national health issue. Research indicates that adult obesity treatment has relatively poor long-term effects. Consequently increased attention has been given to addressing obesity in children. This focus is particularly vital due to the fact that children are one of the populations showing the most rapid increase in overweight and obesity. The purpose of this study was to identify the some of the unique psychological characteristics of obese children by analyzing the intake data of children and adolescents attending the Loma Linda University Growing Fit …

The Radiopacity Of 29 Resin Restorative Materials Compared To Enamel, Peter Lubisich Iv

Loma Linda University Electronic Theses, Dissertations & Projects

The radiodensity of a composite resin can greatly influence a clinician's ability to analyze a restoration for failure. The purpose of this study was to determine the relative radiopacities of 29 resin restorative materials compared to enamel, and to verify the comparable radiopacity between shades of 23 of the restorative composite resin materials. Composite resin discs measuring >5mm in diameter and 2 mm in thickness were made and polymerized. Radiographs were taken of the composite resin discs, which included an aluminum step wedge for calibration and a densitometer was used to measure the radiodensity of the composite resin samples. The …

Quantification Of Bone Deposition In Onlay Bone Grafting, Gregory Alan Kammeyer

Loma Linda University Electronic Theses, Dissertations & Projects

The purpose of this study was to quantitatively evaluate; 1) mandibular posterior mono-cortical block onlay bone graft success rates 2) the amount of grafted bone retained after initial resorption in radiographic projections. Seven consecutive patients with 14 mandibular posterior implant sites requiring bone augmentation for implants of a standard 3.8mm x 10mm dimension were sought. Mono-cortical block bone grafts from the ipsilateral ascending ramus were placed then observed with two standardized CT scans, dividing them into two groups: Group A had the first CT at 4 weeks after graft placement and a second taken before Prosthodontic reconstruction, 36.7 weeks later. …

Cross-Calibration Of Stroke Disability Measures: Bayesian Analysis Of Longitudinal Ordinal Categorical Data Using Negative Dependence, Giovanni Parmigiani, Heidi W. Ashih, Gregory P. Samsa, Pamela W. Duncan, Sue Min Lai, David B. Matchar

Johns Hopkins University, Dept. of Biostatistics Working Papers

It is common to assess disability of stroke patients using standardized scales, such as the Rankin Stroke Outcome Scale (RS) and the Barthel Index (BI). The Rankin Scale, which was designed for applications to stroke, is based on assessing directly the global conditions of a patient. The Barthel Index, which was designed for general applications, is based on a series of questions about the patient’s ability to carry out 10 basis activities of daily living. As both scales are commonly used, but few studies use both, translating between scales is important in gaining an overall understanding of the efficacy of …

An Extended General Location Model For Causal Inference From Data Subject To Noncompliance And Missing Values, Yahong Peng, Rod Little, Trivellore E. Raghuanthan

The University of Michigan Department of Biostatistics Working Paper Series

Noncompliance is a common problem in experiments involving randomized assignment of treatments, and standard analyses based on intention-to treat or treatment received have limitations. An attractive alternative is to estimate the Complier-Average Causal Effect (CACE), which is the average treatment effect for the subpopulation of subjects who would comply under either treatment (Angrist, Imbens and Rubin, 1996, henceforth AIR). We propose an Extended General Location Model to estimate the CACE from data with non-compliance and missing data in the outcome and in baseline covariates. Models for both continuous and categorical outcomes and ignorable and latent ignorable (Frangakis and Rubin, 1999) …

Adjusting For Non-Ignorable Verification Bias In Clinical Studies For Alzheimer’S Disease, Xiao-Hua Zhou, Pete Castelluccio

UW Biostatistics Working Paper Series

A common problem for comparing the relative accuracy of two screening tests for Alzheimer’s disease (D) in a two-stage design study is verification bias. If the verification bias can be assumed to be ignorable, Zhou and Higgs (2000) have proposed a maximum likelihood approach to compare the relative accuracy of screening tests in a two-stage design study. However, if the verification mechanism also depends on the unobserved disease status, the ignorable assumption does not hold. In this paper, we discuss how to use a profile likelihood approach to compare the relative accuracy of two screening tests for AD without assuming …

A New Confidence Interval For The Difference Between Two Binomial Proportions Of Paired Data, Xiao-Hua Zhou, Gengsheng Qin

UW Biostatistics Working Paper Series

Motivated by a study on comparing sensitivities and specificities of two diagnostic tests in a paired design when the sample size is small, we first derived an Edgeworth expansion for the studentized difference between two binomial proportions of paired data. The Edgeworth expansion can help us understand why the usual Wald interval for the difference has poor coverage performance in the small sample size. Based on the Edgeworth expansion, we then derived a transformation based confidence interval for the difference. The new interval removes the skewness in the Edgeworth expansion; the new interval is easy to compute, and its coverage …

The Effect Of Ultrasonics On Fibroblast Cells, Sheila A. Harris

Loma Linda University Electronic Theses, Dissertations & Projects

The field of dentistry uses the debriding properties of ultrasonic vibration. It is unknown if this property is detrimental to the periodontal ligament (PDL) of an avulsed tooth. This information is important when a clinician is faced with a debris-covered avulsed tooth following a traumatic event. The purpose of this study was to determine the effect that ultrasonic vibration has on those PDL cells most numerous and most vital to a successful tooth replantation, namely, the fibroblasts.

Several fibroblast cell sources were exposed to varying ultrasonic times. These included a commercially available cell line of human foreskin fibroblast (HFF), a …

Cultural And Psychological Influences On Diabetic Adherence, Keikilani Mcmillin-Williams

Loma Linda University Electronic Theses, Dissertations & Projects

Diabetes mellitus is a serious disease that poses a particular healthcare challenge because progression is considered controllable (Cox, et al, 1985; Vinicor, et al, 1996) yet treatment adherence, and thus outcome, is very poor (Gonder-Frederick, Cox, & Ritterband, 2002; Goodall, 1991). Culture is a lethal risk factor for diabetic contraction and treatment maintenance. Latinos within the United States are two-to-three times more likely to develop complications and die than non-Latinos (Haffner et al, 1996; Rubin, Peyrot, & Saudek, 1991) and are less likely to adhere to treatment (Lipton, Losey, Giachello, Mendez, & Girotti, 1998). Efforts to eliminate health disparities have …

Improved Confidence Intervals For The Sensitivity At A Fixed Level Of Specificity Of A Continuous-Scale Diagnostic Test, Xiao-Hua Zhou, Gengsheng Qin

UW Biostatistics Working Paper Series

For a continuous-scale test, it is an interest to construct a confidence interval for the sensitivity of the diagnostic test at the cut-off that yields a predetermined level of its specificity (eg. 80%, 90%, or 95%). IN this paper we proposed two new intervals for the sensitivity of a continuous-scale diagnostic test at a fixed level of specificity. We then conducted simulation studies to compare the relative performance of these two intervals with the best existing BCa bootstrap interval, proposed by Platt et al. (2000). Our simulation results showed that the newly proposed intervals are better than the BCa bootstrap …