In tsahota/NMproject: Script Based 'NONMEM' Model Development

## DO NOT MODIFY THIS BLOCK (unless you know what you're doing)
library(knitr)
library(rprojroot)
opts_knit$set(root.dir = find_root(
  has_file(".Rprofile") | 
    is_rstudio_project | 
    is_r_package | 
    is_git_root
))
opts_chunk$set(echo = FALSE)
opts_chunk$set(message = FALSE)
if(!is.null(knitr::opts_knit$get('rmarkdown.pandoc.to'))){
  .m <- params$m  
}

library(NMproject)
library(dplyr)

model object

#wait_finish(.m)

.m %>% status_table()

m_finished <- .m %>%
  subset(status(.) %in% "finished")

m_parent <- parent_run(m_finished[1])

parent results for context

This can be made better later

m_parent %>% rr()

results

m_finished %>%
  coef() %>%
  bind_rows() %>%
  group_by(parameter) %>%
  summarise(low=quantile(FINAL,probs=0.025),
            med=median(FINAL),
            upp=quantile(FINAL,probs=0.975))

Appendix

original control file

m_parent %>% text() %>% dplyr::first()

control file diff

m_parent %>% nm_diff()

tsahota/NMproject documentation built on Oct. 1, 2022, 11:51 a.m.