p_t <- function(sig_sq, betahat, X, y) {
(-1/(2*sig_sq)) * t(y - (X %*% betahat)) %*% (y - (X %*% betahat))
} # full conditional for the knot locations, t
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.