tests/testthat/test_NMreadSim.R

In NMsim: Seamless 'Nonmem' Simulation Platform

## library(devtools)

#### need a function to drop NMsimVersion and NMsimTime from table
fix.time <- function(x){
    meta.x <- attr(x,"NMsimModTab")
    ## meta.x$time.call <- as.POSIXct("2020-02-01 00:01:01",tz="UTC")
    meta.x$NMsimVersion <- NULL
    meta.x$NMsimTime <- NULL

    meta.x$path.lst.read <- NULL
    
    setattr(x,"NMsimModTab",meta.x)
    invisible(x)
}

library(NMdata)
NMdataConf(reset=TRUE)
library(data.table)
data.table::setDTthreads(1)




context("NMreadSim")

## prepare a sim to read
dt.amt <- data.table(DOSE=c(100,400))
dt.amt[,AMT:=DOSE*1000]
dt.amt
doses.sd <- NMcreateDoses(TIME=0,AMT=dt.amt,as.fun="data.table")
doses.sd[,dose:=paste(DOSE,"mg")]
doses.sd[,regimen:="SD"]
doses.sd
dat.sim.sd <- addEVID2(doses.sd,time.sim=0:24,CMT=2,as.fun="data.table")
dat.sim <- copy(dat.sim.sd)
dat.sim[,ROW:=.I]
dat.sim[,BBW:=75]

if(F){

    ## testOutput/NMsim_xgxr021_sd1_NMreadSim_paths.rds
    file.mod <- "testData/nonmem/xgxr021.mod"

    sim1 <- NMsim(file.mod=file.mod,
                  data=dat.sim,
                  dir.sims="testOutput",
                  dir.res="testData/simres",
                  name.sim = "sd1_NMreadSim",
                  seed.nm=2342
                  ## ,reuse.results=TRUE
                 ,nmquiet=F)

    ## unlink("testOutput/xgxr021_sd1_NMreadSim",recursive=T)
}



## if(F){
test_that("Basic",{
    fileRef <- "testReference/NMreadSim_01.rds"
    ## ref <- readRDS(fileRef)
    res1 <- NMreadSim("testData/simres/xgxr021_sd1_NMreadSim_MetaData.rds")

    fix.time(res1)
    
    expect_equal_to_reference(res1,fileRef)

    if(F){
        ref <- readRDS(fileRef)
        compareCols(res1,ref)

        compareCols(
            attributes(res1)$NMsimModTab
           ,
            attributes(ref)$NMsimModTab
           ,keep.names=FALSE)
    }

})
## }

## if(F){
test_that("Reading fst directly",{
    ## NMdataConf(as.fun="data.table")
    
    fileRef <- "testReference/NMreadSim_02.rds"
    ## ref <- readRDS(fileRef)
    res1 <- NMreadSim("testData/simres/xgxr021_sd1_NMreadSim_ResultsData.fst")
    ## library(fst)
    ## res1 <- read_fst("testOutput/xgxr021_sd1_NMreadSim_paths_res.fst",as.data.table=T)

    fix.time(res1)
    
    expect_equal_to_reference(res1,fileRef)

    if(F){
        ref <- readRDS(fileRef)
        compareCols(res1,ref)

### attributes(res1)$NMsimModTab does not exist - thats why one should read the rds
        ## compareCols(
        ##     attributes(res1)$NMsimModTab
        ##    ,
        ##             attributes(ref)$NMsimModTab
        ##            ,keep.names=FALSE
        ##             )
    }

})

## }

## if(F){
test_that("From different wd",{
##    setwd("..")
    fileRef <- "testReference/NMreadSim_03.rds"
    ## ref <- readRDS(fileRef)
    res1 <- NMreadSim("testData/simres/xgxr021_sd1_NMreadSim_MetaData.rds")

    ## setwd("testthat")
    
    fix.time(res1)
    
    ## this should compare to _01 results instead
    expect_equal_to_reference(res1,fileRef)

    
    if(F){
        ref <- readRDS(fileRef)
        compareCols(res1,ref)

        compareCols(
            attributes(res1)$NMsimModTab
           ,
            attributes(ref)$NMsimModTab
           ,keep.names=FALSE)
    }


})
## }