Nothing
loglik.M14 <-
function (pa, z, f)
{
paa <- pa^2 + pa * (1 - pa) * f
pab <- 2 * pa * (1 - pa) * (1 - f)
pbb <- (1 - pa)^2 + pa * (1 - pa) * f
pvec <- c(paa, pab, pbb)
ind <- !(z == 0)
logvec <- log(pvec[ind])
loglik0 <- sum(z[ind] * logvec)
nparam0 <- 2
res <- c(loglik0, nparam0)
return(res)
}
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.