mrgsolve 재현성 테스트

library(mrgsolve)

## 
## Attaching package: 'mrgsolve'

## The following object is masked from 'package:stats':
## 
##     filter

# AIMS - Ryan provided this setting script.
options(nlmixr2.cores = 1)
options(rxode2.numCores = 1)
Sys.setenv(OMP_NUM_THREADS = "1", MKL_NUM_THREADS = "1")
RNGkind(kind = "L'Ecuyer-CMRG")
set.seed(20251023)

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

## Building popex ... done.

out1 <- mod %>% 
  ev(amt = 100) %>% 
  idata_set(expand.idata(WT=runif(10, 40, 140))) %>% 
  zero_re() %>% 
  mrgsim(end = 48, delta = 0.1)
  
out1 %>% plot(DV~time)

비재현성 확인: 시드 설정 없이 재실행

out1_no_seed <- mod %>% 
  ev(amt = 100) %>% 
  idata_set(expand.idata(WT=runif(10, 40, 140))) %>% 
  zero_re() %>% 
  mrgsim(end = 48, delta = 0.1)
  
out1_no_seed %>% plot(DV~time)

identical(as.data.frame(out1), as.data.frame(out1_no_seed))

## [1] FALSE

재현성 확인: 초기 동일 시드로 재실행

# 재현성 확인: 동일 시드로 재실행
set.seed(20251023)
out2 <- mod %>% 
  ev(amt = 100) %>% 
  idata_set(expand.idata(WT=runif(10, 40, 140))) %>% 
  zero_re() %>% 
  mrgsim(end = 48, delta = 0.1)
  
out2 %>% plot(DV~time)

identical(as.data.frame(out1), as.data.frame(out2))

## [1] TRUE