From 16f2fed6ff3351cf6fcf9c1accd79fba13c5445b Mon Sep 17 00:00:00 2001
From: Vincent LAVILLE <vincent.laville@pasteur.fr>
Date: Mon, 19 Mar 2018 15:04:59 +0100
Subject: [PATCH] Fix bug - Variance binomial exposure
---
R/varExp.R | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/R/varExp.R b/R/varExp.R
index bbc6e7e..61c1403 100644
--- a/R/varExp.R
+++ b/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))
}
################################################################################
--
GitLab