# \donttest{
# Data simulation
set.seed(1)
simul <- SimulateGraphical(pk = 20)
# Stability selection
stab <- GraphicalModel(xdata = simul$data)
# Extracting calibrated hyper-parameters
Argmax(stab)
# Extracting calibrated hyper-parameters IDs
ids <- ArgmaxId(stab)
ids
# Relationship between the two functions
stab$Lambda[ids[1], 1]
stab$params$pi_list[ids[2]]
# }
Run the code above in your browser using DataLab