Skip to content

Instantly share code, notes, and snippets.

@tvladeck
Created December 22, 2016 16:12
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save tvladeck/43441e866d8b5b7cecce9bf430cff9c7 to your computer and use it in GitHub Desktop.
Save tvladeck/43441e866d8b5b7cecce9bf430cff9c7 to your computer and use it in GitHub Desktop.
i don't quite get why you run _coda.samples_ twice
# Betas #
samples <- coda.samples(jags, c('beta_age',
'beta_percent_white',
'beta_uninsured',
'beta_unemployed',
'beta_percent_degree',
'beta_income',
'beta_health'), 10000)
samples <- as.data.frame(samples[[1]])
# Other Parameters #
other_parameters <- coda.samples(jags, c('sigma_state',
'sigma_county',
'alpha'), 10000)
other_parameters <- as.data.frame(other_parameters[[1]])
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment