Mapping multiple QTL in experimental crosses
|
|
- Caitlin Neal
- 5 years ago
- Views:
Transcription
1 Human vs mouse Mapping multiple QTL in experimental crosses Karl W Broman Department of Biostatistics & Medical Informatics University of Wisconsin Madison Backcross P P P F BC
2 Intercross Genetic map P P F F Location (cm) 8 F X 5 Phenotype data Genotype data Sugiyama et al. Genomics :-, 5 male mice from the backcross (A B) B Blood pressure after two weeks drinking water with % NaCl X Individuals Blood pressure 5 5 Markers 8
3 Goals Identify quantitative trait loci (QTL) (and interactions among QTL) Interval estimates of QTL location Estimated QTL effects 9 Statistical structure QTL Markers Phenotype Covariates The missing data problem: Markers QTL The model selection problem: QTL, covariates phenotype ANOVA at marker loci Split mice into groups according to genotype at a marker. Do a t-test / ANOVA. Repeat for each marker. 9 Genotype at DMit blood pressure AB Genotype at DMit AB Interval mapping Lander & Botstein (989) Assume a single QTL model. Consider each position in the genome, one at a time, as the location of the putative QTL. Let q = / if the (unobserved) QTL genotype is /AB. (Or // if the QTL genotype is AA/AB/ in an intercross.) Assume y q N(µ q, σ) Calculate p q = Pr(q marker data). y marker data q p q φ(y µ q, σ)
4 LOD scores Permutation test LOD(λ) = log likelihood ratio comparing the hypothesis of a QTL at position λ versus that of no QTL { } Pr(y QTL at λ,ˆµqλ,ˆσ = log λ ) Pr(y no QTL,ˆµ,ˆσ) markers ˆµ qλ, ˆσ λ are the MLEs, assuming a single QTL at position λ. individuals genotype data phenotypes LOD scores maximum LOD score No QTL model: The phenotypes are iid N(µ, σ ). 3 5 LOD curves Permutation results 8 LOD score Genome wide maximum LOD score
5 LOD curves Modeling multiple QTL 8 Reduce residual variation increased power Separate linked QTL Identify interactions among QTL (epistasis) LOD score 5% Estimated effects -dim, -QTL scan 8 cm 3 cm cm Chr cm For each pair of positions, fit the following models: H f : y = µ + β q + β q + γq q + ɛ H a : y = µ + β q + β q + ɛ blood pressure H : y = µ + β q + ɛ H : y = µ + ɛ BA BA BA BA Genotype Genotype Genotype Genotype 9
6 -dim, -QTL scan LOD i and LOD f For each pair of positions, fit the following models: H f : y = µ + β q + β q + γq q + ɛ H a : y = µ + β q + β q + ɛ H : y = µ + β q + ɛ H : y = µ + ɛ LOD f (s, t) = log { L f (s, t) / L } dim, -QTL scan LOD i and LOD f For each pair of positions, fit the following models: H f : y = µ + β q + β q + γq q + ɛ 5 H a : y = µ + β q + β q + ɛ H : y = µ + β q + ɛ H : y = µ + ɛ LOD i (s, t) = log { L f (s, t) / L a (s, t) }
7 LOD i and LOD fv Chr : LOD i and LOD av Position (cm) Position (cm) 8 Estimated effects Hypothesis testing? x x 5 In the past, QTL mapping has been regarded as a task of hypothesis testing. Chr genotype BA Chr genotype BA Is this a QTL? Much of the focus has been on adjusting for test multiplicity. blood pressure 5 5 It is better to view the problem as one of model selection. What set of QTL are well supported? Is there evidence for QTL-QTL interactions? 95 BA 95 BA Model = a defined set of QTL and QTL-QTL interactions (and possibly covariates and QTL-covariate interactions). Chr genotype Chr 5 Genotype 9
8 Model selection What is special here? Class of models Additive models + pairwise interactions + higher-order interactions Regression trees Model fit Maximum likelihood Haley-Knott regression extended Haley-Knott Multiple imputation MCMC Model comparison Estimated prediction error AIC, BIC, penalized likelihood Bayes Model search Forward selection Backward elimination Stepwise selection Randomized algorithms Goal: identify the major players A continuum of ordinal-valued covariates (the genetic loci) Association among the covariates Loci on different chromosomes are independent Along chromosome, a very simple (and known) correlation structure 3 3 Target Exploratory methods Selection of a model includes two types of errors: Miss important terms (QTLs or interactions) Include extraneous terms Unlike in hypothesis testing, we can make both errors at the same time. Identify as many correct terms as possible, while controlling the rate of inclusion of extraneous terms. Condition on a large-effect QTL Reduce residual variation Conditional LOD score: { } Pr(data q, q ) LOD(q q ) = log Pr(data q ) Piece together the putative QTL from the d and d scans Omit loci that no longer look interesting (drop-one-at-a-time analysis) Study potential interactions among the identified loci Scan for additional loci (perhaps allowing interactions), conditional on these 3 33
9 Automation Experience Assistance to the masses Understanding performance Many phenotypes Controls rate of inclusion of extraneous terms Forward selection over-selects Forward selection followed by backward elimination works as well as MCMC Need to define performance criteria Need large-scale simulations Broman & Speed, JRSS B :-5, 3 Additive QTL Epistasis Simple situation: Dense markers Complete genotype data No epistasis y = µ + β j q j + γ jk q j q k + ɛ LOD δɛ (γ) = LOD(γ) T m γ m + T i γ i y = µ + β j q j + ɛ which β j? T m = as chosen previously LOD δ (γ) = LOD(γ) T γ T i =? vs QTL: LOD δ ( ) = LOD δ ({λ}) = LOD({λ}) T 39
10 Idea Models as graphs Imagine there are two additive QTL and consider a d, -QTL scan. T i = 95th percentile of the distribution of max LOD f (s, t) max LOD a (s, t) A C For the mouse genome: T m =.9 (BC) or 3.5 (F ) T H i =. (BC) or.8 (F ) B D Idea Results Imagine there is one QTL and consider a d, -QTL scan. T m + T i = 95th percentile of the distribution of max LOD f (s, t) max LOD (s).3. For the mouse genome: T m =.9 (BC) or 3.5 (F ) T H i =. (BC) or.8 (F ) T L i =.9 (BC) or.9 (F ) T m =.9 T i H =. T i L =.9 T m + T i H = 5.3 T m + T i L =
11 Add an interaction? Add another QTL? T m =.9 T i H =. T i L =.9 T m + T i H = 5.3 T m + T i L = T m =.9 T i H =. T i L =.9 T m + T i H = 5.3 T m + T i L = Add another QTL? Add a pair of QTL? a 3b T m =.9 T i H =. T i L =.9 T m + T i H = 5.3 T m + T i L = T m =.9 T i H =. T i L =.9 T m + T i H = 5.3 T m + T i L =
12 To do Acknowledgments Improve search procedures Study performance (especially relative to other approaches) Measuring model uncertainty Measuring uncertainty in QTL location Ani Manichaikul Gary Churchill Śaunak Sen Terry Speed Brian Yandell Johns Hopkins University Jackson Laboratory University of California, San Francisco University of California, Berkeley University of Wisconsin, Madison Fumihiro Sugiyama Bev Paigen now at University of Tsukuba, Japan Jackson Laboratory Summary QTL mapping is a model selection problem The criterion for comparing models is most important We re focusing on a penalized likelihood method and are close to a practiceable solution
R/qtl workshop. (part 2) Karl Broman. Biostatistics and Medical Informatics University of Wisconsin Madison. kbroman.org
R/qtl workshop (part 2) Karl Broman Biostatistics and Medical Informatics University of Wisconsin Madison kbroman.org github.com/kbroman @kwbroman Example Sugiyama et al. Genomics 71:70-77, 2001 250 male
More informationMapping multiple QTL in experimental crosses
Mapping multiple QTL in experimental crosses Karl W Broman Department of Biostatistics and Medical Informatics University of Wisconsin Madison www.biostat.wisc.edu/~kbroman [ Teaching Miscellaneous lectures]
More informationMultiple QTL mapping
Multiple QTL mapping Karl W Broman Department of Biostatistics Johns Hopkins University www.biostat.jhsph.edu/~kbroman [ Teaching Miscellaneous lectures] 1 Why? Reduce residual variation = increased power
More informationIntroduction to QTL mapping in model organisms
Introduction to QTL mapping in model organisms Karl W Broman Department of Biostatistics and Medical Informatics University of Wisconsin Madison www.biostat.wisc.edu/~kbroman [ Teaching Miscellaneous lectures]
More informationIntroduction to QTL mapping in model organisms
Introduction to QTL mapping in model organisms Karl Broman Biostatistics and Medical Informatics University of Wisconsin Madison kbroman.org github.com/kbroman @kwbroman Backcross P 1 P 2 P 1 F 1 BC 4
More informationIntroduction to QTL mapping in model organisms
Human vs mouse Introduction to QTL mapping in model organisms Karl W Broman Department of Biostatistics Johns Hopkins University www.biostat.jhsph.edu/~kbroman [ Teaching Miscellaneous lectures] www.daviddeen.com
More informationQTL Model Search. Brian S. Yandell, UW-Madison January 2017
QTL Model Search Brian S. Yandell, UW-Madison January 2017 evolution of QTL models original ideas focused on rare & costly markers models & methods refined as technology advanced single marker regression
More informationStatistical issues in QTL mapping in mice
Statistical issues in QTL mapping in mice Karl W Broman Department of Biostatistics Johns Hopkins University http://www.biostat.jhsph.edu/~kbroman Outline Overview of QTL mapping The X chromosome Mapping
More informationIntroduction to QTL mapping in model organisms
Introduction to QTL mapping in model organisms Karl W Broman Department of Biostatistics Johns Hopkins University kbroman@jhsph.edu www.biostat.jhsph.edu/ kbroman Outline Experiments and data Models ANOVA
More informationGene mapping in model organisms
Gene mapping in model organisms Karl W Broman Department of Biostatistics Johns Hopkins University http://www.biostat.jhsph.edu/~kbroman Goal Identify genes that contribute to common human diseases. 2
More informationIntroduction to QTL mapping in model organisms
Introduction to QTL mapping in model organisms Karl W Broman Department of Biostatistics Johns Hopkins University kbroman@jhsph.edu www.biostat.jhsph.edu/ kbroman Outline Experiments and data Models ANOVA
More informationMapping QTL to a phylogenetic tree
Mapping QTL to a phylogenetic tree Karl W Broman Department of Biostatistics & Medical Informatics University of Wisconsin Madison www.biostat.wisc.edu/~kbroman Human vs mouse www.daviddeen.com 3 Intercross
More informationQTL Mapping I: Overview and using Inbred Lines
QTL Mapping I: Overview and using Inbred Lines Key idea: Looking for marker-trait associations in collections of relatives If (say) the mean trait value for marker genotype MM is statisically different
More informationOverview. Background
Overview Implementation of robust methods for locating quantitative trait loci in R Introduction to QTL mapping Andreas Baierl and Andreas Futschik Institute of Statistics and Decision Support Systems
More informationQTL model selection: key players
Bayesian Interval Mapping. Bayesian strategy -9. Markov chain sampling 0-7. sampling genetic architectures 8-5 4. criteria for model selection 6-44 QTL : Bayes Seattle SISG: Yandell 008 QTL model selection:
More informationUse of hidden Markov models for QTL mapping
Use of hidden Markov models for QTL mapping Karl W Broman Department of Biostatistics, Johns Hopkins University December 5, 2006 An important aspect of the QTL mapping problem is the treatment of missing
More informationCausal Model Selection Hypothesis Tests in Systems Genetics
1 Causal Model Selection Hypothesis Tests in Systems Genetics Elias Chaibub Neto and Brian S Yandell SISG 2012 July 13, 2012 2 Correlation and Causation The old view of cause and effect... could only fail;
More informationQTL model selection: key players
QTL Model Selection. Bayesian strategy. Markov chain sampling 3. sampling genetic architectures 4. criteria for model selection Model Selection Seattle SISG: Yandell 0 QTL model selection: key players
More informationQuantile-based permutation thresholds for QTL hotspot analysis: a tutorial
Quantile-based permutation thresholds for QTL hotspot analysis: a tutorial Elias Chaibub Neto and Brian S Yandell September 18, 2013 1 Motivation QTL hotspots, groups of traits co-mapping to the same genomic
More informationInferring Causal Phenotype Networks from Segregating Populat
Inferring Causal Phenotype Networks from Segregating Populations Elias Chaibub Neto chaibub@stat.wisc.edu Statistics Department, University of Wisconsin - Madison July 15, 2008 Overview Introduction Description
More informationCausal Graphical Models in Systems Genetics
1 Causal Graphical Models in Systems Genetics 2013 Network Analysis Short Course - UCLA Human Genetics Elias Chaibub Neto and Brian S Yandell July 17, 2013 Motivation and basic concepts 2 3 Motivation
More informationSelecting explanatory variables with the modified version of Bayesian Information Criterion
Selecting explanatory variables with the modified version of Bayesian Information Criterion Institute of Mathematics and Computer Science, Wrocław University of Technology, Poland in cooperation with J.K.Ghosh,
More informationLocating multiple interacting quantitative trait. loci using rank-based model selection
Genetics: Published Articles Ahead of Print, published on May 16, 2007 as 10.1534/genetics.106.068031 Locating multiple interacting quantitative trait loci using rank-based model selection, 1 Ma lgorzata
More informationLecture 8. QTL Mapping 1: Overview and Using Inbred Lines
Lecture 8 QTL Mapping 1: Overview and Using Inbred Lines Bruce Walsh. jbwalsh@u.arizona.edu. University of Arizona. Notes from a short course taught Jan-Feb 2012 at University of Uppsala While the machinery
More informationBinary trait mapping in experimental crosses with selective genotyping
Genetics: Published Articles Ahead of Print, published on May 4, 2009 as 10.1534/genetics.108.098913 Binary trait mapping in experimental crosses with selective genotyping Ani Manichaikul,1 and Karl W.
More informationTHE problem of identifying the genetic factors un- QTL models because of their ability to separate linked
Copyright 2001 by the Genetics Society of America A Statistical Framework for Quantitative Trait Mapping Śaunak Sen and Gary A. Churchill The Jackson Laboratory, Bar Harbor, Maine 04609 Manuscript received
More informationFast Bayesian Methods for Genetic Mapping Applicable for High-Throughput Datasets
Fast Bayesian Methods for Genetic Mapping Applicable for High-Throughput Datasets Yu-Ling Chang A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial fulfillment
More informationAn applied statistician does probability:
An applied statistician does probability: It s not pretty Karl W Broman Department of Biostatistics & Medical Informatics University of Wisconsin Madison www.biostat.wisc.edu/~kbroman An applied statistician
More informationInferring Genetic Architecture of Complex Biological Processes
Inferring Genetic Architecture of Complex Biological Processes BioPharmaceutical Technology Center Institute (BTCI) Brian S. Yandell University of Wisconsin-Madison http://www.stat.wisc.edu/~yandell/statgen
More informationAnumber of statistical methods are available for map- 1995) in some standard designs. For backcross populations,
Copyright 2004 by the Genetics Society of America DOI: 10.1534/genetics.104.031427 An Efficient Resampling Method for Assessing Genome-Wide Statistical Significance in Mapping Quantitative Trait Loci Fei
More informationBAYESIAN MAPPING OF MULTIPLE QUANTITATIVE TRAIT LOCI
BAYESIAN MAPPING OF MULTIPLE QUANTITATIVE TRAIT LOCI By DÁMARIS SANTANA MORANT A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR
More informationHow the mean changes depends on the other variable. Plots can show what s happening...
Chapter 8 (continued) Section 8.2: Interaction models An interaction model includes one or several cross-product terms. Example: two predictors Y i = β 0 + β 1 x i1 + β 2 x i2 + β 12 x i1 x i2 + ɛ i. How
More informationModel Selection for Multiple QTL
Model Selection for Multiple TL 1. reality of multiple TL 3-8. selecting a class of TL models 9-15 3. comparing TL models 16-4 TL model selection criteria issues of detecting epistasis 4. simulations and
More informationTECHNICAL REPORT NO April 22, Causal Model Selection Tests in Systems Genetics 1
DEPARTMENT OF STATISTICS University of Wisconsin 1300 University Avenue Madison, WI 53706 TECHNICAL REPORT NO. 1157 April 22, 2010 Causal Model Selection Tests in Systems Genetics 1 Elias Chaibub Neto
More informationCausal Model Selection Hypothesis Tests. in Systems Genetics
Causal Model Selection Hypothesis Tests in Systems Genetics Elias Chaibub Neto 1 Aimee T. Broman 2 Mark P Keller 2 Alan D Attie 2 Bin Zhang 1 Jun Zhu 1 Brian S Yandell 3,4 1 Sage Bionetworks, Seattle,
More informationBIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation
BIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation Yujin Chung November 29th, 2016 Fall 2016 Yujin Chung Lec13: MLE Fall 2016 1/24 Previous Parametric tests Mean comparisons (normality assumption)
More informationOne-week Course on Genetic Analysis and Plant Breeding January 2013, CIMMYT, Mexico LOD Threshold and QTL Detection Power Simulation
One-week Course on Genetic Analysis and Plant Breeding 21-2 January 213, CIMMYT, Mexico LOD Threshold and QTL Detection Power Simulation Jiankang Wang, CIMMYT China and CAAS E-mail: jkwang@cgiar.org; wangjiankang@caas.cn
More informationPrediction of the Confidence Interval of Quantitative Trait Loci Location
Behavior Genetics, Vol. 34, No. 4, July 2004 ( 2004) Prediction of the Confidence Interval of Quantitative Trait Loci Location Peter M. Visscher 1,3 and Mike E. Goddard 2 Received 4 Sept. 2003 Final 28
More informationMultiple interval mapping for ordinal traits
Genetics: Published Articles Ahead of Print, published on April 3, 2006 as 10.1534/genetics.105.054619 Multiple interval mapping for ordinal traits Jian Li,,1, Shengchu Wang and Zhao-Bang Zeng,, Bioinformatics
More informationTECHNICAL REPORT NO December 1, 2008
DEPARTMENT OF STATISTICS University of Wisconsin 300 University Avenue Madison, WI 53706 TECHNICAL REPORT NO. 46 December, 2008 Revised on January 27, 2009 Causal Graphical Models in System Genetics: a
More informationA Statistical Framework for Expression Trait Loci (ETL) Mapping. Meng Chen
A Statistical Framework for Expression Trait Loci (ETL) Mapping Meng Chen Prelim Paper in partial fulfillment of the requirements for the Ph.D. program in the Department of Statistics University of Wisconsin-Madison
More informationVariance Component Models for Quantitative Traits. Biostatistics 666
Variance Component Models for Quantitative Traits Biostatistics 666 Today Analysis of quantitative traits Modeling covariance for pairs of individuals estimating heritability Extending the model beyond
More informationQuantile based Permutation Thresholds for QTL Hotspots. Brian S Yandell and Elias Chaibub Neto 17 March 2012
Quantile based Permutation Thresholds for QTL Hotspots Brian S Yandell and Elias Chaibub Neto 17 March 2012 2012 Yandell 1 Fisher on inference We may at once admit that any inference from the particular
More informationBayesian Model Selection for Multiple QTL. Brian S. Yandell
Bayesian Model Selection for Multiple QTL Brian S. Yandell University of Wisconsin-Madison www.stat.wisc.edu/~yandell/statgen Jackson Laboratory, October 005 October 005 Jax Workshop Brian S. Yandell 1
More information. Also, in this case, p i = N1 ) T, (2) where. I γ C N(N 2 2 F + N1 2 Q)
Supplementary information S7 Testing for association at imputed SPs puted SPs Score tests A Score Test needs calculations of the observed data score and information matrix only under the null hypothesis,
More information1. what is the goal of QTL study?
NSF UAB Course 008 Bayesian Interval Mapping Brian S. Yandell, UW-Madison www.stat.wisc.edu/~yandell/statgen overview: multiple QTL approaches Bayesian QTL mapping & model selection data examples in detail
More informationTutorial Session 2. MCMC for the analysis of genetic data on pedigrees:
MCMC for the analysis of genetic data on pedigrees: Tutorial Session 2 Elizabeth Thompson University of Washington Genetic mapping and linkage lod scores Monte Carlo likelihood and likelihood ratio estimation
More informationMethods for QTL analysis
Methods for QTL analysis Julius van der Werf METHODS FOR QTL ANALYSIS... 44 SINGLE VERSUS MULTIPLE MARKERS... 45 DETERMINING ASSOCIATIONS BETWEEN GENETIC MARKERS AND QTL WITH TWO MARKERS... 45 INTERVAL
More informationarxiv: v4 [stat.me] 2 May 2015
Submitted to the Annals of Statistics arxiv: arxiv:0000.0000 BAYESIAN NONPARAMETRIC TESTS VIA SLICED INVERSE MODELING By Bo Jiang, Chao Ye and Jun S. Liu Harvard University, Tsinghua University and Harvard
More informationHierarchical Generalized Linear Models for Multiple QTL Mapping
Genetics: Published Articles Ahead of Print, published on January 1, 009 as 10.1534/genetics.108.099556 Hierarchical Generalized Linear Models for Multiple QTL Mapping Nengun Yi 1,* and Samprit Baneree
More informationLecture 9. QTL Mapping 2: Outbred Populations
Lecture 9 QTL Mapping 2: Outbred Populations Bruce Walsh. Aug 2004. Royal Veterinary and Agricultural University, Denmark The major difference between QTL analysis using inbred-line crosses vs. outbred
More informationLecture 11: Multiple trait models for QTL analysis
Lecture 11: Multiple trait models for QTL analysis Julius van der Werf Multiple trait mapping of QTL...99 Increased power of QTL detection...99 Testing for linked QTL vs pleiotropic QTL...100 Multiple
More informationA new simple method for improving QTL mapping under selective genotyping
Genetics: Early Online, published on September 22, 2014 as 10.1534/genetics.114.168385 A new simple method for improving QTL mapping under selective genotyping Hsin-I Lee a, Hsiang-An Ho a and Chen-Hung
More informationUsers Manual of QTL IciMapping v3.1
8 6 4 2 0 0 0 0 0 0 0 0 0 0 0 0 0 Users Manual of QTL IciMapping v3.1 Jiankang Wang, Huihui Li, Luyan Zhang, Chunhui Li and Lei Meng QZ10 QZ9 QZ8 QZ7 QZ6 QZ5 QZ4 QZ3 QZ2 QZ1 L 120 90 60 30 0 120 90 60
More informationQuantitative Genetics. February 16, 2010
Quantitative Genetics February 16, 2010 A General Model y ij z ij y ij y ij = y i(j) = g j + ε ij 2 1 [zij =j]f ij ( complete likelihood ) j=0 2 ω ij f ij j=0 ( likelihood ) ω i0 = Pr[z ij = j] = q 2 =
More informationLecture 6. QTL Mapping
Lecture 6 QTL Mapping Bruce Walsh. Aug 2003. Nordic Summer Course MAPPING USING INBRED LINE CROSSES We start by considering crosses between inbred lines. The analysis of such crosses illustrates many of
More informationQuantitative genetics
Quantitative genetics Many traits that are important in agriculture, biology and biomedicine are continuous in their phenotypes. For example, Crop Yield Stemwood Volume Plant Disease Resistances Body Weight
More informationBTRY 4830/6830: Quantitative Genomics and Genetics
BTRY 4830/6830: Quantitative Genomics and Genetics Lecture 23: Alternative tests in GWAS / (Brief) Introduction to Bayesian Inference Jason Mezey jgm45@cornell.edu Nov. 13, 2014 (Th) 8:40-9:55 Announcements
More informationData Mining Stat 588
Data Mining Stat 588 Lecture 02: Linear Methods for Regression Department of Statistics & Biostatistics Rutgers University September 13 2011 Regression Problem Quantitative generic output variable Y. Generic
More informationMOST complex traits are influenced by interacting
Copyright Ó 2009 by the Genetics Society of America DOI: 10.1534/genetics.108.099556 Hierarchical Generalized Linear Models for Multiple Quantitative Trait Locus Mapping Nengjun Yi*,1 and Samprit Banerjee
More informationTail probability of linear combinations of chi-square variables and its application to influence analysis in QTL detection
Tail probability of linear combinations of chi-square variables and its application to influence analysis in QTL detection Satoshi Kuriki and Xiaoling Dou (Inst. Statist. Math., Tokyo) ISM Cooperative
More informationStat/F&W Ecol/Hort 572 Review Points Ané, Spring 2010
1 Linear models Y = Xβ + ɛ with ɛ N (0, σ 2 e) or Y N (Xβ, σ 2 e) where the model matrix X contains the information on predictors and β includes all coefficients (intercept, slope(s) etc.). 1. Number of
More informationMULTIPLE-TRAIT MULTIPLE-INTERVAL MAPPING OF QUANTITATIVE-TRAIT LOCI ROBY JOEHANES
MULTIPLE-TRAIT MULTIPLE-INTERVAL MAPPING OF QUANTITATIVE-TRAIT LOCI by ROBY JOEHANES B.S., Universitas Pelita Harapan, Indonesia, 1999 M.S., Kansas State University, 2002 A REPORT submitted in partial
More informationStandard Errors & Confidence Intervals. N(0, I( β) 1 ), I( β) = [ 2 l(β, φ; y) β i β β= β j
Standard Errors & Confidence Intervals β β asy N(0, I( β) 1 ), where I( β) = [ 2 l(β, φ; y) ] β i β β= β j We can obtain asymptotic 100(1 α)% confidence intervals for β j using: β j ± Z 1 α/2 se( β j )
More informationCausal Network Models for Correlated Quantitative Traits. outline
Causal Network Models for Correlated Quantitative Traits Brian S. Yandell UW Madison October 2012 www.stat.wisc.edu/~yandell/statgen Jax SysGen: Yandell 2012 1 outline Correlation and causation Correlatedtraitsinorganized
More informationLecture 8 Genomic Selection
Lecture 8 Genomic Selection Guilherme J. M. Rosa University of Wisconsin-Madison Mixed Models in Quantitative Genetics SISG, Seattle 18 0 Setember 018 OUTLINE Marker Assisted Selection Genomic Selection
More informationExpression QTLs and Mapping of Complex Trait Loci. Paul Schliekelman Statistics Department University of Georgia
Expression QTLs and Mapping of Complex Trait Loci Paul Schliekelman Statistics Department University of Georgia Definitions: Genes, Loci and Alleles A gene codes for a protein. Proteins due everything.
More informationOn the mapping of quantitative trait loci at marker and non-marker locations
Genet. Res., Camb. (2002), 79, pp. 97 106. With 3 figures. 2002 Cambridge University Press DOI: 10.1017 S0016672301005420 Printed in the United Kingdom 97 On the mapping of quantitative trait loci at marker
More informationMLR Model Selection. Author: Nicholas G Reich, Jeff Goldsmith. This material is part of the statsteachr project
MLR Model Selection Author: Nicholas G Reich, Jeff Goldsmith This material is part of the statsteachr project Made available under the Creative Commons Attribution-ShareAlike 3.0 Unported License: http://creativecommons.org/licenses/by-sa/3.0/deed.en
More informationGenome-wide Multiple Loci Mapping in Experimental Crosses by the Iterative Adaptive Penalized Regression
Genetics: Published Articles Ahead of Print, published on February 15, 2010 as 10.1534/genetics.110.114280 Genome-wide Multiple Loci Mapping in Experimental Crosses by the Iterative Adaptive Penalized
More informationMODEL-FREE LINKAGE AND ASSOCIATION MAPPING OF COMPLEX TRAITS USING QUANTITATIVE ENDOPHENOTYPES
MODEL-FREE LINKAGE AND ASSOCIATION MAPPING OF COMPLEX TRAITS USING QUANTITATIVE ENDOPHENOTYPES Saurabh Ghosh Human Genetics Unit Indian Statistical Institute, Kolkata Most common diseases are caused by
More informationMixture Cure Model with an Application to Interval Mapping of Quantitative Trait Loci
Mixture Cure Model with an Application to Interval Mapping of Quantitative Trait Loci Abstract. When censored time-to-event data are used to map quantitative trait loci (QTL), the existence of nonsusceptible
More informationModule 4: Bayesian Methods Lecture 9 A: Default prior selection. Outline
Module 4: Bayesian Methods Lecture 9 A: Default prior selection Peter Ho Departments of Statistics and Biostatistics University of Washington Outline Je reys prior Unit information priors Empirical Bayes
More informationSingle Marker Analysis and Interval Mapping
Single Marker Analysis and Interval Mapping Jiankang Wang, CIMMYT China and CAAS E-mail: jkwang@cgiar.org; wangjiankang@caas.cn Web: http://www.isbreeding.net 1 Comparison of Estimated Recombination Frequency
More informationGenotype Imputation. Biostatistics 666
Genotype Imputation Biostatistics 666 Previously Hidden Markov Models for Relative Pairs Linkage analysis using affected sibling pairs Estimation of pairwise relationships Identity-by-Descent Relatives
More informationHotspots and Causal Inference For Yeast Data
Hotspots and Causal Inference For Yeast Data Elias Chaibub Neto and Brian S Yandell October 24, 2012 Here we reproduce the analysis of the budding yeast genetical genomics data-set presented in Chaibub
More informationwhat is Bayes theorem? posterior = likelihood * prior / C
who was Bayes? Reverend Thomas Bayes (70-76) part-time mathematician buried in Bunhill Cemetary, Moongate, London famous paper in 763 Phil Trans Roy Soc London was Bayes the first with this idea? (Laplace?)
More informationMixture models for heterogeneity in ranked data
Mixture models for heterogeneity in ranked data Brian Francis Lancaster University, UK Regina Dittrich, Reinhold Hatzinger Vienna University of Economics CSDA 2005 Limassol 1 Introduction Social surveys
More informationRegression, Ridge Regression, Lasso
Regression, Ridge Regression, Lasso Fabio G. Cozman - fgcozman@usp.br October 2, 2018 A general definition Regression studies the relationship between a response variable Y and covariates X 1,..., X n.
More informationNIH Public Access Author Manuscript Ann Appl Stat. Author manuscript; available in PMC 2011 January 7.
NIH Public Access Author Manuscript Published in final edited form as: Ann Appl Stat. 2010 March 1; 4(1): 320 339. CAUSAL GRAPHICAL MODELS IN SYSTEMS GENETICS: A UNIFIED FRAMEWORK FOR JOINT INFERENCE OF
More informationAssociation studies and regression
Association studies and regression CM226: Machine Learning for Bioinformatics. Fall 2016 Sriram Sankararaman Acknowledgments: Fei Sha, Ameet Talwalkar Association studies and regression 1 / 104 Administration
More informationMANY important problems in biology and medicine
Copyright Ó 2011 by the Genetics Society of America DOI: 10.1534/genetics.110.122796 A Model Selection Approach for Expression Quantitative Trait Loci (eqtl) Mapping Ping Wang,* John A. Dawson,* Mark P.
More informationQUANTITATIVE trait analysis has many applica- several crosses, since the correlations may not be the
Copyright 001 by the Genetics Society of America Statistical Issues in the Analysis of Quantitative Traits in Combined Crosses Fei Zou, Brian S. Yandell and Jason P. Fine Department of Statistics, University
More informationConsistent high-dimensional Bayesian variable selection via penalized credible regions
Consistent high-dimensional Bayesian variable selection via penalized credible regions Howard Bondell bondell@stat.ncsu.edu Joint work with Brian Reich Howard Bondell p. 1 Outline High-Dimensional Variable
More informationDetection of multiple QTL with epistatic effects under a mixed inheritance model in an outbred population
Genet. Sel. Evol. 36 (2004) 415 433 415 c INRA, EDP Sciences, 2004 DOI: 10.1051/gse:2004009 Original article Detection of multiple QTL with epistatic effects under a mixed inheritance model in an outbred
More informationDimensionality Reduction Techniques (DRT)
Dimensionality Reduction Techniques (DRT) Introduction: Sometimes we have lot of variables in the data for analysis which create multidimensional matrix. To simplify calculation and to get appropriate,
More informationLinear Regression (1/1/17)
STA613/CBB540: Statistical methods in computational biology Linear Regression (1/1/17) Lecturer: Barbara Engelhardt Scribe: Ethan Hada 1. Linear regression 1.1. Linear regression basics. Linear regression
More informationQuantitative Genomics and Genetics BTRY 4830/6830; PBSB
Quantitative Genomics and Genetics BTRY 4830/6830; PBSB.5201.01 Lecture16: Population structure and logistic regression I Jason Mezey jgm45@cornell.edu April 11, 2017 (T) 8:40-9:55 Announcements I April
More informationLinkage analysis and QTL mapping in autotetraploid species. Christine Hackett Biomathematics and Statistics Scotland Dundee DD2 5DA
Linkage analysis and QTL mapping in autotetraploid species Christine Hackett Biomathematics and Statistics Scotland Dundee DD2 5DA Collaborators John Bradshaw Zewei Luo Iain Milne Jim McNicol Data and
More informationTHE basic principle of using genetic markers to study genetic marker map throughout the genome, IM can
Copyright 1999 by the Genetics Society of America Multiple Interval Mapping for Quantitative Trait Loci Chen-Hung Kao,* Zhao-Bang Zeng and Robert D. Teasdale *Institute of Statistical Science, Academia
More informationCalculation of IBD probabilities
Calculation of IBD probabilities David Evans and Stacey Cherny University of Oxford Wellcome Trust Centre for Human Genetics This Session IBD vs IBS Why is IBD important? Calculating IBD probabilities
More informationCalculation of IBD probabilities
Calculation of IBD probabilities David Evans University of Bristol This Session Identity by Descent (IBD) vs Identity by state (IBS) Why is IBD important? Calculating IBD probabilities Lander-Green Algorithm
More informationGWAS. Genotype-Phenotype Association CMSC858P Spring 2012 Hector Corrada Bravo University of Maryland. logistic regression. logistic regression
Genotype-Phenotype Association CMSC858P Spring 202 Hector Corrada Bravo University of Maryland GWAS Genome-wide association studies Scans for SNPs (or other structural variants) that show association with
More informationThe Admixture Model in Linkage Analysis
The Admixture Model in Linkage Analysis Jie Peng D. Siegmund Department of Statistics, Stanford University, Stanford, CA 94305 SUMMARY We study an appropriate version of the score statistic to test the
More informationThe genomes of recombinant inbred lines
The genomes of recombinant inbred lines Karl W Broman Department of Biostatistics Johns Hopkins University http://www.biostat.jhsph.edu/~kbroman C57BL/6 2 1 Recombinant inbred lines (by sibling mating)
More informationAffected Sibling Pairs. Biostatistics 666
Affected Sibling airs Biostatistics 666 Today Discussion of linkage analysis using affected sibling pairs Our exploration will include several components we have seen before: A simple disease model IBD
More informationQUANTITATIVE trait analysis plays an important observed. These assumptions are likely to be false when
Copyright 2004 y the Genetics Society of America DOI: 10.1534/genetics.103.023903 Mapping Quantitative Trait Loci With Censored Oservations Guoqing Diao, D. Y. Lin 1 and Fei Zou Department of Biostatistics,
More informationSTUDIES of genotype phenotype association are central to. Multivariate Analysis of Genotype Phenotype Association HIGHLIGHTED ARTICLE INVESTIGATION
HIGHLIGHTED ARTICLE INVESTIGATION Multivariate Analysis of Genotype Phenotype Association Philipp Mitteroecker,*,1 James M. Cheverud, and Mihaela Pavlicev *Department of Theoretical Biology, University
More informationEmpirical Bayesian LASSO-logistic regression for multiple binary trait locus mapping
Huang et al. BMC Genetics 2013, 14:5 METHODOLOGY ARTICLE Open Access Empirical Bayesian LASSO-logistic regression for multiple binary trait locus mapping Anhui Huang 1, Shizhong Xu 2 and Xiaodong Cai 1*
More informationMANY complex human diseases and traits of bio- tive corrections for multiple testing. Non-Bayesian model
Copyright 2005 by the Genetics Society of America DOI: 10.1534/genetics.104.040386 Bayesian Model Selection for Genome-Wide Epistatic Quantitative Trait Loci Analysis Nengjun Yi,*,,1 Brian S. Yandell,
More information