tests/testthat/test_edgewiseID.R

In SEMID: Identifiability of Linear Structural Equation Models

library(SEMID)
context("Components related to edgewise identification.")

source("helperFunctions.R")

test_that("Edgewise identification does not identify edges erroneously.", {
    # Random test
    set.seed(2323)
    ps <- c(0.3)
    sims <- 30
    ns <- c(5, 6)
    for (p in ps) {
        for (n in ns) {
            for (i in 1:sims) {
                L <- rDirectedAdjMatrix(n, p)
                O <- rUndirectedAdjMat(n, p)
                g <- MixedGraph(L, O)
                gid <- htcID(g)
                eid <- edgewiseID(g)

                expect_true(all(sapply(1:n, FUN = function(x) {
                  all(gid$solvedParents[[x]] %in% eid$solvedParents[[x]])
                })))

                randomIdentificationTest(eid$identifier, L, O, eid$solvedParents)
            }
        }
    }
})



test_that("TO DELETE", {
  set.seed(2323)
  ps <- c(0.3)
  sims <- 30
  ns <- 3
  for (p in ps) {
    for (n in ns) {
      for (i in 1:sims) {
        L <- rDirectedAdjMatrix(n, p)
        O <- rUndirectedAdjMat(n, p)
        g <- MixedGraph(L, O)
        gid <- htcID(g, tianDecompose = F)

        randomIdentificationTest(gid$identifier, L, O, gid$solvedParents)
      }
    }
  }
})