Nothing
normal_update <- function(prior, data, teach=FALSE){
prior <- unlist(prior)
prior_mean <- prior[1]
prior_sd <- prior[2]
data_mean <- data[1]
data_se <- data[2]
prior_precision <- 1 / prior_sd ^ 2
data_precision <- 1 / data_se ^ 2
post_precision <- prior_precision + data_precision
post_sd <- sqrt(1 / post_precision)
post_mean <- weighted.mean(c(prior_mean, data_mean),
c(prior_precision, data_precision))
if(teach==TRUE)
data.frame(Type=c("Prior", "Data", "Posterior"),
Mean=c(prior_mean, data_mean, post_mean),
Precision=c(prior_precision, data_precision,
post_precision),
Stand_Dev=c(prior_sd, data_se, post_sd)) else
c(post_mean, post_sd)
}
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.