knitr::opts_chunk$set(echo = TRUE)
#source("Code/LearnBayes-discrete.R")
n <- 50 b0.true <- 7 b1.true <- 2 tau.true <- .01 x <- 1:n y <- rnorm(n,b0.true+b1.true*x,sd=sqrt(1/tau.true)) plot(x,y) print(lm(y~x))
$$ f(y_i|\lambda) = \frac{\lambda_i^{y_i}}{y_i!}\exp{-\lambda_i}$$ $$ \lambda_i=\exp(\alpha+\beta x_i) $$ $$ \alpha \sim{\rm N}(0,v_\alpha^2)$$ $$ \beta \sim{\rm N}(0,v_\beta^2)$$ $$ p(\alpha,\beta|x,y) \propto \prod\exp(\alpha+\beta x_i)^{y_i}\exp(-\exp(\alpha+\beta x_i)) \exp(-0.5v_a^{-2}\alpha^2) \exp(-0.5v_b^{-2}\beta^2) $$
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.