Statistics and Probability Commons^™

Effect Of Dexamethasone Prodrug On Inflamed Temporomandibular Joints In Juvenile Rats, Mitchell Knudsen, Matthew Bury, Callie Holwegner, Adam L. Reinhardt, Fang Yuan, Yijia Zhang, Peter Giannini, D. B. Marx, Dong Wang, Richard A. Reinhardt

Department of Statistics: Faculty Publications

Introduction: Juvenile idiopathic arthritis (JIA) often causes inflammation of the temporomandibular joint (TMJ) and has been treated with both systemic and intra-articular steroids, with concerns about effects on growing bones. In this study, we evaluated the impact of a macromolecular prodrug of dexamethasone (P-DEX) with inflammation-targeting potential applied systemically or directly to the TMJ.

Methods: Joint inflammation was initiated by injecting two doses of complete Freund’s adjuvant (CFA) at 1-month intervals into the right TMJs of 24 growing Sprague–Dawley male rats (controls on left side). Four additional rats were not manipulated. With the second CFA injection, animals received (1) 5 …

Threshold Models For Genome-Enabled Prediction Of Ordinal Categorical Traits In Plant Breeding, Osval A. Montesinos-López, Abelardo Montesinos-López, Paulino Pérez-Rodríguez, Gustavo De Los Campos, Kent M. Eskridge, José Crossa

Department of Statistics: Faculty Publications

Categorical scores for disease susceptibility or resistance often are recorded in plant breeding. The aim of this study was to introduce genomic models for analyzing ordinal characters and to assess the predictive ability of genomic predictions for ordered categorical phenotypes using a threshold model counterpart of the Genomic Best Linear Unbiased Predictor (i.e., TGBLUP). The threshold model was used to relate a hypothetical underlying scale to the outward categorical response. We present an empirical application where a total of nine models, five without interaction and four with genomic x environment interaction (G·E) and genomic additive x additive x environment interaction …

Measuring Peer Socialization For Adolescent Substance Use:A Comparison Of Perceived And Actual Friends’ Substance Use Effects, Arielle R. Deutsch, Pavel Chernyavskiy, Douglas Steinley, Wendy S. Slutske

Department of Statistics: Faculty Publications

Objective: There has been an increase in the use of social network analysis in studies of peer socialization effects on adolescent substance use. Some researchers argue that social network analyses provide more accurate measures of peer substance use, that the alternate strategy of assessing perceptions of friends’ drug use is biased, and that perceptions of peer use and actual peer use represent different constructs. However, there has been little research directly comparing the two effects, and little is known about the extent to which the measures differ in the magnitude of their influence on adolescent substance use, as well as …

A Copula Based Approach For Design Of Multivariate Random Forests For Drug Sensitivity Prediction, Saad Haider, Raziur Rahman, Souparno Ghosh, Ranadip Pal

Department of Statistics: Faculty Publications

Modeling sensitivity to drugs based on genetic characterizations is a significant challenge in the area of systems medicine. Ensemble based approaches such as Random Forests have been shown to perform well in both individual sensitivity prediction studies and team science based prediction challenges. However, Random Forests generate a deterministic predictive model for each drug based on the genetic characterization of the cell lines and ignores the relationship between different drug sensitivities during model generation. This application motivates the need for generation of multivariate ensemble learning techniques that can increase prediction accuracy and improve variable importance ranking by incorporating the relationships …

Design Of Probabilistic Random Forests With Applications To Anticancer Drug Sensitivity Prediction, Raziur Rahman, Saad Haider, Souparno Gosh, Ranadip Pal

Department of Statistics: Faculty Publications

Random forests consisting of an ensemble of regression trees with equal weights are frequently used for design of predictive models. In this article, we consider an extension of the methodology by representing the regression trees in the form of probabilistic trees and analyzing the nature of heteroscedasticity. The probabilistic tree representation allows for analytical computation of confidence intervals (CIs), and the tree weight optimization is expected to provide stricter CIs with comparable performance in mean error. We approached the ensemble of probabilistic trees’ prediction from the perspectives of a mixture distribution and as a weighted sum of correlated random variables. …

S1: Supplementary Information For Article: A Copula Based Approach For Design Of Multivariate Random Forests For Drug Sensitivity Prediction, Saad Haider, Raziur Rahman, Souparno Ghosh, Ranadip Pal

Department of Statistics: Faculty Publications

Changes in performance with prior feature selection

Random forest (RF) is designed to create uncorrelated trees using random subsets of features in each node of each tree. RF by itself is a great tool for feature selection from a high dimensional set of features. But we observed that the prediction accuracy is improved when a prior feature selection (RELIEFF) [1] approach is implemented. Table A shows the performance of RF, VMRF and CMRF with and without RELIEFF feature selection in 2 drug sets of GDSC.

Performance Analysis for drugsets consisting of more 8 than two drugs

We have generated empirical …

Voc Emissions From Beef Feedlot Pen Surfaces As Affected By Within-Pen Location, Moisture And Temperature, Bryan L. Woodbury, John E. Gilley, David B. Parker, David B. Marx, Roger A. Eigenberg

Department of Statistics: Faculty Publications

A laboratory study was conducted to evaluate the effects of pen location, moisture, and temperature on emissions of volatile organic compounds (VOC) from surface materials obtained from feedlot pens where beef cattle were fed a diet containing 30% wet distillers grain plus solubles. Surface materials were collected from the feed trough (bunk), drainage, and raised areas (mounds) within three feedlot pens. The surface materials were mixed with water to represent dry, wet, or saturated conditions and then incubated at temperatures of 5, 15, 25 and 35 C. A wind tunnel and gas chromatograph-mass spectrometer were used to collect and quantify …

Confidence Intervals For Ranks Of Age-Adjusted Rates Across States Or Counties, Shunpu Zhang, Jun Luo, Li Zhu, David G. Stinchcomb, Dave Campbell, Ginger Carter, Scott Gilkeson, Eric J. Feuer

Department of Statistics: Faculty Publications

Health indices provide information to the general public on the health condition of the community. They can also be used to inform the government’s policy making, to evaluate the effect of a current policy or healthcare program, or for program planning and priority setting. It is a common practice that the health indices across different geographic units are ranked and the ranks are reported as fixed values. We argue that the ranks should be viewed as random and hence should be accompanied by an indication of precision (i.e., the confidence intervals). A technical difficulty in doing so is how to …

Parametric And Nonparametric Statistical Methods For Genomic Selection Of Traits With Additive And Epistatic Genetic Architectures, Reka Howard, Alicia L. Carriquiry, William D. Beavis

Department of Statistics: Faculty Publications

Parametric and nonparametric methods have been developed for purposes of predicting phenotypes. These methods are based on retrospective analyses of empirical data consisting of genotypic and phenotypic scores. Recent reports have indicated that parametric methods are unable to predict phenotypes of traits with known epistatic genetic architectures. Herein, we review parametric methods including least squares regression, ridge regression, Bayesian ridge regression, least absolute shrinkage and selection operator (LASSO), Bayesian LASSO, best linear unbiased prediction (BLUP), Bayes A, Bayes B, Bayes C, and Bayes Cπ. We also review nonparametric methods including Nadaraya-Watson estimator, reproducing kernel Hilbert space, support vector machine regression, …

Practical Guidelines For The Comprehensive Analysis Of Chip-Seq Data, Timonthy Bailey, Pawel Krajewski, Istvan Ladunga, Celine Lefebvre, Qunhua Li, Tao Liu, Pedro Madrigal, Cenny Taslim, Jie Zhang

Department of Statistics: Faculty Publications

Mapping the chromosomal locations of transcription factors, nucleosomes, histone modifications, chromatin remodeling enzymes, chaperones, and polymerases is one of the key tasks of modern biology, as evidenced by the Encyclopedia of DNA Elements (ENCODE) Project. To this end, chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) is the standard methodology. Mapping such protein-DNA interactions in vivo using ChIP-seq presents multiple challenges not only in sample preparation and sequencing but also for computational analysis. Here, we present step-by-step guidelines for the computational analysis of ChIP-seq data. We address all the major steps in the analysis of ChIP-seq data: sequencing depth selection, quality …

Drosophila Melanogaster Selection For Survival After Infection With Bacillus Cereus Spores: Evolutionary Genetic And Phenotypic Investigations Of Respiration And Movement, Junjie Ma, Andrew K. Benson, Stephen D. Kachman, Deidra J. Jacobsen, Lawrence G. Harshman

Department of Statistics: Faculty Publications

Laboratory populations of D.melanogaster have been subjected to selection for survival after live spores of B. cereus were introduced as a pathogenic agent. The present study was designed to investigate correlated traits: respiration as ametabolic trait and movement as a behavioral trait.An underlying hypothesiswas that the evolution of increased survival after B. cereus infection exerts ametabolic cost associated with elevated immunity and this would be detected by increased respiration rates. There was support for this hypothesis in the male response to selection, but not for selected-line females. Two phenotypic effects were also observed in the study. Females especially showed a …

The Change In Differing Leukocyte Populations During Vaccination To Bovine Respiratory Disease And Their Correlations With Lung Scores, Health Records, And Average Daily Gain, R. J. Leach, C. G. Chitko-Mckown, G. L. Bennett, S. A. Jones, Stephen D. Kachman, J. W. Keele, K. A. Leymaster, R. M. Thallman, L. A. Kuehn

Department of Statistics: Faculty Publications

Bovine respiratory disease (BRD) is the most economically important disease in U.S. feedlots. Infection can result in morbidity, mortality, and reduced average daily gain. Cheap and reliable genetic methods of prediction and protection from BRD would be highly advantageous to the industry. The immune response may correlate with BRD incidence. Cattle (n = 2,182) were vaccinated against common viral and bacterial pathogens of BRD. Two blood samples were collected, one during booster vaccination and one 21d later, enabling 3 phenotypes for each trait [prebooster (pre), postbooster (post), and delta (post minus pre)]. From the blood samples innate and adaptive responses …

Prediction In M-Complete Problems With Limited Sample Size, Jennifer Lynn Clarke, Bertrand Clarke, Chi-Wai Yu

Department of Statistics: Faculty Publications

We define a new Bayesian predictor called the posterior weighted median (PWM) and compare its performance to several other predictors including the Bayes model average under squared error loss, the Barbieri-Berger median model predictor, the stacking predictor, and the model average predictor based on Akaike's information criterion. We argue that PWM generally gives better performance than other predictors over a range of M-complete problems. This range is between the M-closed-M-complete boundary and the M-complete- M-open boundary. Indeed, as a problem gets closer to M-open, it seems that M-complete predictive methods begin to break down. Our comparisons rest on extensive simulations …

Comparison Of Molecular Breeding Values Based On Within- And Across-Breed Training In Beef Cattle, Stephen D. Kachman, Matthew L. Spangler, Gary L. Bennett, Kathryn J. Hanford, Larry A. Kuehn, Warren M. Snelling, Mark Thallman, Mahdi Saatchi, Dorian J. Garrick, Robert D. Schnabel, Jeremy F. Taylor, E John Pollak

Department of Statistics: Faculty Publications

Background: Although the efficacy of genomic predictors based on within-breed training looks promising, it is necessary to develop and evaluate across-breed predictors for the technology to be fully applied in the beef industry. The efficacies of genomic predictors trained in one breed and utilized to predict genetic merit in differing breeds based on simulation studies have been reported, as have the efficacies of predictors trained using data from multiple breeds to predict the genetic merit of purebreds. However, comparable studies using beef cattle field data have not been reported.

Methods: Molecular breeding values for weaning and yearling weight …

Murine Gut Microbiota Is Defined By Host Genetics And Modulates Variation Of Metabolic Traits, Autumn M. Mcknite, Maria Elisa Perez-Munoz, Lu Lu, Evan G. Williams, Simon Brewer, Penelope A. Andreux, John W. M. Bastiaansen, Xusheng Wang, Stephen D. Kachman, Johan Auwerx, Robert W. Williams, Andrew K. Benson, Daniel A. Peterson, Daniel C. Ciobanu

Department of Statistics: Faculty Publications

The gastrointestinal tract harbors a complex and diverse microbiota that has an important role in host metabolism. Microbial diversity is influenced by a combination of environmental and host genetic factors and is associated with several polygenic diseases. In this study we combined next-generation sequencing, genetic mapping, and a set of physiological traits of the BXD mouse population to explore genetic factors that explain differences in gut microbiota and its impact on metabolic traits. Molecular profiling of the gut microbiota revealed important quantitative differences in microbial composition among BXD strains. These differences in gut microbial composition are influenced by host-genetics, which …

Transcriptome-Wide Changes In Chlamydomonas Reinhardtii Gene Expression Regulated By Carbon Dioxide And The Co2-Concentrating Mechanism Regulator Cia5/Ccm1, Wei Fang, Yaqing Si, Stephen Douglass, David Casero, Sabeeha S. Merchant, Matteo Pellegrini, Istvan Ladunga, Peng Liu, Martin H. Spalding

Department of Statistics: Faculty Publications

Copyright (c) 2012 American Society of Plant Biologists. Used by permission.

Anticancer Peptidylarginine Deiminase (Pad) Inhibitors Regulate The Autophagy Flux And The Mammalian Target Of Rapamycin Complex 1 Activity, Yuji Wang, Pingxin Li, Shu Wang, Jing Hu, Xiangyun Amy Chen, Jianhui Wu, Megan Fisher, Kira Oshaben, Na Zhao, Ying Gu, Dong Wang, Gong Chen, Yanming Wang

Department of Statistics: Faculty Publications

Tumor suppressor genes are frequently silenced in cancer cells by enzymes catalyzing epigenetic histone modifications. The peptidylarginine deiminase family member PAD4 (also called PADI4) is markedly overexpressed in a majority of human cancers, suggesting that PAD4 is a putative target for cancer treatment. Here, we have generated novel PAD inhibitors with low micromolar IC₅₀ in PAD activity and cancer cell growth inhibition. The lead compound YW3-56 alters the expression of genes controlling the cell cycle and cell death, including SESN2 that encodes an upstream inhibitor of the mammalian target of rapamycin complex 1 (mTORC1) signaling pathway. Guided by the …

Sample Size Under Inverse Negative Binomial Group Testing For Accuracy In Parameter Estimation, Osval Antonio Montesinos-López, Abelardo Montesinos-López, Jose Crossa, Kent M. Eskridge

Department of Statistics: Faculty Publications

Background:The group testing method has been proposed for the detection and estimation of genetically modified plants (adventitious presence of unwanted transgenic plants, AP). For binary response variables (presence or absence), group testing is efficient when the prevalence is low, so that estimation, detection, and sample size methods have been developed under the binomial model. However, when the event is rare (low prevalence

Methodology/Principal Findings: This research proposes three sample size procedures (two computational and one analytic) for estimating prevalence using group testing under inverse (negative) binomial sampling. These methods provide the required number of positive pools (r_m) …

Muscle Organization In Individuals With And Without Pain And Joint Dysfunction, J. C. Nickel, Y. M. Gonzalez, W. D. Mccall, R. Ohrbach, D. B. Marx, H. Liu, L. R. Iwasaki

Department of Statistics: Faculty Publications

Central nervous system organization of masticatory muscles determines the magnitude of joint and muscle forces. Validated computer-assisted models of neuromuscular organization during biting were used to determine organization in individuals with and without temporomandibular disorders (TMD). Ninety-one individuals (47 women, 44 men) were assigned to one of four diagnostic groups based on the presence (+) or absence (-) of pain (P) and bilateral temporomandibular joint disc displacement (DD). Electromyography and bite-forces were measured during right and left incisor and molar biting. Two three-dimensional models employing neuromuscular objectives of minimization of joint loads (MJL) or muscle effort (MME) simulated biting tasks. …

Learning Interdisciplinary Pedagogies, Alison J. Friedow, Erin E. Blankenship, Jennifer L. Green, Walter Stroup

Department of Statistics: Faculty Publications

Advocates of interdisciplinary teaching and learning in higher education suggest that interdisciplinary courses “promise a wide range of desirable educational outcomes for students” (Newell 1994: 35). These outcomes include enhanced affective and cognitive abilities, increased understanding of multiple perspectives, greater appreciation for ambiguity, and superior capacities for creative thinking, among others (35). Despite claims about the possibilities interdisciplinary learning offers, we have few examples of how faculty from different disciplines work together to create interdisciplinary classroom environments where such outcomes can occur. In short, more examples of how faculty from different disciplines actually develop, engage, and revise interdisciplinary pedagogies with …

Drosophila Melanogaster Selection For Survival Of Bacillus Cereus Infection: Life History Trait Indirect Responses, Junjie Ma, Andrew K. Benson, Stephen D. Kachman, Zhen Hu, Lawrence G. Harshman

Department of Statistics: Faculty Publications

To study evolved resistance/tolerance in an insect model, we carried out an experimental evolution study using D. melanogaster and the opportunistic pathogen B. cereus as the agent of selection. The selected lines evolved a 3.0- to 3.3-log increase in the concentration of spores required for 50% mortality after 18–24 generations of selection. In the absence of any treatment, selected lines evolved an increase in egg production and delayed development time. The latter response could be interpreted as a cost of evolution. Alternatively, delayed development might have been a target of selection resulting in increased adult fat body function including production …

Prediction In Several Conventional Contexts, Bertrand Clarke, Jennifer Clarke

Department of Statistics: Faculty Publications

We review predictive techniques from several traditional branches of statistics. Starting with prediction based on the normal model and on the empirical distribution function, we proceed to techniques for various forms of regression and classification. Then, we turn to time series, longitudinal data, and survival analysis. Our focus throughout is on the mechanics of prediction more than on the properties of predictors.

Comment On Article By Sancetta, Bertrand Clarke

Department of Statistics: Faculty Publications

This paper makes a landmark contribution in three senses.

First, it provides many results that are fundamentally important in their own right. I refer specifically to Theorems 3 and 8. Theorem 3 treats arbitrary loss functions by breaking the integral into two terms, one It, where a difference of losses is bounded and another, IIt, where a bound on the moments of a difference of losses must be used. (All notation here is the same as the author's unless noted otherwise.) The treatment of these two terms reveals the role of the relative entropy and how the tails of the …

Effects Of Animal Diet, Manure Application Rate, And Tillage On Transport Of Microorganisms From Manure-Amended Fields, Lisa M. Durso, John Gilley, David B. Marx, Bryan Woodbury

Department of Statistics: Faculty Publications

Manure from cattle fed distillers’ grain or corn diets was applied to fields, and the fields were subjected to rainfall simulation tests. Manure was added at three rates on till and no-till plots. Correlations between microbial transport and runoff characteristics were identified. Results indicate that diet affects phage but not bacterial transport from manure-amended fields.

Variation In Sclerotinia Sclerotiorum Bean Isolates From Multisite Resistance Screening Locations, Lindsey Otto-Hanson, James R, Steadman, Rebecca Higgins, Kent M. Eskridge

Department of Statistics: Faculty Publications

There is no complete resistance to Sclerotinia sclerotiorum, cause of white mold in dry bean (Phaseolus vulgaris). Variable resistance expression is one problem in screening for improved white mold resistance. With no previous information in the literature, pathogen variation in multisite screening nurseries was evaluated as one cause of diverse resistance expression. In all, 10 isolates of S. sclerotiorum used in greenhouse screening and 146 isolates collected from nine white mold field screening nurseries in major bean production areas in the United States were compared using mycelial compatibility groupings (MCGs) and an aggressiveness test. These 10 greenhouse …

Development Of An Internet Based System For Modeling Biotin Metabolism Using Bayesian Networks, Jinglei Zhou, Dong Wang, Vicki L. Schlegel, Janos Zempleni

Department of Statistics: Faculty Publications

Biotin is an essential water-soluble vitamin crucial for maintaining normal body functions. The importance of biotin for human health has been under-appreciated but there is plenty of opportunity for future research with great importance for human health. Currently, carrying out predictions of biotin metabolism involves tedious manual manipulations. In this paper, we report the development of BiotinNet, an internet based program that uses Bayesian networks to integrate published data on various aspects of biotin metabolism. Users can provide a combination of values on the levels of biotin related metabolites to obtain the predictions on other metabolites that are not specified. …

Reference Priors For Exponential Families With Increasing Dimension, Bertrand Clarke, Subhashis Ghosal

Department of Statistics: Faculty Publications

In this article, we establish the asymptotic normality of the posterior distribution for the natural parameter in an exponential family based on independent and identically distributed data. The mode of convergence is expected Kullback-Leibler distance and the number of parameters p is increasing with the sample size n. Using this, we give an asymptotic expansion of the Shannon mutual information valid when p = pm increases at a sufficiently slow rate. The second term in the asymptotic expansion is the largest term that depends on the prior and can be optimized to give Jeffrey's prior as the reference prior in …

Weakly Positioned Nucleosomes Enhance The Transcriptional Competency Of Chromatin, Yaakov Belch, Jingyi Yang, Yang Liu, Sridhar A. Malkaram, Rong Liu, Jean-Jack M. Riethoven, Istvan Ladunga

Department of Statistics: Faculty Publications

Background: Transcription is affected by nucleosomal resistance against polymerase passage. In turn, nucleosomal resistance is determined by DNA sequence, histone chaperones and remodeling enzymes. The contributions of these factors are widely debated: one recent title claims ‘‘… DNA-encoded nucleosome organization…’’ while another title states that ‘‘histone-DNA interactions are not the major determinant of nucleosome positions.’’ These opposing conclusions were drawn from similar experiments analyzed by idealized methods. We attempt to resolve this controversy to reveal nucleosomal competency for transcription.

Methodology/Principal Findings: To this end, we analyzed 26 in vivo, nonlinked, and in vitro genome-wide nucleosome maps/replicates by new, rigorous …

Principal Component Analysis And Biochemical Characterization Of Protein And Starch Reveal Primary Targets For Improving Sorghum Grain, Joshua H. Wong, D. B. Marx, Jeff D. Wilson, Bob B. Buchanan, Peggy G. Lemaux, Jeffrey F. Pedersen

Department of Statistics: Faculty Publications

Limited progress has been made on genetic improvement of the digestibility of sorghum grain because of variability among different varieties. In this study, we applied multiple techniques to assess digestibility of grain from 18 sorghum lines to identify major components responsible for variability. We also identified storage proteins and enzymes as potential targets for genetic modification to improve digestibility. Results from principal component analysis revealed that content of amylose and total starch, together with protein digestibility (PD), accounted for 94% of variation in digestibility. Control of amylose content is understood and manageable. Up-regulation of genes associated with starch accumulation is …

Reference Priors For Exponential Families With Increasing Dimension, Bertrand S. Clarke, Subhashis Ghosal

Department of Statistics: Faculty Publications

In this article, we establish the asymptotic normality of the posterior distribution for the natural parameter in an exponential family based on independent and identically distributed data. The mode of convergence is expected Kullback-Leibler distance and the number of parameters p is increasing with the sample size n. Using this, we give an asymptotic expansion of the Shannon mutual information valid when p = pn increases at a sufficiently slow rate. The second term in the asymptotic expansion is the largest term that depends on the prior and can be optimized to give Jeffreys’ prior as the reference prior in …