Comparing heritability estimates for twin studies + : & Mary Ellen Koran. Tricia Thornton-Wells. Bennett Landman
|
|
- Neil Mitchell
- 6 years ago
- Views:
Transcription
1 Comparing heritability estimates for twin studies + : & Mary Ellen Koran Tricia Thornton-Wells Bennett Landman January 20, 2014
2 Outline Motivation Software for performing heritability analysis Simulations Twin Study (SOLAR & OpenMX) Family Study (SOLAR) Live Demo Conclusion January 20, 2014 IIGC 2014 SOLAR Workshop 2
3 Outline Motivation Software for performing heritability analysis Simulations Twin Study (SOLAR & OpenMX) Family Study (SOLAR) Live Demo Conclusion January 20, 2014 IIGC 2014 SOLAR Workshop 3
4 Heritability (h 2 ) crash course A first question to ask before embarking on a genetic study: is a trait heritable? BEFORE GENOTYPING h 2 = degree of genetic association As h 2 -> 1, a trait increasingly heritable January 20, 2014 IIGC 2014 SOLAR Workshop 4
5 Cortical gray matter density is heritable January 20, 2014 IIGC 2014 SOLAR Workshop 5
6 Cortical area thickness is heritable January 20, 2014 IIGC 2014 SOLAR Workshop 6
7 Resting state fmri connectivity is heritable January 20, 2014 IIGC 2014 SOLAR Workshop 7
8 White matter microstructure is heritable January 20, 2014 IIGC 2014 SOLAR Workshop 8
9 We are just getting started Software for estimating heritability are freely available and functional on commodity hardware. Heritabilities can be estimated using diverse family structures (not just twins). We are just beginning to understand the role of brain phenotypes with substantial degrees of heritability. Missing heritability January 20, 2014 IIGC 2014 SOLAR Workshop 9
10 Choosing study designs Twin studies: h 2 = 2(r MZ r DZ ) Family studies: variance components method January 20, 2014 IIGC 2014 SOLAR Workshop 10
11 Estimating h 2 with variance components Family studies: σ T 2 = σ G 2 + σ E 2 σ G 2 = σ a 2 + σ d 2 Trait variation = genetic + environment Genetic variation = additive + dominant σ E 2 = σ c 2 + σ e 2 Environmental variation = familial/household + random/individual ACE model commonly used σ T 2 σ a 2 + σ c 2 + σ e 2 January 20, 2014 IIGC 2014 SOLAR Workshop 11
12 To wrap up. ACE model A σ a 2 σ a 2 + σ c 2 + σ e 2 C σ c 2 σ a 2 + σ c 2 + σ e 2 E σ e 2 σ a 2 + σ c 2 + σ e 2 Additive genetics Common environmental effects Individual environmental effects or Error h 2 = A = σ a 2 / σ T 2 January 20, 2014 IIGC 2014 SOLAR Workshop 12
13 Outline Motivation Software for performing heritability analysis Simulations Results Twin Study (SOLAR & OpenMX) Family Study (SOLAR) Live Demo Conclusion January 20, 2014 IIGC 2014 SOLAR Workshop 13
14 Estimating Heritability Sequential oligogenic linkage analysis routines Multipoint quantitative-trait linkage analysis in general pedigrees. Almasy, Blangero. Am J Hum Genet Tcl-based Structural equation modeling OpenMx: An Open Source Extended Structural Equation Modeling Framework. Boker, Neale, et.al. Psychometrika R-based January 20, 2014 IIGC 2014 SOLAR Workshop 14
15 SOLAR twins]$ solar SOLAR version (Experimental), last updated on December 31, 2013 Copyright (c) Texas Biomedical Research Institute Enter help for help, exit to exit, doc to browse documentation. solar> load pedigree ped2.ped solar> load phenotypes gen1.csv solar> model new solar> Unloading current pedigree data... solar> trait kindgen1 solar> covar age solar> outdir /data/solar_mek/ solar> polygenic Loading pedigree data from the file /data/solar_mek/10_18_2013_openmx_solar/twins/pair2/ped2.ped... solar> /data/solar_mek/10_18_2013_openmx_solar/twins/pair2/a_0/a_0_c_0/gen1.csv: ID age kindgen1 FAMID solar> solar> solar> solar> solar> ********************************************************************** * Maximize sporadic model * January 20, 2014 IIGC 2014 SOLAR Workshop 15
16 Inputs to SOLAR Pedigree File Phenotype File January 20, 2014 IIGC 2014 SOLAR Workshop 16
17 Output of SOLAR Pedigree: /data/solar_mek/livedemo/twins4/pair20/ped20.ped Phenotypes: /data/solar_mek/livedemo/twins4/pair20/a_50/a_50_c_0/gen10.csv Trait: kindgen10 Individuals: 80 H2r is p = (Significant) H2r Std. Error: Proportion of Variance Due to All Final Covariates Is Loglikelihoods and chi's are in /data/solar_mek/livedemo/twins4/pair20/a_50/a_50_c_0/gen10out//polyg enic.logs.out Best model is named poly and null0 Final models are named poly, spor, nocovar Residual Kurtosis is , within normal range January 20, 2014 IIGC 2014 SOLAR Workshop 17
18 OpenMX January 20, 2014 IIGC 2014 SOLAR Workshop 18
19 Inputs to OpenMX Pedigree Coded in SEM Phenotype Variables January 20, 2014 IIGC 2014 SOLAR Workshop 19
20 Outputs of OpenMX Results stored into R variables January 20, 2014 IIGC 2014 SOLAR Workshop 20
21 Cross Platform Comparison Closed source/black box Easily extendible to complex families Easier to specify complicated pedigrees Single threaded execution Complete list of command descriptions Multi-platform (tcl-based) Open source Difficult to extend to complex families Requires user input for starting values of A, C, E Tends to not converge more often than SOLAR Integrates with R s cluster processing Well managed, active forum, lots of help documents on a wiki January 20, 2014 IIGC 2014 SOLAR Workshop 21
22 Outline Motivation Software for performing heritability analysis Simulations Twin Study (SOLAR & OpenMX) Family Study (SOLAR) Live Demo Conclusion January 20, 2014 IIGC 2014 SOLAR Workshop 22
23 Why use simulation? Debug Performance assessment Method comparison January 20, 2014 IIGC 2014 SOLAR Workshop 23
24 Why simulate? : Debug How do we know if the software is properly Installed (documentation), Configured (system / platform), Used (user error) Isolate errors, crashes, convergence issues Freely sharable design Reproducible results without large data exchange Commoditize heritability analysis January 20, 2014 IIGC 2014 SOLAR Workshop 24
25 Why simulate? : Performance New studies operate on edge cases that may have not been expected by the designers. Are there coding bugs? Are there theoretical bugs? Are crashes due to poorly formatted data? Confirm that a generative model for the study is well understood. Produce empirical power curves for study planning. Assess estimator performance January 20, 2014 IIGC 2014 SOLAR Workshop 25
26 Alternative for Power Estimation I h2power Purpose: Perform heritability power calculations This command performs a power calculation for the currently loaded pedigree, with the following default assumptions: (1) the trait to be studied is either quantitative or dichotomous (e.g. affected/unaffected) (2) the trait to be studied is influenced by additive genetics (3) all pedigree members will be phenotyped for the trait to be studied (unless the -data option is used to exclude those individuals who will not have phenotypic data; see the description of this option below) Also see power. January 20, 2014 IIGC 2014 SOLAR Workshop 26
27 Alternative for Power Estimation II Theoretical power models are available using spectral properties of the pedigree structure. John Blangero, et al. Chapter One - A Kernel of Truth : Statistical Advances in Polygenic Variance Component Models for Complex Human Pedigrees. Advances in Genetics Volume January 20, 2014 IIGC 2014 SOLAR Workshop 27
28 Why simulate? : Comparison Simulations provide scalable access to many (virtual) subjects with known genetic association structure Measure accuracy (bias) of heritability estimators with ground truth Measure precision (variance) of heritability estimators Mean Squared Error = Bias 2 + Variance. Quantify absolute and relative accuracy January 20, 2014 IIGC 2014 SOLAR Workshop 28
29 Modeling family structures MZ and DZ Twins Nuclear Families (Quartets) Grand- Nuclear Families (Octets) January 20, 2014 IIGC 2014 SOLAR Workshop 29
30 Quantitative Trait Simulation Pedigrees Created subjects A simulated C simulated Phenotype files created y = X β + N(0, 2ΦA + γc + 1 E) y = simulated phenotype X = simulated covariate ( age ) β= arbitrary coefficient (.005) N= noise dependent on: A, C, E = specified by user Φ = expected fraction of genome shared between subjects (MZ =.5, DZ =.25, parent-child=.25) γ = assumed common environment between relationship pairs (DZ, MZ, siblings = 1) January 20, 2014 IIGC 2014 SOLAR Workshop 30
31 Example Pedigree Simulation Simulating family structures and phenotypes Twin 1 Twin 2 MZ twins Twin 1 A + C + E A + C Twin 2 A + C A + C + E January 20, 2014 IIGC 2014 SOLAR Workshop 31
32 Example Pedigree Simulation Simulating family structures and phenotypes Twin 1 Twin 2 MZ twins Twin 1 A + C + E A + C Twin 2 A + C A + C + E Twin 1 Twin 2 DZ twins Twin 1 A + C + E.5*A + C Twin 2.5*A + C A + C + E January 20, 2014 IIGC 2014 SOLAR Workshop 32
33 Outline Motivation Software for performing heritability analysis Simulations Twin Study (SOLAR & OpenMX) Family Study (SOLAR) Live Demo Conclusion January 20, 2014 IIGC 2014 SOLAR Workshop 33
34 Twin Study: Cross Platform Comparison January 20, 2014 IIGC 2014 SOLAR Workshop 34
35 Cross Platform Comparison Simulated quantitative phenotype with the ACE Model and covariate of simulated age in twins Pedigrees Created subjects MZ : DZ = 1:1 A simulated C simulated Phenotype files created January 20, 2014 IIGC 2014 SOLAR Workshop 35
36 Family Study of Heritability Estimates with SOLAR Heritability (A) Estimates Simulated Heritability (A) % Number of Subjects Pedigrees Created subjects A simulated C simulated Phenotype files created January 20, 2014 IIGC 2014 SOLAR Workshop 36
37 Cross Platform Comparison-Twin Study C = 0 Simulated Heritability (A) % Bias in Heritability Estimates (A) January 20, 2014 IIGC 2014 SOLAR Workshop 37
38 Cross Platform Comparison-Twin Study C = 0 Simulated Heritability (A) % Bias in Heritability Estimates (A) January 20, 2014 IIGC 2014 SOLAR Workshop 38
39 Cross Platform Comparison-Twin Study C = 0 Heritability (A) % = abs Bias OpenMX abs Bias SOLAR Blue = OpenMX more biased Red = SOLAR more biased January 20, 2014 IIGC 2014 SOLAR Workshop 39
40 Cross Platform Comparison-Twin Study Heritability (A) Estimates C = 0 Bias in A Simulated Heritability (A) % Variance in A January 20, 2014 IIGC 2014 SOLAR Workshop 40
41 Cross Platform Comparison-Twin Study Heritability (A) Estimates C = 0.30 Bias in A Simulated Heritability (A) % Variance in A January 20, 2014 IIGC 2014 SOLAR Workshop 41
42 Cross Platform Comparison-Twin Study Heritability (A) Estimates C = 0.50 Bias in A Simulated Heritability (A) % Variance in A January 20, 2014 IIGC 2014 SOLAR Workshop 42
43 Cross Platform Comparison-Twin Study Heritability (A) Estimates C = 0.70 Bias in A Simulated Heritability (A) % Variance in A January 20, 2014 IIGC 2014 SOLAR Workshop 43
44 Cross Platform Comparison-Twin Study C = 0 C = 0.30 C = 0.50 C = 0.70 Simulated Heritability (A) % = abs Bias OpenMX abs Bias SOLAR Blue = OpenMX more biased Red = SOLAR more biased January 20, 2014 IIGC 2014 SOLAR Workshop 44
45 Outline Motivation Software for performing heritability analysis Simulations Twin Study (SOLAR & OpenMX) Family Study (SOLAR) Live Demo Conclusion January 20, 2014 IIGC 2014 SOLAR Workshop 45
46 Family Study of Heritability Estimates with SOLAR January 20, 2014 IIGC 2014 SOLAR Workshop 46
47 Family Study of Heritability Estimates with SOLAR Heritability (A) Estimates C = 0 Bias in A Simulated Heritability (A) % Variance in A January 20, 2014 IIGC 2014 SOLAR Workshop 47
48 Family Study of Heritability Estimates with SOLAR Heritability (A) Estimates C = 0.30 Bias in A Simulated Heritability (A) % Variance in A January 20, 2014 IIGC 2014 SOLAR Workshop 48
49 Family Study of Heritability Estimates with SOLAR Heritability (A) Estimates C = 0.50 Bias in A Simulated Heritability (A) % Variance in A January 20, 2014 IIGC 2014 SOLAR Workshop 49
50 Family Study of Heritability Estimates with SOLAR Heritability (A) Estimates C = 0.70 Bias in A Simulated Heritability (A) % Variance in A January 20, 2014 IIGC 2014 SOLAR Workshop 50
51 Heritability Estimates with Family Data If you don t know common environment contribution (C) nuclear, grand-nuclear family study in SOLAR twin study in OpenMX If you know C = 0 twin family study in SOLAR If you know C > 0 nuclear, grand-nuclear family study in SOLAR twin family study in OpenMX January 20, 2014 IIGC 2014 SOLAR Workshop 51
52 Outline Motivation Software for performing heritability analysis Simulations Twin Study (SOLAR & OpenMX) Family Study (SOLAR) Live Demo Conclusion January 20, 2014 IIGC 2014 SOLAR Workshop 52
53 Platform Scientific Linux bit OS SOLAR (Experimental), last updated on December 31, 2013 R **NOT 3.X.X** OpenMX January 20, 2014 IIGC 2014 SOLAR Workshop 53
54 Simulation Run Times Simulation Single Core 1.6 GHz OpenMX (including simulations) Scientific Linux 6.4 with Kernel MacBook Pro with VMWare Fusion 7 SOLAR AE Twins 155 min 266 min ACE Twins 242 min 382 min Nuclear Families Grand Nuclear Families min 1803 min January 20, 2014 IIGC 2014 SOLAR Workshop 54
55 January 20, 2014 IIGC 2014 SOLAR Workshop 55
56 Outline Motivation Software for performing heritability analysis Simulations Twin Study (SOLAR & OpenMX) Family Study (SOLAR) Live Demo Conclusion January 20, 2014 IIGC 2014 SOLAR Workshop 56
57 Does any of this matter? h ±0.11 January 20, 2014 IIGC 2014 SOLAR Workshop 57
58 Cross Platform Comparison-Twin Study C = 0 Simulated Heritability (A) % Bias in Heritability Estimates (A) January 20, 2014 IIGC 2014 SOLAR Workshop 58
59 Cross Platform Comparison-Twin Study C = 0 Heritability (A) % = abs Bias OpenMX abs Bias SOLAR Blue = OpenMX more biased Red = SOLAR more biased January 20, 2014 IIGC 2014 SOLAR Workshop 59
60 Side by Side Comparison SOLAR OpenMX January 20, 2014 IIGC 2014 SOLAR Workshop 60
61 Side by Side Comparison SOLAR OpenMX January 20, 2014 IIGC 2014 SOLAR Workshop 61
62 Commonality of the common effect Phenotype A C White matter tracts (Whole brain FA) White Matter (neonatal regional metrics) ~0.5 Low Kochunov, et.al., Neuroimage,2010 Cortical Thickness ~0.7 <0.10 Kremen, et.al. Neuroimage, 2010 Brain Volume ~0.7 <0.05 Kremen, et.al. Neuroimage, 2010 ~0-0.8 ~ Geng, et al., Twin Res Hum Genet. 2012
63 e-science Objectives Provide a freely available platform for learning Virtual Machine available on NITRC Quantify differences between OpenMX and SOLAR Studies both AE and ACE models Illustrate integration with other analysis environments Currently, we are using R. In continuing work, integrate SOLAR s nifti support with pipeline platforms January 20, 2014 IIGC 2014 SOLAR Workshop 63
64 Thank you. Neda Jahanshad Peter Kochunov Tom Nichols Paul Thompson John Blangero David C. Glahn And the ENIGMA DTI working group. MASI Lab Fall 2013 NIH/NIBIB R01 EB January 20, 2014 IIGC 2014 SOLAR Workshop 64
Imaging Genetics: Heritability, Linkage & Association
Imaging Genetics: Heritability, Linkage & Association David C. Glahn, PhD Olin Neuropsychiatry Research Center & Department of Psychiatry, Yale University July 17, 2011 Memory Activation & APOE ε4 Risk
More informationMendelian & Complex Traits. Quantitative Imaging Genomics. Genetics Terminology 2. Genetics Terminology 1. Human Genome. Genetics Terminology 3
Mendelian & Complex Traits Quantitative Imaging Genomics David C. Glahn, PhD Olin Neuropsychiatry Research Center & Department of Psychiatry, Yale University July, 010 Mendelian Trait A trait influenced
More informationNon-parametric methods for linkage analysis
BIOSTT516 Statistical Methods in Genetic Epidemiology utumn 005 Non-parametric methods for linkage analysis To this point, we have discussed model-based linkage analyses. These require one to specify a
More informationWe expand our previous deterministic power
Power of the Classical Twin Design Revisited: II Detection of Common Environmental Variance Peter M. Visscher, 1 Scott Gordon, 1 and Michael C. Neale 1 Genetic Epidemiology, Queensland Institute of Medical
More informationResemblance between Relatives (Part 2) Resemblance Between Relatives (Part 2)
Resemblance Between Relatives (Part 2) Resemblance of Full-Siblings Additive variance components can be estimated using the covariances of the trait values for relatives that do not have dominance effects.
More informationIt is shown that maximum likelihood estimation of
The Use of Linear Mixed Models to Estimate Variance Components from Data on Twin Pairs by Maximum Likelihood Peter M. Visscher, Beben Benyamin, and Ian White Institute of Evolutionary Biology, School of
More informationEstimating genetic variation within families
Estimating genetic variation within families Peter M. Visscher Queensland Institute of Medical Research Brisbane, Australia peter.visscher@qimr.edu.au 1 Overview Estimation of genetic parameters Variation
More informationQuantitative Trait Analysis in Sibling Pairs. Biostatistics 666
Quantitative Trait Analsis in Sibling Pairs Biostatistics 666 Outline Likelihood function for bivariate data Incorporate genetic kinship coefficients Incorporate IBD probabilities The data Pairs of measurements
More informationDiscontinuous Traits. Chapter 22. Quantitative Traits. Types of Quantitative Traits. Few, distinct phenotypes. Also called discrete characters
Discontinuous Traits Few, distinct phenotypes Chapter 22 Also called discrete characters Quantitative Genetics Examples: Pea shape, eye color in Drosophila, Flower color Quantitative Traits Phenotype is
More informationBivariate Genetic Analysis
Bivariate Genetic Analysis Boulder Workshop 2018 Hermine H. Maes, Elizabeth Prom-Wormley, Lucia Colondro Conde, et. al. Questions Univariate Analysis: What are the contributions of additive genetic, dominance/shared
More informationSEMs for Genetic Analysis
SEMs for Genetic Analysis Session 8, Lecture 4 2/3/06 Basic Genetic Analysis of Twin Samples Using SEMs SEMs and variance explained Motivating example Simple correlation test The heritability concept Structural
More informationChapter 2. Linkage Analysis. JenniferH.BarrettandM.DawnTeare. Abstract. 1. Introduction
Chapter 2 Linkage Analysis JenniferH.BarrettandM.DawnTeare Abstract Linkage analysis is used to map genetic loci using observations on relatives. It can be applied to both major gene disorders (parametric
More informationInteraction of Genes and the Environment
Some Traits Are Controlled by Two or More Genes! Phenotypes can be discontinuous or continuous Interaction of Genes and the Environment Chapter 5! Discontinuous variation Phenotypes that fall into two
More informationOUTLIER SUBJECTS PROTOCOL (art_groupoutlier)
OUTLIER SUBJECTS PROTOCOL (art_groupoutlier) Paul K. Mazaika 2/23/2009 Outlier subjects are a problem in fmri data sets for clinical populations. This protocol and program are a method to identify outlier
More informationPerforming. linkage analysis using MERLIN
Performing linkage analysis using MERLIN David Duffy Queensland Institute of Medical Research Brisbane, Australia Overview MERLIN and associated programs Error checking Parametric linkage analysis Nonparametric
More informationImaging Genetics Center, Laboratory of Neuro Imaging, Department of Neurology, UCLA School of Medicine, Los Angeles, CA
Bivariate Genome-Wide Association Study of Genetically Correlated Neuroimaging Phenotypes from DTI and MRI through a Seemingly Unrelated Regression Model Neda Jahanshad 1,**, Priya Bhatt 1,**, Derrek P.
More informationGenomewide Linkage of Forced Mid-Expiratory Flow in Chronic Obstructive Pulmonary Disease
ONLINE DATA SUPPLEMENT Genomewide Linkage of Forced Mid-Expiratory Flow in Chronic Obstructive Pulmonary Disease Dawn L. DeMeo, M.D., M.P.H.,Juan C. Celedón, M.D., Dr.P.H., Christoph Lange, John J. Reilly,
More informationFor more information about how to cite these materials visit
Author(s): Kerby Shedden, Ph.D., 2010 License: Unless otherwise noted, this material is made available under the terms of the Creative Commons Attribution Share Alike 3.0 License: http://creativecommons.org/licenses/by-sa/3.0/
More informationOutline. Introduction GitHub and installation Worked example Stata wishes Discussion. mrrobust: a Stata package for MR-Egger regression type analyses
mrrobust: a Stata package for MR-Egger regression type analyses London Stata User Group Meeting 2017 8 th September 2017 Tom Palmer Wesley Spiller Neil Davies Outline Introduction GitHub and installation
More informationHeritability of surface area and cortical thickness: a comparison between the Human Connectome Project and the UK Biobank dataset
Heritability of surface area and cortical thickness: a comparison between the Human Connectome Project and the UK Biobank dataset Yann Le Guen, Slim Karkar, Antoine Grigis, Cathy Philippe, Jean-François
More informationPath Analysis, SEM and the Classical Twin Model
Path Analysis, SEM and the Classical Twin Model Michael Neale & Frühling Rijsdijk 2 Virginia Institute for Psychiatric and Behavioral Genetics Virginia Commonwealth University 2 MRC SGDP Centre, Institute
More informationIntroduction of Genome wide Complex Trait Analysis (GCTA) Presenter: Yue Ming Chen Location: Stat Gen Workshop Date: 6/7/2013
Introduction of Genome wide Complex Trait Analysis (GCTA) resenter: ue Ming Chen Location: Stat Gen Workshop Date: 6/7/013 Outline Brief review of quantitative genetics Overview of GCTA Ideas Main functions
More informationNew Enhancements: GWAS Workflows with SVS
New Enhancements: GWAS Workflows with SVS August 9 th, 2017 Gabe Rudy VP Product & Engineering 20 most promising Biotech Technology Providers Top 10 Analytics Solution Providers Hype Cycle for Life sciences
More informationToday s Topics. Cracking the Genetic Code. The Process of Genetic Transmission. The Process of Genetic Transmission. Genes
Today s Topics Mechanisms of Heredity Biology of Heredity Genetic Disorders Research Methods in Behavioral Genetics Gene x Environment Interactions The Process of Genetic Transmission Genes: segments of
More informationBenchmark Dose Modeling Cancer Models. Allen Davis, MSPH Jeff Gift, Ph.D. Jay Zhao, Ph.D. National Center for Environmental Assessment, U.S.
Benchmark Dose Modeling Cancer Models Allen Davis, MSPH Jeff Gift, Ph.D. Jay Zhao, Ph.D. National Center for Environmental Assessment, U.S. EPA Disclaimer The views expressed in this presentation are those
More informationNonparametric Linkage Analysis. Nonparametric Linkage Analysis
Limitations of Parametric Linkage Analysis We previously discued parametric linkage analysis Genetic model for the disease must be specified: allele frequency parameters and penetrance parameters Lod scores
More informationHeritability. The Extended Liability-Threshold Model. Polygenic model for continuous trait. Polygenic model for continuous trait.
The Extended Liability-Threshold Model Klaus K. Holst Thomas Scheike, Jacob Hjelmborg 014-05-1 Heritability Twin studies Include both monozygotic (MZ) and dizygotic (DZ) twin
More informationInteraction of Genes and the Environment
Some Traits Are Controlled by Two or More Genes! Phenotypes can be discontinuous or continuous Interaction of Genes and the Environment Chapter 5! Discontinuous variation Phenotypes that fall into two
More informationI conducted an internship at the Texas Biomedical Research Institute in the
Ashley D. Falk ashleydfalk@gmail.com My Internship at Texas Biomedical Research Institute I conducted an internship at the Texas Biomedical Research Institute in the spring of 2013. In this report I will
More informationGraphical Modeling Approaches for Estimating Brain Networks
Graphical Modeling Approaches for Estimating Brain Networks BIOS 516 Suprateek Kundu Department of Biostatistics Emory University. September 28, 2017 Introduction My research focuses on understanding how
More informationLiability Threshold Models
Liability Threshold Models Frühling Rijsdijk & Kate Morley Twin Workshop, Boulder Tuesday March 4 th 2008 Aims Introduce model fitting to categorical data Define liability and describe assumptions of the
More informationThe Late Pretest Problem in Randomized Control Trials of Education Interventions
The Late Pretest Problem in Randomized Control Trials of Education Interventions Peter Z. Schochet ACF Methods Conference, September 2012 In Journal of Educational and Behavioral Statistics, August 2010,
More informationIntroduction to linkage and family based designs to study the genetic epidemiology of complex traits. Harold Snieder
Introduction to linkage and family based designs to study the genetic epidemiology of complex traits Harold Snieder Overview of presentation Designs: population vs. family based Mendelian vs. complex diseases/traits
More informationHeritability. The concept
Heritability The concept What is the Point of Heritability? Is a trait due to nature or nurture? (Genes or environment?) You and I think this is a good point to address, but it is not addressed! What is
More informationImpacts of genomic. pre-selection on classical genetic evaluations. Vincent Ducrocq, Clotilde Patry. Vincent Ducrocq & Clotilde Patry
Impacts of genomic pre-selection on classical genetic evaluations Vincent Ducrocq, Clotilde Patry 1 Overview Theoretical considerations suggest that genomic pre-selection of progeny-tested bulls leads
More informationPenny Rheingans. Evaluation of Visualization
Evaluation CMSC 436/636 Penny Rheingans University of Maryland, Baltimore County Evaluation of Visualization Is this technique useful? For what? To whom? Why is this technique useful? What other techniques
More informationAnalysis of bivariate binomial data: Twin analysis
Analysis of bivariate binomial data: Twin analysis Klaus Holst & Thomas Scheike March 30, 2017 Overview When looking at bivariate binomial data with the aim of learning about the dependence that is present,
More informationThe Inheritance of Complex Traits
The Inheritance of Complex Traits Differences Among Siblings Is due to both Genetic and Environmental Factors VIDEO: Designer Babies Traits Controlled by Two or More Genes Many phenotypes are influenced
More informationSUPPLEMENTAL MATERIAL
1 SUPPLEMENTAL MATERIAL Response time and signal detection time distributions SM Fig. 1. Correct response time (thick solid green curve) and error response time densities (dashed red curve), averaged across
More informationReview. Imagine the following table being obtained as a random. Decision Test Diseased Not Diseased Positive TP FP Negative FN TN
Outline 1. Review sensitivity and specificity 2. Define an ROC curve 3. Define AUC 4. Non-parametric tests for whether or not the test is informative 5. Introduce the binormal ROC model 6. Discuss non-parametric
More informationImplementation Guide for the DoseControl Dosimetry System
1.0 SCOPE The GEX DoseControl System is intended for use as a single dosimetry system solution, specifically developed to satisfy the generic dosimetry needs of radiation sterilization application users
More informationResearch Methods in Forest Sciences: Learning Diary. Yoko Lu December Research process
Research Methods in Forest Sciences: Learning Diary Yoko Lu 285122 9 December 2016 1. Research process It is important to pursue and apply knowledge and understand the world under both natural and social
More informationUnivariate modeling. Sarah Medland
Univariate modeling Sarah Medland Starting at the beginning Data preparation The algebra style used in Mx expects line per case/family (Almost) limitless number of families and variables Missing data Default
More informationCHAPTER 6. Conclusions and Perspectives
CHAPTER 6 Conclusions and Perspectives In Chapter 2 of this thesis, similarities and differences among members of (mainly MZ) twin families in their blood plasma lipidomics profiles were investigated.
More informationAn Introduction to Quantitative Genetics
An Introduction to Quantitative Genetics Mohammad Keramatipour MD, PhD Keramatipour@tums.ac.ir ac ir 1 Mendel s work Laws of inheritance Basic Concepts Applications Predicting outcome of crosses Phenotype
More informationLecture 6 Practice of Linkage Analysis
Lecture 6 Practice of Linkage Analysis Jurg Ott http://lab.rockefeller.edu/ott/ http://www.jurgott.org/pekingu/ The LINKAGE Programs http://www.jurgott.org/linkage/linkagepc.html Input: pedfile Fam ID
More informationSmall Sample Bayesian Factor Analysis. PhUSE 2014 Paper SP03 Dirk Heerwegh
Small Sample Bayesian Factor Analysis PhUSE 2014 Paper SP03 Dirk Heerwegh Overview Factor analysis Maximum likelihood Bayes Simulation Studies Design Results Conclusions Factor Analysis (FA) Explain correlation
More informationMEA DISCUSSION PAPERS
Inference Problems under a Special Form of Heteroskedasticity Helmut Farbmacher, Heinrich Kögel 03-2015 MEA DISCUSSION PAPERS mea Amalienstr. 33_D-80799 Munich_Phone+49 89 38602-355_Fax +49 89 38602-390_www.mea.mpisoc.mpg.de
More informationQuantitative Genetics. Statistics Overview: Mean. Statistics Overview: Variance. Statistics Overview: Distributions. Chapter 22
Quantitative Genetics Chapter Statistics Overview: Distributions Phenotypes on X axis, Frequencies on Y axis Statistics Overview: Mean Measure of central tendency (average) of a group of measurements X
More informationGenetics. The study of heredity. Father of Genetics: Gregor Mendel (mid 1800 s) Developed set of laws that explain how heredity works
Genetics The study of heredity Father of Genetics: Gregor Mendel (mid 1800 s) Developed set of laws that explain how heredity works Father of Genetics: Gregor Mendel original pea plant (input) offspring
More informationCombined Linkage and Association in Mx. Hermine Maes Kate Morley Dorret Boomsma Nick Martin Meike Bartels
Combined Linkage and Association in Mx Hermine Maes Kate Morley Dorret Boomsma Nick Martin Meike Bartels Boulder 2009 Outline Intro to Genetic Epidemiology Progression to Linkage via Path Models Linkage
More informationOBSERVATIONAL MEDICAL OUTCOMES PARTNERSHIP
An empirical approach to measuring and calibrating for error in observational analyses Patrick Ryan on behalf of the OMOP research team 25 April 2013 Consider a typical observational database study: Exploring
More informationWDHS Curriculum Map Probability and Statistics. What is Statistics and how does it relate to you?
WDHS Curriculum Map Probability and Statistics Time Interval/ Unit 1: Introduction to Statistics 1.1-1.3 2 weeks S-IC-1: Understand statistics as a process for making inferences about population parameters
More informationReview and Wrap-up! ESP 178 Applied Research Methods Calvin Thigpen 3/14/17 Adapted from presentation by Prof. Susan Handy
Review and Wrap-up! ESP 178 Applied Research Methods Calvin Thigpen 3/14/17 Adapted from presentation by Prof. Susan Handy Final Proposals Read instructions carefully! Check Canvas for our comments on
More informationChapter 5 INTERACTIONS OF GENES AND THE ENVIRONMENT
Chapter 5 INTERACTIONS OF GENES AND THE ENVIRONMENT Chapter Summary Up to this point, the traits you have been studying have all been controlled by one pair of genes. However, many traits, including some
More informationANALYSIS AND DETECTION OF BRAIN TUMOUR USING IMAGE PROCESSING TECHNIQUES
ANALYSIS AND DETECTION OF BRAIN TUMOUR USING IMAGE PROCESSING TECHNIQUES P.V.Rohini 1, Dr.M.Pushparani 2 1 M.Phil Scholar, Department of Computer Science, Mother Teresa women s university, (India) 2 Professor
More informationDeveloping and evaluating polygenic risk prediction models for stratified disease prevention
Developing and evaluating polygenic risk prediction models for stratified disease prevention Nilanjan Chatterjee 1 3, Jianxin Shi 3 and Montserrat García-Closas 3 Abstract Knowledge of genetics and its
More informationCitation for published version (APA): Ebbes, P. (2004). Latent instrumental variables: a new approach to solve for endogeneity s.n.
University of Groningen Latent instrumental variables Ebbes, P. IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check the document
More informationPopulation Inference post Model Selection in Neuroscience
Population Inference post Model Selection in Neuroscience Genevera I. Allen Dobelman Family Junior Chair, Department of Statistics and Electrical and Computer Engineering, Rice University, Department of
More informationAn Introduction to Quantitative Genetics I. Heather A Lawson Advanced Genetics Spring2018
An Introduction to Quantitative Genetics I Heather A Lawson Advanced Genetics Spring2018 Outline What is Quantitative Genetics? Genotypic Values and Genetic Effects Heritability Linkage Disequilibrium
More informationRussian Journal of Agricultural and Socio-Economic Sciences, 3(15)
ON THE COMPARISON OF BAYESIAN INFORMATION CRITERION AND DRAPER S INFORMATION CRITERION IN SELECTION OF AN ASYMMETRIC PRICE RELATIONSHIP: BOOTSTRAP SIMULATION RESULTS Henry de-graft Acquah, Senior Lecturer
More informationJSM Survey Research Methods Section
Methods and Issues in Trimming Extreme Weights in Sample Surveys Frank Potter and Yuhong Zheng Mathematica Policy Research, P.O. Box 393, Princeton, NJ 08543 Abstract In survey sampling practice, unequal
More information10 Intraclass Correlations under the Mixed Factorial Design
CHAPTER 1 Intraclass Correlations under the Mixed Factorial Design OBJECTIVE This chapter aims at presenting methods for analyzing intraclass correlation coefficients for reliability studies based on a
More informationMULTIFACTORIAL DISEASES. MG L-10 July 7 th 2014
MULTIFACTORIAL DISEASES MG L-10 July 7 th 2014 Genetic Diseases Unifactorial Chromosomal Multifactorial AD Numerical AR Structural X-linked Microdeletions Mitochondrial Spectrum of Alterations in DNA Sequence
More informationSome Thoughts on Calibrating LaModel
Some Thoughts on Calibrating LaModel Dr. Keith A. Heasley Professor Department of Mining Engineering West Virginia University Introduction Recent mine collapses and pillar failures have highlighted hli
More informationScore Tests of Normality in Bivariate Probit Models
Score Tests of Normality in Bivariate Probit Models Anthony Murphy Nuffield College, Oxford OX1 1NF, UK Abstract: A relatively simple and convenient score test of normality in the bivariate probit model
More informationPanel: Towards reproducible brain imaging research
WIN Conference, 24th January 2017 Panel: Towards reproducible brain imaging research Camille Maumet, Alexander Bowring, Thomas Nichols Neuroimaging Statistics, IDH, University of Warwick, UK WIN Conference,
More informationGrowing Flight Software Hands-On Experience: cfs and Educational Outreach
Growing Flight Software Hands-On Experience: cfs and Educational Outreach Allen Brown, Mathew Benson Flight Software California Institute of Technology Pasadena, CA Company Background Odyssey Space Research
More informationStatistical Genetics : Gene Mappin g through Linkag e and Associatio n
Statistical Genetics : Gene Mappin g through Linkag e and Associatio n Benjamin M Neale Manuel AR Ferreira Sarah E Medlan d Danielle Posthuma About the editors List of contributors Preface Acknowledgements
More informationCatherine A. Welch 1*, Séverine Sabia 1,2, Eric Brunner 1, Mika Kivimäki 1 and Martin J. Shipley 1
Welch et al. BMC Medical Research Methodology (2018) 18:89 https://doi.org/10.1186/s12874-018-0548-0 RESEARCH ARTICLE Open Access Does pattern mixture modelling reduce bias due to informative attrition
More informationDan Koller, Ph.D. Medical and Molecular Genetics
Design of Genetic Studies Dan Koller, Ph.D. Research Assistant Professor Medical and Molecular Genetics Genetics and Medicine Over the past decade, advances from genetics have permeated medicine Identification
More informationWhere does "analysis" enter the experimental process?
Lecture Topic : ntroduction to the Principles of Experimental Design Experiment: An exercise designed to determine the effects of one or more variables (treatments) on one or more characteristics (response
More informationDoing Quantitative Research 26E02900, 6 ECTS Lecture 6: Structural Equations Modeling. Olli-Pekka Kauppila Daria Kautto
Doing Quantitative Research 26E02900, 6 ECTS Lecture 6: Structural Equations Modeling Olli-Pekka Kauppila Daria Kautto Session VI, September 20 2017 Learning objectives 1. Get familiar with the basic idea
More informationSample Size Considerations. Todd Alonzo, PhD
Sample Size Considerations Todd Alonzo, PhD 1 Thanks to Nancy Obuchowski for the original version of this presentation. 2 Why do Sample Size Calculations? 1. To minimize the risk of making the wrong conclusion
More informationMultifactorial Inheritance. Prof. Dr. Nedime Serakinci
Multifactorial Inheritance Prof. Dr. Nedime Serakinci GENETICS I. Importance of genetics. Genetic terminology. I. Mendelian Genetics, Mendel s Laws (Law of Segregation, Law of Independent Assortment).
More informationBenchmark: Talend Open Studio vs Pentaho Data Integrator (aka Kettle) V0.23
Benchmark: Talend Open Studio vs Pentaho Data Integrator (aka Kettle) V0.23 MarcRussel@gmail.com Last modified: 2007-07-31 Table of contents Environment... 2 Test 1: Text Input file > Text Output file...
More informationMotivation: Fraud Detection
Outlier Detection Motivation: Fraud Detection http://i.imgur.com/ckkoaop.gif Jian Pei: CMPT 741/459 Data Mining -- Outlier Detection (1) 2 Techniques: Fraud Detection Features Dissimilarity Groups and
More informationThe Effects of Autocorrelated Noise and Biased HRF in fmri Analysis Error Rates
The Effects of Autocorrelated Noise and Biased HRF in fmri Analysis Error Rates Ariana Anderson University of California, Los Angeles Departments of Psychiatry and Behavioral Sciences David Geffen School
More informationProblem 3: Simulated Rheumatoid Arthritis Data
Problem 3: Simulated Rheumatoid Arthritis Data Michael B Miller Michael Li Gregg Lind Soon-Young Jang The plan
More informationNonparametric Estimates of Gene Environment Interaction Using Local Structural. Equation Modeling
April 18, 2015 Nonparametric Estimates of Gene Environment Interaction Using Local Structural Equation Modeling Daniel A. Briley 1, K. Paige Harden 1, Timothy C. Bates 2, Elliot M. Tucker-Drob 1 1 Department
More informationGraph Theory. Steffie Tomson UCLA NITP 2013
Graph Theory Steffie Tomson UCLA NITP 2013 What is graph theory? Mathematical descriptions of relationships between regions Node Opsahl et al. 2010 Why use graphs? Nodal Measures: Participation Coefficient
More informationThe Pretest! Pretest! Pretest! Assignment (Example 2)
The Pretest! Pretest! Pretest! Assignment (Example 2) May 19, 2003 1 Statement of Purpose and Description of Pretest Procedure When one designs a Math 10 exam one hopes to measure whether a student s ability
More informationLec 02: Estimation & Hypothesis Testing in Animal Ecology
Lec 02: Estimation & Hypothesis Testing in Animal Ecology Parameter Estimation from Samples Samples We typically observe systems incompletely, i.e., we sample according to a designed protocol. We then
More informationIndependent Variables Variables (factors) that are manipulated to measure their effect Typically select specific levels of each variable to test
Controlled Experiments experimental investigation of a testable hypothesis, in which conditions are set up to isolate the variables of interest ("independent variables") and test how they affect certain
More informationHeritability estimates on resting state fmri data using ENIGMA analysis pipeline
Heritability estimates on resting state fmri data using ENIGMA analysis pipeline Bhim M. Adhikari Maryland Psychiatry Research Center, Department of Psychiatry, University of Maryland School of Medicine,
More informationSupporting Information Identification of Amino Acids with Sensitive Nanoporous MoS 2 : Towards Machine Learning-Based Prediction
Supporting Information Identification of Amino Acids with Sensitive Nanoporous MoS 2 : Towards Machine Learning-Based Prediction Amir Barati Farimani, Mohammad Heiranian, Narayana R. Aluru 1 Department
More informationENIGMA for Neurorehabilitation: A Large-Scale Meta-Analysis Approach to Modeling Neuroimaging, Genetics, and Behavior
ENIGMA for Neurorehabilitation: A Large-Scale Meta-Analysis Approach to Modeling Neuroimaging, Genetics, and Behavior Sook-Lei Liew, PhD, OTR/L Director, Neural Plasticity and Neurorehabilitation Laboratory
More informationInternational Journal of Engineering Trends and Applications (IJETA) Volume 4 Issue 2, Mar-Apr 2017
RESEARCH ARTICLE OPEN ACCESS Knowledge Based Brain Tumor Segmentation using Local Maxima and Local Minima T. Kalaiselvi [1], P. Sriramakrishnan [2] Department of Computer Science and Applications The Gandhigram
More informationNeuroImage 81 (2013) Contents lists available at SciVerse ScienceDirect. NeuroImage. journal homepage:
NeuroImage 81 (2013) 455 469 Contents lists available at SciVerse ScienceDirect NeuroImage journal homepage: www.elsevier.com/locate/ynimg Multi-site genetic analysis of diffusion images and voxelwise
More informationStep 3 Tutorial #3: Obtaining equations for scoring new cases in an advanced example with quadratic term
Step 3 Tutorial #3: Obtaining equations for scoring new cases in an advanced example with quadratic term DemoData = diabetes.lgf, diabetes.dat, data5.dat We begin by opening a saved 3-class latent class
More informationEmpirical assessment of univariate and bivariate meta-analyses for comparing the accuracy of diagnostic tests
Empirical assessment of univariate and bivariate meta-analyses for comparing the accuracy of diagnostic tests Yemisi Takwoingi, Richard Riley and Jon Deeks Outline Rationale Methods Findings Summary Motivating
More informationStructured models for dengue epidemiology
Structured models for dengue epidemiology submitted by Hannah Woodall for the degree of Doctor of Philosophy of the University of Bath Department of Mathematical Sciences September 24 COPYRIGHT Attention
More informationEmpowered by Psychometrics The Fundamentals of Psychometrics. Jim Wollack University of Wisconsin Madison
Empowered by Psychometrics The Fundamentals of Psychometrics Jim Wollack University of Wisconsin Madison Psycho-what? Psychometrics is the field of study concerned with the measurement of mental and psychological
More informationPredicting Energy Balance Status of Holstein Cows using Mid-Infrared Spectral Data
Predicting Energy Balance Status of Holstein Cows using Mid-Infrared Spectral Data Sinéad Mc Parland, G.Banos, E.Wall, M.P.Coffey, H.Soyeurt, R.F.Veerkamp & D.P.Berry Introduction Energy balance (output-input)
More informationModeling Sentiment with Ridge Regression
Modeling Sentiment with Ridge Regression Luke Segars 2/20/2012 The goal of this project was to generate a linear sentiment model for classifying Amazon book reviews according to their star rank. More generally,
More informationGeneralized Estimating Equations for Depression Dose Regimes
Generalized Estimating Equations for Depression Dose Regimes Karen Walker, Walker Consulting LLC, Menifee CA Generalized Estimating Equations on the average produce consistent estimates of the regression
More informationHow to analyze correlated and longitudinal data?
How to analyze correlated and longitudinal data? Niloofar Ramezani, University of Northern Colorado, Greeley, Colorado ABSTRACT Longitudinal and correlated data are extensively used across disciplines
More informationQuality of Life. The assessment, analysis and reporting of patient-reported outcomes. Third Edition
Quality of Life The assessment, analysis and reporting of patient-reported outcomes Third Edition PETER M. FAYERS Institute of Applied Health Sciences, University ofaberdeen School of Medicine and Dentistry,
More informationPopulation-adjusted treatment comparisons Overview and recommendations from the NICE Decision Support Unit
Population-adjusted treatment comparisons Overview and recommendations from the NICE Decision Support Unit David M Phillippo, University of Bristol 3 Available from www.nicedsu.org.uk Outline Background
More informationSum of Neurally Distinct Stimulus- and Task-Related Components.
SUPPLEMENTARY MATERIAL for Cardoso et al. 22 The Neuroimaging Signal is a Linear Sum of Neurally Distinct Stimulus- and Task-Related Components. : Appendix: Homogeneous Linear ( Null ) and Modified Linear
More information