tests/testthat/test-standino_mg.R

In JihongZ/StanDCM: Diagnostic Classification Model using Stan

library(testthat)
context("Test multigroup models")

test_that("multigoup DINO run",{
  # save.name="DINO_uninf_multiG"
  # Qmatrix=Qmatrix;response.matrix=respMatrix;GroupID=c(rep(1,500),rep(2,500));
  # fixeditem.vector=c(1,3);class.equal=F;chain.num=3;init.list='random';
  # save.path=getwd();
  # warmup=10;control.list<-list(adapt_delta=0.82);iter=200
  options(mc.cores = parallel::detectCores())
  StanDINO_mG.script(Qmatrix,2,c(1,2),class.equal=F)
  StanDINO_mG.run(Qmatrix = Qmatrix,
                  response.matrix= respMatrix,
                  GroupID = c(rep(1,500),rep(2,500)),
                  fixeditem.vector = c(1,3),
                  class.equal = F,
                  script.path=NA,
                  save.path=getwd(),
                  save.name = "DINO_uninf_multiG",
                  iter = 200,
                  warmup = 10,
                  chain.num = 3,
                  init.list = 'random',
                  control.list = list(adapt_delta=0.82))

  })

test_that("multigoup DINA run",{
  # save.name="DINO_uninf_multiG"
  # Qmatrix=Qmatrix;response.matrix=respMatrix;GroupID=c(rep(1,500),rep(2,500));
  # fixeditem.vector=c(1,3);class.equal=F;chain.num=3;init.list='random';
  # save.path=getwd();
  # warmup=10;control.list<-list(adapt_delta=0.82);iter=200
  options(mc.cores = parallel::detectCores())
  StanDINA_mG.script(Qmatrix,2,c(1,2),class.equal=F)
  StanDINA_mG.run(Qmatrix = Qmatrix,
                  response.matrix= respMatrix,
                  GroupID = c(rep(1,500),rep(2,500)),
                  fixeditem.vector = c(1,3),
                  class.equal = F,
                  script.path=NA,
                  save.path=getwd(),
                  save.name = "DINA_uninf_multiG",
                  iter = 200,
                  warmup = 10,
                  chain.num = 3,
                  init.list = 'random',
                  control.list = list(adapt_delta=0.82))

})

test_that("multigoup LCDM run",{
  # save.name="DINO_uninf_multiG"
  # Qmatrix=Qmatrix;response.matrix=respMatrix;GroupID=c(rep(1,500),rep(2,500));
  # fixeditem.vector=c(1,3);class.equal=F;chain.num=3;init.list='random';
  # save.path=getwd();
  # warmup=10;control.list<-list(adapt_delta=0.82);iter=200
  options(mc.cores = parallel::detectCores())
  StanLCDM_mG.script(Qmatrix,2,c(1,2),class.equal=F)
  StanLCDM_mG.run(Qmatrix = Qmatrix,
                  response.matrix= respMatrix,
                  GroupID = c(rep(1,500),rep(2,500)),
                  fixeditem.vector = c(1,3),
                  class.equal = F,
                  script.path=NA,
                  save.path=getwd(),
                  save.name = "LCDM_uninf_multiG",
                  iter = 200,
                  warmup = 10,
                  chain.num = 3,
                  init.list = 'random',
                  control.list = list(adapt_delta=0.82))

})

test_that("multigoup NCRUM run",{
  # save.name="DINO_uninf_multiG"
  # Qmatrix=Qmatrix;response.matrix=respMatrix;GroupID=c(rep(1,500),rep(2,500));
  # fixeditem.vector=c(1,3);class.equal=F;chain.num=3;init.list='random';
  # save.path=getwd();
  # warmup=10;control.list<-list(adapt_delta=0.82);iter=200
  options(mc.cores = parallel::detectCores())
  StanNCRUM_mG.script(Qmatrix,2,c(1,2),class.equal=F)
  StanNCRUM_mG.run(Qmatrix = Qmatrix,
                  response.matrix= respMatrix,
                  GroupID = c(rep(1,500),rep(2,500)),
                  fixeditem.vector = c(1,3),
                  class.equal = F,
                  script.path=NA,
                  save.path=getwd(),
                  save.name = "NCRUM_uninf_multiG",
                  iter = 200,
                  warmup = 10,
                  chain.num = 3,
                  init.list = 'random',
                  control.list = list(adapt_delta=0.82))

})

test_that("multigoup CRUM run",{
  # save.name="DINO_uninf_multiG"
  # Qmatrix=Qmatrix;response.matrix=respMatrix;GroupID=c(rep(1,500),rep(2,500));
  # fixeditem.vector=c(1,3);class.equal=F;chain.num=3;init.list='random';
  # save.path=getwd();
  # warmup=10;control.list<-list(adapt_delta=0.82);iter=200
  options(mc.cores = parallel::detectCores())
  StanCRUM_mG.script(Qmatrix,2,c(1,2),class.equal=F)
  StanCRUM_mG.run(Qmatrix = Qmatrix,
                  response.matrix= respMatrix,
                  GroupID = c(rep(1,500),rep(2,500)),
                  fixeditem.vector = c(1,3),
                  class.equal = F,
                  script.path=NA,
                  save.path=getwd(),
                  save.name = "CRUM_uninf_multiG",
                  iter = 200,
                  warmup = 10,
                  chain.num = 3,
                  init.list = 'random',
                  control.list = list(adapt_delta=0.82))

})