log.kern.iwish <-
function(W, v, S)
{
k <- dim(S)[[1]]
out <- - ((v+k+1)/2)*log(det(W)) - (1/2)*trace(S%*%solve(W))
return(out)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.