In mrgsolve/mrgsim: Scripted simulation with 'mrgsolve'

Scripted simulation runs with mrgsolve

library(knitr)
knitr::opts_chunk$set(comment='.', message=FALSE)

library(mrgsim)
library(mrgsolve)
library(dplyr)

Load a model

mod <- mread("popex", modlib())

param(mod)

Load a scripted simulation run

run <- load_run("inst/yaml/foo.yaml")

Simulate

out <- sim_run(mod,run)

runtime <- system.time(out <- sim_run(mod,run))

runtime

dim(out)

out %>%
  group_by(arm) %>% 
  filter(time==4032) %>%
  summarise(Mean = mean(DV), Min = min(DV), Max = max(DV), N=n())

The simulation run specification

yam <- readLines("inst/yaml/foo.yaml")
cat(yam, sep="\n")

Run specification for pediatric / adult pk

yam <- readLines("inst/yaml/ped_adult.yaml")
cat(yam, sep="\n")

Run specification details

The run has arms

arms(run)

treatment sequences

sequences(run)

treatment periods

run$periods

covariate sets

covsets(run)

sampling times

run$designs

mrgsolve/mrgsim documentation built on May 30, 2019, 8:07 p.m.