Fixed reverse GP function

jameshay218 · May 17, 2021 · eea4390 · eea4390
1 parent dfdaff3
commit eea4390
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/R/incidence_functions.R b/R/incidence_functions.R
@@ -105,6 +105,7 @@ reverse_gp_model <- function(desired_probs, pars, times){
 
 
   ps <- sum(desired_probs)*desired_probs + 0.000000001
+  ps <- pars["overall_prob"]*ps/sum(ps)
   k1 <- -log((1/ps) -1)
 
   k_solve <- (k1 %*% solve(t(L_K)))[1,]