rejectsampling=function (logf, tpar, dmax, n, data)
{
d = length(tpar$m)
theta = rmt(n, mean = c(tpar$m), S = tpar$var, df = tpar$df)
lf = matrix(0, c(dim(theta)[1], 1))
for (j in 1:dim(theta)[1]) lf[j] = logf(theta[j, ], data)
lg = dmt(theta, mean = c(tpar$m), S = tpar$var, df = tpar$df,
log = TRUE)
if (d == 1) {
prob = exp(c(lf) - lg - dmax)
return(theta[runif(n) < prob])
}
else {
prob = exp(lf - lg - dmax)
return(theta[runif(n) < prob, ])
}
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.