sd0 <- 1.1
phi0 <- sd0^2
theta0 <- 82.4
sdobs <- 1.7
phi <- sdobs^2
x <- 82.1
phi1 <- (phi0^(-1) + phi^(-1))^(-1)
sd1 <- sqrt(phi1)
theta1 <- phi1*(theta0/phi0 + x/phi)
hypvalue <- 83.0
pi0 <- pnorm(hypvalue,theta0,sd0)
priorodds <- pi0/(1-pi0)
p0 <- pnorm(hypvalue,theta1,sd1)
postodds <- p0/(1-p0)
B <- postodds/priorodds
cat("Prior prob",pi0,"and prior odds",priorodds,"\n")
cat("Prior prob",p0,"and prior odds",postodds,"\n")
cat("Bayes factor",B,"\n")