In AstraZeneca/blueprint: create models easily

library(blueprint)
library(tidyverse)

blueprint <- Blueprint$new("nonmem")

templates <- load_templates("nonmem")
models <- available_models("nonmem")

map2(templates, names(templates), ~ {
  cat(glue::glue("--------  {.y} -------- \n\n"))
  cat(.x)
})

one_cmpt_iv <- blueprint %>% 
  parameters(CL = parameter(5,  
                    lower_bound = 0.01, 
                    upper_bound = 10),
         V = parameter(34, lower_bound = 0.1)) 

one_cmpt_iv$get_all_elements()

one_cmpt_iv <- one_cmpt_iv %>%
  constants(STD_WT = 70) %>%
  residual_error(prop = sigma_param(0.04, "PROP", comment = "Proportional Error")) %>%
  with_data(head(Theoph)) %>%
  from_path("path/to/my_data.csv") 

one_cmpt_iv %>%
  hierarchies(CL = 0.04, V = 0.09) %>% 
  model_type(models$one_cmt_iv) %>%
  use_template(templates$compartmental)%>% 
  render() %>% cat()

one_cmpt_iv %>%
  hierarchies(b1 = block(0.04, 0.01, 0.09, param_names = c("CL", "V"))) %>% 
  model_type(models$one_cmt_iv) %>%
  use_template(templates$compartmental) %>% 
  render() %>% cat()

one_cmpt_iv %>%
  hierarchies(bl = block(0.04, 0.06, 0.09, params = c("CL", "V"))) %>% 
  use_template(templates$one_cmt_iv) %>% render() %>% cat()

AstraZeneca/blueprint documentation built on May 17, 2019, 11:32 a.m.