Fix bug - Variance binomial exposure

R/varExp.R

@@ -73,7 +73,7 @@ calculateParamsFromIndParams <- function(N, m, v) {
 #'
 calculateParamsFromCounts <- function(N, n_expo) {
     meanval <- sum(n_expo) / sum(N)
-    c(meanval, meanval * (1 - meanval))
+    c(meanval, sum(N) / (sum(N) - 1) * meanval * (1 - meanval))
 }
 
 ################################################################################