fixing bug in R version of eppasm.R

mrc-ide · Aug 9, 2024 · 9b1da23 · 9b1da23
1 parent c56bc27
commit 9b1da23
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/eppasm.R b/R/eppasm.R
@@ -250,7 +250,7 @@ simmod <- function(fp, VERSION = "C") {
           testneg_infections_ha <- infections_ha * (testnegpop[,,hivn.idx,i] / hivn_pop_ha)
 
           grad_tn[ , , hivn.idx] <- grad_tn[ , , hivn.idx] - testneg_infections_ha
-          grad_tn[ , , hivp.idx] <- grad_tn[ , , hivp.idx] - testneg_infections_ha
+          grad_tn[ , , hivp.idx] <- grad_tn[ , , hivp.idx] + testneg_infections_ha
         }
 
         ## Do new diagnoses