GenCB 511 Coarse Notes Population Genetics NONRANDOM MATING & GENETIC DRIFT

Size: px
Start display at page:

Download "GenCB 511 Coarse Notes Population Genetics NONRANDOM MATING & GENETIC DRIFT"

Transcription

1 NONRANDOM MATING & GENETIC DRIFT NONRANDOM MATING/INBREEDING READING: Hartl & Clark, Wll dstngush two tyes of nonrandom matng: (1) Assortatve matng: matng between ndvduals wth smlar henotyes or among ndvduals that occur n a artcular locaton. () Inbreedng: matng between related ndvduals. Both tyes of nonrandom matng may have smlar consequences snce ndvduals wth smlar henotyes often have smlar genotyes. It s often dffcult to searate cause from effect. E.g., ndvduals wth smlar henotyes may mate because a) henotyc assortatve matng occurs; b) matng wth relatves s referred; c) matngs are rmarly based on roxmty. Poulaton subdvson: The Wahlund Effect It turns out that oulaton subdvson er se can effect the dstrbuton of genotyes n the entre oulaton. Consder a locus wth alleles (A and a) and a collecton of solated suboulatons numbered 1,, 3,... Let the frequences of A and a n suboulaton be and q. Assumng random matng wthn each (solated) suboulaton: Freq(AA) n suboulaton Freq(Aa) n suboulaton q Freq(aa) n suboulaton q Let Avg ( ) average freq. of A across all suboulatons. ( ) 1 - Lkewse, let q Avg q What s the average frequency of each genotye over all the suboulatons? Consder AA homozygotes frst: From Fun Facts: Var( X) E X So, Avg ( ) ( ) [ E( X) ] so E( X ) [ E( X) ] + Var X ( ) + Var ( ) II-1

2 Smlarly, for aa homozygotes: Avg q ( ) q q + Var q Fnally, for Aa heterozygotes: Avg q A Thought Exerment: ( ) q + Var( ) snce Var(q) Var(1 ) Var(). ( ) 1 q 1 q Var ( ) q Var( ). Suose genotyes were randomly samled from a oulaton whose substructure was unknown. The frequences of A and a n samle would be and q. Wth random matng, would then exect to fnd genotyes n roortons AA : Aa :aa : q : q. But, the genotye frequences observed would be AA : Aa :aa : q : q + Var( ): q Var( ) : q + Var( ). I.e., would fnd an excess of homozygotes and a defct of heterozygotes, comared to exectatons. Why? Smly because of oulaton subdvson and, n artcular, varance n allele frequences across suboulatons. Gven across-suboulatons dfferences n allele frequences, the aarent excess n homozygotes and defct of heterozygotes from what s exected were the entre oulaton to mate at random defnes what s called the Wahlund Effect. The Wahlund effect s a common cause of non-conformty to Hardy-Wenberg exectatons n oulaton samles. INBREEDING Wll now consder the genetc consequences of matng between relatves: nbreedng. In 190's, Sewall Wrght nvented an ngenous aroach to trackng genotyes through edgrees based on the robablty that allele coes are dentcal by descent. Wll follow the French genetcst Malecot's reworkng of Wrght's method here. e.g., B C D E G II-

3 Identty by descent (IBD): Two alleles are dentcal by descent f (1) both are descended from the same allele n a common ancestor or () one allele s descended from the other. Wll mean IBD relatve to a secfc base oulaton (whose alleles are deemed to be not IBD). Defnton: The nbreedng coeffcent, f, of an ndvdual s the robablty that ts two gene coes at a locus are dentcal by descent. Once f s known, t's not hard to fnd the robabltes that s AA, Aa, or aa: Consder a randomly chosen ndvdual: Wth robablty f, both gene coes n that ndvdual are IBD. Then both wll be A f the allele they were coed from n the base oulaton were A. But, A occurs wth frequency n the base oulaton, so the robablty of beng an AA gven both genes are IBD s. the robablty of gettng two A alleles that are IBD s f. Lkewse, the robablty of beng aa wth both genes are IBD s f ( 1 ). Wth robablty 1 f, the two genes n an ndvdual wll not be IBD. Must have descended from dfferent allele coes n the base oulaton. Assumng the coes are made ndeendently, then wth robablty, the coed alleles are both A (genotye AA), etc. Puttng ths all together, have: ( 1 f ) + f P Aa ( 1 f ) ( 1 ) ( )( 1 ) + f ( 1 ) P aa 1 f Great! So now only need to determne f. Thanks to Wrght, ths s very easy to do. C E.g., Let's fnd f n the edgree above. c c' Need only concentrate on the central art of the edgree: f Prob e c ( ) Prob( c c ) Prob ( c g) E G e g II-3

4 So the exected genotye frequences for ths edgree are: ( 7 8) P Aa ( 7 8)( 1 ) 7 8 ( )( 1 ) + 1 8( 1 ) ASIDE: What s the average frequency of A among ndvduals wth nbreedng coeffcent f? Freq A [ ] + 1 [( 1 f ) ( 1 ) ] [ ] + f ( 1 f ) + f ( ) + 1 P 1 f Aa ( ) + f ( 1 f ) + Inbreedng does not affect allele frequences on average, but does affect the robabltes that A's or a's co-occur n an ndvdual. Comutng Inbreedng Coeffcents (n general) Say we want to fnd f I n the edgree at rght: Rules: B (1) Enumerate each loo () Each loo must C D a) go through each ndvdual no more than once b) only change from u to down once (3) Multly by 1 for each assage through an ndvdual If the assage through an ndvdual nvolves a change of drecton (u/down) multly by 1 + f ndvdual. ( ) nstead of 1/, where f s the nbreedng coeffcent for that (4) Add the robabltes of each loo. For the above examle: f I G E C G E D B ( ) 1 D 1 E G ( ) 1 E ( ) 1 ( ) G 1+ fe I II-4

5 Some forbdden loos: IGECBDEGI (goes through G twce) IDBCEGI(already counted) IGECBDEI (loo ECBDE already accounted for n f E ) Evolutonary Alcaton: Kn Selecton Probablty that two ndvduals share an allele descended from a common ancestor s called the knsh coeffcent or coeffcent of consangunty. Knsh coeffcent between ndvduals A and B s denoted F AB. What s F G n the last examle? Clearly, t s the nbreedng coeffcent of ther offsrng I, f I 7 3. The connecton between f and F: The knsh coeffcent of two ndvduals s equal to the nbreedng coeffcent of ther (erhas hyothetcal) offsrng e.g., F mother,daughter 1/4 (assumng unrelated, non-nbred arents). Knsh coeffcents useful n studyng evoluton of socal (esecally "altrustc") trats. In 1964, W. D. Hamlton roosed a rule-of-thumb to determne whether a rare allele wll be favored by selecton: A rare mutaton whch affects the ftness of ts carrer and others wll sread f br > c where b "beneft" ncrease n ftness to recents of acton c "cost" loss n ftness to actor r "degree of relatedness". can use above rules to fnd r: r F Another vew of f : f s the roorton by whch heterozygosty s reduced relatve to a random matng grou wth the same allele frequences. If f 0, the oulaton s n Hardy-Wenberg roortons. If f 1, all ndvduals are homozygous. In general: Het f ( 1 f )Het 0 where Het 0 ( 1 ) q Consder, e.g., the Wahlund effect: Overall frequences of A and a are and q. If ndvduals from all suboulatons mate at random, exect to fnd q heterozygotes. Would actually observe q Var( ) heterozygotes. Thus, n ths case q Var( ) ( 1 f ) q. In order for ths to be true, f Var( ) q. [Suggested exercse: show ths.] Note: f 0 even though there s no defct of heterozygotes wthn suboulatons and no obvous nbreedng! II-5

Kinship Coefficients. Biostatistics 666

Kinship Coefficients. Biostatistics 666 Knshp Coeffcents Bostatstcs 666 Today s Lecture Genetc analyses requre relatonshps to be specfed Msspecfed relatonshps lead to tests of napproprate sze Inflated Type I error Decreased power Knshp Coeffcents

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Exerments-I MODULE III LECTURE - 2 EXPERIMENTAL DESIGN MODELS Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 2 We consder the models

More information

Quantitative Genetic Models Least Squares Genetic Model. Hardy-Weinberg (1908) Principle. change of allele & genotype frequency over generations

Quantitative Genetic Models Least Squares Genetic Model. Hardy-Weinberg (1908) Principle. change of allele & genotype frequency over generations Quanttatve Genetc Models Least Squares Genetc Model Hardy-Wenberg (1908) Prncple partton of effects P = G + E + G E P s phenotypc effect G s genetc effect E s envronmental effect G E s nteracton effect

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Recall that quantitative genetics is based on the extension of Mendelian principles to polygenic traits.

Recall that quantitative genetics is based on the extension of Mendelian principles to polygenic traits. BIOSTT/STT551, Statstcal enetcs II: Quanttatve Trats Wnter 004 Sources of varaton for multlocus trats and Handout Readng: Chapter 5 and 6. Extensons to Multlocus trats Recall that quanttatve genetcs s

More information

On New Selection Procedures for Unequal Probability Sampling

On New Selection Procedures for Unequal Probability Sampling Int. J. Oen Problems Comt. Math., Vol. 4, o. 1, March 011 ISS 1998-66; Coyrght ICSRS Publcaton, 011 www.-csrs.org On ew Selecton Procedures for Unequal Probablty Samlng Muhammad Qaser Shahbaz, Saman Shahbaz

More information

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1 Random varables Measure of central tendences and varablty (means and varances) Jont densty functons and ndependence Measures of assocaton (covarance and correlaton) Interestng result Condtonal dstrbutons

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Exerments-I MODULE II LECTURE - GENERAL LINEAR HYPOTHESIS AND ANALYSIS OF VARIANCE Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 3.

More information

LINEAR REGRESSION ANALYSIS. MODULE VIII Lecture Indicator Variables

LINEAR REGRESSION ANALYSIS. MODULE VIII Lecture Indicator Variables LINEAR REGRESSION ANALYSIS MODULE VIII Lecture - 7 Indcator Varables Dr. Shalabh Department of Maematcs and Statstcs Indan Insttute of Technology Kanpur Indcator varables versus quanttatve explanatory

More information

Maximum Likelihood Estimation and Binary Dependent Variables

Maximum Likelihood Estimation and Binary Dependent Variables MLE and Bnary Deendent Varables Maxmum Lkelhood Estmaton and Bnary Deendent Varables. Startng wth a Smle Examle: Bernoull Trals Lets start wth a smle examle: Teams A and B lay one another 0 tmes; A wns

More information

Naïve Bayes Classifier

Naïve Bayes Classifier 9/8/07 MIST.6060 Busness Intellgence and Data Mnng Naïve Bayes Classfer Termnology Predctors: the attrbutes (varables) whose values are used for redcton and classfcaton. Predctors are also called nut varables,

More information

+, where 0 x N - n. k k

+, where 0 x N - n. k k CO 745, Mdterm Len Cabrera. A multle choce eam has questons, each of whch has ossble answers. A student nows the correct answer to n of these questons. For the remanng - n questons, he checs the answers

More information

Michael Batty. Alan Wilson Plenary Session Entropy, Complexity, & Information in Spatial Analysis

Michael Batty. Alan Wilson Plenary Session Entropy, Complexity, & Information in Spatial Analysis Alan Wlson Plenary Sesson Entroy, Comlexty, & Informaton n Satal Analyss Mchael Batty m.batty@ucl.ac.uk @jmchaelbatty htt://www.comlexcty.nfo/ htt://www.satalcomlexty.nfo/ for Advanced Satal Analyss CentreCentre

More information

THERMODYNAMICS. Temperature

THERMODYNAMICS. Temperature HERMODYNMICS hermodynamcs s the henomenologcal scence whch descrbes the behavor of macroscoc objects n terms of a small number of macroscoc arameters. s an examle, to descrbe a gas n terms of volume ressure

More information

= z 20 z n. (k 20) + 4 z k = 4

= z 20 z n. (k 20) + 4 z k = 4 Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5

More information

Managing Capacity Through Reward Programs. on-line companion page. Byung-Do Kim Seoul National University College of Business Administration

Managing Capacity Through Reward Programs. on-line companion page. Byung-Do Kim Seoul National University College of Business Administration Managng Caacty Through eward Programs on-lne comanon age Byung-Do Km Seoul Natonal Unversty College of Busness Admnstraton Mengze Sh Unversty of Toronto otman School of Management Toronto ON M5S E6 Canada

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Linear system of the Schrödinger equation Notes on Quantum Mechanics

Linear system of the Schrödinger equation Notes on Quantum Mechanics Lnear sstem of the Schrödnger equaton Notes on Quantum Mechancs htt://quantum.bu.edu/notes/quantummechancs/lnearsstems.df Last udated Wednesda, October 9, 003 :0:08 Corght 003 Dan Dll (dan@bu.edu) Deartment

More information

A General Class of Selection Procedures and Modified Murthy Estimator

A General Class of Selection Procedures and Modified Murthy Estimator ISS 684-8403 Journal of Statstcs Volume 4, 007,. 3-9 A General Class of Selecton Procedures and Modfed Murthy Estmator Abdul Bast and Muhammad Qasar Shahbaz Abstract A new selecton rocedure for unequal

More information

Markov Chain Monte Carlo Lecture 6

Markov Chain Monte Carlo Lecture 6 where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways

More information

Power-sum problem, Bernoulli Numbers and Bernoulli Polynomials.

Power-sum problem, Bernoulli Numbers and Bernoulli Polynomials. Power-sum roblem, Bernoull Numbers and Bernoull Polynomals. Arady M. Alt Defnton 1 Power um Problem Fnd the sum n : 1... n where, n N or, usng sum notaton, n n n closed form. Recurrence for n Exercse Usng

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Expected Value and Variance

Expected Value and Variance MATH 38 Expected Value and Varance Dr. Neal, WKU We now shall dscuss how to fnd the average and standard devaton of a random varable X. Expected Value Defnton. The expected value (or average value, or

More information

Checking Pairwise Relationships. Lecture 19 Biostatistics 666

Checking Pairwise Relationships. Lecture 19 Biostatistics 666 Checkng Parwse Relatonshps Lecture 19 Bostatstcs 666 Last Lecture: Markov Model for Multpont Analyss X X X 1 3 X M P X 1 I P X I P X 3 I P X M I 1 3 M I 1 I I 3 I M P I I P I 3 I P... 1 IBD states along

More information

Lecture 10: May 6, 2013

Lecture 10: May 6, 2013 TTIC/CMSC 31150 Mathematcal Toolkt Sprng 013 Madhur Tulsan Lecture 10: May 6, 013 Scrbe: Wenje Luo In today s lecture, we manly talked about random walk on graphs and ntroduce the concept of graph expander,

More information

ISQS 6348 Final Open notes, no books. Points out of 100 in parentheses. Y 1 ε 2

ISQS 6348 Final Open notes, no books. Points out of 100 in parentheses. Y 1 ε 2 ISQS 6348 Fnal Open notes, no books. Ponts out of 100 n parentheses. 1. The followng path dagram s gven: ε 1 Y 1 ε F Y 1.A. (10) Wrte down the usual model and assumptons that are mpled by ths dagram. Soluton:

More information

Conservation of Angular Momentum = "Spin"

Conservation of Angular Momentum = Spin Page 1 of 6 Conservaton of Angular Momentum = "Spn" We can assgn a drecton to the angular velocty: drecton of = drecton of axs + rght hand rule (wth rght hand, curl fngers n drecton of rotaton, thumb ponts

More information

Algorithms for factoring

Algorithms for factoring CSA E0 235: Crytograhy Arl 9,2015 Instructor: Arta Patra Algorthms for factorng Submtted by: Jay Oza, Nranjan Sngh Introducton Factorsaton of large ntegers has been a wdely studed toc manly because of

More information

6. Hamilton s Equations

6. Hamilton s Equations 6. Hamlton s Equatons Mchael Fowler A Dynamcal System s Path n Confguraton Sace and n State Sace The story so far: For a mechancal system wth n degrees of freedom, the satal confguraton at some nstant

More information

Multipoint Analysis for Sibling Pairs. Biostatistics 666 Lecture 18

Multipoint Analysis for Sibling Pairs. Biostatistics 666 Lecture 18 Multpont Analyss for Sblng ars Bostatstcs 666 Lecture 8 revously Lnkage analyss wth pars of ndvduals Non-paraetrc BS Methods Maxu Lkelhood BD Based Method ossble Trangle Constrant AS Methods Covered So

More information

Linear Regression Analysis: Terminology and Notation

Linear Regression Analysis: Terminology and Notation ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented

More information

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA 4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected

More information

Statistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics )

Statistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics ) Ismor Fscher, 8//008 Stat 54 / -8.3 Summary Statstcs Measures of Center and Spread Dstrbuton of dscrete contnuous POPULATION Random Varable, numercal True center =??? True spread =???? parameters ( populaton

More information

Lecture 3: Probability Distributions

Lecture 3: Probability Distributions Lecture 3: Probablty Dstrbutons Random Varables Let us begn by defnng a sample space as a set of outcomes from an experment. We denote ths by S. A random varable s a functon whch maps outcomes nto the

More information

Solutions (mostly for odd-numbered exercises)

Solutions (mostly for odd-numbered exercises) Solutons (mostly for odd-numbered exercses) c 005 A. Coln Cameron and Pravn K. Trved "Mcroeconometrcs: Methods and Alcatons" 1. Chater 1: Introducton o exercses.. Chater : Causal and oncausal Models o

More information

First Year Examination Department of Statistics, University of Florida

First Year Examination Department of Statistics, University of Florida Frst Year Examnaton Department of Statstcs, Unversty of Florda May 7, 010, 8:00 am - 1:00 noon Instructons: 1. You have four hours to answer questons n ths examnaton.. You must show your work to receve

More information

A simple model of group selection that cannot be analyzed with inclusive fitness

A simple model of group selection that cannot be analyzed with inclusive fitness A smple model of group selecton that cannot be analyzed wth nclusve ftness 1.0 Hamlton s rule (bology) 0.5 0.0 0.0 0.5 1.0 Folk theorem (economcs) Ant-group selecton team Pro-group selecton team Jerry

More information

Title: Bounds and normalization of the composite linkage disequilibrium coefficient.

Title: Bounds and normalization of the composite linkage disequilibrium coefficient. Ttle: ounds and normalzaton of the composte lnkage dsequlbrum coeffcent. (Genetc Epdemology 004 7:5-57 uthor: Dmtr V. Zaykn * * GlaxoSmthKlne Inc., Research Trangle Park, NC Contact detals: Dmtr Zaykn

More information

Homework Notes Week 7

Homework Notes Week 7 Homework Notes Week 7 Math 4 Sprng 4 #4 (a Complete the proof n example 5 that s an nner product (the Frobenus nner product on M n n (F In the example propertes (a and (d have already been verfed so we

More information

Comparing two Quantiles: the Burr Type X and Weibull Cases

Comparing two Quantiles: the Burr Type X and Weibull Cases IOSR Journal of Mathematcs (IOSR-JM) e-issn: 78-578, -ISSN: 39-765X. Volume, Issue 5 Ver. VII (Se. - Oct.06), PP 8-40 www.osrjournals.org Comarng two Quantles: the Burr Tye X and Webull Cases Mohammed

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y) Secton 1.5 Correlaton In the prevous sectons, we looked at regresson and the value r was a measurement of how much of the varaton n y can be attrbuted to the lnear relatonshp between y and x. In ths secton,

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Lecture 3 Stat102, Spring 2007

Lecture 3 Stat102, Spring 2007 Lecture 3 Stat0, Sprng 007 Chapter 3. 3.: Introducton to regresson analyss Lnear regresson as a descrptve technque The least-squares equatons Chapter 3.3 Samplng dstrbuton of b 0, b. Contnued n net lecture

More information

Chapter 8 Indicator Variables

Chapter 8 Indicator Variables Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n

More information

Notes on Frequency Estimation in Data Streams

Notes on Frequency Estimation in Data Streams Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands Content. Inference on Regresson Parameters a. Fndng Mean, s.d and covarance amongst estmates.. Confdence Intervals and Workng Hotellng Bands 3. Cochran s Theorem 4. General Lnear Testng 5. Measures of

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

Allele Frequency Estimation

Allele Frequency Estimation Allele Frequency Estimation Examle: ABO blood tyes ABO genetic locus exhibits three alleles: A, B, and O Four henotyes: A, B, AB, and O Genotye A/A A/O A/B B/B B/O O/O Phenotye A A AB B B O Data: Observed

More information

ENG 8801/ Special Topics in Computer Engineering: Pattern Recognition. Memorial University of Newfoundland Pattern Recognition

ENG 8801/ Special Topics in Computer Engineering: Pattern Recognition. Memorial University of Newfoundland Pattern Recognition EG 880/988 - Specal opcs n Computer Engneerng: Pattern Recognton Memoral Unversty of ewfoundland Pattern Recognton Lecture 7 May 3, 006 http://wwwengrmunca/~charlesr Offce Hours: uesdays hursdays 8:30-9:30

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs hyscs 151 Lecture Canoncal Transformatons (Chater 9) What We Dd Last Tme Drect Condtons Q j Q j = = j, Q, j, Q, Necessary and suffcent j j for Canoncal Transf. = = j Q, Q, j Q, Q, Infntesmal CT

More information

On Unequal Probability Sampling Without Replacement Sample Size 2

On Unequal Probability Sampling Without Replacement Sample Size 2 Int J Oen Problems Com Math, Vol, o, March 009 On Unequal Probablt Samlng Wthout Relacement Samle Sze aser A Alodat Deartment of Mathematcs, Irbd atonal Unverst, Jordan e-mal: n_odat@ahoocom Communcated

More information

Logistic regression with one predictor. STK4900/ Lecture 7. Program

Logistic regression with one predictor. STK4900/ Lecture 7. Program Logstc regresson wth one redctor STK49/99 - Lecture 7 Program. Logstc regresson wth one redctor 2. Maxmum lkelhood estmaton 3. Logstc regresson wth several redctors 4. Devance and lkelhood rato tests 5.

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

More information

Hidden Markov Model Cheat Sheet

Hidden Markov Model Cheat Sheet Hdden Markov Model Cheat Sheet (GIT ID: dc2f391536d67ed5847290d5250d4baae103487e) Ths document s a cheat sheet on Hdden Markov Models (HMMs). It resembles lecture notes, excet that t cuts to the chase

More information

Introduction to Econometrics (3 rd Updated Edition, Global Edition) Solutions to Odd-Numbered End-of-Chapter Exercises: Chapter 13

Introduction to Econometrics (3 rd Updated Edition, Global Edition) Solutions to Odd-Numbered End-of-Chapter Exercises: Chapter 13 Introducton to Econometrcs (3 rd Updated Edton, Global Edton by James H. Stock and Mark W. Watson Solutons to Odd-Numbered End-of-Chapter Exercses: Chapter 13 (Ths verson August 17, 014 Stock/Watson -

More information

NEWTON S LAWS. These laws only apply when viewed from an inertial coordinate system (unaccelerated system).

NEWTON S LAWS. These laws only apply when viewed from an inertial coordinate system (unaccelerated system). EWTO S LAWS Consder two partcles. 1 1. If 1 0 then 0 wth p 1 m1v. 1 1 2. 1.. 3. 11 These laws only apply when vewed from an nertal coordnate system (unaccelerated system). consder a collecton of partcles

More information

Pattern Recognition. Approximating class densities, Bayesian classifier, Errors in Biometric Systems

Pattern Recognition. Approximating class densities, Bayesian classifier, Errors in Biometric Systems htt://.cubs.buffalo.edu attern Recognton Aromatng class denstes, Bayesan classfer, Errors n Bometrc Systems B. W. Slverman, Densty estmaton for statstcs and data analyss. London: Chaman and Hall, 986.

More information

Computational Biology Lecture 8: Substitution matrices Saad Mneimneh

Computational Biology Lecture 8: Substitution matrices Saad Mneimneh Computatonal Bology Lecture 8: Substtuton matrces Saad Mnemneh As we have ntroduced last tme, smple scorng schemes lke + or a match, - or a msmatch and -2 or a gap are not justable bologcally, especally

More information

Ph 219a/CS 219a. Exercises Due: Wednesday 23 October 2013

Ph 219a/CS 219a. Exercises Due: Wednesday 23 October 2013 1 Ph 219a/CS 219a Exercses Due: Wednesday 23 October 2013 1.1 How far apart are two quantum states? Consder two quantum states descrbed by densty operators ρ and ρ n an N-dmensonal Hlbert space, and consder

More information

What Independencies does a Bayes Net Model? Bayesian Networks: Independencies and Inference. Quick proof that independence is symmetric

What Independencies does a Bayes Net Model? Bayesian Networks: Independencies and Inference. Quick proof that independence is symmetric Bayesan Networks: Indeendences and Inference Scott Daves and ndrew Moore Note to other teachers and users of these sldes. ndrew and Scott would be delghted f you found ths source materal useful n gvng

More information

find (x): given element x, return the canonical element of the set containing x;

find (x): given element x, return the canonical element of the set containing x; COS 43 Sprng, 009 Dsjont Set Unon Problem: Mantan a collecton of dsjont sets. Two operatons: fnd the set contanng a gven element; unte two sets nto one (destructvely). Approach: Canoncal element method:

More information

SMALL AREA ESTIMATES FROM THE AMERICAN COMMUNITY SURVEY USING A HOUSING UNIT MODEL

SMALL AREA ESTIMATES FROM THE AMERICAN COMMUNITY SURVEY USING A HOUSING UNIT MODEL SMALL AREA ESTIMATES FROM THE AMERICAN COMMUNITY SURVEY USING A HOUSING UNIT MODEL Nanak Chand and Donald Malec U.S. Bureau of the Census Abstract The Amercan Communty Survey (ACS) s desgned to, ultmately,

More information

Ken G Dodds AgResearch, Invermay Agricultural Centre, Mosgiel, New Zealand

Ken G Dodds AgResearch, Invermay Agricultural Centre, Mosgiel, New Zealand STATISTICAL GENETICS Ken G Dodds AgResearch, Invermay Agrcultural Centre, Mosgel, New Zealand Keywords: Quanttatve Genetcs, Relatedness, Inbreedng, Breedng, Hertablty, Selecton, Lnkage, Mappng, Genetc

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

Lecture 08 Robust Design & Time symmetrical models

Lecture 08 Robust Design & Time symmetrical models WILD 75 - Analyss of Wldlfe Poulatons of 3 www. auburn.edu/~grandjb/wldo Lecture 8 Robust Desgn & Tme symmetrcal models Pollock s Robust desgn. Readngs endall, W. L., J. D. chols, and J. E. Hnes. 997.

More information

MATH 5707 HOMEWORK 4 SOLUTIONS 2. 2 i 2p i E(X i ) + E(Xi 2 ) ä i=1. i=1

MATH 5707 HOMEWORK 4 SOLUTIONS 2. 2 i 2p i E(X i ) + E(Xi 2 ) ä i=1. i=1 MATH 5707 HOMEWORK 4 SOLUTIONS CİHAN BAHRAN 1. Let v 1,..., v n R m, all lengths v are not larger than 1. Let p 1,..., p n [0, 1] be arbtrary and set w = p 1 v 1 + + p n v n. Then there exst ε 1,..., ε

More information

Machine Learning. Classification. Theory of Classification and Nonparametric Classifier. Representing data: Hypothesis (classifier) Eric Xing

Machine Learning. Classification. Theory of Classification and Nonparametric Classifier. Representing data: Hypothesis (classifier) Eric Xing Machne Learnng 0-70/5 70/5-78, 78, Fall 008 Theory of Classfcaton and Nonarametrc Classfer Erc ng Lecture, Setember 0, 008 Readng: Cha.,5 CB and handouts Classfcaton Reresentng data: M K Hyothess classfer

More information

Effects of Ignoring Correlations When Computing Sample Chi-Square. John W. Fowler February 26, 2012

Effects of Ignoring Correlations When Computing Sample Chi-Square. John W. Fowler February 26, 2012 Effects of Ignorng Correlatons When Computng Sample Ch-Square John W. Fowler February 6, 0 It can happen that ch-square must be computed for a sample whose elements are correlated to an unknown extent.

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

Using Genetic Algorithms in System Identification

Using Genetic Algorithms in System Identification Usng Genetc Algorthms n System Identfcaton Ecaterna Vladu Deartment of Electrcal Engneerng and Informaton Technology, Unversty of Oradea, Unverstat, 410087 Oradea, Româna Phone: +40259408435, Fax: +40259408408,

More information

Cathy Walker March 5, 2010

Cathy Walker March 5, 2010 Cathy Walker March 5, 010 Part : Problem Set 1. What s the level of measurement for the followng varables? a) SAT scores b) Number of tests or quzzes n statstcal course c) Acres of land devoted to corn

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore Sesson Outlne Introducton to classfcaton problems and dscrete choce models. Introducton to Logstcs Regresson. Logstc functon and Logt functon. Maxmum Lkelhood Estmator (MLE) for estmaton of LR parameters.

More information

BIO Lab 2: TWO-LEVEL NORMAL MODELS with school children popularity data

BIO Lab 2: TWO-LEVEL NORMAL MODELS with school children popularity data Lab : TWO-LEVEL NORMAL MODELS wth school chldren popularty data Purpose: Introduce basc two-level models for normally dstrbuted responses usng STATA. In partcular, we dscuss Random ntercept models wthout

More information

Pre-Talbot ANSS. Michael Andrews Department of Mathematics MIT. April 2, 2013

Pre-Talbot ANSS. Michael Andrews Department of Mathematics MIT. April 2, 2013 Pre-Talbot ANSS Mchael Andrews Deartment of Mathematcs MIT Arl 2, 203 The mage of J We have an unbased ma SO = colm n SO(n) colm n Ω n S n = QS 0 mang nto the -comonent of QS 0. The ma nduced by SO QS

More information

REDUCTION MODULO p. We will prove the reduction modulo p theorem in the general form as given by exercise 4.12, p. 143, of [1].

REDUCTION MODULO p. We will prove the reduction modulo p theorem in the general form as given by exercise 4.12, p. 143, of [1]. REDUCTION MODULO p. IAN KIMING We wll prove the reducton modulo p theorem n the general form as gven by exercse 4.12, p. 143, of [1]. We consder an ellptc curve E defned over Q and gven by a Weerstraß

More information

Credit Card Pricing and Impact of Adverse Selection

Credit Card Pricing and Impact of Adverse Selection Credt Card Prcng and Impact of Adverse Selecton Bo Huang and Lyn C. Thomas Unversty of Southampton Contents Background Aucton model of credt card solctaton - Errors n probablty of beng Good - Errors n

More information

STAT 3014/3914. Semester 2 Applied Statistics Solution to Tutorial 13

STAT 3014/3914. Semester 2 Applied Statistics Solution to Tutorial 13 STAT 304/394 Semester Appled Statstcs 05 Soluton to Tutoral 3. Note that s the total mleage for branch. a) -stage cluster sample Cluster branches N ; n 4) Element cars M 80; m 40) Populaton mean no. of

More information

Differentiating Gaussian Processes

Differentiating Gaussian Processes Dfferentatng Gaussan Processes Andrew McHutchon Aprl 17, 013 1 Frst Order Dervatve of the Posteror Mean The posteror mean of a GP s gven by, f = x, X KX, X 1 y x, X α 1 Only the x, X term depends on the

More information

HMMT February 2016 February 20, 2016

HMMT February 2016 February 20, 2016 HMMT February 016 February 0, 016 Combnatorcs 1. For postve ntegers n, let S n be the set of ntegers x such that n dstnct lnes, no three concurrent, can dvde a plane nto x regons (for example, S = {3,

More information

Statistics Chapter 4

Statistics Chapter 4 Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment

More information

Mixture of Gaussians Expectation Maximization (EM) Part 2

Mixture of Gaussians Expectation Maximization (EM) Part 2 Mture of Gaussans Eectaton Mamaton EM Part 2 Most of the sldes are due to Chrstoher Bsho BCS Summer School Eeter 2003. The rest of the sldes are based on lecture notes by A. Ng Lmtatons of K-means Hard

More information

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1 Problem Set 4 Suggested Solutons Problem (A) The market demand functon s the soluton to the followng utlty-maxmzaton roblem (UMP): The Lagrangean: ( x, x, x ) = + max U x, x, x x x x st.. x + x + x y x,

More information

Solutions for Tutorial 1

Solutions for Tutorial 1 Toc 1: Sem-drect roducts Solutons for Tutoral 1 1. Show that the tetrahedral grou s somorhc to the sem-drect roduct of the Klen four grou and a cyclc grou of order three: T = K 4 (Z/3Z). 2. Show further

More information

Application of artificial intelligence in earthquake forecasting

Application of artificial intelligence in earthquake forecasting Alcaton of artfcal ntellgence n earthquae forecastng Zhou Shengu, Wang Chengmn and Ma L Center for Analyss and Predcton of CSB, 63 Fuxng Road, Bejng 00036 P.R.Chna (e-mal zhou@ca.ac.cn; hone: 86 0 6827

More information

19. Genetic Drift. The biological context. There are four basic consequences of genetic drift:

19. Genetic Drift. The biological context. There are four basic consequences of genetic drift: 9. Genetic Drift Genetic drift is the alteration of gene frequencies due to sampling variation from one generation to the next. It operates to some degree in all finite populations, but can be significant

More information

18.1 Introduction and Recap

18.1 Introduction and Recap CS787: Advanced Algorthms Scrbe: Pryananda Shenoy and Shjn Kong Lecturer: Shuch Chawla Topc: Streamng Algorthmscontnued) Date: 0/26/2007 We contnue talng about streamng algorthms n ths lecture, ncludng

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 151 Lecture 22 Canoncal Transformatons (Chater 9) What We Dd Last Tme Drect Condtons Q j Q j = = j P, Q, P j, P Q, P Necessary and suffcent P j P j for Canoncal Transf. = = j Q, Q, P j

More information

SELECTED PROOFS. DeMorgan s formulas: The first one is clear from Venn diagram, or the following truth table:

SELECTED PROOFS. DeMorgan s formulas: The first one is clear from Venn diagram, or the following truth table: SELECTED PROOFS DeMorgan s formulas: The frst one s clear from Venn dagram, or the followng truth table: A B A B A B Ā B Ā B T T T F F F F T F T F F T F F T T F T F F F F F T T T T The second one can be

More information

MA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials

MA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have

More information

Lecture 3 Examples and Problems

Lecture 3 Examples and Problems Lecture 3 Examles and Problems Mechancs & thermodynamcs Equartton Frst Law of Thermodynamcs Ideal gases Isothermal and adabatc rocesses Readng: Elements Ch. 1-3 Lecture 3, 1 Wllam Thomson (1824 1907) a.k.a.

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

LECTURE # How does one test whether a population is in the HW equilibrium? (i) try the following example: Genotype Observed AA 50 Aa 0 aa 50

LECTURE # How does one test whether a population is in the HW equilibrium? (i) try the following example: Genotype Observed AA 50 Aa 0 aa 50 LECTURE #10 A. The Hardy-Weinberg Equilibrium 1. From the definitions of p and q, and of p 2, 2pq, and q 2, an equilibrium is indicated (p + q) 2 = p 2 + 2pq + q 2 : if p and q remain constant, and if

More information

MAGNETISM MAGNETIC DIPOLES

MAGNETISM MAGNETIC DIPOLES MAGNETISM We now turn to magnetsm. Ths has actually been used for longer than electrcty. People were usng compasses to sal around the Medterranean Sea several hundred years BC. However t was not understood

More information

Processes of Evolution

Processes of Evolution 15 Processes of Evolution Forces of Evolution Concept 15.4 Selection Can Be Stabilizing, Directional, or Disruptive Natural selection can act on quantitative traits in three ways: Stabilizing selection

More information

Definition. Measures of Dispersion. Measures of Dispersion. Definition. The Range. Measures of Dispersion 3/24/2014

Definition. Measures of Dispersion. Measures of Dispersion. Definition. The Range. Measures of Dispersion 3/24/2014 Measures of Dsperson Defenton Range Interquartle Range Varance and Standard Devaton Defnton Measures of dsperson are descrptve statstcs that descrbe how smlar a set of scores are to each other The more

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Vapnik-Chervonenkis theory

Vapnik-Chervonenkis theory Vapnk-Chervonenks theory Rs Kondor June 13, 2008 For the purposes of ths lecture, we restrct ourselves to the bnary supervsed batch learnng settng. We assume that we have an nput space X, and an unknown

More information