Calculate Omega matrices
calc_omega(eval, D_l)
| eval | vector of eigenvalues from decomposition of relatedness matrix |
|---|---|
| D_l | vector of length d_size |
list of length 2. First entry in the list is the symmetric matrix OmegaU. Second entry in the list is the symmetric matrix OmegaE.
#> [[1]] #> [,1] [,2] [,3] [,4] [,5] [,6] #> [1,] 0.01814778 0.01848321 0.01883127 0.01919269 0.01956826 0.01995882 #> [2,] 0.01960022 0.01999206 0.02039990 0.02082472 0.02126762 0.02172976 #> [,7] [,8] [,9] [,10] [,11] [,12] #> [1,] 0.02036529 0.02078865 0.02122999 0.02169048 0.02217139 0.02267411 #> [2,] 0.02221243 0.02271703 0.02324509 0.02379828 0.02437844 0.02498760 #> [,13] [,14] [,15] [,16] [,17] [,18] #> [1,] 0.02320015 0.02375118 0.02432903 0.02493569 0.02557338 0.02624454 #> [2,] 0.02562798 0.02630205 0.02701254 0.02776247 0.02855524 0.02939461 #> [,19] [,20] [,21] [,22] [,23] [,24] #> [1,] 0.02695188 0.02769840 0.02848746 0.02932279 0.03020859 0.03114957 #> [2,] 0.03028482 0.03123063 0.03223743 0.03331130 0.03445918 0.03568899 #> [,25] [,26] [,27] [,28] [,29] [,30] #> [1,] 0.03215107 0.0332191 0.03436052 0.03558318 0.03689606 0.03830953 #> [2,] 0.03700983 0.0384322 0.03996827 0.04163224 0.04344078 0.04541358 #> [,31] [,32] [,33] [,34] [,35] [,36] #> [1,] 0.03983561 0.04148832 0.04328411 0.04524238 0.04738625 0.04974341 #> [2,] 0.04757409 0.04995043 0.05257666 0.05549437 0.05875494 0.06242257 #> [,37] [,38] [,39] [,40] [,41] [,42] #> [1,] 0.05234734 0.05523895 0.05846870 0.06209958 0.06621128 0.07090606 #> [2,] 0.06657858 0.07132746 0.07680582 0.08319574 0.09074537 0.09980193 #> [,43] [,44] [,45] [,46] [,47] [,48] [,49] #> [1,] 0.07631743 0.0826230 0.09006439 0.09897886 0.1098519 0.1234085 0.1407823 #> [2,] 0.11086663 0.1246907 0.14245326 0.16611719 0.1992093 0.2487657 0.3311427 #> [,50] #> [1,] 0.1638493 #> [2,] 0.4950871 #> #> [[2]] #> [,1] [,2] [,3] [,4] [,5] [,6] [,7] #> [1,] 0.9073889 0.9056771 0.9039009 0.9020565 0.9001399 0.8981468 0.8960725 #> [2,] 0.9800108 0.9796112 0.9791952 0.9787620 0.9783103 0.9778390 0.9773467 #> [,8] [,9] [,10] [,11] [,12] [,13] [,14] #> [1,] 0.8939120 0.8916598 0.8893098 0.8868557 0.8842902 0.8816057 0.8787937 #> [2,] 0.9768321 0.9762936 0.9757294 0.9751377 0.9745165 0.9738634 0.9731759 #> [,15] [,16] [,17] [,18] [,19] [,20] [,21] #> [1,] 0.8758449 0.8727490 0.8694948 0.8660697 0.8624601 0.8586504 0.8546237 #> [2,] 0.9724514 0.9716865 0.9708780 0.9700220 0.9691141 0.9681495 0.9671228 #> [,22] [,23] [,24] [,25] [,26] [,27] [,28] #> [1,] 0.8503609 0.8458405 0.8410385 0.8359277 0.8304774 0.8246525 0.8184131 #> [2,] 0.9660276 0.9648569 0.9636027 0.9622556 0.9608051 0.9592385 0.9575415 #> [,29] [,30] [,31] [,32] [,33] [,34] [,35] #> [1,] 0.8117132 0.8045001 0.7967122 0.7882781 0.7791140 0.7691205 0.758180 #> [2,] 0.9556971 0.9536851 0.9514817 0.9490582 0.9463799 0.9434043 0.940079 #> [,36] [,37] [,38] [,39] [,40] [,41] [,42] #> [1,] 0.7461511 0.7328628 0.7181064 0.7016244 0.6830954 0.6621128 0.6381545 #> [2,] 0.9363386 0.9321001 0.9272569 0.9216698 0.9151531 0.9074537 0.8982173 #> [,43] [,44] [,45] [,46] [,47] [,48] [,49] #> [1,] 0.6105394 0.5783610 0.5403863 0.4948943 0.4394075 0.3702256 0.2815645 #> [2,] 0.8869330 0.8728347 0.8547195 0.8305860 0.7968371 0.7462971 0.6622853 #> [,50] #> [1,] 0.1638493 #> [2,] 0.4950871 #>