Nothing
winit <-
function(x,n,G,mug,sigmainv,vg,p,sg,zmat){
w <- matrix(0,n,G)
delta <- matrix(0,n,G)
for(g in 1:G){
delta[,g] <- mahalanobis(x, mug[g,], sigmainv[,,g], inverted=TRUE)
w[,g] <- (vg[g]+p)/(vg[g]+delta[,g])
}
w
}
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.