Nothing
crit <-
function(theta, S,n, lam1, lam2, penalize.diagonal) # theta = list of pXp matrices, length k
{
p = dim(S[[1]])[1]
K = length(S)
lam1 = penalty.as.matrix(lam1,p,penalize.diagonal=penalize.diagonal)
lam2 = penalty.as.matrix(lam2,p,penalize.diagonal=TRUE)
crit = 0
for(k in 1:length(theta))
{
# add log det that was entered as an argument, or else calculate it
crit = crit+n[k]*log(det(theta[[k]]))-n[k]*sum(S[[k]]*theta[[k]])-sum(lam1*abs(theta[[k]]))
for(kp in k:length(theta))
{
crit = crit - sum(lam2*abs(theta[[k]]-theta[[kp]]))
}
}
return(crit)
}
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.