tests/test-medicin-macro.R

In tagteam/heaven: Data Preparation Routines for Medical Registry Data

### test-medicin-macro.R --- 
##----------------------------------------------------------------------
## Author: Thomas Alexander Gerds
## Created: Mar 11 2019 (12:57) 
## Version: 
## Last-Updated: Jan 17 2020 (07:52) 
##           By: Thomas Alexander Gerds
##     Update #: 12
##----------------------------------------------------------------------
## 
### Commentary: 
## 
### Change Log:
##----------------------------------------------------------------------
## 
### Code:
library(testthat)
library(data.table)
library(heaven)

testmed <- data.table(pnr=1,apk=1L,atc="C03",eksd=as.Date(c("2010-03-18","2010-05-25","2010-08-10","2010-10-05","2010-11-24","2011-01-13","2011-03-07","2011-05-03","2011-06-06","2011-08-02","2011-09-27","2012-02-27","2012-04-23","2012-08-06","2012-11-05","2013-01-15","2013-05-14","2013-08-05","2013-10-28","2014-02-04")),packsize=100,strnum=40)
testlpr <- simAdmissionData(1)
testlpr[,pattype:=0]
drug1=list(atc=c("C03"),
           maxdepot=4000,
           period=as.Date(c("2000-01-01", "2016-12-31")),
           prescriptionwindow=2,
           doses=list(value=c(40, 400, 500,750),min = c(20, 200, 250,750),max = c(160, 800, 1000,750),def = c(40, 400, 500,750)))

test_that("minimal example",{
    a <- medicinMacro(drugs=list("drug1"=drug1),
                      drugdb=testmed,
                      admdb=testlpr,
                      admdb.datevars=c("inddto","uddto"))
    expect_equal(names(a[[1]]),c("pnr","dose","firstday","lastday","exposure.days"))
})
##
test_that("inconsistent pnr format",{
    testmed1 <- copy(testmed)
    testmed1[,pnr:=as.character(pnr)]
    expect_error(medicinMacro(drugs=list("drug1"=drug1),
                              drugdb=testmed1,
                              admdb=testlpr,
                              admdb.datevars=c("inddto","uddto")))
})

test_that("pnr as character",{
    a <- medicinMacro(drugs=list("drug1"=drug1),drugdb=testmed,admdb=testlpr,admdb.datevars=c("inddto","uddto"))
    testmed[,pnr:=as.character(pnr)]
    testlpr[,pnr:=as.character(pnr)]
    u <- medicinMacro(drugs=list("drug1"=drug1),drugdb=testmed,admdb=testlpr,admdb.datevars=c("inddto","uddto"))
    a[[1]][,pnr:=as.character(pnr)]
    expect_equal(a,u)
})

test_that("rounding issue",{
    testmed2 <- testmed[1:3]
    testmed2[,apk:=as.numeric(apk)]
    testmed2[1,apk:=0.00999998]
    testlpr <- getAdmLimits(testlpr,collapse=TRUE)
    b <- medicinMacro(drugs=list("drug1"=drug1),drugdb=testmed2[1],admdb=testlpr)
    b
})

######################################################################
### test-medicin-macro.R ends here