R Under development (unstable) (2024-02-05 r85863 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(pcalg)
> suppressWarnings(RNGversion("3.5.0"))
> 
> ## bugfix 02.05.2020 ####
> set.seed(542) ## seed=542 und p=2; seed=20499 und p=3
> p <- 2 ## 20
> rDAG <- randomDAG(p, prob = 0.2, lB=-1, uB=1)
> B <- as(rDAG, "matrix")
> cov.th <- t(solve(diag(p) - B)) %*% solve(diag(p) - B)
> data <- rmvDAG(1e4, rDAG, errDist = "normal")
> stopifnot(cov(data)[2,2] > 1.6)
> ##
> 
> set.seed(100)
> 
> wmat <- rbind(c(0,1,0,0,0),
+               c(0,0,0,1,0),
+               c(0,0,0,1,0),
+               c(0,0,0,0,1),
+               c(0,0,0,0,0))
> colnames(wmat) <- rownames(wmat) <- c("1","2","3","4","5")
> print.table(wmat, zero.print=".")
  1 2 3 4 5
1 . 1 . . .
2 . . . 1 .
3 . . . 1 .
4 . . . . 1
5 . . . . .
> 
> g <- as(wmat,"graphNEL")
> 
> e.true <- 0
> var.true <- 5
> 
> dat <- rmvDAG(1000,g)
> x5 <- dat[,5]
> 
> ## test mean
> if (t.test(x5,alternative="two.sided")$p.value<0.05) {
+   stop("Test of rmvDAG: Mean not correct!")
+ }
> 
> ## test variance
> if (var.test(x5,rnorm(1000,0,sqrt(5)),ratio=1,
+              alternative="two.sided")$p.value<0.05) {
+   stop("Test of rmvDAG: Variance not correct!")
+ }
> 
> ###----- Check  gmG  generation: ---> ../man/gmG.Rd
> 
> ## Used to generate "gmG"
> set.seed(40)
> p <- 8
> n <- 5000
> ## true DAG:
> vars <- c("Author", "Bar", "Ctrl", "Goal", paste0("V",5:8))
> gGtrue <- randomDAG(p, prob = 0.3, V = vars)
> x <- rmvDAG(n, gGtrue, back.compatible=TRUE)
> 
> data(gmG)
> 
> ## gmG, gmI were produced on 64-bit -> very small difference even in weights:
> stopifnot(all.equal(gGtrue, gmG$g,  tol=6e-16),
+           all.equal(x,      gmG$ x, tol=1e-15))
> 
> ###----- Check  gmI  generation: ---> ../man/gmI.Rd
> 
> ## Used to generate "gmI"
> set.seed(123)
> p <- 7
> n <- 10000
> myDAG <- randomDAG(p, prob = 0.2) ## true DAG
> datI <- rmvDAG(n, myDAG, back.compatible=TRUE)
> 
> data(gmI)
> stopifnot(all.equal(myDAG, gmI$ g, tol=6e-16),# for 32-bit
+           all.equal(datI,  gmI$ x, tol=1e-15))
> 
> 
> proc.time()
   user  system elapsed 
   1.01    0.20    1.21