Linear Modeling with Bayesian Statistics
|
|
- Archibald Todd
- 5 years ago
- Views:
Transcription
1 Linear Modeling with Bayesian Statistics Bayesian Approach I I I I I Estimate probability of a parameter State degree of believe in specific parameter values Evaluate probability of hypothesis given the data Incorporate prior knowledge Recognizes that Data is one realization of some parameter distribution
2 Today Markov Chain Monte-Carlo Sampling Applying MCMC & Bayesian Inference to Multiparameter Models Evaluating convergence and fit of MCMC fit models Modifying our Prior Bayes Theorem Expanded p(θ X) = p(θ X) = p(x θ)p (θ) j i=0 p(x θ i )p(θ i ) - Algebraically Solvable p(x θ)p (θ) p(x θ)p(θ)dθ - Analytically Solveable for Conjugate Priors p(x θ)p (θ η)p(η)dη p(θ X) = p(x θ)p(θ)dθdη - Hierarchical Model: need numerical integration approach with random hyperparameters
3 Trace of (Intercept) Density of (Intercept) Iterations N = 1000 Bandwidth = If we cannot analytically solve a distribution, we can still simulate from it: Chose a set of starting values X at t=0 Chose a random set of parameters, Y, based on X Calculate an acceptance ratio, α, based on P(Y)/P(X) If α 1 X(t+1) = Y Otherwise, select a uniorm random number between 0 and 1, U If U α, X(t+1) = Y. Otherwise, X(t+1) = X. Rinse and repeat (Note, this is the Metropolis-Hastings Algorithm - there are others)
4 This is a time series. To use it for inference to sample from the final stationary distribution: Discard a burn in set of samples Thin your chain to reduce temporal autocorrelation Examine chain for convergence on your posterior distribution Evaluate multiple chains to ensure convergence to a single distribution Many different samplers using different decision rules for f. We use the Gibbs Sampler commonly. Value Iteration
5 Value Iteration Value Iteration
6 Value Iteration Value Iteration
7 Value Iteration Value Iteration
8 Thinned Value Iteration Linear Modeling with MCMC
9 Software Options for MCMC WinBUGS OpenBUGS JAGS STAN MCMCglmm in R MCMCpack in R BUGS code for a Simple Linear Regression model { Prior alpha dnorm(0,0.001) beta dnorm(0,0.001) sigma dunif(0,100) Likelihood for (i in 1:n){ y[i] dnorm(mu[i],tau) mu[i] <- alpha + beta*x[i] } }
10 Back to the Wolves with MCMCglmm & coda pups inbreeding.coefficient MCMCglmm Powerful package that fits bayesian models using Metropolis-Hastings MCMC Available tools to tweak priors Flexible in error distribution Can accomodate random effects Built to deal with phylogenetic autocorrelation Can handle other autocorrelation (temporal, spatial) as well
11 MCMCglmm v. LM wolf_lm <- lm(pups inbreeding.coefficient, data=wolves) wolf_mcmc <- MCMCglmm(pups inbreeding.coefficient, data=wolves, verbose=false) Coefficients are Very Similar summary(wolf_mcmc)$solutions post.mean l-95% CI u-95% CI eff.samp (Intercept) inbreeding.coefficient pmcmc (Intercept) inbreeding.coefficient summary(wolf_lm)$coef Estimate Std. Error t value (Intercept) inbreeding.coefficient Pr(> t ) (Intercept) 3.143e-08 inbreeding.coefficient 1.633e-03
12 All of the assumptions of a linear model with a normal error distribution still apply! Make a QQ-Plot! wolf_resid <- wolves$pups - predict(wolf_mcmc) qqnorm(wolf_resid) qqline(wolf_resid) Normal Q Q Plot Sample Quantiles
13 Check Normality of Residuals Chisq - problems with small values q <- qqnorm(wolf_resid) Normal Q Q Plot Sample Quantiles Theoretical Quantiles chisq.test(q$y, q$x) Fitted v. Residuals! Pearson's Chi-squared test plot(fitted(wolf_mcmc), wolf_resid) abline(a=0, b=0, lwd=1, lty=2) data: q$y and q$x X-squared = 432, df = 414, p-value = KolmogorovâĂŞSmirnov Problems with Ties ks.test(wolf_resid) Error: argument "y" is missing, with no default wolf_resid Shapiro-Wilks Test Very Sensitive shapiro.test(wolf_resid) Shapiro-Wilk normality test data: wolf_resid fitted(wolf_mcmc)
14 MCMC has its own assumptions! No Inferences Without Convergence!!! coda package Analysis of MCMC chains Variety of tools for convergence Works with single or multiple chains
15 Does this look Converged? plot(wolf_mcmc, trace=true, smooth=t) Trace of (Intercept) Density of (Intercept) Iterations N = 1000 Bandwidth = Trace of inbreeding.coefficient Density of inbreeding.coefficient Iterations N = 1000 Bandwidth = Trace of units Density of units Multiple Chains Every MCMC run (chain) starts with different start values Each run SHOULD converge on the same answer Chains should overlap make a list of MCMC chains using Iterations Sol which contains parameters N = 1000 Bandwidth = wolfchains2 <- as.mcmc.list(lapply(1:2, function(i) MCMCglmm(pups inbreeding.coefficient, data=wolves, verbose=false)$sol ))
16 Multiple Chains plot(wolfchains2) Trace of (Intercept) Density of (Intercept) Iterations N = 1000 Bandwidth = Trace of inbreeding.coefficient Density of inbreeding.coefficient Iterations N = 1000 Bandwidth = Multiple Chains Trace of (Intercept) Density of (Intercept) Iterations N = 1000 Bandwidth = Trace of inbreeding.coefficient Density of inbreeding.coefficient Iterations N = 1000 Bandwidth = 0.766
17 New Information from Many Chains Iterations = 3001:12991 Thinning interval = 10 Number of chains = 2 Sample size per chain = Empirical mean and standard deviation for each variable, plus standard error of the mean: Mean SD Naive SE Time-series SE (Intercept) inbreeding.coefficient Quantiles for each variable: 2.5% 25% 50% 75% 97.5% (Intercept) inbreeding.coefficient But am I REALLY Converging? Gelman-Rubin Diagnostic Averaging across chains should shrink estimates to a mean How much shrinkage is happening after 10, 20, 30, 40,... iterations? If you have converged, the scale reduction should be 1 - i.e., none. Other diagnostics, too - but this is fairly standard
18 Gelman-Rubin Plots gelman.plot(wolfchains4) (Intercept) inbreeding.coefficient shrink factor median 97.5% shrink factor median 97.5% last iteration in chain last iteration in chain Gelman-Rubin By the Numbers gelman.diag(wolfchains4) Potential scale reduction factors: Point est. Upper C.I. (Intercept) inbreeding.coefficient Multivariate psrf 1
19 So what about that Prior thing? Priors The Default Prior for an Intercept-Slope Model: MV Norm ( µ = ( ) 0, V = 0 Intercept = row 1, Slope = Row2 ( ) 1e10 0 ) 0 1e10
20 What if we have a prior? MV Norm ( µ = ( ) 0, V = 10 ( ) 1e10 0 ) 0 1 strongprior<-list(b=list(mu=c(0,-5), V=diag(c(1e10, 1)))) Put a Prior On It wolf_mcmc_prior <- MCMCglmm(pups inbreeding.coefficient, data=wolves, verbose=false, prior=strongprior)
21 Influence of a Prior on this Small Dataset summary(wolf_mcmc)$solutions post.mean l-95% CI u-95% CI eff.samp (Intercept) inbreeding.coefficient pmcmc (Intercept) inbreeding.coefficient summary(wolf_mcmc_prior)$solutions post.mean l-95% CI u-95% CI eff.samp (Intercept) inbreeding.coefficient pmcmc (Intercept) inbreeding.coefficient 0.001
Markov Chain Monte Carlo (part 1)
Markov Chain Monte Carlo (part 1) Edps 590BAY Carolyn J. Anderson Department of Educational Psychology c Board of Trustees, University of Illinois Spring 2018 Depending on the book that you select for
More informationMCMC Diagnostics. Yingbo Li MATH Clemson University. Yingbo Li (Clemson) MCMC Diagnostics MATH / 24
MCMC Diagnostics Yingbo Li Clemson University MATH 9810 Yingbo Li (Clemson) MCMC Diagnostics MATH 9810 1 / 24 Convergence to Posterior Distribution Theory proves that if a Gibbs sampler iterates enough,
More informationBayesian Modelling with JAGS and R
Bayesian Modelling with JAGS and R Martyn Plummer International Agency for Research on Cancer Rencontres R, 3 July 2012 CRAN Task View Bayesian Inference The CRAN Task View Bayesian Inference is maintained
More informationIssues in MCMC use for Bayesian model fitting. Practical Considerations for WinBUGS Users
Practical Considerations for WinBUGS Users Kate Cowles, Ph.D. Department of Statistics and Actuarial Science University of Iowa 22S:138 Lecture 12 Oct. 3, 2003 Issues in MCMC use for Bayesian model fitting
More informationA Basic Example of ANOVA in JAGS Joel S Steele
A Basic Example of ANOVA in JAGS Joel S Steele The purpose This demonstration is intended to show how a simple one-way ANOVA can be coded and run in the JAGS framework. This is by no means an exhaustive
More informationwinbugs and openbugs
Eric F. Lock UMN Division of Biostatistics, SPH elock@umn.edu 04/19/2017 Bayesian estimation software Several stand-alone applications and add-ons to estimate Bayesian models Stand-alone applications:
More informationBAYESIAN OUTPUT ANALYSIS PROGRAM (BOA) VERSION 1.0 USER S MANUAL
BAYESIAN OUTPUT ANALYSIS PROGRAM (BOA) VERSION 1.0 USER S MANUAL Brian J. Smith January 8, 2003 Contents 1 Getting Started 4 1.1 Hardware/Software Requirements.................... 4 1.2 Obtaining BOA..............................
More informationST440/540: Applied Bayesian Analysis. (5) Multi-parameter models - Initial values and convergence diagn
(5) Multi-parameter models - Initial values and convergence diagnostics Tuning the MCMC algoritm MCMC is beautiful because it can handle virtually any statistical model and it is usually pretty easy to
More informationJournal of Statistical Software
JSS Journal of Statistical Software December 2007, Volume 23, Issue 9. http://www.jstatsoft.org/ WinBUGSio: A SAS Macro for the Remote Execution of WinBUGS Michael K. Smith Pfizer Global Research and Development
More informationIntroduction to WinBUGS
Introduction to WinBUGS Joon Jin Song Department of Statistics Texas A&M University Introduction: BUGS BUGS: Bayesian inference Using Gibbs Sampling Bayesian Analysis of Complex Statistical Models using
More informationBayesian Statistics Group 8th March Slice samplers. (A very brief introduction) The basic idea
Bayesian Statistics Group 8th March 2000 Slice samplers (A very brief introduction) The basic idea lacements To sample from a distribution, simply sample uniformly from the region under the density function
More informationBayesian data analysis using R
Bayesian data analysis using R BAYESIAN DATA ANALYSIS USING R Jouni Kerman, Samantha Cook, and Andrew Gelman Introduction Bayesian data analysis includes but is not limited to Bayesian inference (Gelman
More informationBUGS: Language, engines, and interfaces
BUGS: Language, engines, and interfaces Patrick Breheny January 17 Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/18 The BUGS framework The BUGS project (Bayesian inference using Gibbs Sampling)
More informationBayesian Computation with JAGS
JAGS is Just Another Gibbs Sampler Cross-platform Accessible from within R Bayesian Computation with JAGS What I did Downloaded and installed JAGS. In the R package installer, downloaded rjags and dependencies.
More information1 Methods for Posterior Simulation
1 Methods for Posterior Simulation Let p(θ y) be the posterior. simulation. Koop presents four methods for (posterior) 1. Monte Carlo integration: draw from p(θ y). 2. Gibbs sampler: sequentially drawing
More informationJAGS for Rats. Jonathan Rougier School of Mathematics University of Bristol UK. Version 2.0, compiled February 23, 2017
JAGS for Rats Jonathan Rougier School of Mathematics University of Bristol UK Version 2.0, compiled February 23, 2017 1 Introduction This is a worksheet about JAGS, using rjags. These must both be installed
More informationMCMC Methods for data modeling
MCMC Methods for data modeling Kenneth Scerri Department of Automatic Control and Systems Engineering Introduction 1. Symposium on Data Modelling 2. Outline: a. Definition and uses of MCMC b. MCMC algorithms
More informationParameterization Issues and Diagnostics in MCMC
Parameterization Issues and Diagnostics in MCMC Gill Chapter 10 & 12 November 10, 2008 Convergence to Posterior Distribution Theory tells us that if we run the Gibbs sampler long enough the samples we
More informationAn Introduction to Using WinBUGS for Cost-Effectiveness Analyses in Health Economics
Practical 1: Getting started in OpenBUGS Slide 1 An Introduction to Using WinBUGS for Cost-Effectiveness Analyses in Health Economics Dr. Christian Asseburg Centre for Health Economics Practical 1 Getting
More informationClustering Relational Data using the Infinite Relational Model
Clustering Relational Data using the Infinite Relational Model Ana Daglis Supervised by: Matthew Ludkin September 4, 2015 Ana Daglis Clustering Data using the Infinite Relational Model September 4, 2015
More informationOverview. Monte Carlo Methods. Statistics & Bayesian Inference Lecture 3. Situation At End Of Last Week
Statistics & Bayesian Inference Lecture 3 Joe Zuntz Overview Overview & Motivation Metropolis Hastings Monte Carlo Methods Importance sampling Direct sampling Gibbs sampling Monte-Carlo Markov Chains Emcee
More informationFrom Bayesian Analysis of Item Response Theory Models Using SAS. Full book available for purchase here.
From Bayesian Analysis of Item Response Theory Models Using SAS. Full book available for purchase here. Contents About this Book...ix About the Authors... xiii Acknowledgments... xv Chapter 1: Item Response
More informationR package mcll for Monte Carlo local likelihood estimation
R package mcll for Monte Carlo local likelihood estimation Minjeong Jeon University of California, Berkeley Sophia Rabe-Hesketh University of California, Berkeley February 4, 2013 Cari Kaufman University
More information1. Start WinBUGS by double clicking on the WinBUGS icon (or double click on the file WinBUGS14.exe in the WinBUGS14 directory in C:\Program Files).
Hints on using WinBUGS 1 Running a model in WinBUGS 1. Start WinBUGS by double clicking on the WinBUGS icon (or double click on the file WinBUGS14.exe in the WinBUGS14 directory in C:\Program Files). 2.
More informationA GENERAL GIBBS SAMPLING ALGORITHM FOR ANALYZING LINEAR MODELS USING THE SAS SYSTEM
A GENERAL GIBBS SAMPLING ALGORITHM FOR ANALYZING LINEAR MODELS USING THE SAS SYSTEM Jayawant Mandrekar, Daniel J. Sargent, Paul J. Novotny, Jeff A. Sloan Mayo Clinic, Rochester, MN 55905 ABSTRACT A general
More informationISyE8843A, Brani Vidakovic Handout 14
ISyE8843A, Brani Vidakovic Handout 4 BUGS BUGS is freely available software for constructing Bayesian statistical models and evaluating them using MCMC methodology. BUGS and WINBUGS are distributed freely
More informationMarkov chain Monte Carlo methods
Markov chain Monte Carlo methods (supplementary material) see also the applet http://www.lbreyer.com/classic.html February 9 6 Independent Hastings Metropolis Sampler Outline Independent Hastings Metropolis
More informationThe boa Package. May 3, 2005
The boa Package May 3, 2005 Version 1.1.5-2 Date 2005-05-02 Title Bayesian Output Analysis Program (BOA) for MCMC Author Maintainer Depends R (>= 1.7) A menu-driven program and
More informationStatistical Matching using Fractional Imputation
Statistical Matching using Fractional Imputation Jae-Kwang Kim 1 Iowa State University 1 Joint work with Emily Berg and Taesung Park 1 Introduction 2 Classical Approaches 3 Proposed method 4 Application:
More informationBayesian Estimation for Skew Normal Distributions Using Data Augmentation
The Korean Communications in Statistics Vol. 12 No. 2, 2005 pp. 323-333 Bayesian Estimation for Skew Normal Distributions Using Data Augmentation Hea-Jung Kim 1) Abstract In this paper, we develop a MCMC
More informationPackage beast. March 16, 2018
Type Package Package beast March 16, 2018 Title Bayesian Estimation of Change-Points in the Slope of Multivariate Time-Series Version 1.1 Date 2018-03-16 Author Maintainer Assume that
More informationSemiparametric Mixed Effecs with Hierarchical DP Mixture
Semiparametric Mixed Effecs with Hierarchical DP Mixture R topics documented: April 21, 2007 hdpm-package........................................ 1 hdpm............................................ 2 hdpmfitsetup........................................
More informationAn Ecological Modeler s Primer on JAGS
1 An Ecological Modeler s Primer on JAGS 2 3 N. Thompson Hobbs May 19, 2014 4 5 6 Natural Resource Ecology Laboratory, Department of Ecosystem Science and Sustainability, and Graduate Degree Program in
More informationTutorial using BEAST v2.4.1 Troubleshooting David A. Rasmussen
Tutorial using BEAST v2.4.1 Troubleshooting David A. Rasmussen 1 Background The primary goal of most phylogenetic analyses in BEAST is to infer the posterior distribution of trees and associated model
More informationNested Sampling: Introduction and Implementation
UNIVERSITY OF TEXAS AT SAN ANTONIO Nested Sampling: Introduction and Implementation Liang Jing May 2009 1 1 ABSTRACT Nested Sampling is a new technique to calculate the evidence, Z = P(D M) = p(d θ, M)p(θ
More informationAMCMC: An R interface for adaptive MCMC
AMCMC: An R interface for adaptive MCMC by Jeffrey S. Rosenthal * (February 2007) Abstract. We describe AMCMC, a software package for running adaptive MCMC algorithms on user-supplied density functions.
More informationIntroduction to Applied Bayesian Modeling A brief JAGS and R2jags tutorial
Introduction to Applied Bayesian Modeling A brief JAGS and R2jags tutorial Johannes Karreth University of Georgia jkarreth@uga.edu ICPSR Summer Program 2011 Last updated on July 7, 2011 1 What are JAGS,
More informationPackage bsam. July 1, 2017
Type Package Package bsam July 1, 2017 Title Bayesian State-Space Models for Animal Movement Version 1.1.2 Depends R (>= 3.3.0), rjags (>= 4-6) Imports coda (>= 0.18-1), dplyr (>= 0.5.0), ggplot2 (>= 2.1.0),
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Lecture 17 EM CS/CNS/EE 155 Andreas Krause Announcements Project poster session on Thursday Dec 3, 4-6pm in Annenberg 2 nd floor atrium! Easels, poster boards and cookies
More informationBayesian Analysis of Extended Lomax Distribution
Bayesian Analysis of Extended Lomax Distribution Shankar Kumar Shrestha and Vijay Kumar 2 Public Youth Campus, Tribhuvan University, Nepal 2 Department of Mathematics and Statistics DDU Gorakhpur University,
More informationYEVHEN YANKOVSKYY. Specialist, Ukrainian Marine State Technical University, 1999 M.A., Kyiv-Mohyla Academy, 2003 M.S., Iowa State University, 2006
APPLICATION OF A GIBBS SAMPLER TO ESTIMATING PARAMETERS OF A HIERARCHICAL NORMAL MODEL WITH A TIME TREND AND TESTING FOR EXISTENCE OF THE GLOBAL WARMING by YEVHEN YANKOVSKYY Specialist, Ukrainian Marine
More informationUsing the package glmbfp: a binary regression example.
Using the package glmbfp: a binary regression example. Daniel Sabanés Bové 3rd August 2017 This short vignette shall introduce into the usage of the package glmbfp. For more information on the methodology,
More informationCHAPTER 1 INTRODUCTION
Introduction CHAPTER 1 INTRODUCTION Mplus is a statistical modeling program that provides researchers with a flexible tool to analyze their data. Mplus offers researchers a wide choice of models, estimators,
More informationA quick introduction to First Bayes
A quick introduction to First Bayes Lawrence Joseph October 1, 2003 1 Introduction This document very briefly reviews the main features of the First Bayes statistical teaching package. For full details,
More informationCS281 Section 9: Graph Models and Practical MCMC
CS281 Section 9: Graph Models and Practical MCMC Scott Linderman November 11, 213 Now that we have a few MCMC inference algorithms in our toolbox, let s try them out on some random graph models. Graphs
More informationFitting Bayesian item response models in Stata and Stan
The Stata Journal (2017) 17, Number 2, pp. 343 357 Fitting Bayesian item response models in Stata and Stan Robert L. Grant BayesCamp Croydon, UK robert@bayescamp.com Daniel C. Furr University of California
More informationChapter 1. Introduction
Chapter 1 Introduction A Monte Carlo method is a compuational method that uses random numbers to compute (estimate) some quantity of interest. Very often the quantity we want to compute is the mean of
More informationBART STAT8810, Fall 2017
BART STAT8810, Fall 2017 M.T. Pratola November 1, 2017 Today BART: Bayesian Additive Regression Trees BART: Bayesian Additive Regression Trees Additive model generalizes the single-tree regression model:
More informationRunning WinBUGS from within R
Applied Bayesian Inference A Running WinBUGS from within R, KIT WS 2010/11 1 Running WinBUGS from within R 1 Batch Mode Although WinBUGS includes methods to analyze the output of the Markov chains like
More informationPackage mfa. R topics documented: July 11, 2018
Package mfa July 11, 2018 Title Bayesian hierarchical mixture of factor analyzers for modelling genomic bifurcations Version 1.2.0 MFA models genomic bifurcations using a Bayesian hierarchical mixture
More informationIntroduction to Bayesian Analysis in Stata
tools Introduction to Bayesian Analysis in Gustavo Sánchez Corp LLC September 15, 2017 Porto, Portugal tools 1 Bayesian analysis: 2 Basic Concepts The tools 14: The command 15: The bayes prefix Postestimation
More informationAnalysis of Incomplete Multivariate Data
Analysis of Incomplete Multivariate Data J. L. Schafer Department of Statistics The Pennsylvania State University USA CHAPMAN & HALL/CRC A CR.C Press Company Boca Raton London New York Washington, D.C.
More informationRegression Analysis and Linear Regression Models
Regression Analysis and Linear Regression Models University of Trento - FBK 2 March, 2015 (UNITN-FBK) Regression Analysis and Linear Regression Models 2 March, 2015 1 / 33 Relationship between numerical
More informationMissing Data Analysis for the Employee Dataset
Missing Data Analysis for the Employee Dataset 67% of the observations have missing values! Modeling Setup Random Variables: Y i =(Y i1,...,y ip ) 0 =(Y i,obs, Y i,miss ) 0 R i =(R i1,...,r ip ) 0 ( 1
More informationCase Study IV: Bayesian clustering of Alzheimer patients
Case Study IV: Bayesian clustering of Alzheimer patients Mike Wiper and Conchi Ausín Department of Statistics Universidad Carlos III de Madrid Advanced Statistics and Data Mining Summer School 2nd - 6th
More informationCorrelation and Regression
How are X and Y Related Correlation and Regression Causation (regression) p(y X=) Correlation p(y=, X=) http://kcd.com/552/ Correlation Can be Induced b Man Mechanisms # of Pups 1 2 3 4 5 6 7 8 Inbreeding
More informationApproximate Bayesian Computation. Alireza Shafaei - April 2016
Approximate Bayesian Computation Alireza Shafaei - April 2016 The Problem Given a dataset, we are interested in. The Problem Given a dataset, we are interested in. The Problem Given a dataset, we are interested
More informationbrms: An R Package for Bayesian Multilevel Models using Stan
brms: An R Package for Bayesian Multilevel Models using Stan Paul Bürkner Institut für Psychologie Westfälische Wilhelms-Universität Münster 26.02.2016 Agenda Agenda 1 Short introduction to Stan 2 The
More informationBayesian Spatiotemporal Modeling with Hierarchical Spatial Priors for fmri
Bayesian Spatiotemporal Modeling with Hierarchical Spatial Priors for fmri Galin L. Jones 1 School of Statistics University of Minnesota March 2015 1 Joint with Martin Bezener and John Hughes Experiment
More informationINTRO TO THE METROPOLIS ALGORITHM
INTRO TO THE METROPOLIS ALGORITHM A famous reliability experiment was performed where n = 23 ball bearings were tested and the number of revolutions were recorded. The observations in ballbearing2.dat
More informationPackage ggmcmc. August 29, 2016
Package ggmcmc August 29, 2016 Title Tools for Analyzing MCMC Simulations from Bayesian Inference Tools for assessing and diagnosing convergence of Markov Chain Monte Carlo simulations, as well as for
More informationStatistical techniques for data analysis in Cosmology
Statistical techniques for data analysis in Cosmology arxiv:0712.3028; arxiv:0911.3105 Numerical recipes (the bible ) Licia Verde ICREA & ICC UB-IEEC http://icc.ub.edu/~liciaverde outline Lecture 1: Introduction
More informationA noninformative Bayesian approach to small area estimation
A noninformative Bayesian approach to small area estimation Glen Meeden School of Statistics University of Minnesota Minneapolis, MN 55455 glen@stat.umn.edu September 2001 Revised May 2002 Research supported
More informationA Short History of Markov Chain Monte Carlo
A Short History of Markov Chain Monte Carlo Christian Robert and George Casella 2010 Introduction Lack of computing machinery, or background on Markov chains, or hesitation to trust in the practicality
More informationWinBUGS Tutorial Outline
WinBUGS Tutorial Outline August 4, 2005 Authors: Ed Merkle & Trisha Van Zandt WinBUGS overview Description of Example Data Example 1: Confidence in Signal Detection Example 2: Mean RT Model Example 3:
More informationVariability in Annual Temperature Profiles
Variability in Annual Temperature Profiles A Multivariate Spatial Analysis of Regional Climate Model Output Tamara Greasby, Stephan Sain Institute for Mathematics Applied to Geosciences, National Center
More informationPackage atmcmc. February 19, 2015
Type Package Package atmcmc February 19, 2015 Title Automatically Tuned Markov Chain Monte Carlo Version 1.0 Date 2014-09-16 Author Jinyoung Yang Maintainer Jinyoung Yang
More informationQuantitative Biology II!
Quantitative Biology II! Lecture 3: Markov Chain Monte Carlo! March 9, 2015! 2! Plan for Today!! Introduction to Sampling!! Introduction to MCMC!! Metropolis Algorithm!! Metropolis-Hastings Algorithm!!
More informationGraphical Models, Bayesian Method, Sampling, and Variational Inference
Graphical Models, Bayesian Method, Sampling, and Variational Inference With Application in Function MRI Analysis and Other Imaging Problems Wei Liu Scientific Computing and Imaging Institute University
More informationRecap: The E-M algorithm. Biostatistics 615/815 Lecture 22: Gibbs Sampling. Recap - Local minimization methods
Recap: The E-M algorithm Biostatistics 615/815 Lecture 22: Gibbs Sampling Expectation step (E-step) Given the current estimates of parameters λ (t), calculate the conditional distribution of latent variable
More informationWinBUGS A Bayesian modelling framework: Concepts, structure, and extensibility
Statistics and Computing (2000) 10, 325 337 WinBUGS A Bayesian modelling framework: Concepts, structure, and extensibility DAVID J. LUNN, ANDREW THOMAS, NICKY BEST and DAVID SPIEGELHALTER Department of
More informationNONPARAMETRIC REGRESSION WIT MEASUREMENT ERROR: SOME RECENT PR David Ruppert Cornell University
NONPARAMETRIC REGRESSION WIT MEASUREMENT ERROR: SOME RECENT PR David Ruppert Cornell University www.orie.cornell.edu/ davidr (These transparencies, preprints, and references a link to Recent Talks and
More informationCanopy Light: Synthesizing multiple data sources
Canopy Light: Synthesizing multiple data sources Tree growth depends upon light (previous example, lab 7) Hard to measure how much light an ADULT tree receives Multiple sources of proxy data Exposed Canopy
More informationA Comparison of Two MCMC Algorithms for Hierarchical Mixture Models
A Comparison of Two MCMC Algorithms for Hierarchical Mixture Models Russell G. Almond Florida State University Abstract Mixture models form an important class of models for unsupervised learning, allowing
More informationTest Run to Check the Installation of JAGS & Rjags
John Miyamoto File = D:\bugs\test.jags.install.docm 1 Test Run to Check the Installation of JAGS & Rjags The following annotated code is extracted from John Kruschke's R scripts, "E:\btut\r\BernBetaBugsFull.R"
More informationMODELING AND OPTIMIZATION OF A ROBUST GAS SENSOR
MODELING AND OPTIMIZATION OF A ROBUST GAS SENSOR Margarita A. Rebolledo C., Sebastian Krey, Thomas Bartz-Beielstein, Oliver Flasch, Andreas Fischbach, Jörg Stork SPOTSeven Lab, TH Köln, Gummersbach, Germany
More informationSTATISTICS (STAT) Statistics (STAT) 1
Statistics (STAT) 1 STATISTICS (STAT) STAT 2013 Elementary Statistics (A) Prerequisites: MATH 1483 or MATH 1513, each with a grade of "C" or better; or an acceptable placement score (see placement.okstate.edu).
More informationMissing Data Analysis for the Employee Dataset
Missing Data Analysis for the Employee Dataset 67% of the observations have missing values! Modeling Setup For our analysis goals we would like to do: Y X N (X, 2 I) and then interpret the coefficients
More informationEstimation of Item Response Models
Estimation of Item Response Models Lecture #5 ICPSR Item Response Theory Workshop Lecture #5: 1of 39 The Big Picture of Estimation ESTIMATOR = Maximum Likelihood; Mplus Any questions? answers Lecture #5:
More informationJAGS Version user manual. Martyn Plummer
JAGS Version 4.3.0 user manual Martyn Plummer 28 June 2017 Contents 1 Introduction 4 1.1 Downloading JAGS................................. 4 1.2 Getting help..................................... 4 1.3
More informationTHE UNIVERSITY OF BRITISH COLUMBIA FORESTRY 430 and 533. Time: 50 minutes 40 Marks FRST Marks FRST 533 (extra questions)
THE UNIVERSITY OF BRITISH COLUMBIA FORESTRY 430 and 533 MIDTERM EXAMINATION: October 14, 2005 Instructor: Val LeMay Time: 50 minutes 40 Marks FRST 430 50 Marks FRST 533 (extra questions) This examination
More informationDemo 09-1a: One Sample (Single Group) Models
March 8, 2017 File = D:\P548\demo.09-1a.p548.w17.docm 1 John Miyamoto (email: jmiyamot@uw.edu) Psych 548: Bayesian Statistics, Modeling & Reasoning Winter 2017 Course website: http://faculty.washington.edu/jmiyamot/p548/p548-set.htm
More informationBayesian Robust Inference of Differential Gene Expression The bridge package
Bayesian Robust Inference of Differential Gene Expression The bridge package Raphael Gottardo October 30, 2017 Contents Department Statistics, University of Washington http://www.rglab.org raph@stat.washington.edu
More informationConvergence and Efficiency of Adaptive MCMC. Jinyoung Yang
Convergence and Efficiency of Adaptive MCMC by Jinyoung Yang A thesis submitted in conformity with the requirements for the degree of Doctor of Philosophy Graduate Department of Statistical Sciences University
More informationMCMC GGUM v1.2 User s Guide
MCMC GGUM v1.2 User s Guide Wei Wang University of Central Florida Jimmy de la Torre Rutgers, The State University of New Jersey Fritz Drasgow University of Illinois at Urbana-Champaign Travis Meade and
More informationBayesian parameter estimation in Ecolego using an adaptive Metropolis-Hastings-within-Gibbs algorithm
IT 16 052 Examensarbete 30 hp Juni 2016 Bayesian parameter estimation in Ecolego using an adaptive Metropolis-Hastings-within-Gibbs algorithm Sverrir Þorgeirsson Institutionen för informationsteknologi
More informationMonte Carlo for Spatial Models
Monte Carlo for Spatial Models Murali Haran Department of Statistics Penn State University Penn State Computational Science Lectures April 2007 Spatial Models Lots of scientific questions involve analyzing
More informationSlice sampler algorithm for generalized Pareto distribution
Slice sampler algorithm for generalized Pareto distribution Mohammad Rostami, Mohd Bakri Adam Yahya, Mohamed Hisham Yahya, Noor Akma Ibrahim Abstract In this paper, we developed the slice sampler algorithm
More informationADAPTIVE METROPOLIS-HASTINGS SAMPLING, OR MONTE CARLO KERNEL ESTIMATION
ADAPTIVE METROPOLIS-HASTINGS SAMPLING, OR MONTE CARLO KERNEL ESTIMATION CHRISTOPHER A. SIMS Abstract. A new algorithm for sampling from an arbitrary pdf. 1. Introduction Consider the standard problem of
More informationMissing Data and Imputation
Missing Data and Imputation NINA ORWITZ OCTOBER 30 TH, 2017 Outline Types of missing data Simple methods for dealing with missing data Single and multiple imputation R example Missing data is a complex
More informationApproaches to Missing Data
Approaches to Missing Data A Presentation by Russell Barbour, Ph.D. Center for Interdisciplinary Research on AIDS (CIRA) and Eugenia Buta, Ph.D. CIRA and The Yale Center of Analytical Studies (YCAS) April
More informationPackage jagsui. December 12, 2017
Version 1.4.9 Date 2017-12-08 Package jagsui December 12, 2017 Title A Wrapper Around 'rjags' to Streamline 'JAGS' Analyses Author Ken Kellner Maintainer Ken Kellner
More informationRJaCGH, a package for analysis of
RJaCGH, a package for analysis of CGH arrays with Reversible Jump MCMC 1. CGH Arrays: Biological problem: Changes in number of DNA copies are associated to cancer activity. Microarray technology: Oscar
More informationMSA101/MVE Lecture 5
MSA101/MVE187 2017 Lecture 5 Petter Mostad Chalmers University September 12, 2017 1 / 15 Importance sampling MC integration computes h(x)f (x) dx where f (x) is a probability density function, by simulating
More informationMarkov Random Fields and Gibbs Sampling for Image Denoising
Markov Random Fields and Gibbs Sampling for Image Denoising Chang Yue Electrical Engineering Stanford University changyue@stanfoed.edu Abstract This project applies Gibbs Sampling based on different Markov
More informationMixed Effects Models. Biljana Jonoska Stojkova Applied Statistics and Data Science Group (ASDa) Department of Statistics, UBC.
Mixed Effects Models Biljana Jonoska Stojkova Applied Statistics and Data Science Group (ASDa) Department of Statistics, UBC March 6, 2018 Resources for statistical assistance Department of Statistics
More informationPackage HKprocess. R topics documented: September 6, Type Package. Title Hurst-Kolmogorov process. Version
Package HKprocess September 6, 2014 Type Package Title Hurst-Kolmogorov process Version 0.0-1 Date 2014-09-06 Author Maintainer Imports MCMCpack (>= 1.3-3), gtools(>= 3.4.1) Depends
More informationImage analysis. Computer Vision and Classification Image Segmentation. 7 Image analysis
7 Computer Vision and Classification 413 / 458 Computer Vision and Classification The k-nearest-neighbor method The k-nearest-neighbor (knn) procedure has been used in data analysis and machine learning
More informationMissing Data and Imputation
Missing Data and Imputation Hoff Chapter 7, GH Chapter 25 April 21, 2017 Bednets and Malaria Y:presence or absence of parasites in a blood smear AGE: age of child BEDNET: bed net use (exposure) GREEN:greenness
More informationUP School of Statistics Student Council Education and Research
w UP School of Statistics Student Council Education and Research erho.weebly.com 0 erhomyhero@gmail.com f /erhoismyhero t @erhomyhero S133_HOA_001 Statistics 133 Bayesian Statistical Inference Use of R
More information