1 / 11

Opinionated

Opinionated. Lessons. in Statistics. by Bill Press. #6 The Towne Family Tree. Neat example (with some biology):. Individual identity, or ancestry, can be determined by “variable length short tandem repeats” (STRs) in the genome.

gareth-boone
Download Presentation

Opinionated

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Opinionated Lessons in Statistics by Bill Press #6 The Towne Family Tree Professor William H. Press, Department of Computer Science, the University of Texas at Austin

  2. Neat example (with some biology): Individual identity, or ancestry, can be determined by “variable length short tandem repeats” (STRs) in the genome. ~0.5% mutation prob per STR per generation (though highly variable) if use Y chromosome only, get paternal ancestry There are companies that sell “certificates” with your genotype. A bit opportunistic, since in a few years your whole genome will be sequenced by your health plan. Professor William H. Press, Department of Computer Science, the University of Texas at Austin

  3. Margaret, my ex-wife, is really into the Towne family.(And, she’s neither a biologist nor a Towne.) Professor William H. Press, Department of Computer Science, the University of Texas at Austin

  4. Here’s data from Margaret on 8 recent Townes (identified only by T code).(We’ll use this data several times in the next few of lectures.) T3 T11 T2 T5 T13 T4 T8 T6 Professor William H. Press, Department of Computer Science, the University of Texas at Austin

  5. Or, just showing the changes from what we impute for William: T3 T11 T2 T5 T13 T4 T8 T6 Professor William H. Press, Department of Computer Science, the University of Texas at Austin

  6. a b c d e Unraveling dependencies Another important idea is “conditional independence” Example: b and e are “conditionally independent given a” while b and d are not conditionally independent given a: ? Professor William H. Press, Department of Computer Science, the University of Texas at Austin

  7. Let’s do a Bayesian estimation of the parameter r, the mutation rate per locus per generation. let’s assume no back mutations!(their effect on this data set would be small) D=0 N=1 D=0 N=3 bin(0,3x37,r) N=6 D=0 bin(0,3x37,r) N=9 D=0 bin(0,6x37,r) D=5 N=10 D=4 (of 12) D=23 D=3 N=11 bin(1,10x37,r) D=0 D=1 D=1 Professor William H. Press, Department of Computer Science, the University of Texas at Austin bin(1,5x37,r) bin(0,5x37,r) bin(1,11x37,r)

  8. 1 0 1 0 r r 1 R R ( j ) ( ) d P d d l P 1 0 t t a a p a r a m e e r s r r r o g = = r r r 1 ( j ) ( j ) ( ) ( j ) d d d P P P P t t t £ £ r a a / a a r r / a a r ( j ) ( ) ( ) ( ) ( ) d b b b b P i i i i 0 3 3 7 0 3 3 7 1 5 3 7 0 5 3 7 t £ £ £ £ a a r n r n r n r n r = r ; ; ; ; ; ; ; ; ( ) ( ) ( ) b b b i i i 0 6 3 7 1 1 1 3 7 3 1 0 3 7 £ £ £ £ n r n r n r ; ; ; ; ; ; So we have a statistical model for the data!Any statistical model is just a way to compute Our model is not “exact”, but statistical models rarely (never?) are. neglects backmutations assumes single probability for all loci etc. The model is: Bayes estimation of the parameter: What kind of prior is this???It is called “log-uniform” The log-uniform prior has equal probability ineach order of magnitude.It is often taken as the non-informative prior when you don’t even know the order of magnitude of the (positive) quantity. It is an “improper prior” since its integral is infinite. This is almost always ok, but it is possible to construct paradoxes with improper priors (e.g., the “marginalization paradox”) Professor William H. Press, Department of Computer Science, the University of Texas at Austin

  9. ( j ) d P t r a a Here is the plot of the (normalized) This is (almost) real biology. We’ve measured the mutation probability, per locus per generation of Y chromosome STRs. This tells us something about the actual DNA replication machinery! Professor William H. Press, Department of Computer Science, the University of Texas at Austin

  10. It really did matter (a bit) that we sorted out the conditional dependencies correctly.Here’s a comparison to doing it wrong by assuming all data independent: The true dependencies allow somewhat larger values of r, because we don’t wrongly count the D=0 branches multiple times We’ll come back to the Towne family for some fancier stuff later! Ignoring conditional dependencies and just multiplying the probabilities of the data as if they were independent is called naïve Bayes. People often do this. It is mathematically incorrect, but sometimes it is all you can do! Professor William H. Press, Department of Computer Science, the University of Texas at Austin

  11. ( j ) [ ] ( ) ( ) d l d d l b b P i i 5 9 3 7 4 1 0 1 2 t £ £ £ a a r o m o e n r n r = ; ; ; ; now include these A small cloud: The way we “trimmed” the data mattered. (And should trouble us a bit!) Here’s the effect of including T11 and T13, both of which seemed to be outliers: Editing outliers is a tricky issue that we will return to when we learn about mixture models. Professor William H. Press, Department of Computer Science, the University of Texas at Austin

More Related