in Computer Simulations for Bioinformatics
|
|
- Baldric Cummings
- 5 years ago
- Views:
Transcription
1 Bi04a_ Unit 04a: Stochastic Processes and their Applications in Computer Simulations for Bioinformatics Stochastic Processes and their Applications in Computer Simulations for Bioinformatics Basic Probability Concepts revisited Crude Monte Carlo Markov-Chain Monte Carlo a) Metropolis-Hastings b) Gibbs-Sampling Real Example: How to find the Consensus Sequence (Unit 4c, Start) Simulated Annealing using Markov-Chain-Monte Carlo Real Example: How to find the Consensus Sequence (Unit 4c, ctd) Genetic Algorithms a) General applications used to optimize structure of multidimensional objects: molecule conformations, sequence alignments, etc. Real example: Multiple sequence alignment by GA Program SAGA b) Genetic Computing: Genetic Algorithms used to optimize computer programs Bi04a_2
2 Basic Probability Concepts revisited Bi04a_3 Discrete distributions absolute frequencies relative frequencies limiting behaviours as N probabilites & cumulative probabilities Continuous distributions absolute frequencies relative frequencies limiting behaviours as N probabilites & cumulative probabilities Basic Concepts revisited Bi04a_4. Discrete Distributions random variable takes only a few possible values example: ξ {, 2, 3, 4, 5, 6} for (normal) dice N-times throwing a dice: {, 2, 3, 4, 5, 6, } 6 i = i= H H H H H H absolute Frequences H 3 -times dice shows 3 H N Summation Condition for absolute Frequences 2 6 h h2 h6 relative frequences 6 i= H H H =, =,... = N N N h = i Normalization Condition for relative Frequences " irgendein Wert muß ja schließlich auftreten..." 2
3 Basic Concepts revisited Bi04a_5 Limiting behaviour as N {h,...h 6 } {p...p 6 } Hi p Probability is limit of i = lim = lim hi N N N relative frequencies as N 6 6 hi i= i= = p = i Normalization condition to Generalize: a) Consider an arbitrary number of values (L instead of 6) b) Consider unequal probabilities p i p j for i j L L hi i= i= = p = i Basic Concepts revisited Bi04a_6 Example for limiting behaviour of relative frequencies as N {h, h 2,...h 6 } N {p, p 2,...p 6 } small N: {h, h 2,...h 6 } we obtain very different results for each lap large N: {h, h 2,...h 6 } we obtain more similar results for each lap 3
4 Basic Concepts revisited Bi04a_7 p(i) = Pr(ξ=i) probability distribution F(i) = Pr(ξ i) () p() i F i = i l= cumulative probability Basic Concepts revisited Bi04a_8 A simple example: normal dice Mafia dice p F / p F /6 L l= () = () = p() l F i p l i l=
5 Basic Concepts revisited Bi04a_9 2. Continuous Distributions random variable takes arbitrary (real) values within given intervall ex.: ξ, 0 ξ finite intervals -2 ξ 5 0 ξ infinite intervals - ξ Basic Concepts revisited Bi04a_0 Continuous Distributions, ctd. p(x) uniform distribution x x+dx p(x) is defined such that Pr(x ξ x + dx) = p(x)dx dx... infinitesimally small interval shaded area ( ) b Pr a ξ b = p ( ξ ) dξ a ξ Since the random variable can take any (in between) value, the concept of discrete probability values must be generalized to a continuous probability distribution function ( ) ( ) Pr ξ = p ξ dξ = [a,b]... real intervall Normalization for probability density function (p.d.f.) irgendwo muß ξ ja liegen 5
6 Basic Concepts revisited Bi04a_ Continuous Distributions, ctd. specifically consider interval [-,x] p p(x) F(x) if this holds, we know from calculus: x P(- ξ x ) = x ξ ( ξ) ξ = ( ) = Pr ( ξ ) p d F x x ( ) F( x) p x = d dx F(- ) = 0 F(+ ) = F(x) is increasing since p(x) 0 probability density function (p.d.f.) cumulative distribution function Basic Concepts revisited Bi04a_2 Example: linearly increasing p.d.f. on [0,]: p(x) = 2x p(x) F(x) Normalization: Cumulative distribution function: 2 2x p( x) dx= = 0= x 2 x 2 p( ξ) dξ = ξ = x 0 0 6
7 Basic Concepts revisited Bi04a_3 from Kalos (986) Basic Concepts revisited Bi04a_4 how to get a random variable ξ? physically: throw 0 sided prism to construct real number as 0.n n 2 n 3 n 4...n L [0,) from Morgan (984) on a computer: use random number generator! (see expendix) 7
8 Expendix : How to generate random numbers Bi04a_5 Real random numbers (from nature, i.e. physical process) from Morgan (984) Pseudo random numbers (PRNs) (from computer programs) x n+ =a x n +b (mod m) congruential method e.g. a = 573, b = 9, m = 000, x 0 =89 x ξ = n n [0,] m seed uniform distributed Really random versus Pseudo random. Bi04a_6 Real random numbers (from nature, i.e. physical process) no prediction possible randomness (of e.g. dice) due to uncertaincy principle of quantum physics statistical tests will not reveal systematic predictability awkward for simulation Pseudo random numbers (PRNs) (from computer programs) strictly predictable randomness due to algorithm, its parameters & seed statistical tests will not easily reveal systematic predictability however: beware of loops handy for simulation 8
9 Expendix : PRNs Quality Bi04a_7 Bad Generator from Morgan (984) Improved Generator Expendix : PRNs from most common distributions Bi04a_8 ξ [0,] uniformly distributed is available in (almost) every programing language (library) ξ N(µ,σ) normal (Gaussian) distribution mean standard deviation is available in many programing language (libraries) from Kalos (986) from Kalos (986) 9
10 Expendix : Ho to get PRNs from other distributions Bi04a_9 statistical and mathematical program libraries deriving desired distribution from uniform PRNs on [0,]: inversion method other, even more sophisticated methods... rejection method (works for arbitrary distributions wanted!) table lookup-method Bi04a_20 Expendix : Linear transformations of PRNs uniform PRNs uniform PRNs Gaussian PRNs Gaussian PRNs arbitrary distributed PRNs (? not generally predictable) 0 N 0 - trafo 0 trafo + ξ N. ξ [0,N] CAUTION:! what happens at start & end of interval! ξ 2 ξ- [-,+] 0 A B- trafo ξ A+(B-A) ξ [A,B] 0 0
11 Expendix : General transformations of PRNs Bi04a_2 from Morgan (984) consider cumulative distribution function for argumentation And for details: see Morgan (984), p.29 or Kalos (986), p.40 y ( ) = Pr ( ) = ( ) = ( ) F y Y y X x F X ( ) dfx x dx fy ( y) = dx dy dx fy( y) = fx( x) dy dy = f ( x) dx x d dy Apply differentiatiation operator to both sides of equation! Use chain-rule for differentiation taking absolute values will make formula valid for decreasing and increasing transformation functions! Expendix : General transformations of PRNs Bi04a_22 An illustration of the transformation y = x relating the densities f X (x) = e -x, f Y (y) = 2ye -y2. The shaded regions have equal areas. from Morgan (984)
12 Expendix : Random numbers from libraries Bi04a_23 SAS NAG IMSL NORMAL generates a normally distributed G05CAF uniform over (0,) nextbeta beta distribution pseudo-random variate RANBIN generates an observation from a G05DAF uniform over (a,b) nextbinomial binomial distribution binomial distribution RANCAU generates a Cauchy deviate G05DBF exponential nextcauchy Cauchy distribution RANEXP generates an exponential deviate G05DDF Normal nextchisquared Chi-squared distribution RANGAM generates an observation from a gamma distribution G05DYF discrete uniform nextexponential standard exponential distribution RANNOR generates a normal deviate G05DRF Poisson nextexponentialmix mixture of two exponential distributions RANPOI generates an observation from a Poisson distribution G05FEF Beta distribution (multiple) nextgamma standard gamma distribution RANTBL generates deviates from a tabled G05DFF Cauchy distribution nextgeometric geometric distribution probability mass function RANTRI generates an observation from a triangular distribution G05DHF Chi-square distribution nexthypergeometric hypergeometric distribution RANUNI generates a uniform deviate G05DKF F-distribution nextlogarithmic logarithmic distribution UNIFORM generates a pseudo-random variate G05FFF Gamma distribution (multiple) nextlognormal lognormal distribution uniformly distributed on the interval (0,) G05DCF Logistic distribution nextmultivariatenormal multivariate normal distribution G05DEF Lognormal distribution nextnegativebinomial negative binomial distribution G05DJF Student's t-distribution nextnormal standard normal distribution using an inverse CDF method G05FSF von Mises distribution nextnormalar standard normal distribution using an acceptance/rejection method G05DPF Weibull distribution nextpoisson Poisson distribution nextstudentst Student's t distribution next Triangular triangular distribution nextvonmises von Mises distribution nextweibull Weibull distribution Expendix : Rejection method Bi04a_24 p wanted distribution (for some peculiar reason...) a recipe: generate ξ [a,b], uniform generate η [0, max(p)], uniform b ξ Rejection method is simple! works for the most fancy distributions normalizes itself! if η<p(ξ) take ξ as next random number ( Points) otherwise discard (ξ, η) ( Points) 2
13 Expendix : Rejection method, ctd. Bi04a_25 rejection method is simple but: may be inefficient: Simple but dead slow! Expendix : Select randomly from a group of n individuals according to weights (probabilities) ( table lookup method for a discrete random variable ) Bi04a_26 ) Case : equal probabilities (=weights) ξ individual # 2 3 n equidistant intervals 0 /n 2/n 3/n (n-)/n n/n = recipe: draw pseudo random number ξ from uniform distribution over [0,] ξ individual selected: i = int +, = min (, ) / i i n n 3
14 Expendix : Select randomly from a group of n individuals according to weights (probabilities) ( table lookup method for a discrete random variable ) Bi04a_27 2) Case 2: different probabilities (=weights), e.g. p i proportional e -E(s i )/τ ξ individual # n equidistant intervals 0 n ( )/ n Esi τ Es ( i) / τ pi = e e Normalization : p i= i= i= recipe: draw pseudo random number ξ from uniform distribution over [0,] individual i selected so that: i p < ξ l l= l= i p l Literature on Random Numbers Bi04a_28 Morgan,B.J.T Elements of Simulation. Chapman and Hall, New York. Rosanow,J.A Wahrscheinlichkeitstheorie. Rowohlt Taschenbuch Verlag, Hamburg. Kalos,M.H. and P.A.Whitlock 986. Monte Carlo methods. Wiley, New York. SAS Institute SAS User's Guide: Basics. Cary. NAg Fortran Library Mark 6. Oxford. JMSL Reference Manual Visual Numerics, San Ramon. 4
16 : Markov Chain Monte Carlo (MCMC)
10-708: Probabilistic Graphical Models 10-708, Spring 2014 16 : Markov Chain Monte Carlo MCMC Lecturer: Matthew Gormley Scribes: Yining Wang, Renato Negrinho 1 Sampling from low-dimensional distributions
More informationMalvin H. Kalos, Paula A. Whitlock. Monte Carlo Methods. Second Revised and Enlarged Edition WILEY- BLACKWELL. WILEY-VCH Verlag GmbH & Co.
Malvin H. Kalos, Paula A. Whitlock Monte Carlo Methods Second Revised and Enlarged Edition WILEY- BLACKWELL WILEY-VCH Verlag GmbH & Co. KGaA v I Contents Preface to the Second Edition IX Preface to the
More informationExploring Monte Carlo Methods
Exploring Monte Carlo Methods William L Dunn J. Kenneth Shultis AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO ELSEVIER Academic Press Is an imprint
More information2 Random Variable Generation
2 Random Variable Generation Most Monte Carlo computations require, as a starting point, a sequence of i.i.d. random variables with given marginal distribution. We describe here some of the basic methods
More informationNAG Fortran Library Chapter Introduction. G01 Simple Calculations on Statistical Data
G01 Simple Calculations on Statistical Data Introduction G01 NAG Fortran Library Chapter Introduction G01 Simple Calculations on Statistical Data Contents 1 Scope of the Chapter... 2 2 Background to the
More informationMonte Carlo and cold gases. Lode Pollet.
Monte Carlo and cold gases Lode Pollet lpollet@physics.harvard.edu 1 Outline Classical Monte Carlo The Monte Carlo trick Markov chains Metropolis algorithm Ising model critical slowing down Quantum Monte
More informationIndependent Events. Two events are independent if knowing that one occurs does not change the probability of the other occurring
Independent Events Two events are independent if knowing that one occurs does not change the probability of the other occurring Conditional probability is denoted P(A B), which is defined to be: P(A and
More informationNotation Precedence Diagram
Notation Precedence Diagram xix xxi CHAPTER 1 Introduction 1 1.1. Systems, Models, and Simulation 1 1.2. Verification, Approximation, and Validation 8 1.2.1. Verifying a Program 9 1.2.2. Approximation
More informationMonte Carlo Methods. Handbook of. University ofqueensland. Thomas Taimre. Zdravko I. Botev. Dirk P. Kroese. Universite de Montreal
Handbook of Monte Carlo Methods Dirk P. Kroese University ofqueensland Thomas Taimre University ofqueensland Zdravko I. Botev Universite de Montreal A JOHN WILEY & SONS, INC., PUBLICATION Preface Acknowledgments
More informationToday: Fundamentals of Monte Carlo
Today: Fundamentals of Monte Carlo What is Monte Carlo? Named at Los Alamos in 940 s after the casino. Any method which uses (pseudo)random numbers as an essential part of the algorithm. Stochastic - not
More informationSampling Distributions Allen & Tildesley pgs and Numerical Recipes on random numbers
Sampling Distributions Allen & Tildesley pgs. 345-351 and Numerical Recipes on random numbers Today I explain how to generate a non-uniform probability distributions. These are used in importance sampling
More information19 : Slice Sampling and HMC
10-708: Probabilistic Graphical Models 10-708, Spring 2018 19 : Slice Sampling and HMC Lecturer: Kayhan Batmanghelich Scribes: Boxiang Lyu 1 MCMC (Auxiliary Variables Methods) In inference, we are often
More informationData Analysis I. Dr Martin Hendry, Dept of Physics and Astronomy University of Glasgow, UK. 10 lectures, beginning October 2006
Astronomical p( y x, I) p( x, I) p ( x y, I) = p( y, I) Data Analysis I Dr Martin Hendry, Dept of Physics and Astronomy University of Glasgow, UK 10 lectures, beginning October 2006 4. Monte Carlo Methods
More informationMonte Carlo Simulations
Monte Carlo Simulations What are Monte Carlo Simulations and why ones them? Pseudo Random Number generators Creating a realization of a general PDF The Bootstrap approach A real life example: LOFAR simulations
More informationToday: Fundamentals of Monte Carlo
Today: Fundamentals of Monte Carlo What is Monte Carlo? Named at Los Alamos in 1940 s after the casino. Any method which uses (pseudo)random numbers as an essential part of the algorithm. Stochastic -
More informationA Probability Primer. A random walk down a probabilistic path leading to some stochastic thoughts on chance events and uncertain outcomes.
A Probability Primer A random walk down a probabilistic path leading to some stochastic thoughts on chance events and uncertain outcomes. Are you holding all the cards?? Random Events A random event, E,
More informationSampling Distributions Allen & Tildesley pgs and Numerical Recipes on random numbers
Sampling Distributions Allen & Tildesley pgs. 345-351 and Numerical Recipes on random numbers Today I explain how to generate a non-uniform probability distributions. These are used in importance sampling
More information( ) ( ) Monte Carlo Methods Interested in. E f X = f x d x. Examples:
Monte Carlo Methods Interested in Examples: µ E f X = f x d x Type I error rate of a hypothesis test Mean width of a confidence interval procedure Evaluating a likelihood Finding posterior mean and variance
More informationChapter Learning Objectives. Probability Distributions and Probability Density Functions. Continuous Random Variables
Chapter 4: Continuous Random Variables and Probability s 4-1 Continuous Random Variables 4-2 Probability s and Probability Density Functions 4-3 Cumulative Functions 4-4 Mean and Variance of a Continuous
More informationProbability Distributions Columns (a) through (d)
Discrete Probability Distributions Columns (a) through (d) Probability Mass Distribution Description Notes Notation or Density Function --------------------(PMF or PDF)-------------------- (a) (b) (c)
More informationPattern Recognition and Machine Learning. Bishop Chapter 11: Sampling Methods
Pattern Recognition and Machine Learning Chapter 11: Sampling Methods Elise Arnaud Jakob Verbeek May 22, 2008 Outline of the chapter 11.1 Basic Sampling Algorithms 11.2 Markov Chain Monte Carlo 11.3 Gibbs
More information6 Markov Chain Monte Carlo (MCMC)
6 Markov Chain Monte Carlo (MCMC) The underlying idea in MCMC is to replace the iid samples of basic MC methods, with dependent samples from an ergodic Markov chain, whose limiting (stationary) distribution
More informationComputer Vision Group Prof. Daniel Cremers. 11. Sampling Methods
Prof. Daniel Cremers 11. Sampling Methods Sampling Methods Sampling Methods are widely used in Computer Science as an approximation of a deterministic algorithm to represent uncertainty without a parametric
More informationGeneration from simple discrete distributions
S-38.3148 Simulation of data networks / Generation of random variables 1(18) Generation from simple discrete distributions Note! This is just a more clear and readable version of the same slide that was
More informationSlides 5: Random Number Extensions
Slides 5: Random Number Extensions We previously considered a few examples of simulating real processes. In order to mimic real randomness of events such as arrival times we considered the use of random
More informationMarkov Chain Monte Carlo (MCMC)
Markov Chain Monte Carlo (MCMC Dependent Sampling Suppose we wish to sample from a density π, and we can evaluate π as a function but have no means to directly generate a sample. Rejection sampling can
More informationSimulating Random Variables
Simulating Random Variables Timothy Hanson Department of Statistics, University of South Carolina Stat 740: Statistical Computing 1 / 23 R has many built-in random number generators... Beta, gamma (also
More informationToday: Fundamentals of Monte Carlo
Today: Fundamentals of Monte Carlo What is Monte Carlo? Named at Los Alamos in 1940 s after the casino. Any method which uses (pseudo)random numbers as an essential part of the algorithm. Stochastic -
More informationIrr. Statistical Methods in Experimental Physics. 2nd Edition. Frederick James. World Scientific. CERN, Switzerland
Frederick James CERN, Switzerland Statistical Methods in Experimental Physics 2nd Edition r i Irr 1- r ri Ibn World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI HONG KONG TAIPEI CHENNAI CONTENTS
More informationCIVL Continuous Distributions
CIVL 3103 Continuous Distributions Learning Objectives - Continuous Distributions Define continuous distributions, and identify common distributions applicable to engineering problems. Identify the appropriate
More informationfunctions Poisson distribution Normal distribution Arbitrary functions
Physics 433: Computational Physics Lecture 6 Random number distributions Generation of random numbers of various distribuition functions Normal distribution Poisson distribution Arbitrary functions Random
More informationCheng Soon Ong & Christian Walder. Canberra February June 2018
Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 2018 Outlines Overview Introduction Linear Algebra Probability Linear Regression
More informationRandom processes and probability distributions. Phys 420/580 Lecture 20
Random processes and probability distributions Phys 420/580 Lecture 20 Random processes Many physical processes are random in character: e.g., nuclear decay (Poisson distributed event count) P (k, τ) =
More informationNumerical Analysis for Statisticians
Kenneth Lange Numerical Analysis for Statisticians Springer Contents Preface v 1 Recurrence Relations 1 1.1 Introduction 1 1.2 Binomial CoefRcients 1 1.3 Number of Partitions of a Set 2 1.4 Horner's Method
More informationS6880 #7. Generate Non-uniform Random Number #1
S6880 #7 Generate Non-uniform Random Number #1 Outline 1 Inversion Method Inversion Method Examples Application to Discrete Distributions Using Inversion Method 2 Composition Method Composition Method
More informationStatistical Methods in Particle Physics
Statistical Methods in Particle Physics Lecture 3 October 29, 2012 Silvia Masciocchi, GSI Darmstadt s.masciocchi@gsi.de Winter Semester 2012 / 13 Outline Reminder: Probability density function Cumulative
More informationContinuous Random Variables and Continuous Distributions
Continuous Random Variables and Continuous Distributions Continuous Random Variables and Continuous Distributions Expectation & Variance of Continuous Random Variables ( 5.2) The Uniform Random Variable
More informationMAS1302 Computational Probability and Statistics
MAS1302 Computational Probability and Statistics April 23, 2008 3. Simulating continuous random behaviour 3.1 The Continuous Uniform U(0,1) Distribution We have already used this random variable a great
More informationCh3. Generating Random Variates with Non-Uniform Distributions
ST4231, Semester I, 2003-2004 Ch3. Generating Random Variates with Non-Uniform Distributions This chapter mainly focuses on methods for generating non-uniform random numbers based on the built-in standard
More informationGeneral Principles in Random Variates Generation
General Principles in Random Variates Generation E. Moulines and G. Fort Telecom ParisTech June 2015 Bibliography : Luc Devroye, Non-Uniform Random Variate Generator, Springer-Verlag (1986) available on
More informationContinuous Random Variables
Continuous Random Variables Recall: For discrete random variables, only a finite or countably infinite number of possible values with positive probability. Often, there is interest in random variables
More information2905 Queueing Theory and Simulation PART IV: SIMULATION
2905 Queueing Theory and Simulation PART IV: SIMULATION 22 Random Numbers A fundamental step in a simulation study is the generation of random numbers, where a random number represents the value of a random
More informationMarkov Chain Monte Carlo
1 Motivation 1.1 Bayesian Learning Markov Chain Monte Carlo Yale Chang In Bayesian learning, given data X, we make assumptions on the generative process of X by introducing hidden variables Z: p(z): prior
More informationAppendix F. Computational Statistics Toolbox. The Computational Statistics Toolbox can be downloaded from:
Appendix F Computational Statistics Toolbox The Computational Statistics Toolbox can be downloaded from: http://www.infinityassociates.com http://lib.stat.cmu.edu. Please review the readme file for installation
More information2. Variance and Covariance: We will now derive some classic properties of variance and covariance. Assume real-valued random variables X and Y.
CS450 Final Review Problems Fall 08 Solutions or worked answers provided Problems -6 are based on the midterm review Identical problems are marked recap] Please consult previous recitations and textbook
More informationWeek 1 Quantitative Analysis of Financial Markets Distributions A
Week 1 Quantitative Analysis of Financial Markets Distributions A Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October
More informationComputer Vision Group Prof. Daniel Cremers. 11. Sampling Methods: Markov Chain Monte Carlo
Group Prof. Daniel Cremers 11. Sampling Methods: Markov Chain Monte Carlo Markov Chain Monte Carlo In high-dimensional spaces, rejection sampling and importance sampling are very inefficient An alternative
More informationComputer Science, Informatik 4 Communication and Distributed Systems. Simulation. Discrete-Event System Simulation. Dr.
Simulation Discrete-Event System Simulation Chapter 6 andom-variate Generation Purpose & Overview Develop understanding of generating samples from a specified distribution as input to a simulation model.
More informationLecture 6: Monte-Carlo methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 6: Monte-Carlo methods Readings Recommended: handout on Classes site, from notes by Weinan E, Tiejun Li, and Eric Vanden-Eijnden Optional:
More informationReview of Statistical Terminology
Review of Statistical Terminology An experiment is a process whose outcome is not known with certainty. The experiment s sample space S is the set of all possible outcomes. A random variable is a function
More informationMethods of Data Analysis Random numbers, Monte Carlo integration, and Stochastic Simulation Algorithm (SSA / Gillespie)
Methods of Data Analysis Random numbers, Monte Carlo integration, and Stochastic Simulation Algorithm (SSA / Gillespie) Week 1 1 Motivation Random numbers (RNs) are of course only pseudo-random when generated
More information57:022 Principles of Design II Final Exam Solutions - Spring 1997
57:022 Principles of Design II Final Exam Solutions - Spring 1997 Part: I II III IV V VI Total Possible Pts: 52 10 12 16 13 12 115 PART ONE Indicate "+" if True and "o" if False: + a. If a component's
More informationCommon ontinuous random variables
Common ontinuous random variables CE 311S Earlier, we saw a number of distribution families Binomial Negative binomial Hypergeometric Poisson These were useful because they represented common situations:
More informationDistribution Fitting (Censored Data)
Distribution Fitting (Censored Data) Summary... 1 Data Input... 2 Analysis Summary... 3 Analysis Options... 4 Goodness-of-Fit Tests... 6 Frequency Histogram... 8 Comparison of Alternative Distributions...
More informationMonte Carlo Methods. PHY 688: Numerical Methods for (Astro)Physics
Monte Carlo Methods Random Numbers How random is random? From Pang: we want Long period before sequence repeats Little correlation between numbers (plot ri+1 vs ri should fill the plane) Fast Typical random
More informationProbability Distributions for Continuous Variables. Probability Distributions for Continuous Variables
Probability Distributions for Continuous Variables Probability Distributions for Continuous Variables Let X = lake depth at a randomly chosen point on lake surface If we draw the histogram so that the
More informationLecture 7 and 8: Markov Chain Monte Carlo
Lecture 7 and 8: Markov Chain Monte Carlo 4F13: Machine Learning Zoubin Ghahramani and Carl Edward Rasmussen Department of Engineering University of Cambridge http://mlg.eng.cam.ac.uk/teaching/4f13/ Ghahramani
More informationTABLE OF CONTENTS CHAPTER 1 COMBINATORIAL PROBABILITY 1
TABLE OF CONTENTS CHAPTER 1 COMBINATORIAL PROBABILITY 1 1.1 The Probability Model...1 1.2 Finite Discrete Models with Equally Likely Outcomes...5 1.2.1 Tree Diagrams...6 1.2.2 The Multiplication Principle...8
More informationIntroduction to Machine Learning CMU-10701
Introduction to Machine Learning CMU-10701 Markov Chain Monte Carlo Methods Barnabás Póczos & Aarti Singh Contents Markov Chain Monte Carlo Methods Goal & Motivation Sampling Rejection Importance Markov
More information2008 Winton. Review of Statistical Terminology
1 Review of Statistical Terminology 2 Formal Terminology An experiment is a process whose outcome is not known with certainty The experiment s sample space S is the set of all possible outcomes. A random
More informationChapter 5: Generating Random Numbers from Distributions
Chapter 5: Generating Random Numbers from Distributions See Reading Assignment OR441-DrKhalid Nowibet 1 Review 1 Inverse Transform Generate a number u i between 0 and 1 (one U-axis) and then find the corresponding
More informationGenerating Random Variables
Generating Random Variables These slides are created by Dr. Yih Huang of George Mason University. Students registered in Dr. Huang's courses at GMU can make a single machine-readable copy and print a single
More informationChapter 4: Monte Carlo Methods. Paisan Nakmahachalasint
Chapter 4: Monte Carlo Methods Paisan Nakmahachalasint Introduction Monte Carlo Methods are a class of computational algorithms that rely on repeated random sampling to compute their results. Monte Carlo
More informationSTAT 3610: Review of Probability Distributions
STAT 3610: Review of Probability Distributions Mark Carpenter Professor of Statistics Department of Mathematics and Statistics August 25, 2015 Support of a Random Variable Definition The support of a random
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 20, 2012 Introduction Introduction The world of the model-builder
More informationBayesian Inference and MCMC
Bayesian Inference and MCMC Aryan Arbabi Partly based on MCMC slides from CSC412 Fall 2018 1 / 18 Bayesian Inference - Motivation Consider we have a data set D = {x 1,..., x n }. E.g each x i can be the
More informationA quick introduction to Markov chains and Markov chain Monte Carlo (revised version)
A quick introduction to Markov chains and Markov chain Monte Carlo (revised version) Rasmus Waagepetersen Institute of Mathematical Sciences Aalborg University 1 Introduction These notes are intended to
More informationComputer Vision Group Prof. Daniel Cremers. 14. Sampling Methods
Prof. Daniel Cremers 14. Sampling Methods Sampling Methods Sampling Methods are widely used in Computer Science as an approximation of a deterministic algorithm to represent uncertainty without a parametric
More informationMonte Carlo Techniques
Physics 75.502 Part III: Monte Carlo Methods 40 Monte Carlo Techniques Monte Carlo refers to any procedure that makes use of random numbers. Monte Carlo methods are used in: Simulation of natural phenomena
More informationLecture 6: Markov Chain Monte Carlo
Lecture 6: Markov Chain Monte Carlo D. Jason Koskinen koskinen@nbi.ku.dk Photo by Howard Jackman University of Copenhagen Advanced Methods in Applied Statistics Feb - Apr 2016 Niels Bohr Institute 2 Outline
More informationIn manycomputationaleconomicapplications, one must compute thede nite n
Chapter 6 Numerical Integration In manycomputationaleconomicapplications, one must compute thede nite n integral of a real-valued function f de ned on some interval I of
More informationDifferentiation and Integration
Differentiation and Integration (Lectures on Numerical Analysis for Economists II) Jesús Fernández-Villaverde 1 and Pablo Guerrón 2 February 12, 2018 1 University of Pennsylvania 2 Boston College Motivation
More informationISyE 3044 Fall 2017 Test #1a Solutions
1 NAME ISyE 344 Fall 217 Test #1a Solutions This test is 75 minutes. You re allowed one cheat sheet. Good luck! 1. Suppose X has p.d.f. f(x) = 4x 3, < x < 1. Find E[ 2 X 2 3]. Solution: By LOTUS, we have
More informationProbability and Distributions
Probability and Distributions What is a statistical model? A statistical model is a set of assumptions by which the hypothetical population distribution of data is inferred. It is typically postulated
More informationContinuous Random Variables. and Probability Distributions. Continuous Random Variables and Probability Distributions ( ) ( ) Chapter 4 4.
UCLA STAT 11 A Applied Probability & Statistics for Engineers Instructor: Ivo Dinov, Asst. Prof. In Statistics and Neurology Teaching Assistant: Christopher Barr University of California, Los Angeles,
More informationRandom variables. DS GA 1002 Probability and Statistics for Data Science.
Random variables DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Motivation Random variables model numerical quantities
More information9.07 Introduction to Probability and Statistics for Brain and Cognitive Sciences Emery N. Brown
9.07 Introduction to Probability and Statistics for Brain and Cognitive Sciences Emery N. Brown I. Objectives Lecture 5: Conditional Distributions and Functions of Jointly Distributed Random Variables
More informationStatistics and data analyses
Statistics and data analyses Designing experiments Measuring time Instrumental quality Precision Standard deviation depends on Number of measurements Detection quality Systematics and methology σ tot =
More informationNonlife Actuarial Models. Chapter 14 Basic Monte Carlo Methods
Nonlife Actuarial Models Chapter 14 Basic Monte Carlo Methods Learning Objectives 1. Generation of uniform random numbers, mixed congruential method 2. Low discrepancy sequence 3. Inversion transformation
More informationMathematical Statistics 1 Math A 6330
Mathematical Statistics 1 Math A 6330 Chapter 2 Transformations and Expectations Mohamed I. Riffi Department of Mathematics Islamic University of Gaza September 14, 2015 Outline 1 Distributions of Functions
More informationStat 451 Lecture Notes Simulating Random Variables
Stat 451 Lecture Notes 05 12 Simulating Random Variables Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapter 6 in Givens & Hoeting, Chapter 22 in Lange, and Chapter 2 in Robert & Casella 2 Updated:
More informationHow does the computer generate observations from various distributions specified after input analysis?
1 How does the computer generate observations from various distributions specified after input analysis? There are two main components to the generation of observations from probability distributions.
More informationChing-Han Hsu, BMES, National Tsing Hua University c 2015 by Ching-Han Hsu, Ph.D., BMIR Lab. = a + b 2. b a. x a b a = 12
Lecture 5 Continuous Random Variables BMIR Lecture Series in Probability and Statistics Ching-Han Hsu, BMES, National Tsing Hua University c 215 by Ching-Han Hsu, Ph.D., BMIR Lab 5.1 1 Uniform Distribution
More informationTransformations and Expectations
Transformations and Expectations 1 Distributions of Functions of a Random Variable If is a random variable with cdf F (x), then any function of, say g(), is also a random variable. Sine Y = g() is a function
More informationContinuous random variables
Continuous random variables CE 311S What was the difference between discrete and continuous random variables? The possible outcomes of a discrete random variable (finite or infinite) can be listed out;
More informationWhy the Summation Test Results in a Benford, and not a Uniform Distribution, for Data that Conforms to a Log Normal Distribution. R C Hall, MSEE, BSEE
Why the Summation Test Results in a Benford, and not a Uniform Distribution, for Data that Conforms to a Log Normal Distribution R C Hall, MSEE, BSEE e-mail: rhall2448@aol.com Abstract The Summation test
More informationStatistics, Data Analysis, and Simulation SS 2013
Statistics, Data Analysis, and Simulation SS 213 8.128.73 Statistik, Datenanalyse und Simulation Dr. Michael O. Distler Mainz, 23. April 213 What we ve learned so far Fundamental
More informationFrom Random Numbers to Monte Carlo. Random Numbers, Random Walks, Diffusion, Monte Carlo integration, and all that
From Random Numbers to Monte Carlo Random Numbers, Random Walks, Diffusion, Monte Carlo integration, and all that Random Walk Through Life Random Walk Through Life If you flip the coin 5 times you will
More informationThe random variable 1
The random variable 1 Contents 1. Definition 2. Distribution and density function 3. Specific random variables 4. Functions of one random variable 5. Mean and variance 2 The random variable A random variable
More informationMonte Carlo Integration. Computer Graphics CMU /15-662, Fall 2016
Monte Carlo Integration Computer Graphics CMU 15-462/15-662, Fall 2016 Talk Announcement Jovan Popovic, Senior Principal Scientist at Adobe Research will be giving a seminar on Character Animator -- Monday
More informationF denotes cumulative density. denotes probability density function; (.)
BAYESIAN ANALYSIS: FOREWORDS Notation. System means the real thing and a model is an assumed mathematical form for the system.. he probability model class M contains the set of the all admissible models
More informationMCMC: Markov Chain Monte Carlo
I529: Machine Learning in Bioinformatics (Spring 2013) MCMC: Markov Chain Monte Carlo Yuzhen Ye School of Informatics and Computing Indiana University, Bloomington Spring 2013 Contents Review of Markov
More informationExperimental Design and Statistics - AGA47A
Experimental Design and Statistics - AGA47A Czech University of Life Sciences in Prague Department of Genetics and Breeding Fall/Winter 2014/2015 Matúš Maciak (@ A 211) Office Hours: M 14:00 15:30 W 15:30
More informationMachine Learning using Bayesian Approaches
Machine Learning using Bayesian Approaches Sargur N. Srihari University at Buffalo, State University of New York 1 Outline 1. Progress in ML and PR 2. Fully Bayesian Approach 1. Probability theory Bayes
More informationHANDBOOK OF APPLICABLE MATHEMATICS
HANDBOOK OF APPLICABLE MATHEMATICS Chief Editor: Walter Ledermann Volume II: Probability Emlyn Lloyd University oflancaster A Wiley-Interscience Publication JOHN WILEY & SONS Chichester - New York - Brisbane
More informationGOV 2001/ 1002/ Stat E-200 Section 1 Probability Review
GOV 2001/ 1002/ Stat E-200 Section 1 Probability Review Solé Prillaman Harvard University January 28, 2015 1 / 54 LOGISTICS Course Website: j.mp/g2001 lecture notes, videos, announcements Canvas: problem
More informationLecture Notes 1 Probability and Random Variables. Conditional Probability and Independence. Functions of a Random Variable
Lecture Notes 1 Probability and Random Variables Probability Spaces Conditional Probability and Independence Random Variables Functions of a Random Variable Generation of a Random Variable Jointly Distributed
More informationDistributions of Functions of Random Variables. 5.1 Functions of One Random Variable
Distributions of Functions of Random Variables 5.1 Functions of One Random Variable 5.2 Transformations of Two Random Variables 5.3 Several Random Variables 5.4 The Moment-Generating Function Technique
More informationStat 516, Homework 1
Stat 516, Homework 1 Due date: October 7 1. Consider an urn with n distinct balls numbered 1,..., n. We sample balls from the urn with replacement. Let N be the number of draws until we encounter a ball
More informationMonte Carlo Methods in High Energy Physics I
Helmholtz International Workshop -- CALC 2009, July 10--20, Dubna Monte Carlo Methods in High Energy Physics CALC2009 - July 20 10, Dubna 2 Contents 3 Introduction Simple definition: A Monte Carlo technique
More information