Fix zero dose + evid=2 to not be counted as dose

nlmixr2 · Feb 7, 2024 · 3e9b16b · 3e9b16b
1 parent 8249e25
commit 3e9b16b
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 2 deletions.
diff --git a/inst/include/rxode2parseHandleEvid.h b/inst/include/rxode2parseHandleEvid.h
@@ -229,6 +229,13 @@ static inline int handleTlastInlineUpateDosingInformation(rx_solving_options_ind
     }
     break;
   }
+  if (curDose[0] == 0.0 &&
+      getEvidP1(ind, ind->idx) == 2) {
+    // evid=2 can add zero dose to turn on a compartment;
+    // therefore if the current dose is zero and the next evid = 2
+    // don't treat it as a new dose.
+    return 0;
+  }
   return 1;
 }
 

diff --git a/src/etTran.cpp b/src/etTran.cpp
@@ -1220,7 +1220,7 @@ List etTransParse(List inData, List mv, bool addCmt=false,
     else caddl = inAddl[i];
     if (IntegerVector::is_na(caddl)) caddl = 0;
     // EVID flag
-    if (evidCol == -1){
+    if (evidCol == -1) {
       // Missing EVID
       if (rateI == 0 && (ISNA(camt) || camt == 0.0)){
         cevid = 0;
@@ -1252,7 +1252,7 @@ List etTransParse(List inData, List mv, bool addCmt=false,
         cevid=0;
       }
     }
-    switch(cevid){
+    switch(cevid) {
     case 0:
       // Observation
       nobs++;