Advanced Data Analysis - Lecture Notes

Lecture notes for Advanced Data Analysis (ADA1 Stat 427/527 and ADA2 Stat 428/528), Department of Mathematics and Statistics, University of New Mexico, Fall 2016-Spring 2017. Additional material including RMarkdown templates for in-class and homework exercises, datasets, R code, and video lectures are available on the course websites: https://statacumen.com/teaching/ada1 and https://statacumen.com/teaching/ada2 . Contents I ADA1: Software 0 Introduction to R, Rstudio, and ggplot II ADA1: Summaries and displays, and one-, two-, and many-way tests of means 1 Summarizing and Displaying Data 2 Estimation in One-Sample Problems 3 Two-Sample Inferences 4 Checking Assumptions 5 One-Way Analysis of Variance III ADA1: Nonparametric, categorical, and regression methods 6 Nonparametric Methods 7 Categorical Data Analysis 8 Correlation and Regression IV ADA1: Additional topics 9 Introduction to the Bootstrap 10 Power and Sample size 11 Data Cleaning V ADA2: Review of ADA1 1 R statistical software and review VI ADA2: Introduction to multiple regression and model selection 2 Introduction to Multiple Linear Regression 3 A Taste of Model Selection for Multiple Regression VII ADA2: Experimental design and observational studies 4 One Factor Designs and Extensions 5 Paired Experiments and Randomized Block Experiments 6 A Short Discussion of Observational Studies VIII ADA2: ANCOVA and logistic regression 7 Analysis of Covariance: Comparing Regression Lines 8 Polynomial Regression 9 Discussion of Response Models with Factors and Predictors 10 Automated Model Selection for Multiple Regression 11 Logistic Regression IX ADA2: Multivariate Methods 12 An Introduction to Multivariate Methods 13 Principal Component Analysis 14 Cluster Analysis 15 Multivariate Analysis of Variance 16 Discriminant Analysis 17 Classificationhttps://digitalrepository.unm.edu/unm_oer/1002/thumbnail.jp

Individual differences in toddlers' social understanding and prosocial behavior: Disposition or socialization?

We examined how individual differences in social understanding contribute to variability in early-appearing prosocial behavior. Moreover, potential sources of variability in social understanding were explored and examined as additional possible predictors of prosocial behavior. Using a multi-method approach with both observed and parent-report measures, 325 children aged 18-30 months were administered measures of social understanding (e.g., use of emotion words; self-understanding), prosocial behavior (in separate tasks measuring instrumental helping, empathic helping, and sharing, as well as parent-reported prosociality at home), temperament (fearfulness, shyness, and social fear), and parental socialization of prosocial behavior in the family. Individual differences in social understanding predicted variability in empathic helping and parent-reported prosociality, but not instrumental helping or sharing. Parental socialization of prosocial behavior was positively associated with toddlers' social understanding, prosocial behavior at home, and instrumental helping in the lab, and negatively associated with sharing (possibly reflecting parents' increased efforts to encourage children who were less likely to share). Further, socialization moderated the association between social understanding and prosocial behavior, such that social understanding was less predictive of prosocial behavior among children whose parents took a more active role in socializing their prosociality. None of the dimensions of temperament was associated with either social understanding or prosocial behavior. Parental socialization of prosocial behavior is thus an important source of variability in children's early prosociality, acting in concert with early differences in social understanding, with different patterns of influence for different subtypes of prosocial behavior

Comparison of Temperature Measurements in the Middle Atmosphere by Satellite with Profiles Obtained by Meteorological Rockets

Measurements using the inflatable falling sphere technique have occasionally been used to obtain temperature results from density data and thereby provide comparison with temperature profiles obtained by satellite sounders in the mesosphere and stratosphere. To insure density measurements within narrow time frames and close in space, the inflatable falling sphere is launched within seconds of the nearly overhead satellite pass. Sphere measurements can be used to validate remotely measured temperatures but also have the advantage of measuring small-scale atmospheric features. Even so, with the dearth of remaining falling spheres available (the manufacture of these systems has been discontinued), it may be time to consider whether the remote measurements are mature enough to stand alone. Three field studies are considered, one in 2003 from Northern Sweden, and two in 2010 from the vicinity of Kwajalein Atoll in the South Pacific and from Barking Sands, Hawaii. All three sites are used to compare temperature retrievals between satellite and in situ falling spheres. The major satellite instruments employed are SABER, MLS, and AIRS. The comparisons indicate that remotely measured temperatures mimic the sphere temperature measurements quite well. The data also confirm that satellite retrievals, while not always at the exact location required for detailed studies in space and time, compare sufficiently well to be highly useful. Although the falling sphere will provide a measurement at a specific location and time, satellites only pass a given location daily or less frequently. This report reveals that averaged satellite measurements can provide temperatures and densities comparable to those obtained from the falling sphere, thereby providing a reliable measure of global temperatur

Multi-center implementation of rapid whole genome sequencing provides additional evidence of its utility in the pediatric inpatient setting

ObjectiveMulti-center implementation of rapid whole genome sequencing with assessment of the clinical utility of rapid whole genome sequencing (rWGS), including positive, negative and uncertain results, in admitted infants with a suspected genetic disease.Study designrWGS tests were ordered at eight hospitals between November 2017 and April 2020. Investigators completed a survey of demographic data, Human Phenotype Ontology (HPO) terms, test results and impacts of results on clinical care.ResultsA total of 188 patients, on general hospital floors and intensive care unit (ICU) settings, underwent rWGS testing. Racial and ethnic characteristics of the tested infants were broadly representative of births in the country at large. 35% of infants received a diagnostic result in a median of 6 days. The most common HPO terms for tested infants indicated an abnormality of the nervous system, followed by the cardiovascular system, the digestive system, the respiratory system and the head and neck. Providers indicated a major change in clinical management because of rWGS for 32% of infants tested overall and 70% of those with a diagnostic result. Also, 7% of infants with a negative rWGS result and 23% with a variant of unknown significance (VUS) had a major change in management due to testing.ConclusionsOur study demonstrates that the implementation of rWGS is feasible across diverse institutions, and provides additional evidence to support the clinical utility of rWGS in a demographically representative sample of admitted infants and includes assessment of the clinical impact of uncertain rWGS results in addition to both positive and negative results

A Bayesian Approach to Estimate OJD Prevalence From Pooled Fecal Samples of Variable Pool Size

This paper describes a Bayesian approach to prevalence estimation based on pooled samples that accommodates variation in pool size and adjusts for test imperfection. A logistic model was developed for pooled fecal culture (PFC) sensitivity as a function of pool size and a logistic mixed model for ovine Johne’s disease (OJD) prevalence as a function of covariates that were found significant in a recent OJD risk factor study conducted in Australia. Available data on these factors and prior information about prevalence and sensitivity were incorporated into a Bayesian model to estimate OJD prevalence from PFC data. Overall, posterior cohort OJD prevalence was estimated to be 0.16 (range of prevalences across cohorts 0.002 to 0.72). The average prevalence was higher in wethers than ewes. PFC sensitivities for pool sizes 10, 30 and 50 were estimated to be 0.91 (95% probability intervals 0.80, 0.96), 0.85 (0.80, 0.90) and 0.77 (0.65, 0.88), respectively. Posterior specificity of PFC was almost perfect though based primarily on the prior. Results suggest the Bayesian model successfully estimated the animal-level prevalence after accounting for variable pool size and imperfect test parameters. The method can be easily adapted for other conditions and diseases where pooled samples are collected. WinBugs code for the article is available online. Key Words: Diagnostic test Mycobacterium Ovine Johne’s disease Paratuberculosis Prevalence Sensitivit

Phenotypic Plasticity and Contemporary Evolution in Introduced Populations: Evidence from Translocated Populations of White Sands Pupfish (Cyrpinodon tularosa)

Contemporary evolution has been shown in a few studies to be an important component of colonization ability, but seldom have researchers considered whether phenotypic plasticity facilitates directional evolution from the invasion event. In the current study, we evaluated body shape divergence of the New Mexico State-threatened White Sands pupfish (Cyprinodon tularosa) that were introduced to brackish, lacustrine habitats at two different time in the recent past (approximately 30 years and 1 year previously) from the same source population (saline river environment). Pupfish body shape is correlated with environmental salinity: fish from saline habitats are characterized by slender body shapes, whereas fish from fresher, yet brackish springs are deep-bodied. In this study, lacustrine populations consisted of an approximately 30-year old population and several 1-year old populations, all introduced from the same source. The body shape divergence of the 30-year old population was significant and greater than any of the divergences of the 1-year old populations (which were for the most part not significant). Nonetheless, all body shape changes exhibited body deepening in less saline environments. We conclude that phenotypic plasticity potentially facilitates directional evolution of body deepening for introduced pupfish populations

Single Molecule Analysis Research Tool (SMART): An Integrated Approach for Analyzing Single Molecule Data

Single molecule studies have expanded rapidly over the past decade and have the ability to provide an unprecedented level of understanding of biological systems. A common challenge upon introduction of novel, data-rich approaches is the management, processing, and analysis of the complex data sets that are generated. We provide a standardized approach for analyzing these data in the freely available software package SMART: Single Molecule Analysis Research Tool. SMART provides a format for organizing and easily accessing single molecule data, a general hidden Markov modeling algorithm for fitting an array of possible models specified by the user, a standardized data structure and graphical user interfaces to streamline the analysis and visualization of data. This approach guides experimental design, facilitating acquisition of the maximal information from single molecule experiments. SMART also provides a standardized format to allow dissemination of single molecule data and transparency in the analysis of reported data