### R code from vignette source 'case-con-NN-sol.rnw'

###################################################
### code chunk number 1: case-con-NN-sol.rnw:15-18
###################################################
EOR.crude <- (8/(1000-8)) / (2/(5000-2))
EOR.c <- (8/1)/(2/5) 
cbind( EOR.crude, EOR.c )


###################################################
### code chunk number 2: case-con-NN-sol.rnw:26-29
###################################################
SE.crude <- sqrt( 1/8 + 1/(1000-8) + 1/2 + 1/(5000-2) )
SE.c <- sqrt( 1/8 + 1/2 )  
cbind( SE.crude, SE.c )


###################################################
### code chunk number 3: case-con-NN-sol.rnw:33-36
###################################################
EOR.c95low <- EOR.c / exp(1.96*SE.c) 
EOR.c95up  <- EOR.c * exp(1.96*SE.c)
round( cbind( EOR.c, SE.c, EOR.c95low, EOR.c95up ), 3 )