Nothing
EMGaussMax <-
function(Z, W, yw, ptype, mutype, sdtype, sd.fixed=0.05, p) {
ZW <- W %*% Z
mu <- EMMax.mu(yw, ZW, mutype)
sigma <- EMMax.sd(yw, ZW, mu, sdtype, sd.fixed)
p <- EMMax.p(Z, diag(W), p, ptype)
list(mu=mu,sigma=sigma,p=p)
#so mu is optimized first, these are used to optimize sigma
#the new p are derived from the current p's of the samples
}
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.