sc <- read.table(file('stdin')) | |
scratio <- sc$V2 / sc$V3 | |
scratio <- scratio[scratio > 0] | |
# Take the log of the ratio. Our null hypothesis is a normal distribution | |
# around zero. | |
tt <- t.test(log(scratio)) | |
tt | |
# This gives us the geo-mean as we are taking the exponent of the linear mean | |
# of logarithms. | |
1 - 1/exp(tt$estimate) | |
# Likewise for the confidence interval. | |
1 - 1/exp(tt$conf.int) |