Skip Navigation

Your Environment. Your Health.

Shyamal D. Peddada, Ph.D.

Biostatistics & Computational Biology Branch

Shyamal D. Peddada, Ph.D.
Shyamal D. Peddada, Ph.D.
Acting Branch Chief, Biostatistics & Computational Biology Branch and Principal Investigator
Tel 919-541-1122
111 T W Alexander Dr
Rall Building
Research Triangle Park, NC 27709

Research Summary

Shyamal Peddada, Ph.D., is the Acting Branch Chief and Principal Investigator of the Biostatistics and Computational Biology Branch and holds a secondary appointment in the NIEHS Epidemiology Branch.

We develop broadly applicable statistical methods that are motivated by applications in environmental sciences. Constraints arise naturally in many scientific problems. For example, in dose-response studies the mean response is often expected to be monotonic. Similarly, mean expression of differentially expressed genes may not vary arbitrarily across ordered conditions (e.g. stages of cancer) but may have a systematic pattern of expression over the conditions. These patterns can be nonparametrically represented using mathematical inequalities. In some cases, the parameters (and the data) may be constrained by a pyramid (or “simplex”), as in the case of microbiome data, or a circle, as in the case of circadian clock or cell-cycle gene expression data. Ignoring such additional information may not only result in inefficiency (i.e. loss of power or increase in sample size requirements), but more importantly, may potentially result in misinterpretation of the scientific data. We develop statistical methods that exploit the underlying constraints on the biological parameter and consequently our methods are not only more efficient but they also provide better interpretation of the data. Our group also develops methods that are suitable for analyzing high dimensional data obtained from genomics (e.g. microarray, RNA-Seq, etc.), quantitative high-throughput screening (qHTS) assays, microbiome, etc. Software developed by the group are freely accessible through links provided below.


The following software was developed in this research program and is freely available to download.

  • Order Restricted Inference for Oscillatory Systems (ORIOS) for Detecting Rhythmic Signals
    (R code developed by Yolanda Larriba, University of Valladolid, Valladolid, Spain).
    ORIOS is a model free order restricted inference based algorithm that detects rhythmic components (e.g. transcripts or genes) s participating in oscillatory systems such as the circadian clock. Although this software can be used for any oscillatory data, for simplicity of description, throughout this file we shall use the term “circadian clock data” rather than “oscillatory data” and “genes” in place of “components” of an oscillatory system.  The strength of model free methodology such as the order restricted inference is that, instead of using a mathematical model to describe the shape or pattern of expression, it uses mathematical inequalities to describe patterns.  Thus it is robust and is not limited by any mathematical model that may be rigid and restricted. ORIOS not only identifies rhythmic genes, it also classifies them into four typical classes of genes, called cyclical, quasi cyclical, non-flat and non-periodic, and flat, according to its signal shape.  Cyclical and quasi cyclical genes are declared as rhythmic, while non-flat and non-periodic, and flat are declared as non-rhythmic genes. Compared to some commonly used rhythmicity detection algorithms in, ORIOS has substantially higher power to detect true rhythmic genes, while also declaring substantially fewer non-rhythmic genes as rhythmic.

  • Analysis of Compositional Microbiomes (ANCOM) data
    (R-code developed by Dr. Siddhartha Mandal, Norwegian Institute of Public Health and the Shiny app developed by Dr. Casey Jelsema, Biostatistics and Computational Biology Branch, NIEHS)
    This R package is designed for comparing the abundance of individual taxa in two populations using log-ratios of abundance. This software is based on the ANCOM methodology developed in Mandal et al. (2015)
    • Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, and Peddada SD (2015). Analysis of composition of microbiomes: a novel method for studying microbial composition. Microbial Ecology in Health and Disease, 26, 1 – 7.

  • R Code for Estimating of Global Relative Order of Peak Expression Satisfied by a Set of Oscillatory Genes(28KB)
    (Programmed by Ms. Sandra Barragán, University of Valladolid, Spain):
    For a given collection of oscillatory genes (e.g. cell-cycle genes or circadian clock genes) with phase angles estimated from multiple experiments, in this software we estimate the relative order of peak expression among the genes. It contains 2 functions written in R, called Aggregation of Circular Orders (ACO), which is based on a solution to the traveling salesman problem, and Circular Local Minimization (CLM) algorithm which is used to smooth the solution obtained from ACO. To run these programs the user should first download the companion R package called \emph{isocir} from CRAN .

  • Circular FSA("/Rhythmyx/assembler/render?sys_contentid=51492&sys_revision=2&sys_variantid=639&sys_context=0&sys_authtype=0&sys_siteid=&sys_folderid=" sys_dependentvariantid="639" sys_dependentid="51492" inlinetype="rxhyperlink" rxinlineslot="103" sys_dependentid="51492" sys_siteid="" sys_folderid="")
    (Programmed by Prof. Miguel Fernandez and Ms. Sandra Barragán, University of Valladolid, Spain):
    For a given set of estimates of angular parameters, this software can be used for testing whether the corresponding angular parameters satisfy a pre-specified order around a unit circle.
    • Fernandez M, Rueda C, Peddada SD* (2012) Identification of a core set of signature cell-cycle genes whose relative order of time to peak expression is conserved across species. Nucleic Acids Research, 40,2823-32. Epub 2011/12/03. doi: 10.1093/nar/gkr1077.
    • Rueda C, Fernandez M, Peddada SD* (2009). Estimation of parameters subject to order restriction on a circle with application to estimation of phase angles of cell-cycle genes. J. Amer. Statist. Assoc., 104, 338-347.

  • R code for fitting Random Periods Model ("/Rhythmyx/assembler/render?sys_contentid=51500&sys_revision=2&sys_variantid=639&sys_context=0&sys_authtype=0&sys_siteid=&sys_folderid=" sys_dependentvariantid="639" sys_dependentid="51500" inlinetype="rxhyperlink" rxinlineslot="103" sys_dependentid="51500" sys_siteid="" sys_folderid="")
    (Programmed by Mr. Sai Cheemalapati, High School Summer Intern):
    For a given periodic time-course data (e.g. expression of cell-cycle gene or a circadian clock gene), this program fits the nonlinear random periods model (RPM) and estimates all parameters of the model, namely, the intercept, slope, amplitude, phase, period and attenuation. It also provides the information matrix associated with the estimates.
    • Liu D, Umbach D, Peddada SD, Li L, Crockett P, Weinberg C (2004). A Random-Periods Model for Expression of Cell-Cycle Genes. Proceedings of National Academy of Sciences, 101, No. 19, 7240-7245.

  • ORIOGEN 4.01 - Order Restricted Inference for Ordered Gene Expression and Multiple Pairwise Comparisons("/Rhythmyx/assembler/render?sys_contentid=35065&sys_revision=8&sys_variantid=639&sys_context=0&sys_authtype=0&sys_siteid=&sys_folderid=" sys_dependentvariantid="639" sys_dependentid="35065" inlinetype="rxhyperlink" rxinlineslot="103" sys_dependentid="35065" sys_siteid="" sys_folderid=""):
    This software is designed for comparing two or more experimental groups. There are two options available within software, with one used for analyzing ordered experimental conditions (e.g. time, dose, tumor stages, etc.). Under this option, the software can handle an independent sample case, as well as a dependent sample case (e.g. repeated measurements). The residual bootstrap methodology used in this software is robust to any underlying dependence structure. The method controls the FDR at the desired level. The second option is suitable for pairwise comparisons and is not limited to ordered experimental conditions. Thus, for any given design, the second option allows one to make all desired pairwise comparisons among the experimental groups. In addition it allows one to make directional inferences (such as up or down regulated genes etc.). The method controls for the overall mixed directional false discovery rates (mdFDR).
    • Guo W, Peddada SD* (2008). Adaptive Choice of the Number of Bootstrap Samples in Large Scale Multiple Testing. Statistical Applications in Genetics and Molecular Biology, 7 (1), Art. 13.
    • Peddada SD*, Harris S, Zajd J, Harvey E (2005). ORIOGEN: Order Restricted Inference for Ordered Gene Expression data. Bioinformatics, 21, 3933-3934.
    • Peddada SD*, Lobenhofer L, Li L, Afshari C, Weinberg C, Umbach D (2003). Gene selection and clustering for time-course and dose-response microarray experiments using order-restricted inference. Bioinformatics, 19, 834-841.

  • Constrained Linear Mixed Effects (CLME) for analyzing mixed and fixed models under inequality constraints.
    (Programmed by Dr. Casey M. Jelsema, Research Fellow, Biostatistics Branch, NIEHS):
    In many applications, such as in dose-response studies or time-course experiments, researchers are interested in testing for specific inequality constraints or patterns among the means of experimental groups. This R package is designed to test for such inequality patterns using a robust residual bootstrap based methodology which does not require the data to be normally distributed. Furthermore, this software can also handle the situation when covariates and/or random effects are present. Thus, for example, this package can be used in the context of repeated measurement designs with covariates. This package comes with a user friendly graphical interface so no programming is necessary to run this package. All the user needs to do is to provide input source of the data and select options from the interface.

Selected Publications

  1. Kaul A, Davidov O, Peddada Sd. Structural zeros in high-dimensional data with applications to microbiome studies. Biostatistics (Oxford, England) 2017 ():-. [Abstract]
  2. Hill-Burns EM , Debelius JW, Morton JT, Wissemann WT, Lewis MR , Wallen ZD, Peddada SD , Factor SA, Molho E, Zabetian CP, Knight R, Payami H. Parkinson's disease and Parkinson's disease medications have distinct signatures of the gut microbiome. Movement disorders: official journal of the Movement Disorder Society 2017 ():-. [Abstract]
  3. Mandal M , Godfrey K , McDonald D , van Treuren W , Bjørnholt JV , Midvedt T, Moen B , Rudi K , Knight R , Peddada SD , Eggesbø M. Fat and vitamin intakes during pregnancy have stronger relations with a pro-inflammatory maternal microbiota than does carbohydrate intake. Microbiome 2016 4(1):55-. [Abstract]
  4. Larriba Y, Rueda C, Fernandez MA, Peddada SD. Order restricted inference for oscillatory systems for detecting rhythmic signals. Nucleic acids research 2016 ():-. [Abstract]
  5. Hayes S, Pandiri AR, Ton TV, Hong HH, Clayton NP, Shockley KR, Peddada SD, Gerrish K, Wyde M, Sills RC, Hoenerhoff MJ. (2016) Renal Cell Carcinomas in Vinylidene Chloride-exposed Male B6C3F1 Mice Are Characterized by Oxidative Stress and TP53 Pathway Dysregulation. Toxicologic pathology 44(1):71-87. [Abstract]
  6. Yamashita, H, Hoenerhoff, MJ, Peddada, SD, Sills, RC, Pandiri, AR. Chemical Exacerbation of Light-induced Retinal Degeneration in F344/N Rats in National Toxicology Program Rodent Bioassays. Toxicologic pathology 2016 44(6):892-903. [Abstract]
  7. Reese SE, Zhao S, Wu M, Joubert BR, Parr CH, Haberg S, Ueland PM, Nilsen RM, Midttun O, Vollset SE, Peddada S, Nystad W, London SJ. DNA Methylation Score as a Biomarker in Newborns for Sustained Maternal Smoking during Pregnancy. Environmental health perspectives 2016 ():-. [Abstract]
  8. Harvey JB, Hong HH, Bhusari S, Ton TV, Wang Y, Foley JF, Peddada SD, Hooth M, DeVito M, Nyska A, Pandiri AR, Hoenerhoff MJ. F344/NTac Rats Chronically Exposed to Bromodichloroacetic Acid Develop Mammary Adenocarcinomas With Mixed Luminal/Basal Phenotype and Tgfβ Dysregulation. Veterinary pathology 2016 53(1):170-181. [Abstract]
  9. Rueda C, Fernandez MA, Barragan S, Mardia KV, Peddada SD. Circular piecewise regression with applications to cell-cycle data. Biometrics 2016 ():-. [Abstract]
  10. Jelsema C, Peddada SD (2016). An R Package for Linear Mixed Effects Models under Inequality Constraints. Journal of Statistical Software (In press).
  11. Rebera SO, Siebler PH, Donner NC, Morton JT, Smith DG, Kopelman JM, Lowe KR, Campbell K, Fox JH, Hassell JE, Greenwood BN, Janscha C, Lechner A, Uschold-Schmidt N, Füchsl AM, Langgartner D, Walker FR, Hale MW, Perez GL, Treuren WV, González A, Halweg-Edwards AL, Fleshner M, Raison CL, Rook GA, Peddada SD, Knight R, Lowry CA (2016). Immunization with a heat-killed preparation of the environmental bacterium Mycobacterium vaccae promotes stress resilience in mice. Proc. National Acad. Sci. (In press).
  12. Joubert B, den Dekker H, Felix J, Bohlin J, Ligthart S, Beckett E, Tiemeier H, van Meurs J, Uitterlinden A, Hofman A, Haberg S, Reese S, Peters M, Andreassen B, Steegers E, Nilsen R, Vollset S, Midttun O, Ueland P, Franco O, Dehghan A, de Jongste J, Wu M, Wang T, Peddada S, Jaddoe V, Nystad W, Duijts L, London S. Maternal plasma folate impacts differential DNA methylation in an epigenome-wide meta-analysis of newborns. Nat Commun. 2016; Feb 7:10577. Epub 11 Feb 2016. PMID: 26861414. PMCID: PMC4749955. [Abstract]
  13. Grandhi A, Guo W, Peddada SD. A multiple testing procedure for multi-dimensional pairwise comparisons with application to gene expression studies. BMC bioinformatics 2016 (17)104. [Full Text]
  14. Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, Peddada SD (2015). Analysis of composition of microbiomes: a novel method for studying microbial composition. Microbial Ecology in Health and Disease, 26, 1 – 7. [Abstract]
  15. Davidov O, Peddada S. (2013). Testing for the Multivariate Stochastic Order among Ordered Experimental Groups with Application to Dose-Response Studies. Biometrics 69(4):982-990. [Abstract]
  16. White RA, Bjørnholt J, Baird DD, Midtvedt T, Harris JR, Pagano M, Hide W, Rudi K, Moen B, Iszatt N, Peddada SD, Eggesbø M (2013). Novel Developmental Analyses Identify Longitudinal Patterns of Early Gut Microbiota that Affect Infant Growth. PLoS Computational Biology 9(5):e1003042. [Abstract]
  17. Davidov O, Peddada SD (2013). The linear stochastic order and directed inference for multivariate ordered distributions. Annals of Statistics 41(1):1-40.  [Abstract]
  18. Lim C, Sen PK, Peddada SD. (2013). Robust Analysis of High Throughput Screening (HTS) Assay Data. Technometrics 55(2):150-160. [Abstract]
  19. Fernandez M, Rueda C, Peddada SD (2012). Identification of a core set of signature cell-cycle genes whose relative order of time to peak expression is conserved across species. Nucleic Acids Research, 40,2823-32, doi: 10.1093/nar/gkr1077. [Abstract]
  20. Davidov O, Peddada SD (2011). Order restricted inference for multivariate binary data with application to toxicology. J. Amer. Statist. Assoc., 106, 1394-1404.
  21. Guo W, Sarkar SK, Peddada, SD (2010). Controlling False Discoveries in Multidimensional Directional Decisions, with Applications to Gene Expression Data on Ordered Categories. Biometrics, 66, 485 - 492. [Abstract]
  22. Peddada, SD, Dinse, G and Kissling, G (2007). Incorporating Historical Control Data When Comparing Tumor Incidence Rates. J. Amer. Statist. Assoc., 102, 1212-1220.
  23. Peddada SD, Lobenhofer L, Li L, Afshari C, Weinberg C and Umbach D. (2003). Gene Selection and Clustering for Time-course and Dose-response Microarray Experiments using Order-restricted Inference. Bioinformatics, 19, 834-841. [Abstract]
  24. Hwang JTG and Peddada SD (1994). Confidence Interval Estimation Subject to Order Restrictions. Annals of Statistics, 22, 67-93.