Build mr.mash prior covariance matrices
Source:R/regularizedRegressionWrappers.R
buildMrmashPriorMatrices.RdShared helper used by both mrmashWrapper (individual-level) and
mrmashRssWeights (summary statistics). Constructs the S0
list of prior covariance matrices via the canonical mixture
(mr.mashr::compute_canonical_covs) and optional data-driven matrices,
expanded over a scaling grid via mr.mashr::expand_covs. The prior grid
is derived from Bhat and Shat via computeGrid when not
supplied.
Usage
buildMrmashPriorMatrices(
Bhat,
Shat,
K = NULL,
dataDrivenPriorMatrices = NULL,
canonicalPriorMatrices = TRUE,
priorGrid = NULL,
canonicalCovs = MrmashCanonicalCovsOptions(),
expandCovs = MrmashExpandCovsOptions()
)Arguments
- Bhat
Numeric matrix of effect-size estimates (variants x conditions).
- Shat
Numeric matrix of standard errors (variants x conditions).
- K
Number of conditions. When NULL, inferred from
ncol(Bhat).- dataDrivenPriorMatrices
Optional list with element
U(list of raw covariance matrices) computed e.g. bycomputeCovDiag, or bymashr::cov_flash, which returns such a list directly.- canonicalPriorMatrices
Logical. When TRUE (default for RSS), include the standard canonical mixture from
mr.mashr::compute_canonical_covs(). When FALSE,dataDrivenPriorMatricesmust be supplied.- priorGrid
Optional pre-computed scaling grid (numeric vector). When NULL, derived from
Bhat,ShatviacomputeGrid().- canonicalCovs
Arguments for
mr.mashr::compute_canonical_covs()—singletons,hetgrid— built withMrmashCanonicalCovsOptions. Empty by default, which leaves that function's own defaults in place.- expandCovs
Arguments for
mr.mashr::expand_covs()—zeromat— built withMrmashExpandCovsOptions. Empty by default, which leaves that function's own defaults in place.
Value
A list with components S0 (the expanded list of prior
covariance matrices) and prior_grid (the scaling grid that was
used).
Examples
Bhat <- matrix(rnorm(15), 5, 3)
Shat <- matrix(abs(rnorm(15)) + 0.1, 5, 3)
buildMrmashPriorMatrices(Bhat = Bhat, Shat = Shat)
#> $S0
#> $S0$null
#> [,1] [,2] [,3]
#> [1,] 0 0 0
#> [2,] 0 0 0
#> [3,] 0 0 0
#>
#> $S0$singleton1_grid1
#> [,1] [,2] [,3]
#> [1,] 0.0001357805 0 0
#> [2,] 0.0000000000 0 0
#> [3,] 0.0000000000 0 0
#>
#> $S0$singleton1_grid2
#> [,1] [,2] [,3]
#> [1,] 0.0002715609 0 0
#> [2,] 0.0000000000 0 0
#> [3,] 0.0000000000 0 0
#>
#> $S0$singleton1_grid3
#> [,1] [,2] [,3]
#> [1,] 0.0005431219 0 0
#> [2,] 0.0000000000 0 0
#> [3,] 0.0000000000 0 0
#>
#> $S0$singleton1_grid4
#> [,1] [,2] [,3]
#> [1,] 0.001086244 0 0
#> [2,] 0.000000000 0 0
#> [3,] 0.000000000 0 0
#>
#> $S0$singleton1_grid5
#> [,1] [,2] [,3]
#> [1,] 0.002172487 0 0
#> [2,] 0.000000000 0 0
#> [3,] 0.000000000 0 0
#>
#> $S0$singleton1_grid6
#> [,1] [,2] [,3]
#> [1,] 0.004344975 0 0
#> [2,] 0.000000000 0 0
#> [3,] 0.000000000 0 0
#>
#> $S0$singleton1_grid7
#> [,1] [,2] [,3]
#> [1,] 0.00868995 0 0
#> [2,] 0.00000000 0 0
#> [3,] 0.00000000 0 0
#>
#> $S0$singleton1_grid8
#> [,1] [,2] [,3]
#> [1,] 0.0173799 0 0
#> [2,] 0.0000000 0 0
#> [3,] 0.0000000 0 0
#>
#> $S0$singleton1_grid9
#> [,1] [,2] [,3]
#> [1,] 0.0347598 0 0
#> [2,] 0.0000000 0 0
#> [3,] 0.0000000 0 0
#>
#> $S0$singleton1_grid10
#> [,1] [,2] [,3]
#> [1,] 0.0695196 0 0
#> [2,] 0.0000000 0 0
#> [3,] 0.0000000 0 0
#>
#> $S0$singleton1_grid11
#> [,1] [,2] [,3]
#> [1,] 0.1390392 0 0
#> [2,] 0.0000000 0 0
#> [3,] 0.0000000 0 0
#>
#> $S0$singleton1_grid12
#> [,1] [,2] [,3]
#> [1,] 0.2780784 0 0
#> [2,] 0.0000000 0 0
#> [3,] 0.0000000 0 0
#>
#> $S0$singleton1_grid13
#> [,1] [,2] [,3]
#> [1,] 0.5561568 0 0
#> [2,] 0.0000000 0 0
#> [3,] 0.0000000 0 0
#>
#> $S0$singleton1_grid14
#> [,1] [,2] [,3]
#> [1,] 1.112314 0 0
#> [2,] 0.000000 0 0
#> [3,] 0.000000 0 0
#>
#> $S0$singleton1_grid15
#> [,1] [,2] [,3]
#> [1,] 2.224627 0 0
#> [2,] 0.000000 0 0
#> [3,] 0.000000 0 0
#>
#> $S0$singleton1_grid16
#> [,1] [,2] [,3]
#> [1,] 4.449254 0 0
#> [2,] 0.000000 0 0
#> [3,] 0.000000 0 0
#>
#> $S0$singleton1_grid17
#> [,1] [,2] [,3]
#> [1,] 8.898509 0 0
#> [2,] 0.000000 0 0
#> [3,] 0.000000 0 0
#>
#> $S0$singleton1_grid18
#> [,1] [,2] [,3]
#> [1,] 17.79702 0 0
#> [2,] 0.00000 0 0
#> [3,] 0.00000 0 0
#>
#> $S0$singleton2_grid1
#> [,1] [,2] [,3]
#> [1,] 0 0.0000000000 0
#> [2,] 0 0.0001357805 0
#> [3,] 0 0.0000000000 0
#>
#> $S0$singleton2_grid2
#> [,1] [,2] [,3]
#> [1,] 0 0.0000000000 0
#> [2,] 0 0.0002715609 0
#> [3,] 0 0.0000000000 0
#>
#> $S0$singleton2_grid3
#> [,1] [,2] [,3]
#> [1,] 0 0.0000000000 0
#> [2,] 0 0.0005431219 0
#> [3,] 0 0.0000000000 0
#>
#> $S0$singleton2_grid4
#> [,1] [,2] [,3]
#> [1,] 0 0.000000000 0
#> [2,] 0 0.001086244 0
#> [3,] 0 0.000000000 0
#>
#> $S0$singleton2_grid5
#> [,1] [,2] [,3]
#> [1,] 0 0.000000000 0
#> [2,] 0 0.002172487 0
#> [3,] 0 0.000000000 0
#>
#> $S0$singleton2_grid6
#> [,1] [,2] [,3]
#> [1,] 0 0.000000000 0
#> [2,] 0 0.004344975 0
#> [3,] 0 0.000000000 0
#>
#> $S0$singleton2_grid7
#> [,1] [,2] [,3]
#> [1,] 0 0.00000000 0
#> [2,] 0 0.00868995 0
#> [3,] 0 0.00000000 0
#>
#> $S0$singleton2_grid8
#> [,1] [,2] [,3]
#> [1,] 0 0.0000000 0
#> [2,] 0 0.0173799 0
#> [3,] 0 0.0000000 0
#>
#> $S0$singleton2_grid9
#> [,1] [,2] [,3]
#> [1,] 0 0.0000000 0
#> [2,] 0 0.0347598 0
#> [3,] 0 0.0000000 0
#>
#> $S0$singleton2_grid10
#> [,1] [,2] [,3]
#> [1,] 0 0.0000000 0
#> [2,] 0 0.0695196 0
#> [3,] 0 0.0000000 0
#>
#> $S0$singleton2_grid11
#> [,1] [,2] [,3]
#> [1,] 0 0.0000000 0
#> [2,] 0 0.1390392 0
#> [3,] 0 0.0000000 0
#>
#> $S0$singleton2_grid12
#> [,1] [,2] [,3]
#> [1,] 0 0.0000000 0
#> [2,] 0 0.2780784 0
#> [3,] 0 0.0000000 0
#>
#> $S0$singleton2_grid13
#> [,1] [,2] [,3]
#> [1,] 0 0.0000000 0
#> [2,] 0 0.5561568 0
#> [3,] 0 0.0000000 0
#>
#> $S0$singleton2_grid14
#> [,1] [,2] [,3]
#> [1,] 0 0.000000 0
#> [2,] 0 1.112314 0
#> [3,] 0 0.000000 0
#>
#> $S0$singleton2_grid15
#> [,1] [,2] [,3]
#> [1,] 0 0.000000 0
#> [2,] 0 2.224627 0
#> [3,] 0 0.000000 0
#>
#> $S0$singleton2_grid16
#> [,1] [,2] [,3]
#> [1,] 0 0.000000 0
#> [2,] 0 4.449254 0
#> [3,] 0 0.000000 0
#>
#> $S0$singleton2_grid17
#> [,1] [,2] [,3]
#> [1,] 0 0.000000 0
#> [2,] 0 8.898509 0
#> [3,] 0 0.000000 0
#>
#> $S0$singleton2_grid18
#> [,1] [,2] [,3]
#> [1,] 0 0.00000 0
#> [2,] 0 17.79702 0
#> [3,] 0 0.00000 0
#>
#> $S0$singleton3_grid1
#> [,1] [,2] [,3]
#> [1,] 0 0 0.0000000000
#> [2,] 0 0 0.0000000000
#> [3,] 0 0 0.0001357805
#>
#> $S0$singleton3_grid2
#> [,1] [,2] [,3]
#> [1,] 0 0 0.0000000000
#> [2,] 0 0 0.0000000000
#> [3,] 0 0 0.0002715609
#>
#> $S0$singleton3_grid3
#> [,1] [,2] [,3]
#> [1,] 0 0 0.0000000000
#> [2,] 0 0 0.0000000000
#> [3,] 0 0 0.0005431219
#>
#> $S0$singleton3_grid4
#> [,1] [,2] [,3]
#> [1,] 0 0 0.000000000
#> [2,] 0 0 0.000000000
#> [3,] 0 0 0.001086244
#>
#> $S0$singleton3_grid5
#> [,1] [,2] [,3]
#> [1,] 0 0 0.000000000
#> [2,] 0 0 0.000000000
#> [3,] 0 0 0.002172487
#>
#> $S0$singleton3_grid6
#> [,1] [,2] [,3]
#> [1,] 0 0 0.000000000
#> [2,] 0 0 0.000000000
#> [3,] 0 0 0.004344975
#>
#> $S0$singleton3_grid7
#> [,1] [,2] [,3]
#> [1,] 0 0 0.00000000
#> [2,] 0 0 0.00000000
#> [3,] 0 0 0.00868995
#>
#> $S0$singleton3_grid8
#> [,1] [,2] [,3]
#> [1,] 0 0 0.0000000
#> [2,] 0 0 0.0000000
#> [3,] 0 0 0.0173799
#>
#> $S0$singleton3_grid9
#> [,1] [,2] [,3]
#> [1,] 0 0 0.0000000
#> [2,] 0 0 0.0000000
#> [3,] 0 0 0.0347598
#>
#> $S0$singleton3_grid10
#> [,1] [,2] [,3]
#> [1,] 0 0 0.0000000
#> [2,] 0 0 0.0000000
#> [3,] 0 0 0.0695196
#>
#> $S0$singleton3_grid11
#> [,1] [,2] [,3]
#> [1,] 0 0 0.0000000
#> [2,] 0 0 0.0000000
#> [3,] 0 0 0.1390392
#>
#> $S0$singleton3_grid12
#> [,1] [,2] [,3]
#> [1,] 0 0 0.0000000
#> [2,] 0 0 0.0000000
#> [3,] 0 0 0.2780784
#>
#> $S0$singleton3_grid13
#> [,1] [,2] [,3]
#> [1,] 0 0 0.0000000
#> [2,] 0 0 0.0000000
#> [3,] 0 0 0.5561568
#>
#> $S0$singleton3_grid14
#> [,1] [,2] [,3]
#> [1,] 0 0 0.000000
#> [2,] 0 0 0.000000
#> [3,] 0 0 1.112314
#>
#> $S0$singleton3_grid15
#> [,1] [,2] [,3]
#> [1,] 0 0 0.000000
#> [2,] 0 0 0.000000
#> [3,] 0 0 2.224627
#>
#> $S0$singleton3_grid16
#> [,1] [,2] [,3]
#> [1,] 0 0 0.000000
#> [2,] 0 0 0.000000
#> [3,] 0 0 4.449254
#>
#> $S0$singleton3_grid17
#> [,1] [,2] [,3]
#> [1,] 0 0 0.000000
#> [2,] 0 0 0.000000
#> [3,] 0 0 8.898509
#>
#> $S0$singleton3_grid18
#> [,1] [,2] [,3]
#> [1,] 0 0 0.00000
#> [2,] 0 0 0.00000
#> [3,] 0 0 17.79702
#>
#> $S0$independent_grid1
#> [,1] [,2] [,3]
#> [1,] 0.0001357805 0.0000000000 0.0000000000
#> [2,] 0.0000000000 0.0001357805 0.0000000000
#> [3,] 0.0000000000 0.0000000000 0.0001357805
#>
#> $S0$independent_grid2
#> [,1] [,2] [,3]
#> [1,] 0.0002715609 0.0000000000 0.0000000000
#> [2,] 0.0000000000 0.0002715609 0.0000000000
#> [3,] 0.0000000000 0.0000000000 0.0002715609
#>
#> $S0$independent_grid3
#> [,1] [,2] [,3]
#> [1,] 0.0005431219 0.0000000000 0.0000000000
#> [2,] 0.0000000000 0.0005431219 0.0000000000
#> [3,] 0.0000000000 0.0000000000 0.0005431219
#>
#> $S0$independent_grid4
#> [,1] [,2] [,3]
#> [1,] 0.001086244 0.000000000 0.000000000
#> [2,] 0.000000000 0.001086244 0.000000000
#> [3,] 0.000000000 0.000000000 0.001086244
#>
#> $S0$independent_grid5
#> [,1] [,2] [,3]
#> [1,] 0.002172487 0.000000000 0.000000000
#> [2,] 0.000000000 0.002172487 0.000000000
#> [3,] 0.000000000 0.000000000 0.002172487
#>
#> $S0$independent_grid6
#> [,1] [,2] [,3]
#> [1,] 0.004344975 0.000000000 0.000000000
#> [2,] 0.000000000 0.004344975 0.000000000
#> [3,] 0.000000000 0.000000000 0.004344975
#>
#> $S0$independent_grid7
#> [,1] [,2] [,3]
#> [1,] 0.00868995 0.00000000 0.00000000
#> [2,] 0.00000000 0.00868995 0.00000000
#> [3,] 0.00000000 0.00000000 0.00868995
#>
#> $S0$independent_grid8
#> [,1] [,2] [,3]
#> [1,] 0.0173799 0.0000000 0.0000000
#> [2,] 0.0000000 0.0173799 0.0000000
#> [3,] 0.0000000 0.0000000 0.0173799
#>
#> $S0$independent_grid9
#> [,1] [,2] [,3]
#> [1,] 0.0347598 0.0000000 0.0000000
#> [2,] 0.0000000 0.0347598 0.0000000
#> [3,] 0.0000000 0.0000000 0.0347598
#>
#> $S0$independent_grid10
#> [,1] [,2] [,3]
#> [1,] 0.0695196 0.0000000 0.0000000
#> [2,] 0.0000000 0.0695196 0.0000000
#> [3,] 0.0000000 0.0000000 0.0695196
#>
#> $S0$independent_grid11
#> [,1] [,2] [,3]
#> [1,] 0.1390392 0.0000000 0.0000000
#> [2,] 0.0000000 0.1390392 0.0000000
#> [3,] 0.0000000 0.0000000 0.1390392
#>
#> $S0$independent_grid12
#> [,1] [,2] [,3]
#> [1,] 0.2780784 0.0000000 0.0000000
#> [2,] 0.0000000 0.2780784 0.0000000
#> [3,] 0.0000000 0.0000000 0.2780784
#>
#> $S0$independent_grid13
#> [,1] [,2] [,3]
#> [1,] 0.5561568 0.0000000 0.0000000
#> [2,] 0.0000000 0.5561568 0.0000000
#> [3,] 0.0000000 0.0000000 0.5561568
#>
#> $S0$independent_grid14
#> [,1] [,2] [,3]
#> [1,] 1.112314 0.000000 0.000000
#> [2,] 0.000000 1.112314 0.000000
#> [3,] 0.000000 0.000000 1.112314
#>
#> $S0$independent_grid15
#> [,1] [,2] [,3]
#> [1,] 2.224627 0.000000 0.000000
#> [2,] 0.000000 2.224627 0.000000
#> [3,] 0.000000 0.000000 2.224627
#>
#> $S0$independent_grid16
#> [,1] [,2] [,3]
#> [1,] 4.449254 0.000000 0.000000
#> [2,] 0.000000 4.449254 0.000000
#> [3,] 0.000000 0.000000 4.449254
#>
#> $S0$independent_grid17
#> [,1] [,2] [,3]
#> [1,] 8.898509 0.000000 0.000000
#> [2,] 0.000000 8.898509 0.000000
#> [3,] 0.000000 0.000000 8.898509
#>
#> $S0$independent_grid18
#> [,1] [,2] [,3]
#> [1,] 17.79702 0.00000 0.00000
#> [2,] 0.00000 17.79702 0.00000
#> [3,] 0.00000 0.00000 17.79702
#>
#> $S0$shared0.25_grid1
#> [,1] [,2] [,3]
#> [1,] 1.357805e-04 3.394512e-05 3.394512e-05
#> [2,] 3.394512e-05 1.357805e-04 3.394512e-05
#> [3,] 3.394512e-05 3.394512e-05 1.357805e-04
#>
#> $S0$shared0.25_grid2
#> [,1] [,2] [,3]
#> [1,] 2.715609e-04 6.789023e-05 6.789023e-05
#> [2,] 6.789023e-05 2.715609e-04 6.789023e-05
#> [3,] 6.789023e-05 6.789023e-05 2.715609e-04
#>
#> $S0$shared0.25_grid3
#> [,1] [,2] [,3]
#> [1,] 0.0005431219 0.0001357805 0.0001357805
#> [2,] 0.0001357805 0.0005431219 0.0001357805
#> [3,] 0.0001357805 0.0001357805 0.0005431219
#>
#> $S0$shared0.25_grid4
#> [,1] [,2] [,3]
#> [1,] 0.0010862437 0.0002715609 0.0002715609
#> [2,] 0.0002715609 0.0010862437 0.0002715609
#> [3,] 0.0002715609 0.0002715609 0.0010862437
#>
#> $S0$shared0.25_grid5
#> [,1] [,2] [,3]
#> [1,] 0.0021724875 0.0005431219 0.0005431219
#> [2,] 0.0005431219 0.0021724875 0.0005431219
#> [3,] 0.0005431219 0.0005431219 0.0021724875
#>
#> $S0$shared0.25_grid6
#> [,1] [,2] [,3]
#> [1,] 0.004344975 0.001086244 0.001086244
#> [2,] 0.001086244 0.004344975 0.001086244
#> [3,] 0.001086244 0.001086244 0.004344975
#>
#> $S0$shared0.25_grid7
#> [,1] [,2] [,3]
#> [1,] 0.008689950 0.002172487 0.002172487
#> [2,] 0.002172487 0.008689950 0.002172487
#> [3,] 0.002172487 0.002172487 0.008689950
#>
#> $S0$shared0.25_grid8
#> [,1] [,2] [,3]
#> [1,] 0.017379900 0.004344975 0.004344975
#> [2,] 0.004344975 0.017379900 0.004344975
#> [3,] 0.004344975 0.004344975 0.017379900
#>
#> $S0$shared0.25_grid9
#> [,1] [,2] [,3]
#> [1,] 0.03475980 0.00868995 0.00868995
#> [2,] 0.00868995 0.03475980 0.00868995
#> [3,] 0.00868995 0.00868995 0.03475980
#>
#> $S0$shared0.25_grid10
#> [,1] [,2] [,3]
#> [1,] 0.0695196 0.0173799 0.0173799
#> [2,] 0.0173799 0.0695196 0.0173799
#> [3,] 0.0173799 0.0173799 0.0695196
#>
#> $S0$shared0.25_grid11
#> [,1] [,2] [,3]
#> [1,] 0.1390392 0.0347598 0.0347598
#> [2,] 0.0347598 0.1390392 0.0347598
#> [3,] 0.0347598 0.0347598 0.1390392
#>
#> $S0$shared0.25_grid12
#> [,1] [,2] [,3]
#> [1,] 0.2780784 0.0695196 0.0695196
#> [2,] 0.0695196 0.2780784 0.0695196
#> [3,] 0.0695196 0.0695196 0.2780784
#>
#> $S0$shared0.25_grid13
#> [,1] [,2] [,3]
#> [1,] 0.5561568 0.1390392 0.1390392
#> [2,] 0.1390392 0.5561568 0.1390392
#> [3,] 0.1390392 0.1390392 0.5561568
#>
#> $S0$shared0.25_grid14
#> [,1] [,2] [,3]
#> [1,] 1.1123136 0.2780784 0.2780784
#> [2,] 0.2780784 1.1123136 0.2780784
#> [3,] 0.2780784 0.2780784 1.1123136
#>
#> $S0$shared0.25_grid15
#> [,1] [,2] [,3]
#> [1,] 2.2246272 0.5561568 0.5561568
#> [2,] 0.5561568 2.2246272 0.5561568
#> [3,] 0.5561568 0.5561568 2.2246272
#>
#> $S0$shared0.25_grid16
#> [,1] [,2] [,3]
#> [1,] 4.449254 1.112314 1.112314
#> [2,] 1.112314 4.449254 1.112314
#> [3,] 1.112314 1.112314 4.449254
#>
#> $S0$shared0.25_grid17
#> [,1] [,2] [,3]
#> [1,] 8.898509 2.224627 2.224627
#> [2,] 2.224627 8.898509 2.224627
#> [3,] 2.224627 2.224627 8.898509
#>
#> $S0$shared0.25_grid18
#> [,1] [,2] [,3]
#> [1,] 17.797017 4.449254 4.449254
#> [2,] 4.449254 17.797017 4.449254
#> [3,] 4.449254 4.449254 17.797017
#>
#> $S0$shared0.5_grid1
#> [,1] [,2] [,3]
#> [1,] 1.357805e-04 6.789023e-05 6.789023e-05
#> [2,] 6.789023e-05 1.357805e-04 6.789023e-05
#> [3,] 6.789023e-05 6.789023e-05 1.357805e-04
#>
#> $S0$shared0.5_grid2
#> [,1] [,2] [,3]
#> [1,] 0.0002715609 0.0001357805 0.0001357805
#> [2,] 0.0001357805 0.0002715609 0.0001357805
#> [3,] 0.0001357805 0.0001357805 0.0002715609
#>
#> $S0$shared0.5_grid3
#> [,1] [,2] [,3]
#> [1,] 0.0005431219 0.0002715609 0.0002715609
#> [2,] 0.0002715609 0.0005431219 0.0002715609
#> [3,] 0.0002715609 0.0002715609 0.0005431219
#>
#> $S0$shared0.5_grid4
#> [,1] [,2] [,3]
#> [1,] 0.0010862437 0.0005431219 0.0005431219
#> [2,] 0.0005431219 0.0010862437 0.0005431219
#> [3,] 0.0005431219 0.0005431219 0.0010862437
#>
#> $S0$shared0.5_grid5
#> [,1] [,2] [,3]
#> [1,] 0.002172487 0.001086244 0.001086244
#> [2,] 0.001086244 0.002172487 0.001086244
#> [3,] 0.001086244 0.001086244 0.002172487
#>
#> $S0$shared0.5_grid6
#> [,1] [,2] [,3]
#> [1,] 0.004344975 0.002172487 0.002172487
#> [2,] 0.002172487 0.004344975 0.002172487
#> [3,] 0.002172487 0.002172487 0.004344975
#>
#> $S0$shared0.5_grid7
#> [,1] [,2] [,3]
#> [1,] 0.008689950 0.004344975 0.004344975
#> [2,] 0.004344975 0.008689950 0.004344975
#> [3,] 0.004344975 0.004344975 0.008689950
#>
#> $S0$shared0.5_grid8
#> [,1] [,2] [,3]
#> [1,] 0.01737990 0.00868995 0.00868995
#> [2,] 0.00868995 0.01737990 0.00868995
#> [3,] 0.00868995 0.00868995 0.01737990
#>
#> $S0$shared0.5_grid9
#> [,1] [,2] [,3]
#> [1,] 0.0347598 0.0173799 0.0173799
#> [2,] 0.0173799 0.0347598 0.0173799
#> [3,] 0.0173799 0.0173799 0.0347598
#>
#> $S0$shared0.5_grid10
#> [,1] [,2] [,3]
#> [1,] 0.0695196 0.0347598 0.0347598
#> [2,] 0.0347598 0.0695196 0.0347598
#> [3,] 0.0347598 0.0347598 0.0695196
#>
#> $S0$shared0.5_grid11
#> [,1] [,2] [,3]
#> [1,] 0.1390392 0.0695196 0.0695196
#> [2,] 0.0695196 0.1390392 0.0695196
#> [3,] 0.0695196 0.0695196 0.1390392
#>
#> $S0$shared0.5_grid12
#> [,1] [,2] [,3]
#> [1,] 0.2780784 0.1390392 0.1390392
#> [2,] 0.1390392 0.2780784 0.1390392
#> [3,] 0.1390392 0.1390392 0.2780784
#>
#> $S0$shared0.5_grid13
#> [,1] [,2] [,3]
#> [1,] 0.5561568 0.2780784 0.2780784
#> [2,] 0.2780784 0.5561568 0.2780784
#> [3,] 0.2780784 0.2780784 0.5561568
#>
#> $S0$shared0.5_grid14
#> [,1] [,2] [,3]
#> [1,] 1.1123136 0.5561568 0.5561568
#> [2,] 0.5561568 1.1123136 0.5561568
#> [3,] 0.5561568 0.5561568 1.1123136
#>
#> $S0$shared0.5_grid15
#> [,1] [,2] [,3]
#> [1,] 2.224627 1.112314 1.112314
#> [2,] 1.112314 2.224627 1.112314
#> [3,] 1.112314 1.112314 2.224627
#>
#> $S0$shared0.5_grid16
#> [,1] [,2] [,3]
#> [1,] 4.449254 2.224627 2.224627
#> [2,] 2.224627 4.449254 2.224627
#> [3,] 2.224627 2.224627 4.449254
#>
#> $S0$shared0.5_grid17
#> [,1] [,2] [,3]
#> [1,] 8.898509 4.449254 4.449254
#> [2,] 4.449254 8.898509 4.449254
#> [3,] 4.449254 4.449254 8.898509
#>
#> $S0$shared0.5_grid18
#> [,1] [,2] [,3]
#> [1,] 17.797017 8.898509 8.898509
#> [2,] 8.898509 17.797017 8.898509
#> [3,] 8.898509 8.898509 17.797017
#>
#> $S0$shared0.75_grid1
#> [,1] [,2] [,3]
#> [1,] 0.0001357805 0.0001018353 0.0001018353
#> [2,] 0.0001018353 0.0001357805 0.0001018353
#> [3,] 0.0001018353 0.0001018353 0.0001357805
#>
#> $S0$shared0.75_grid2
#> [,1] [,2] [,3]
#> [1,] 0.0002715609 0.0002036707 0.0002036707
#> [2,] 0.0002036707 0.0002715609 0.0002036707
#> [3,] 0.0002036707 0.0002036707 0.0002715609
#>
#> $S0$shared0.75_grid3
#> [,1] [,2] [,3]
#> [1,] 0.0005431219 0.0004073414 0.0004073414
#> [2,] 0.0004073414 0.0005431219 0.0004073414
#> [3,] 0.0004073414 0.0004073414 0.0005431219
#>
#> $S0$shared0.75_grid4
#> [,1] [,2] [,3]
#> [1,] 0.0010862437 0.0008146828 0.0008146828
#> [2,] 0.0008146828 0.0010862437 0.0008146828
#> [3,] 0.0008146828 0.0008146828 0.0010862437
#>
#> $S0$shared0.75_grid5
#> [,1] [,2] [,3]
#> [1,] 0.002172487 0.001629366 0.001629366
#> [2,] 0.001629366 0.002172487 0.001629366
#> [3,] 0.001629366 0.001629366 0.002172487
#>
#> $S0$shared0.75_grid6
#> [,1] [,2] [,3]
#> [1,] 0.004344975 0.003258731 0.003258731
#> [2,] 0.003258731 0.004344975 0.003258731
#> [3,] 0.003258731 0.003258731 0.004344975
#>
#> $S0$shared0.75_grid7
#> [,1] [,2] [,3]
#> [1,] 0.008689950 0.006517462 0.006517462
#> [2,] 0.006517462 0.008689950 0.006517462
#> [3,] 0.006517462 0.006517462 0.008689950
#>
#> $S0$shared0.75_grid8
#> [,1] [,2] [,3]
#> [1,] 0.01737990 0.01303492 0.01303492
#> [2,] 0.01303492 0.01737990 0.01303492
#> [3,] 0.01303492 0.01303492 0.01737990
#>
#> $S0$shared0.75_grid9
#> [,1] [,2] [,3]
#> [1,] 0.03475980 0.02606985 0.02606985
#> [2,] 0.02606985 0.03475980 0.02606985
#> [3,] 0.02606985 0.02606985 0.03475980
#>
#> $S0$shared0.75_grid10
#> [,1] [,2] [,3]
#> [1,] 0.0695196 0.0521397 0.0521397
#> [2,] 0.0521397 0.0695196 0.0521397
#> [3,] 0.0521397 0.0521397 0.0695196
#>
#> $S0$shared0.75_grid11
#> [,1] [,2] [,3]
#> [1,] 0.1390392 0.1042794 0.1042794
#> [2,] 0.1042794 0.1390392 0.1042794
#> [3,] 0.1042794 0.1042794 0.1390392
#>
#> $S0$shared0.75_grid12
#> [,1] [,2] [,3]
#> [1,] 0.2780784 0.2085588 0.2085588
#> [2,] 0.2085588 0.2780784 0.2085588
#> [3,] 0.2085588 0.2085588 0.2780784
#>
#> $S0$shared0.75_grid13
#> [,1] [,2] [,3]
#> [1,] 0.5561568 0.4171176 0.4171176
#> [2,] 0.4171176 0.5561568 0.4171176
#> [3,] 0.4171176 0.4171176 0.5561568
#>
#> $S0$shared0.75_grid14
#> [,1] [,2] [,3]
#> [1,] 1.1123136 0.8342352 0.8342352
#> [2,] 0.8342352 1.1123136 0.8342352
#> [3,] 0.8342352 0.8342352 1.1123136
#>
#> $S0$shared0.75_grid15
#> [,1] [,2] [,3]
#> [1,] 2.224627 1.668470 1.668470
#> [2,] 1.668470 2.224627 1.668470
#> [3,] 1.668470 1.668470 2.224627
#>
#> $S0$shared0.75_grid16
#> [,1] [,2] [,3]
#> [1,] 4.449254 3.336941 3.336941
#> [2,] 3.336941 4.449254 3.336941
#> [3,] 3.336941 3.336941 4.449254
#>
#> $S0$shared0.75_grid17
#> [,1] [,2] [,3]
#> [1,] 8.898509 6.673881 6.673881
#> [2,] 6.673881 8.898509 6.673881
#> [3,] 6.673881 6.673881 8.898509
#>
#> $S0$shared0.75_grid18
#> [,1] [,2] [,3]
#> [1,] 17.79702 13.34776 13.34776
#> [2,] 13.34776 17.79702 13.34776
#> [3,] 13.34776 13.34776 17.79702
#>
#> $S0$shared1_grid1
#> [,1] [,2] [,3]
#> [1,] 0.0001357805 0.0001357805 0.0001357805
#> [2,] 0.0001357805 0.0001357805 0.0001357805
#> [3,] 0.0001357805 0.0001357805 0.0001357805
#>
#> $S0$shared1_grid2
#> [,1] [,2] [,3]
#> [1,] 0.0002715609 0.0002715609 0.0002715609
#> [2,] 0.0002715609 0.0002715609 0.0002715609
#> [3,] 0.0002715609 0.0002715609 0.0002715609
#>
#> $S0$shared1_grid3
#> [,1] [,2] [,3]
#> [1,] 0.0005431219 0.0005431219 0.0005431219
#> [2,] 0.0005431219 0.0005431219 0.0005431219
#> [3,] 0.0005431219 0.0005431219 0.0005431219
#>
#> $S0$shared1_grid4
#> [,1] [,2] [,3]
#> [1,] 0.001086244 0.001086244 0.001086244
#> [2,] 0.001086244 0.001086244 0.001086244
#> [3,] 0.001086244 0.001086244 0.001086244
#>
#> $S0$shared1_grid5
#> [,1] [,2] [,3]
#> [1,] 0.002172487 0.002172487 0.002172487
#> [2,] 0.002172487 0.002172487 0.002172487
#> [3,] 0.002172487 0.002172487 0.002172487
#>
#> $S0$shared1_grid6
#> [,1] [,2] [,3]
#> [1,] 0.004344975 0.004344975 0.004344975
#> [2,] 0.004344975 0.004344975 0.004344975
#> [3,] 0.004344975 0.004344975 0.004344975
#>
#> $S0$shared1_grid7
#> [,1] [,2] [,3]
#> [1,] 0.00868995 0.00868995 0.00868995
#> [2,] 0.00868995 0.00868995 0.00868995
#> [3,] 0.00868995 0.00868995 0.00868995
#>
#> $S0$shared1_grid8
#> [,1] [,2] [,3]
#> [1,] 0.0173799 0.0173799 0.0173799
#> [2,] 0.0173799 0.0173799 0.0173799
#> [3,] 0.0173799 0.0173799 0.0173799
#>
#> $S0$shared1_grid9
#> [,1] [,2] [,3]
#> [1,] 0.0347598 0.0347598 0.0347598
#> [2,] 0.0347598 0.0347598 0.0347598
#> [3,] 0.0347598 0.0347598 0.0347598
#>
#> $S0$shared1_grid10
#> [,1] [,2] [,3]
#> [1,] 0.0695196 0.0695196 0.0695196
#> [2,] 0.0695196 0.0695196 0.0695196
#> [3,] 0.0695196 0.0695196 0.0695196
#>
#> $S0$shared1_grid11
#> [,1] [,2] [,3]
#> [1,] 0.1390392 0.1390392 0.1390392
#> [2,] 0.1390392 0.1390392 0.1390392
#> [3,] 0.1390392 0.1390392 0.1390392
#>
#> $S0$shared1_grid12
#> [,1] [,2] [,3]
#> [1,] 0.2780784 0.2780784 0.2780784
#> [2,] 0.2780784 0.2780784 0.2780784
#> [3,] 0.2780784 0.2780784 0.2780784
#>
#> $S0$shared1_grid13
#> [,1] [,2] [,3]
#> [1,] 0.5561568 0.5561568 0.5561568
#> [2,] 0.5561568 0.5561568 0.5561568
#> [3,] 0.5561568 0.5561568 0.5561568
#>
#> $S0$shared1_grid14
#> [,1] [,2] [,3]
#> [1,] 1.112314 1.112314 1.112314
#> [2,] 1.112314 1.112314 1.112314
#> [3,] 1.112314 1.112314 1.112314
#>
#> $S0$shared1_grid15
#> [,1] [,2] [,3]
#> [1,] 2.224627 2.224627 2.224627
#> [2,] 2.224627 2.224627 2.224627
#> [3,] 2.224627 2.224627 2.224627
#>
#> $S0$shared1_grid16
#> [,1] [,2] [,3]
#> [1,] 4.449254 4.449254 4.449254
#> [2,] 4.449254 4.449254 4.449254
#> [3,] 4.449254 4.449254 4.449254
#>
#> $S0$shared1_grid17
#> [,1] [,2] [,3]
#> [1,] 8.898509 8.898509 8.898509
#> [2,] 8.898509 8.898509 8.898509
#> [3,] 8.898509 8.898509 8.898509
#>
#> $S0$shared1_grid18
#> [,1] [,2] [,3]
#> [1,] 17.79702 17.79702 17.79702
#> [2,] 17.79702 17.79702 17.79702
#> [3,] 17.79702 17.79702 17.79702
#>
#>
#> $priorGrid
#> [1] 1.357805e-04 2.715609e-04 5.431219e-04 1.086244e-03 2.172487e-03
#> [6] 4.344975e-03 8.689950e-03 1.737990e-02 3.475980e-02 6.951960e-02
#> [11] 1.390392e-01 2.780784e-01 5.561568e-01 1.112314e+00 2.224627e+00
#> [16] 4.449254e+00 8.898509e+00 1.779702e+01
#>