in progress

Author: d-morrison

R/antibody_decay_function.R

@@ -25,7 +25,7 @@ antibody_decay_curve = function(
     alpha = 0.00002192627,
     rho = 2
     # rho = 2 # exponential decay?
-    )
+)
 {
 
   t1 = t1f(
@@ -41,10 +41,18 @@ antibody_decay_curve = function(
     mu_y = mu_y,
     t1 = t1)
 
+  phase2_rho1 <- y1 * exp(-alpha * (t - t1))
+  phase2_default <-
+    y1 * (1+ (rho-1) * (y1 ^ (rho - 1)) * alpha * (t - t1))^(-1/(rho - 1))
   yt = ifelse(
     t < t1,
     y0 * exp(mu_y * t),
-    y1 * (1+ (rho-1) * (y1 ^ (rho - 1)) * alpha * (t - t1))^(-1/(rho - 1))
+    ifelse(
+      rho == 1,
+      phase2_rho1,
+      phase2_default
+
+    )
   )
   return(yt)
 }