R/p_t.R

Defines functions p_t

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
aashen12/BARS documentation built on Aug. 21, 2020, 2:44 p.m.