calc_omega: Calculate Omega matrices

In gemma2: GEMMA Multivariate Linear Mixed Model

Description

Calculate Omega matrices

Usage

calc_omega(eval, D_l)

Arguments

eval	vector of eigenvalues from decomposition of relatedness matrix
D_l	vector of length d_size

Value

list of length 2. First entry in the list is the symmetric matrix OmegaU. Second entry in the list is the symmetric matrix OmegaE.

Examples

calc_omega(eval = 50:1, D_l = runif(2))

Example output

[[1]]
           [,1]       [,2]       [,3]       [,4]       [,5]       [,6]
[1,] 0.01952614 0.01991500 0.02031967 0.02074112 0.02118042 0.02163874
[2,] 0.01950279 0.01989072 0.02029439 0.02071478 0.02115296 0.02161008
           [,7]       [,8]       [,9]      [,10]      [,11]      [,12]
[1,] 0.02211733 0.02261757 0.02314097 0.02368916 0.02426395 0.02486733
[2,] 0.02208739 0.02258626 0.02310819 0.02365481 0.02422791 0.02482948
          [,13]      [,14]      [,15]      [,16]      [,17]      [,18]
[1,] 0.02550148 0.02616883 0.02687203 0.02761408 0.02839827 0.02922831
[2,] 0.02546168 0.02612692 0.02682784 0.02756742 0.02834893 0.02917604
          [,19]      [,20]      [,21]      [,22]      [,23]      [,24]
[1,] 0.03010832 0.03104297 0.03203751 0.03309789 0.03423086 0.03544414
[2,] 0.03005286 0.03098402 0.03197472 0.03303088 0.03415918 0.03536730
          [,25]      [,26]      [,27]      [,28]      [,29]      [,30]
[1,] 0.03674659 0.03814842 0.03966144 0.04129943 0.04307855 0.04501786
[2,] 0.03666401 0.03805942 0.03956525 0.04119515 0.04296510 0.04489397
          [,31]      [,32]      [,33]      [,34]      [,35]      [,36]
[1,] 0.04714000 0.04947211 0.05204699 0.05490461 0.05809425 0.06167735
[2,] 0.04700418 0.04932254 0.05188147 0.05472045 0.05788811 0.06144505
          [,37]      [,38]      [,39]      [,40]      [,41]      [,42]
[1,] 0.06573149 0.07035610 0.07568070 0.08187723 0.08917895 0.09791051
[2,] 0.06546771 0.07005399 0.07533124 0.08146835 0.08869412 0.09732640
         [,43]     [,44]     [,45]     [,46]     [,47]     [,48]    [,49]
[1,] 0.1085375 0.1217521 0.1386307 0.1609423 0.1918131 0.2373375 0.311196
[2,] 0.1078201 0.1208502 0.1374626 0.1593700 0.1895841 0.2339343 0.305371
         [,50]
[1,] 0.4517919
[2,] 0.4396175

[[2]]
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.9763068 0.9758350 0.9753439 0.9748326 0.9742995 0.9737434 0.9731626
[2,] 0.9751397 0.9746452 0.9741306 0.9735948 0.9730362 0.9724535 0.9718451
          [,8]      [,9]     [,10]     [,11]     [,12]     [,13]     [,14]
[1,] 0.9725556 0.9719206 0.9712554 0.9705579 0.9698258 0.9690563 0.9682465
[2,] 0.9712092 0.9705439 0.9698471 0.9691166 0.9683497 0.9675439 0.9666959
         [,15]     [,16]     [,17]     [,18]     [,19]     [,20]     [,21]
[1,] 0.9673932 0.9664928 0.9655413 0.9645341 0.9634663 0.9623322 0.9611254
[2,] 0.9658024 0.9648597 0.9638635 0.9628092 0.9616915 0.9605045 0.9592417
         [,22]     [,23]     [,24]     [,25]     [,26]     [,27]     [,28]
[1,] 0.9598387 0.9584640 0.9569918 0.9554114 0.9537104 0.9518745 0.9498869
[2,] 0.9578954 0.9564571 0.9549171 0.9532642 0.9514855 0.9495660 0.9474884
         [,29]     [,30]     [,31]     [,32]     [,33]     [,34]     [,35]
[1,] 0.9477281 0.9453750 0.9427999 0.9399701 0.9368457 0.9333783 0.9295080
[2,] 0.9452322 0.9427734 0.9400835 0.9371283 0.9338664 0.9302476 0.9262097
         [,36]     [,37]     [,38]     [,39]     [,40]     [,41]     [,42]
[1,] 0.9251602 0.9202409 0.9146293 0.9081684 0.9006495 0.8917895 0.8811946
[2,] 0.9216757 0.9165480 0.9107018 0.9039749 0.8961519 0.8869412 0.8759376
         [,43]     [,44]     [,45]     [,46]     [,47]     [,48]     [,49]
[1,] 0.8682997 0.8522649 0.8317843 0.8047113 0.7672524 0.7120126 0.6223921
[2,] 0.8625611 0.8459516 0.8247757 0.7968502 0.7583363 0.7018028 0.6107421
         [,50]
[1,] 0.4517919
[2,] 0.4396175

gemma2 documentation built on Oct. 24, 2020, 5:06 p.m.