Bayesian Computation with JAGS
|
|
- Eustace Leonard
- 5 years ago
- Views:
Transcription
1 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. rm(list=ls()) library(rjags) Loading required package: coda Linked to JAGS Loaded modules: basemod,bugs # Try the coffee taste test. 60 chose the new blend. # Uniform prior, posterior is Beta(61,41) x = 60 # The model specification model_string <- "model{ + X ~ dbinom(theta,100) # p(x theta) + theta ~ dbeta(1,1) # pi(theta) + }" model <- jags.model(textconnection(model_string), + data = list(x=x)) Compiling model graph Resolving undeclared variables Allocating nodes Graph Size: 4 Initializing model update(model, 10000); # Burnin for samples ************************************************** 100% post <- jags.samples(model, + variable.names=c("theta"), + n.iter=20000) ************************************************** 100% summary(post) Length Class Mode theta mcarray numeric postvals = as.numeric(post[[1]]) hist(postvals,probability=t,xlab=expression(theta),xlim=c(0,1), main='posterior distribution') Page 1 of 15
2 Page 2 of 15
3 acf(postvals) # Auto-correlation function Page 3 of 15
4 a = 61; b = 41 a/(a+b) # Expected value of theta given X [1] mean(postvals) [1] a*b/((a+b)^2*(a+b+1)) # Variance of theta given X [1] var(postvals) [1] # It totally worked. # Coda samples come with diagnostics as "methods." # This seems to be more popular. samp <- coda.samples(model, + variable.names=c("theta"), + n.iter=500, progress.bar="none") summary(samp) Iterations = 30001:30500 Thinning interval = 1 Number of chains = 1 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 Quantiles for each variable: 2.5% 25% 50% 75% 97.5% plot(samp) Page 4 of 15
5 Maybe JAGS knows too much. Try normal data with a non-conjugate prior. y i Normal(mu,tau) # Tau is precision = 1/variance mu Normal(0,100) tau LogNormal(0,1) Quit R and re-start. First see what LogNormal(0,1) looks like. rm(list=ls()) tau = seq(from = 0,to=10,length=200) Density = dlnorm(tau, 0,1) plot(tau,density,type='l',xlab=expression(tau),main="log-normal density of the precision") Page 5 of 15
6 rm(list=ls()) set.seed(9999) y <- rnorm(n = 50, mean = 100, sd = 15) # True precision = 1/15^2 = c(mean(y),sqrt(var(y))) [1] library(rjags) Loading required package: coda Linked to JAGS Loaded modules: basemod,bugs # The model specification model_string <- "model{ + for(i in 1:length(y)) { + y[i] ~ dnorm(mu, tau) + } + mu ~ dnorm(0, ) # A diffuse prior, but still proper + tau ~ dlnorm(0, 1) + }" # Running the model model <- jags.model(textconnection(model_string), data = list(y = y), n.chains = 1, n.adapt= 10000) % # n.adapt: "Initial sampling phase during which the samplers adapt their # behaviour to maximize their efficiency" update(model, 10000); # Burnin for samples ************************************************** 100% normal_samples <- coda.samples(model, variable.names=c("mu", "tau"), n.iter=500) ************************************************** 100% # normal_samples is a list of matrices -- one item long. summary(normal_samples); plot(normal_samples) Iterations = 20001:20500 Thinning interval = 1 Number of chains = 1 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 mu 1.006e e e-02 tau 5.611e e e Quantiles for each variable: 2.5% 25% 50% 75% 97.5% mu e tau e Page 6 of 15
7 normal_samples[[1]][1:10,] # First 10 rows mu tau [1,] [2,] [3,] [4,] [5,] [6,] [7,] [8,] [9,] [10,] Page 7 of 15
8 musim = normal_samples[[1]][,1] # First column is.numeric(musim) [1] TRUE summary(musim) Iterations = 20001:20500 Thinning interval = 1 Number of chains = 1 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 Quantiles for each variable: 2.5% 25% 50% 75% 97.5% acf(musim) Page 8 of 15
9 tausim = normal_samples[[1]][,2] # Second column summary(tausim) Iterations = 20001:20500 Thinning interval = 1 Number of chains = 1 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 5.611e e e e Quantiles for each variable: 2.5% 25% 50% 75% 97.5% acf(tausim) Page 9 of 15
10 sigmasim = 1/sqrt(tausim) # Convert to standard deviations -- more familiar summary(sigmasim) Iterations = 20001:20500 Thinning interval = 1 Number of chains = 1 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 Quantiles for each variable: 2.5% 25% 50% 75% 97.5% acf(sigmasim) Page 10 of 15
11 # plot(musim,sigmasim) # Error plot(as.numeric(musim),as.numeric(sigmasim)) cor(musim,sigmasim) [1] Page 11 of 15
12 # Experiment and try to break something rm(list=ls()) set.seed(9999) y <- rnorm(n = 500, mean = 100, sd = 15) # Much bigger sample size this time c(mean(y),sqrt(var(y))) [1] library(rjags) # The model specification model_string <- "model{ + for(i in 1:length(Y)) { + Y[i] ~ dnorm(mu, tau) + } + mu ~ dexp(1) # Standard exponential (Just guessing at the notation) + tau ~ dunif(0, 1) # Uniform prior + }" # Running the model model <- jags.model(textconnection(model_string), data = list(y = y), n.chains = 1, n.adapt= 10000) % update(model, 10000); # Burnin for samples ************************************************** 100% normal_samples <- coda.samples(model, variable.names=c("mu", "tau"), n.iter=20000) ************************************************** 100% summary(normal_samples) Iterations = 20001:40000 Thinning interval = 1 Number of chains = 1 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 mu e tau e Quantiles for each variable: 2.5% 25% 50% 75% 97.5% mu 2.575e tau 5.202e plot(normal_samples) Page 12 of 15
13 musim = normal_samples[[1]][,1] # First column max(musim) [1] sort(musim)[19995:20000] [1] acf(musim) Page 13 of 15
14 Page 14 of 15
15 tausim = normal_samples[[1]][,2] # Second column sigmasim = 1/sqrt(tausim) # Convert to standard deviations -- more familiar summary(sigmasim) Iterations = 20001:40000 Thinning interval = 1 Number of chains = 1 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 Quantiles for each variable: 2.5% 25% 50% 75% 97.5% acf(sigmasim) Page 15 of 15
A 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 informationLinear Modeling with Bayesian Statistics
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
More informationrjags Introduction Parameters: How to determine the parameters of a statistical model
January 25, 207 File = E:\inet\p548\demo.04-2.rjags.intro.docm John Miyamoto (email: jmiyamot@uw.edu) Psych 548: Bayesian Statistics, Modeling & Reasoning Winter 207 Course website: http://faculty.washington.edu/jmiyamot/p548/p548-set.htm
More informationMarkov 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 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 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 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 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 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 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 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 informationTest Run to Check the Installation of OpenBUGS, JAGS, BRugs, R2OpenBUGS, Rjags & R2jags
John Miyamoto File = E:\bugs\test.bugs.install.docm 1 Test Run to Check the Installation of OpenBUGS, JAGS, BRugs, R2OpenBUGS, Rjags & R2jags The following annotated code is extracted from John Kruschke's
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 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 informationIllustration of Adaptive MCMC
Illustration of Adaptive MCMC Andrew O. Finley and Sudipto Banerjee September 5, 2014 1 Overview The adaptmetropgibb function in spbayes generates MCMC samples for a continuous random vector using an adaptive
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 informationPackage rjags. R topics documented: October 19, 2018
Version 4-8 Date 2018-10-19 Title Bayesian Graphical Models using MCMC Depends R (>= 2.14.0), coda (>= 0.13) SystemRequirements JAGS 4.x.y URL http://mcmc-jags.sourceforge.net Suggests tcltk Interface
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 informationPackage lira. October 27, 2017
Type Package Title LInear Regression in Astronomy Version 2.0.0 Date 2017-10-27 Author Mauro Sereno Package lira October 27, 2017 Maintainer Mauro Sereno Description Performs Bayesian
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 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 informationPackage BayesCR. September 11, 2017
Type Package Package BayesCR September 11, 2017 Title Bayesian Analysis of Censored Regression Models Under Scale Mixture of Skew Normal Distributions Version 2.1 Author Aldo M. Garay ,
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 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 informationOrange tree growth: A demonstration and evaluation of nonlinear mixed-effects models in R, ADMB, and BUGS
Orange tree growth: A demonstration and evaluation of nonlinear mixed-effects models in R, ADMB, and BUGS Arni Magnusson, Mark Maunder, and Ben Bolker February 11, 2013 Contents 1 Introduction 1 2 Data
More informationPackage visualize. April 28, 2017
Type Package Package visualize April 28, 2017 Title Graph Probability Distributions with User Supplied Parameters and Statistics Version 4.3.0 Date 2017-04-27 Depends R (>= 3.0.0) Graphs the pdf or pmf
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 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 informationPackage bacon. October 31, 2018
Type Package Package October 31, 2018 Title Controlling bias and inflation in association studies using the empirical null distribution Version 1.10.0 Author Maarten van Iterson [aut, cre], Erik van Zwet
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 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 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 informationBeyond the black box: Using, programming, and sharing hierarchical modeling algorithms such as MCMC and Sequential MC using NIMBLE
Beyond the black box: Using, programming, and sharing hierarchical modeling algorithms such as MCMC and Sequential MC using NIMBLE Christopher Paciorek UC Berkeley Statistics Perry de Valpine (PI) Daniel
More informationPackage bdpopt. March 30, 2016
Version 1.0-1 Date 2016-03-29 Title Optimisation of Bayesian Decision Problems Author Sebastian Jobjörnsson [aut, cre] Package bdpopt March 30, 2016 Maintainer Depends R (>= 3.0.2) Optimisation of the
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 informationPackage TBSSurvival. January 5, 2017
Version 1.3 Date 2017-01-05 Package TBSSurvival January 5, 2017 Title Survival Analysis using a Transform-Both-Sides Model Author Adriano Polpo , Cassio de Campos , D.
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 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 informationPackage sparsereg. R topics documented: March 10, Type Package
Type Package Package sparsereg March 10, 2016 Title Sparse Bayesian Models for Regression, Subgroup Analysis, and Panel Data Version 1.2 Date 2016-03-01 Author Marc Ratkovic and Dustin Tingley Maintainer
More informationApplied Bayesian Modeling Using JAGS and BUGS via R
Applied Bayesian Modeling Using JAGS and BUGS via R Johannes Karreth Ursinus College jkarreth@ursinus.edu ICPSR Summer Program 2017 All code used in this tutorial can also be found on my github page at
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 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 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 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 informationPackage SafeBayes. October 20, 2016
Type Package Package SafeBayes October 20, 2016 Title Generalized and Safe-Bayesian Ridge and Lasso Regression Version 1.1 Date 2016-10-17 Depends R (>= 3.1.2), stats Description Functions for Generalized
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 informationPackage BAMBI. R topics documented: August 28, 2017
Type Package Title Bivariate Angular Mixture Models Version 1.1.1 Date 2017-08-23 Author Saptarshi Chakraborty, Samuel W.K. Wong Package BAMBI August 28, 2017 Maintainer Saptarshi Chakraborty
More informationSimulated example: Estimating diet proportions form fatty acids and stable isotopes
Simulated example: Estimating diet proportions form fatty acids and stable isotopes Philipp Neubauer Dragonfly Science, PO Box 755, Wellington 6, New Zealand June, Preamble DISCLAIMER: This is an evolving
More informationLogistic Regression. (Dichotomous predicted variable) Tim Frasier
Logistic Regression (Dichotomous predicted variable) Tim Frasier Copyright Tim Frasier This work is licensed under the Creative Commons Attribution 4.0 International license. Click here for more information.
More informationPackage TBSSurvival. July 1, 2012
Package TBSSurvival July 1, 2012 Version 1.0 Date 2012-06-30 Title TBS Model R package Author Adriano Polpo , Cassio de Campos , D. Sinha , Stuart
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 information1 Pencil and Paper stuff
Spring 2008 - Stat C141/ Bioeng C141 - Statistics for Bioinformatics Course Website: http://www.stat.berkeley.edu/users/hhuang/141c-2008.html Section Website: http://www.stat.berkeley.edu/users/mgoldman
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 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 informationSimulating from the Polya posterior by Glen Meeden, March 06
1 Introduction Simulating from the Polya posterior by Glen Meeden, glen@stat.umn.edu March 06 The Polya posterior is an objective Bayesian approach to finite population sampling. In its simplest form it
More informationPractical Bayesian Analysis for Failure Time Data
Practical Bayesian Analysis for Failure Time Data Best Practice Authored by: Michael Harman 15 June 2017 The goal of the STAT COE is to assist in developing rigorous, defensible test strategies to more
More informationMCMC for Bayesian Inference Metropolis: Solutions
MCMC for Bayesian Inference Metropolis: Solutions Below are the solutions to these exercises on MCMC for Bayesian Inference Metropolis: Exercises. #### Run the folowing lines before doing the exercises:
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 informationChapter 3. Bootstrap. 3.1 Introduction. 3.2 The general idea
Chapter 3 Bootstrap 3.1 Introduction The estimation of parameters in probability distributions is a basic problem in statistics that one tends to encounter already during the very first course on the subject.
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 informationR Programming Basics - Useful Builtin Functions for Statistics
R Programming Basics - Useful Builtin Functions for Statistics Vectorized Arithmetic - most arthimetic operations in R work on vectors. Here are a few commonly used summary statistics. testvect = c(1,3,5,2,9,10,7,8,6)
More informationNIMBLE User Manual. NIMBLE Development Team. Version 0.1
NIMBLE User Manual NIMBLE Development Team Version 0.1 Contents 1 Welcome to NIMBLE 5 1.1 Why something new?............................... 5 1.2 What does NIMBLE do?............................. 6 1.3
More informationPackage bayescl. April 14, 2017
Package bayescl April 14, 2017 Version 0.0.1 Date 2017-04-10 Title Bayesian Inference on a GPU using OpenCL Author Rok Cesnovar, Erik Strumbelj Maintainer Rok Cesnovar Description
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 informationChapter 6: Simulation Using Spread-Sheets (Excel)
Chapter 6: Simulation Using Spread-Sheets (Excel) Refer to Reading Assignments 1 Simulation Using Spread-Sheets (Excel) OBJECTIVES To be able to Generate random numbers within a spreadsheet environment.
More informationPackage boral. R topics documented: May 16, 2018
Package boral May 16, 2018 Title Bayesian Ordination and Regression AnaLysis Version 1.6.1 Date 2018-06-30 Author Francis K.C. Hui , with contributions from Wade Blanchard
More informationPackage logitnorm. R topics documented: February 20, 2015
Title Functions for the al distribution. Version 0.8.29 Date 2012-08-24 Author Package February 20, 2015 Maintainer Density, distribution, quantile and random generation function
More informationPackage gibbs.met documentation
Version 1.1-2 Package gibbs.met documentation Title Naive Gibbs Sampling with Metropolis Steps Author Longhai Li of March 26, 2008 Maintainer Longhai Li
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 Bergm. R topics documented: September 25, Type Package
Type Package Package Bergm September 25, 2018 Title Bayesian Exponential Random Graph Models Version 4.2.0 Date 2018-09-25 Author Alberto Caimo [aut, cre], Lampros Bouranis [aut], Robert Krause [aut] Nial
More informationStata goes BUGS (via R)
(via R) Department of Political Science I University of Mannheim 31. March 2006 4th German Stata Users Group press ctrl + l to start presentation Problem 1 Problem You are a Stata user... and have a complicated
More informationBayesian inference for psychology, part III: Parameter estimation in nonstandard models
Psychon Bull Rev (2018) 25:77 101 https://doi.org/10.3758/s13423-017-1394-5 Bayesian inference for psychology, part III: Parameter estimation in nonstandard models Dora Matzke 1 Udo Boehm 2 Joachim Vandekerckhove
More informationGetting started with simulating data in R: some helpful functions and how to use them Ariel Muldoon August 28, 2018
Getting started with simulating data in R: some helpful functions and how to use them Ariel Muldoon August 28, 2018 Contents Overview 2 Generating random numbers 2 rnorm() to generate random numbers from
More informationTying Hypothesis Tests, Power, and More Together with Functions
Tying Hypothesis Tests, Power, and More Together with Functions General Homework Notes Read the entire question. Show your answers. Don t give me the code and assume it shows the answers (e.g., plots,
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 informationBUGS language. Collect several definitions by using indexing and looping:
BUGS language Collect several definitions by using indexing and looping: for(i in 1:K){ X[i] ~ dbin(theta[i],n[i]) Give K as fixed value in the data listing. Every definition within for(i in 1:K){ should
More informationPackage TreeBUGS. December 18, 2018
Version 1.4.1 Date 2018-12-18 Package TreeBUGS December 18, 2018 Title Hierarchical Multinomial Processing Tree Modeling Author Daniel W. Heck [aut, cre], Nina R. Arnold [aut, dtc], Denis Arnold [aut],
More information=!"#$%!! 2! 2 (1 +!! )! (1 +!! )!! 2!! 2!!
MCEM is not a good choice when you are able to comput a close form of pdf/cdf. What's worse? If one only knows the kernel of pdf, things getts very very boring. For example, assmume!!!"#!!,!!!! =!"#$%!!
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 informationR Package rube (Really Useful WinBUGS (or JAGS) Enhancer) Version 0.3-5, March 26, 2014 Author: Howard J. Seltman
R Package rube (Really Useful WinBUGS (or JAGS) Enhancer) Version 0.3-5, March 26, 2014 Author: Howard J. Seltman 1 Introduction These notes are on the world wide web at: http://www.stat.cmu.edu/ hseltman/rube/rubemanual.pdf.
More informationPackage LCMCR. R topics documented: July 8, 2017
Package LCMCR July 8, 2017 Type Package Title Bayesian Non-Parametric Latent-Class Capture-Recapture Version 0.4.3 Date 2017-07-07 Author Daniel Manrique-Vallier Bayesian population size estimation using
More informationPackage bmeta. R topics documented: January 8, Type Package
Type Package Package bmeta January 8, 2016 Title Bayesian Meta-Analysis and Meta-Regression Version 0.1.2 Date 2016-01-08 Author Tao Ding, Gianluca Baio Maintainer Gianluca Baio
More informationPackage dclone. February 26, Type Package
Type Package Package dclone February 26, 2018 Title Data Cloning and MCMC Tools for Maximum Likelihood Methods Version 2.2-0 Date 2018-02-26 Author Peter Solymos Maintainer Peter Solymos
More informationPackage pcnetmeta. R topics documented: November 7, 2017
Type Package Title Patient-Centered Network Meta-Analysis Version 2.6 Date 2017-11-07 Author Lifeng Lin, Jing Zhang, and Haitao Chu Maintainer Lifeng Lin Depends R (>= 2.14.0), rjags
More informationPackage mederrrank. R topics documented: July 8, 2015
Package mederrrank July 8, 2015 Title Bayesian Methods for Identifying the Most Harmful Medication Errors Version 0.0.8 Date 2015-07-06 Two distinct but related statistical approaches to the problem of
More informationNIMBLE User Manual. NIMBLE Development Team. Version 0.4-1
NIMBLE User Manual NIMBLE Development Team Version 0.4-1 Contents 1 Welcome to NIMBLE 5 1.1 Why something new?............................... 5 1.2 What does NIMBLE do?............................. 6 1.3
More informationLAB #2: SAMPLING, SAMPLING DISTRIBUTIONS, AND THE CLT
NAVAL POSTGRADUATE SCHOOL LAB #2: SAMPLING, SAMPLING DISTRIBUTIONS, AND THE CLT Statistics (OA3102) Lab #2: Sampling, Sampling Distributions, and the Central Limit Theorem Goal: Use R to demonstrate sampling
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 informationHierarchical modeling of individual subjects nested within treatment conditions
February 22, 2017 File = D:\bcm\ch10.indiv.difs.mem.reten.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 informationPackage MfUSampler. June 13, 2017
Package MfUSampler June 13, 2017 Type Package Title Multivariate-from-Univariate (MfU) MCMC Sampler Version 1.0.4 Date 2017-06-09 Author Alireza S. Mahani, Mansour T.A. Sharabiani Maintainer Alireza S.
More informationEcography. Supplementary material
Ecography ECOG-02306 Yen, J. D. L., Thomson, J. R., Keith, J. M., Paganin, D. M. and Mac Nally, R. 2016. How do different aspects of biodiversity change through time? A case study on an Australian bird
More informationThe GSM Package. August 9, 2007
The GSM Package August 9, 2007 Title Gamma Shape Mixture This package implements a Bayesian approach for estimation of a mixture of gamma distributions in which the mixing occurs over the shape parameter.
More informationPackage CNVrd2. March 18, 2019
Type Package Package CNVrd2 March 18, 2019 Title CNVrd2: a read depth-based method to detect and genotype complex common copy number variants from next generation sequencing data. Version 1.20.0 Date 2014-10-04
More informationrcppbugs native MCMC for R
rcppbugs native MCMC for R Whit Armstrong armstrong.whit@gmail.com KLS Diversified Asset Management May 12, 2012 rcppbugs native MCMC for R 1 / 33 What is this rcppbugs thing? Or as Paul Johnson put it
More informationPackage tsbugs. February 20, 2015
Type Package Title Create time series BUGS models. Version 1.2 Author Guy J. Abel Package tsbugs February 20, 2015 Maintainer ``Guy J. Abel'' Depends R(>= 2.15.2) Suggests R2OpenBUGS,
More informationPackage gibbs.met. February 19, 2015
Version 1.1-3 Title Naive Gibbs Sampling with Metropolis Steps Author Longhai Li Package gibbs.met February 19, 2015 Maintainer Longhai Li Depends R (>=
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 informationBayesian Analysis with Stata
Bayesian Analysis with Stata John Thompson Department of Health Sciences Univeristy of Leicester john.thompson@le.ac.uk March 19, 2014 1 Introduction The manual offers a printed version of the information
More informationReddit Recommendation System Daniel Poon, Yu Wu, David (Qifan) Zhang CS229, Stanford University December 11 th, 2011
Reddit Recommendation System Daniel Poon, Yu Wu, David (Qifan) Zhang CS229, Stanford University December 11 th, 2011 1. Introduction Reddit is one of the most popular online social news websites with millions
More informationPackage BKPC. March 13, 2018
Type Package Title Bayesian Kernel Projection Classifier Version 1.0.1 Date 2018-03-06 Author K. Domijan Maintainer K. Domijan Package BKPC March 13, 2018 Description Bayesian kernel
More information