emp_var_function <- function(mu, scm_list) {
P <- length(scm_list)
P^(-1)*Reduce("+",lapply(scm_list, function(x) RaoDist(x, mu)))
}
eta <- function(sigma) {
-1/(2*sigma^2)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.