x <- c(3,4,2,1,2,3)
n <- length(x)
T <- sum(x)
nuprime <- 2*T + 1
S1 <- 2*n
postmean <- nuprime/S1
postvar <- 2*nuprime/S1^2
cat("Posterior has mean",postmean,"and variance",postvar,"\n")