From 6c2cbf0c68b343601e098f6c4699066c305f4aa8 Mon Sep 17 00:00:00 2001
From: Bill Denney <wdenney@humanpredictions.com>
Date: Sat, 21 Oct 2023 05:06:04 -0400
Subject: [PATCH] Add PKNCA to README.md

---
 README.Rmd | 12 +++++++++++-
 README.md  | 32 +++++++++++++++++++++-----------
 2 files changed, 32 insertions(+), 12 deletions(-)

diff --git a/README.Rmd b/README.Rmd
index e33a8daa..6bfe4004 100644
--- a/README.Rmd
+++ b/README.Rmd
@@ -74,7 +74,17 @@ command line.
 
 ## Example
 
-After installed, if you use the standard interface, you can convert to Monolix with
+After installed, if you use the standard interface, you can obtain new initial estimates with PKNCA:
+
+```r
+mod <-
+  nlmixr2(
+    nlmixrFun, nlmmixrData, est = "pknca",
+    control = pkncaControl(concu = "ng/mL", doseu = "mg", timeu = "hr", volumeu = "L")
+  )
+```
+
+or, you can convert to Monolix with
 
 ```r
 mod <- nlmixr(nlmixrFun, nlmmixrData, est="monolix")
diff --git a/README.md b/README.md
index 2837c896..93e7e8c8 100644
--- a/README.md
+++ b/README.md
@@ -46,17 +46,17 @@ Otherwise you can always install from GitHub:
 
 Babelmixr2 can help you by:
 
-  - Running your nlmixr2 model in a commercial nonlinear mixed effects
-    modeling tool like
-    [`NONMEM`](https://nlmixr2.github.io/babelmixr2/articles/running-nonmem.html)
-    or `Monolix`
+- Running your nlmixr2 model in a commercial nonlinear mixed effects
+  modeling tool like
+  [`NONMEM`](https://nlmixr2.github.io/babelmixr2/articles/running-nonmem.html)
+  or `Monolix`
 
-  - Convert your [`NONMEM` model to a nlmixr2
-    model](https://nlmixr2.github.io/nonmem2rx/articles/convert-nlmixr2.html)
-    (in conjunction with `nonmem2rx`)
+- Convert your [`NONMEM` model to a nlmixr2
+  model](https://nlmixr2.github.io/nonmem2rx/articles/convert-nlmixr2.html)
+  (in conjunction with `nonmem2rx`)
 
-  - Calculate scaling factors and automatically add initial conditions
-    based on non-compartmental analysis (using `PKNCA`)
+- Calculate scaling factors and automatically add initial conditions
+  based on non-compartmental analysis (using `PKNCA`)
 
 ## Monolix Setup
 
@@ -68,8 +68,18 @@ instead of using a command line.
 
 ## Example
 
-After installed, if you use the standard interface, you can convert to
-Monolix with
+After installed, if you use the standard interface, you can obtain new
+initial estimates with PKNCA:
+
+``` r
+mod <-
+  nlmixr2(
+    nlmixrFun, nlmmixrData, est = "pknca",
+    control = pkncaControl(concu = "ng/mL", doseu = "mg", timeu = "hr", volumeu = "L")
+  )
+```
+
+or, you can convert to Monolix with
 
 ``` r
 mod <- nlmixr(nlmixrFun, nlmmixrData, est="monolix")