morie_chol_lower
Examples
set.seed(1)
n <- 12
A <- matrix(rnorm(n * 6), nrow = n)
G <- morie_grm(A) + diag(0.3, n)
rmorie:::morie_chol_lower(G)
#> [,1] [,2] [,3] [,4] [,5] [,6]
#> [1,] 1.15201549 0.00000000 0.00000000 0.00000000 0.000000000 0.00000000
#> [2,] 0.12838618 1.04034457 0.00000000 0.00000000 0.000000000 0.00000000
#> [3,] -0.01283583 -0.36341511 1.01846129 0.00000000 0.000000000 0.00000000
#> [4,] -0.63896371 -0.08711234 -0.20288008 1.03282532 0.000000000 0.00000000
#> [5,] -0.22919210 0.13341694 -0.14076774 -0.07552165 0.655886686 0.00000000
#> [6,] 0.35277750 -0.39106219 0.06254605 -0.26199216 0.137970198 0.90340303
#> [7,] -0.46184761 0.04658816 0.22680111 -0.65419087 -0.068579321 -0.12767445
#> [8,] -0.07904839 -0.02070600 0.21022231 -0.06330491 -0.197044817 -0.63644010
#> [9,] 0.12297173 -0.24743122 -0.42165046 -0.23933193 0.100865298 0.21665318
#> [10,] 0.67486673 -0.43406532 -0.21240411 0.01563658 0.042589084 0.31570296
#> [11,] -0.49790276 0.03005544 -0.21694855 0.48962580 -0.002048928 -0.35316306
#> [12,] -0.25081406 0.54961617 0.06589384 0.30590258 -0.026528804 0.05705384
#> [,7] [,8] [,9] [,10] [,11] [,12]
#> [1,] 0.00000000 0.000000000 0.000000000 0.00000000 0.00000000 0.0000000
#> [2,] 0.00000000 0.000000000 0.000000000 0.00000000 0.00000000 0.0000000
#> [3,] 0.00000000 0.000000000 0.000000000 0.00000000 0.00000000 0.0000000
#> [4,] 0.00000000 0.000000000 0.000000000 0.00000000 0.00000000 0.0000000
#> [5,] 0.00000000 0.000000000 0.000000000 0.00000000 0.00000000 0.0000000
#> [6,] 0.00000000 0.000000000 0.000000000 0.00000000 0.00000000 0.0000000
#> [7,] 1.09272048 0.000000000 0.000000000 0.00000000 0.00000000 0.0000000
#> [8,] 0.09633511 0.834576279 0.000000000 0.00000000 0.00000000 0.0000000
#> [9,] 0.06378849 -0.026367093 0.706601299 0.00000000 0.00000000 0.0000000
#> [10,] -0.39348624 -0.007635104 0.191664941 0.67989568 0.00000000 0.0000000
#> [11,] -0.20176392 0.190378384 0.009946306 -0.01193821 0.65337673 0.0000000
#> [12,] 0.04859924 -0.300125091 -0.265698462 -0.19196973 -0.06950237 0.7085324