tests/testthat/test-PPMC.R

In JihongZ/StanDCM: Diagnostic Classification Model using Stan

library(testthat)
context("Test PPMC function")


test_that("PPMC test",{
  mod1 <- StanDINA.run(Qmatrix,respMatrix, iter=10,init.list='cdm', warmup = 20)
  mod2 <- StanDINO.run(Qmatrix,respMatrix, iter=10,init.list='cdm', warmup = 20)
  mod3 <- StanLCDM.run(Qmatrix,respMatrix, iter=10,init.list='cdm', warmup = 20)
  mod4 <- StanNCRUM.run(Qmatrix,respMatrix, iter=10,init.list='cdm', warmup = 20)
  mod5 <- StanCRUM.run(Qmatrix,respMatrix, iter=10,init.list='cdm', warmup = 20)
  mod6 <- StanORDM.run(Qmatrix,respMatrix, iter=10,init.list='cdm', warmup = 20)
  ppmc.mod1 <- StanDCM.ppmc(mod1, respMatrix, n.sim = 500, n.burnin = 10, plot.option = TRUE)
  ppmc.mod2 <- StanDCM.ppmc(mod2, respMatrix, n.sim = 500, n.burnin = 10, plot.option = TRUE)
  ppmc.mod3 <- StanDCM.ppmc(mod3, respMatrix, n.sim = 500, n.burnin = 10, plot.option = TRUE)
  ppmc.mod4 <- StanDCM.ppmc(mod4, respMatrix, n.sim = 500, n.burnin = 10, plot.option = TRUE)
  ppmc.mod5 <- StanDCM.ppmc(mod5, respMatrix, n.sim = 500, n.burnin = 10, plot.option = TRUE)
  ppmc.mod6 <- StanDCM.ppmc(mod6, respMatrix, n.sim = 500, n.burnin = 10, plot.option = TRUE)
  print(ppmc.mod1)
  print(ppmc.mod2)
  print(ppmc.mod3)
  print(ppmc.mod4)
  print(ppmc.mod5)
  print(ppmc.mod6)
})