Get bootstrap probabilities
Examples
LiNGAM_sample_1000 <- generate_lingam_sample_6()
bs_model <- lingam_direct_bootstrap(LiNGAM_sample_1000$data,
n_sampling = 30L, reg_method = "ols", seed = 42
)
#> Bootstrap: 30 iterations, method=ols (sequential)
#> iteration 1 / 30
#> iteration 10 / 30
#> iteration 20 / 30
#> iteration 30 / 30
#> Completed in 0.2 seconds.
get_probabilities(bs_model)
#> [,1] [,2] [,3] [,4] [,5] [,6]
#> [1,] 0.00000000 0.03333333 0.53333333 0.9666667 0.03333333 0.00000000
#> [2,] 0.96666667 0.00000000 0.96666667 0.9666667 0.60000000 0.63333333
#> [3,] 0.46666667 0.03333333 0.00000000 0.9666667 0.03333333 0.43333333
#> [4,] 0.03333333 0.03333333 0.03333333 0.0000000 0.03333333 0.03333333
#> [5,] 0.96666667 0.40000000 0.96666667 0.9666667 0.00000000 0.56666667
#> [6,] 1.00000000 0.36666667 0.56666667 0.9666667 0.43333333 0.00000000
