analyze.R
392 Bytes
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)