counts <- scan("/home/wkarcher/Desktop/Home/RT08/4/counts.dat")
sigma <- sd(counts)^2
mu <- mean(counts)
iod <- sigma/mu
iod