mattblackwell

## Function to implment the MCMC sampler
polling.gibbs.kal <- function(y, size, dates, states, big, ## data
                              F, X=NULL, XL = NULL, Z=NULL,     ## covariates
                              m.0, C.0, b = NULL, B = NULL,  ## hyper-parameters
                              c.0, d.0, e.0, f.0, g.0, h.0,
                              starts, phi, sigma.a, nu, beta,          ## starting values
                              tune, A, tune.phi,
                              iters = 10000, thin = 1, burnin = 0) { ## mcmc stuff

  aa <- proc.time()

mattblackwell

### Create plots for MO talk

library(tikzDevice)

x <- rnorm(200,0,.4)

y <- x + rnorm(200, 0, .4)

night <- rgb(26,31,30, max=255)
beach <- rgb(227,223,186, max=255)

mattblackwell

## Function to implment the MCMC sampler
polling.gibbs.kal.fixedS <- function(y, size, dates, states, big, ## data
                              F, X=NULL, XL = NULL, Z=NULL,     ## covariates
                              m.0, C.0, b = NULL, B = NULL,  ## hyper-parameters
                              starts, phi, sigma.a, nu, beta, k,         ## starting values
                              tune, A,
                              iters = 10000, thin = 1, burnin = 0) { ## mcmc stuff

  aa <- proc.time()
  ## pre-gibbs cleaning (calculate variances, setup times, etc)

mattblackwell

library(Amelia)
data(africa)
a.out <- amelia(africa, ts = "year", cs = "country")

overimpute(a.out, "gdp_pc")
overimpute(a.out, "gdp_pc", subset = country == "Burkina Faso")
overimpute(a.out, "gdp_pc", subset = gdp_pc >=500)

mattblackwell

#! /bin/sh

if [ -r "$1" ]; then
    mdfile="$1"
    Rscript -e "rmarkdown::render('${mdfile}', 'all')"
else
    exit 0
fi

mattblackwell

mattblackwell

group.mean <- function(x, group) {
  out <- tapply(x, group, mean, na.rm = TRUE)
  out[group]
}

data(state.fips, package = "maps")
state.fips <- unique(state.fips[,c("fips","abb")])
state.fips$abb <- as.character(state.fips$abb)
state.fips <- rbind(state.fips, c(2, "AK"))
state.fips <- rbind(state.fips, c(15, "HI"))

mattblackwell

---
title: 'Beamer Overlay Test'
output:
  beamer_presentation:
    incremental: true
    theme: CambridgeUS
---


```{r, echo = FALSE}