Nothing
extendedBIC <-
function(gamma,omegahat,S,n) {
p = nrow(omegahat)
es = sum(omegahat[upper.tri(omegahat)]!=0)
return(-log(det(omegahat)) + sum(diag(omegahat %*% S)) + es * (log(n)/n) + es * gamma * (4*log(p)/n))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.