R Under development (unstable) (2023-08-30 r85030 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2023 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(GGMselect)
Loading required package: mvtnorm
> p=30
> n=30
> eta=0.13
> dmax=3
> iG = 4
> iS = 20
> set.seed(iG)
> Gr <- simulateGraph(p,eta)
> set.seed(iS*(pi/3.1415)**iG)
> X <- rmvnorm(n, mean=rep(0,p), sigma=Gr$C)
> 
> K=2.5
> mon1GRestALL <- selectFast(X, dmax, K, family="LA")
> resFam <- selectMyFam(X, list(mon1GRestALL$LA$G), K)
> 
> print(all.equal(mon1GRestALL$LA$G, resFam$G))
[1] TRUE
> print(all.equal(mon1GRestALL$LA$Neighb[,1, drop=F], resFam$Neighb))
[1] TRUE
> print(all.equal(mon1GRestALL$LA$crit.min, resFam$crit.min))
[1] TRUE
> 
> K=c(2.5, 1)
> zz=selectFast(X, dmax, K, family="LA")
> resFam <- selectMyFam(X, list(zz$LA$G[,,1], zz$LA$G[,,2]), K=K)
> 
> print(all.equal(zz$LA$G, resFam$G))
[1] TRUE
> print(all.equal(as.vector(zz$LA$Neighb[,1,1]), as.vector(resFam$Neighb[,,1])))
[1] "Numeric: lengths (30, 60) differ"
> print(all.equal(zz$LA$crit.min, resFam$crit.min))
[1] TRUE
> 
> proc.time()
   user  system elapsed 
   0.39    0.10    0.43