Apply a fitted mashr model (from
mashModelFit) to a target SumStats set, returning the
posterior matrices (PosteriorMean, PosteriorSD, lfsr,
...). This is the posterior step of the mash workflow – mash_fit's
second step and mash_posterior's per-analysis-unit computation.
Usage
mashPosterior(
model,
sumStats,
alpha,
vhat = NULL,
excludeCondition = character(0),
outputPosteriorCov = TRUE,
inputScale = c("auto", "beta", "z"),
mashDataArgs = MashDataOptions(),
posteriorArgs = MashPosteriorOptions()
)Arguments
- model
A fitted mashr model (the
mashModelFitoutput).- sumStats
A single
QtlSumStats/GwasSumStats– the target (e.g. strong) set to compute posteriors on.- alpha
mash
alpha(match the fit).- vhat
Residual correlation matrix (
V);NULL-> identity.- excludeCondition
Character vector of condition (column) names to drop from the target AND the model before computing posteriors (the model's covariances are resized via
updateMashModelCov). Default none.- outputPosteriorCov
Return the full posterior covariance array (needed by
fitMashContrast). DefaultTRUE.- inputScale
SumStats -> matrix conversion scale.
- mashDataArgs
Extra arguments for
mashr::mash_set_data(), built withMashDataOptions– for examplezero_Bhat_Shat_resetorzero_Shat_reset.Bhat,Shat,alphaandVare supplied by pecotmr and are refused by the constructor.- posteriorArgs
Extra arguments for
mashr::mash_compute_posterior_matrices(), built withMashPosteriorOptions– for examplealgorithm.version.g,data,A,output_posterior_covandseedare owned by pecotmr and are refused by the constructor.
Value
The mashr::mash_compute_posterior_matrices() result: a list of
PosteriorMean / PosteriorSD / lfsr /
NegativeProb (+ PosteriorCov when outputPosteriorCov =
TRUE).
Examples
data(mashInputExample)
mi <- mashInputExample
mk <- function(b, s) {
qtlSumStatsFromBetaMatrix(as.matrix(mi[[b]]), as.matrix(mi[[s]]),
studyName = "mash")
}
ssl <- list(strong = mk("strong.b", "strong.s"),
random = mk("random.b", "random.s"), null = mk("null.b", "null.s"))
conds <- colnames(mi$strong.b)
vhat <- diag(length(conds))
dimnames(vhat) <- list(conds, conds)
prior <- mashPriorCovariances(ssl, alpha = 0, vhat = vhat,
components = "canonical")
#> - Computing 17 x 313 likelihood matrix.
#> - Likelihood calculations took 0.01 seconds.
#> - Fitting model with 313 mixture components.
#> - Model fitting took 0.03 seconds.
model <- mashModelFit(ssl, alpha = 0, priorCovariances = prior,
vhat = vhat)
#> - Computing 68 x 287 likelihood matrix.
#> - Likelihood calculations took 0.02 seconds.
#> - Fitting model with 287 mixture components.
#> - Model fitting took 0.37 seconds.
#> - Computing posterior matrices.
#> - Computation allocated took 0.00 seconds.
mashPosterior(model, mk("strong.b", "strong.s"), alpha = 0, vhat = vhat)
#> $PosteriorMean
#> ALL Ast End Exc
#> mash::mash::var1 1.732327e-01 1.699445e-01 1.696735e-01 1.721133e-01
#> mash::mash::var2 1.463980e-02 1.386809e-02 1.405698e-02 1.395979e-02
#> mash::mash::var3 1.993913e-01 1.781980e-01 1.826245e-01 2.011939e-01
#> mash::mash::var4 -1.068118e-01 -1.022474e-01 -1.017032e-01 -1.077036e-01
#> mash::mash::var5 -2.053818e-02 -2.239960e-02 -5.349338e-02 -2.064765e-02
#> mash::mash::var6 4.391819e-17 6.058851e-17 7.635642e-17 1.006284e-17
#> mash::mash::var7 -6.703892e-10 -4.050672e-01 -6.765331e-10 -3.401429e-10
#> mash::mash::var8 -3.236640e-07 -6.304466e-03 -7.239392e-07 -3.114503e-07
#> mash::mash::var9 8.611694e-02 6.118350e-02 7.413988e-02 6.629913e-02
#> mash::mash::var10 -8.198549e-04 -1.359032e-04 -2.015088e-03 -1.326624e-03
#> mash::mash::var11 1.783160e-01 1.650132e-01 1.629840e-01 1.761605e-01
#> mash::mash::var12 -1.018259e-01 -4.743449e-02 -5.689337e-02 -1.022507e-01
#> mash::mash::var13 -1.347660e-01 -1.239028e-01 -9.437686e-02 -1.275104e-01
#> mash::mash::var14 -3.706634e-04 -2.857872e-04 -3.348290e-01 -4.238503e-04
#> mash::mash::var15 5.215931e-02 2.479819e-02 2.744594e-02 6.550792e-02
#> mash::mash::var16 4.299903e-06 -3.672818e-07 1.039914e-05 4.609306e-06
#> mash::mash::var17 1.921340e-01 1.631047e-01 1.705078e-01 1.738614e-01
#> Inh Mic OPC Oli
#> mash::mash::var1 1.709787e-01 1.707482e-01 1.703722e-01 1.708895e-01
#> mash::mash::var2 1.452499e-02 1.450407e-02 1.411593e-02 1.511437e-02
#> mash::mash::var3 1.888205e-01 1.813981e-01 1.843592e-01 1.834557e-01
#> mash::mash::var4 -1.039959e-01 -1.022543e-01 -9.974423e-02 -9.963496e-02
#> mash::mash::var5 -2.026891e-02 -2.041690e-02 -2.037971e-02 -2.044313e-02
#> mash::mash::var6 7.428756e-17 7.552491e-17 1.122940e+00 1.100110e-16
#> mash::mash::var7 -6.388503e-10 -7.015579e-10 -9.248674e-10 -9.107671e-10
#> mash::mash::var8 -6.012371e-07 -3.319319e-07 -4.752720e-01 -7.269175e-07
#> mash::mash::var9 4.395225e-02 1.581745e-01 7.702635e-02 6.974477e-02
#> mash::mash::var10 -4.844039e-04 9.275211e-04 3.166563e-02 -2.910397e-04
#> mash::mash::var11 1.690050e-01 1.648797e-01 1.664258e-01 1.634323e-01
#> mash::mash::var12 -3.937860e-02 -1.616142e-02 -4.399621e-02 -5.227220e-02
#> mash::mash::var13 -1.209970e-01 -7.959069e-02 -1.025180e-01 -9.898423e-02
#> mash::mash::var14 -3.521015e-04 -3.311382e-04 -3.078692e-04 -3.222916e-04
#> mash::mash::var15 1.402996e-02 2.848861e-02 4.082076e-02 4.078157e-02
#> mash::mash::var16 3.190104e-06 2.652281e-06 -4.194478e-01 4.744552e-06
#> mash::mash::var17 1.697109e-01 1.735160e-01 1.727000e-01 2.034895e-01
#>
#> $PosteriorSD
#> ALL Ast End Exc
#> mash::mash::var1 1.502447e-02 1.688453e-02 1.740473e-02 1.471528e-02
#> mash::mash::var2 2.606255e-02 2.949004e-02 2.943047e-02 2.558846e-02
#> mash::mash::var3 1.922490e-02 2.891423e-02 2.582288e-02 1.957110e-02
#> mash::mash::var4 1.702712e-02 1.927235e-02 1.974095e-02 1.709753e-02
#> mash::mash::var5 3.273777e-02 4.277619e-02 1.028692e-01 3.281609e-02
#> mash::mash::var6 1.745953e-09 1.954450e-09 2.218536e-09 1.722333e-09
#> mash::mash::var7 6.488618e-06 4.329820e-02 6.766114e-06 3.581778e-06
#> mash::mash::var8 9.906185e-05 3.778567e-02 2.091828e-04 9.415196e-05
#> mash::mash::var9 1.886598e-02 2.139762e-02 2.255807e-02 1.981090e-02
#> mash::mash::var10 1.280080e-02 1.753385e-02 2.300350e-02 1.304867e-02
#> mash::mash::var11 1.593113e-02 2.117360e-02 2.324226e-02 1.538093e-02
#> mash::mash::var12 1.346505e-02 1.974596e-02 2.101680e-02 1.208822e-02
#> mash::mash::var13 1.740902e-02 2.126724e-02 2.342758e-02 1.684578e-02
#> mash::mash::var14 3.980383e-03 3.347386e-03 7.265223e-02 4.504785e-03
#> mash::mash::var15 1.349154e-02 1.893949e-02 2.148635e-02 1.365707e-02
#> mash::mash::var16 5.367714e-04 6.649222e-04 1.211284e-03 5.415268e-04
#> mash::mash::var17 2.338643e-02 3.342581e-02 3.069182e-02 2.919324e-02
#> Inh Mic OPC Oli
#> mash::mash::var1 1.575017e-02 1.601107e-02 1.649348e-02 1.594440e-02
#> mash::mash::var2 2.610719e-02 2.613024e-02 3.192826e-02 2.661704e-02
#> mash::mash::var3 2.117362e-02 2.643011e-02 2.454707e-02 2.495738e-02
#> mash::mash::var4 1.851143e-02 1.915370e-02 2.102297e-02 2.087942e-02
#> mash::mash::var5 3.262486e-02 3.272224e-02 3.824855e-02 3.273288e-02
#> mash::mash::var6 1.986741e-09 1.962058e-09 1.293032e-01 2.442393e-09
#> mash::mash::var7 6.347557e-06 6.984117e-06 8.992471e-06 8.825135e-06
#> mash::mash::var8 1.716949e-04 1.086054e-04 9.633369e-02 2.087332e-04
#> mash::mash::var9 2.222979e-02 2.042649e-02 2.193413e-02 2.282935e-02
#> mash::mash::var10 1.293109e-02 1.362759e-02 1.134100e-01 1.305095e-02
#> mash::mash::var11 1.832016e-02 2.150174e-02 2.041350e-02 2.267848e-02
#> mash::mash::var12 1.619940e-02 1.481524e-02 2.065651e-02 2.056812e-02
#> mash::mash::var13 1.892982e-02 2.582343e-02 2.274536e-02 2.256814e-02
#> mash::mash::var14 3.888967e-03 3.789987e-03 3.575900e-03 3.716320e-03
#> mash::mash::var15 1.937477e-02 2.100275e-02 2.027656e-02 1.916823e-02
#> mash::mash::var16 5.568730e-04 5.696012e-04 9.650108e-02 5.689820e-04
#> mash::mash::var17 3.080708e-02 2.933231e-02 2.936241e-02 2.639379e-02
#>
#> $lfdr
#> ALL Ast End Exc
#> mash::mash::var1 1.538797e-31 1.488460e-31 1.525647e-31 1.538797e-31
#> mash::mash::var2 6.088163e-01 6.015143e-01 5.987778e-01 6.088163e-01
#> mash::mash::var3 2.156670e-30 2.143070e-30 2.133419e-30 2.156670e-30
#> mash::mash::var4 5.756681e-10 5.663776e-10 5.687998e-10 5.756681e-10
#> mash::mash::var5 5.719734e-01 5.613362e-01 4.665297e-01 5.719734e-01
#> mash::mash::var6 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00
#> mash::mash::var7 1.000000e+00 4.715089e-16 1.000000e+00 1.000000e+00
#> mash::mash::var8 9.999861e-01 9.725467e-01 9.999861e-01 9.999861e-01
#> mash::mash::var9 8.418545e-14 8.387032e-14 8.340594e-14 8.418545e-14
#> mash::mash::var10 7.889133e-01 7.816736e-01 7.743506e-01 7.889133e-01
#> mash::mash::var11 7.735517e-37 7.257203e-37 7.679827e-37 7.735517e-37
#> mash::mash::var12 1.013588e-23 1.010251e-23 1.000088e-23 1.013588e-23
#> mash::mash::var13 1.216794e-18 4.690133e-22 1.216789e-18 1.216794e-18
#> mash::mash::var14 9.901091e-01 9.901084e-01 3.291260e-04 9.901091e-01
#> mash::mash::var15 1.023198e-06 1.020683e-06 1.006869e-06 1.023198e-06
#> mash::mash::var16 9.992206e-01 9.991948e-01 9.991462e-01 9.992206e-01
#> mash::mash::var17 1.726319e-17 1.709965e-17 1.693644e-17 1.726319e-17
#> Inh Mic OPC Oli
#> mash::mash::var1 1.538797e-31 1.538797e-31 1.494386e-31 1.538797e-31
#> mash::mash::var2 6.088163e-01 6.088163e-01 6.002560e-01 6.088163e-01
#> mash::mash::var3 2.156670e-30 2.156670e-30 2.087089e-30 2.156670e-30
#> mash::mash::var4 5.756681e-10 5.756681e-10 5.711012e-10 5.756681e-10
#> mash::mash::var5 5.719734e-01 5.719734e-01 5.644135e-01 5.719734e-01
#> mash::mash::var6 1.000000e+00 1.000000e+00 4.510376e-15 1.000000e+00
#> mash::mash::var7 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00
#> mash::mash::var8 9.999861e-01 9.999861e-01 2.743942e-02 9.999861e-01
#> mash::mash::var9 8.418545e-14 8.418545e-14 8.312190e-14 8.418545e-14
#> mash::mash::var10 7.889133e-01 7.889133e-01 7.006890e-01 7.889133e-01
#> mash::mash::var11 7.735517e-37 7.735517e-37 7.290720e-37 7.735517e-37
#> mash::mash::var12 1.013588e-23 1.013588e-23 1.005940e-23 1.013588e-23
#> mash::mash::var13 1.216794e-18 1.216794e-18 1.216791e-18 1.216794e-18
#> mash::mash::var14 9.901091e-01 9.901091e-01 9.901078e-01 9.901091e-01
#> mash::mash::var15 1.023198e-06 1.023198e-06 1.009930e-06 1.023198e-06
#> mash::mash::var16 9.992206e-01 9.992206e-01 4.213281e-03 9.992206e-01
#> mash::mash::var17 1.726319e-17 1.726319e-17 1.686670e-17 1.726319e-17
#>
#> $NegativeProb
#> ALL Ast End Exc
#> mash::mash::var1 6.400173e-23 2.427943e-12 4.069359e-11 1.931396e-21
#> mash::mash::var2 4.068437e-02 4.960566e-02 5.042974e-02 4.577236e-02
#> mash::mash::var3 6.060366e-22 4.266370e-11 8.942074e-12 2.028720e-23
#> mash::mash::var4 1.000000e+00 9.999860e-01 9.999732e-01 1.000000e+00
#> mash::mash::var5 3.928936e-01 4.006612e-01 4.993679e-01 3.934585e-01
#> mash::mash::var6 1.187824e-15 9.808549e-16 8.587202e-16 1.680910e-15
#> mash::mash::var7 1.132368e-08 1.000000e+00 1.128971e-08 1.105312e-08
#> mash::mash::var8 1.343194e-05 2.745334e-02 1.382569e-05 1.353260e-05
#> mash::mash::var9 2.501499e-06 2.122463e-03 5.069906e-04 4.090461e-04
#> mash::mash::var10 1.171562e-01 1.133249e-01 1.250280e-01 1.238671e-01
#> mash::mash::var11 1.087129e-26 3.435832e-12 9.377734e-11 6.814873e-24
#> mash::mash::var12 1.000000e+00 9.918524e-01 9.966057e-01 1.000000e+00
#> mash::mash::var13 1.000000e+00 1.000000e+00 9.999829e-01 1.000000e+00
#> mash::mash::var14 9.823207e-03 9.245603e-03 9.996409e-01 9.878036e-03
#> mash::mash::var15 5.925879e-05 9.652142e-02 1.036959e-01 2.337597e-07
#> mash::mash::var16 2.982419e-04 3.807507e-04 3.222087e-04 2.928254e-04
#> mash::mash::var17 9.884720e-16 1.780601e-08 6.247141e-09 1.335607e-09
#> Inh Mic OPC Oli
#> mash::mash::var1 6.120759e-14 3.491757e-13 4.487130e-12 5.417589e-13
#> mash::mash::var2 4.230731e-02 4.268923e-02 4.930001e-02 3.861405e-02
#> mash::mash::var3 5.351605e-16 5.869271e-12 1.700391e-12 1.356685e-12
#> mash::mash::var4 9.999959e-01 9.999889e-01 9.999229e-01 9.999444e-01
#> mash::mash::var5 3.909899e-01 3.918189e-01 3.950984e-01 3.920063e-01
#> mash::mash::var6 8.023148e-16 7.642709e-16 3.702226e-16 4.957889e-16
#> mash::mash::var7 1.130359e-08 1.129994e-08 1.132303e-08 1.132357e-08
#> mash::mash::var8 1.387073e-05 1.285098e-05 9.725606e-01 1.385167e-05
#> mash::mash::var9 2.401017e-02 5.286388e-15 2.226192e-04 1.125129e-03
#> mash::mash::var10 1.125437e-01 9.467875e-02 1.066051e-01 1.100024e-01
#> mash::mash::var11 5.568018e-14 1.282698e-11 6.305963e-12 2.868757e-11
#> mash::mash::var12 9.924687e-01 8.623340e-01 9.834099e-01 9.944798e-01
#> mash::mash::var13 1.000000e+00 9.996390e-01 9.999965e-01 9.999956e-01
#> mash::mash::var14 9.661778e-03 9.415711e-03 9.299908e-03 9.357904e-03
#> mash::mash::var15 2.367165e-01 9.030640e-02 2.380802e-02 1.816198e-02
#> mash::mash::var16 3.206680e-04 3.301492e-04 9.954149e-01 2.945906e-04
#> mash::mash::var17 5.241641e-09 1.545712e-09 1.128256e-09 1.015330e-15
#>
#> $lfsr
#> ALL Ast End Exc
#> mash::mash::var1 6.400173e-23 2.427943e-12 4.069359e-11 1.931396e-21
#> mash::mash::var2 6.495007e-01 6.511200e-01 6.492075e-01 6.545886e-01
#> mash::mash::var3 6.060366e-22 4.266370e-11 8.942074e-12 2.028720e-23
#> mash::mash::var4 1.044163e-08 1.400697e-05 2.683731e-05 3.350299e-09
#> mash::mash::var5 6.071064e-01 5.993388e-01 5.006321e-01 6.065415e-01
#> mash::mash::var6 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00
#> mash::mash::var7 1.000000e+00 4.440892e-16 1.000000e+00 1.000000e+00
#> mash::mash::var8 9.999866e-01 9.725467e-01 9.999862e-01 9.999865e-01
#> mash::mash::var9 2.501499e-06 2.122463e-03 5.069906e-04 4.090461e-04
#> mash::mash::var10 8.828438e-01 8.866751e-01 8.749720e-01 8.761329e-01
#> mash::mash::var11 1.087129e-26 3.435832e-12 9.377734e-11 6.814873e-24
#> mash::mash::var12 1.965095e-14 8.147553e-03 3.394277e-03 0.000000e+00
#> mash::mash::var13 3.330669e-15 7.524208e-09 1.708834e-05 4.474199e-14
#> mash::mash::var14 9.901768e-01 9.907544e-01 3.591089e-04 9.901220e-01
#> mash::mash::var15 6.028199e-05 9.652244e-02 1.036969e-01 1.256957e-06
#> mash::mash::var16 9.995189e-01 9.995755e-01 9.994684e-01 9.995134e-01
#> mash::mash::var17 1.005735e-15 1.780601e-08 6.247141e-09 1.335607e-09
#> Inh Mic OPC Oli
#> mash::mash::var1 6.120759e-14 3.491757e-13 4.487130e-12 5.417589e-13
#> mash::mash::var2 6.511236e-01 6.515055e-01 6.495560e-01 6.474303e-01
#> mash::mash::var3 5.351605e-16 5.869271e-12 1.700391e-12 1.356685e-12
#> mash::mash::var4 4.053404e-06 1.110639e-05 7.713362e-05 5.556849e-05
#> mash::mash::var5 6.090101e-01 6.081811e-01 6.049016e-01 6.079937e-01
#> mash::mash::var6 1.000000e+00 1.000000e+00 4.880599e-15 1.000000e+00
#> mash::mash::var7 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00
#> mash::mash::var8 9.999861e-01 9.999871e-01 2.743942e-02 9.999861e-01
#> mash::mash::var9 2.401017e-02 8.947184e-14 2.226192e-04 1.125129e-03
#> mash::mash::var10 8.874563e-01 8.835921e-01 8.072940e-01 8.899976e-01
#> mash::mash::var11 5.568018e-14 1.282698e-11 6.305963e-12 2.868757e-11
#> mash::mash::var12 7.531342e-03 1.376660e-01 1.659011e-02 5.520197e-03
#> mash::mash::var13 2.371752e-10 3.610442e-04 3.452905e-06 4.436665e-06
#> mash::mash::var14 9.903382e-01 9.905843e-01 9.907001e-01 9.906421e-01
#> mash::mash::var15 2.367175e-01 9.030743e-02 2.380903e-02 1.816301e-02
#> mash::mash::var16 9.995413e-01 9.995508e-01 4.585106e-03 9.995152e-01
#> mash::mash::var17 5.241641e-09 1.545712e-09 1.128256e-09 1.032594e-15
#>
#> $PosteriorCov
#> , , mash::mash::var1
#>
#> ALL Ast End Exc Inh
#> ALL 0.0002257348 0.0001770547 0.0001757787 0.0001920376 0.0001857486
#> Ast 0.0001770547 0.0002850875 0.0002622572 0.0002016080 0.0002314367
#> End 0.0001757787 0.0002622572 0.0003029247 0.0002031491 0.0002361957
#> Exc 0.0001920376 0.0002016080 0.0002031491 0.0002165395 0.0001993931
#> Inh 0.0001857486 0.0002314367 0.0002361957 0.0001993931 0.0002480678
#> Mic 0.0001841904 0.0002371367 0.0002425230 0.0002002500 0.0002215644
#> OPC 0.0001816364 0.0002464225 0.0002528313 0.0002016371 0.0002271756
#> Oli 0.0001857918 0.0002344683 0.0002395306 0.0002003441 0.0002202686
#> Mic OPC Oli
#> ALL 0.0001841904 0.0001816364 0.0001857918
#> Ast 0.0002371367 0.0002464225 0.0002344683
#> End 0.0002425230 0.0002528313 0.0002395306
#> Exc 0.0002002500 0.0002016371 0.0002003441
#> Inh 0.0002215644 0.0002271756 0.0002202686
#> Mic 0.0002563545 0.0002319806 0.0002239399
#> OPC 0.0002319806 0.0002720350 0.0002299150
#> Oli 0.0002239399 0.0002299150 0.0002542240
#>
#> , , mash::mash::var2
#>
#> ALL Ast End Exc Inh
#> ALL 0.0006792566 0.0006379467 0.0006376914 0.0006261379 0.0006404640
#> Ast 0.0006379467 0.0008696626 0.0006370003 0.0006262528 0.0006388081
#> End 0.0006376914 0.0006370003 0.0008661525 0.0006257951 0.0006385382
#> Exc 0.0006261379 0.0006262528 0.0006257951 0.0006547691 0.0006270826
#> Inh 0.0006404640 0.0006388081 0.0006385382 0.0006270826 0.0006815853
#> Mic 0.0006408105 0.0006392594 0.0006389889 0.0006275360 0.0006414099
#> OPC 0.0006372808 0.0006365637 0.0006362370 0.0006253306 0.0006381144
#> Oli 0.0006529519 0.0006493394 0.0006492476 0.0006365208 0.0006530079
#> Mic OPC Oli
#> ALL 0.0006408105 0.0006372808 0.0006529519
#> Ast 0.0006392594 0.0006365637 0.0006493394
#> End 0.0006389889 0.0006362370 0.0006492476
#> Exc 0.0006275360 0.0006253306 0.0006365208
#> Inh 0.0006414099 0.0006381144 0.0006530079
#> Mic 0.0006827895 0.0006385623 0.0006533118
#> OPC 0.0006385623 0.0010194135 0.0006488700
#> Oli 0.0006533118 0.0006488700 0.0007084670
#>
#> , , mash::mash::var3
#>
#> ALL Ast End Exc Inh
#> ALL 0.0003695966 0.0001733287 0.0001958887 0.0002548082 0.0002098062
#> Ast 0.0001733287 0.0008360328 0.0005922983 0.0001285342 0.0004307384
#> End 0.0001958887 0.0005922983 0.0006668210 0.0001615323 0.0003918824
#> Exc 0.0002548082 0.0001285342 0.0001615323 0.0003830280 0.0001904576
#> Inh 0.0002098062 0.0004307384 0.0003918824 0.0001904576 0.0004483222
#> Mic 0.0001871157 0.0006171327 0.0005353755 0.0001499239 0.0003999350
#> OPC 0.0002016938 0.0005493337 0.0004849505 0.0001714961 0.0003733904
#> Oli 0.0001955853 0.0005680489 0.0004982868 0.0001632909 0.0003797035
#> Mic OPC Oli
#> ALL 0.0001871157 0.0002016938 0.0001955853
#> Ast 0.0006171327 0.0005493337 0.0005680489
#> End 0.0005353755 0.0004849505 0.0004982868
#> Exc 0.0001499239 0.0001714961 0.0001632909
#> Inh 0.0003999350 0.0003733904 0.0003797035
#> Mic 0.0006985508 0.0004997269 0.0005146258
#> OPC 0.0004997269 0.0006025588 0.0004673503
#> Oli 0.0005146258 0.0004673503 0.0006228708
#>
#> , , mash::mash::var4
#>
#> ALL Ast End Exc Inh
#> ALL 0.0002899227 0.0002190907 0.0002197172 0.0002168970 0.0002216277
#> Ast 0.0002190907 0.0003714234 0.0002849519 0.0002094841 0.0002599094
#> End 0.0002197172 0.0002849519 0.0003897051 0.0002088595 0.0002648043
#> Exc 0.0002168970 0.0002094841 0.0002088595 0.0002923256 0.0002159605
#> Inh 0.0002216277 0.0002599094 0.0002648043 0.0002159605 0.0003426729
#> Mic 0.0002180747 0.0002762645 0.0002835297 0.0002085018 0.0002585809
#> OPC 0.0002160385 0.0003037515 0.0003145515 0.0002007822 0.0002749743
#> Oli 0.0002136448 0.0003020217 0.0003128869 0.0001981820 0.0002727931
#> Mic OPC Oli
#> ALL 0.0002180747 0.0002160385 0.0002136448
#> Ast 0.0002762645 0.0003037515 0.0003020217
#> End 0.0002835297 0.0003145515 0.0003128869
#> Exc 0.0002085018 0.0002007822 0.0001981820
#> Inh 0.0002585809 0.0002749743 0.0002727931
#> Mic 0.0003668641 0.0003023002 0.0003006233
#> OPC 0.0003023002 0.0004419654 0.0003445873
#> Oli 0.0003006233 0.0003445873 0.0004359500
#>
#> , , mash::mash::var5
#>
#> ALL Ast End Exc Inh
#> ALL 0.0010717614 0.0009912916 0.0003751202 0.0010334042 0.001026680
#> Ast 0.0009912916 0.0018298027 0.0002738809 0.0009935054 0.000987733
#> End 0.0003751202 0.0002738809 0.0105820678 0.0003744390 0.000378682
#> Exc 0.0010334042 0.0009935054 0.0003744390 0.0010768957 0.001028941
#> Inh 0.0010266800 0.0009877330 0.0003786820 0.0010289410 0.001064381
#> Mic 0.0010300231 0.0009907149 0.0003777502 0.0010324069 0.001026048
#> OPC 0.0010265757 0.0009873989 0.0003751157 0.0010288681 0.001022826
#> Oli 0.0010304237 0.0009910487 0.0003773905 0.0010328290 0.001026393
#> Mic OPC Oli
#> ALL 0.0010300231 0.0010265757 0.0010304237
#> Ast 0.0009907149 0.0009873989 0.0009910487
#> End 0.0003777502 0.0003751157 0.0003773905
#> Exc 0.0010324069 0.0010288681 0.0010328290
#> Inh 0.0010260477 0.0010228257 0.0010263932
#> Mic 0.0010707450 0.0010259138 0.0010296339
#> OPC 0.0010259138 0.0014629518 0.0010262662
#> Oli 0.0010296339 0.0010262662 0.0010714413
#>
#> , , mash::mash::var6
#>
#> ALL Ast End Exc Inh
#> ALL 3.048351e-18 2.347427e-18 2.492419e-18 1.889244e-18 2.433495e-18
#> Ast 2.347427e-18 3.819876e-18 2.881930e-18 1.857758e-18 2.825442e-18
#> End 2.492419e-18 2.881930e-18 4.921900e-18 1.857732e-18 3.101938e-18
#> Exc 1.889244e-18 1.857758e-18 1.857732e-18 2.966432e-18 1.788545e-18
#> Inh 2.433495e-18 2.825442e-18 3.101938e-18 1.788545e-18 3.947141e-18
#> Mic 2.391774e-18 2.790367e-18 3.070225e-18 1.729675e-18 3.021721e-18
#> OPC -4.628917e-17 -6.422233e-17 -8.139281e-17 -9.654809e-18 -7.911124e-17
#> Oli 2.791554e-18 3.420627e-18 3.853016e-18 1.703122e-18 3.807192e-18
#> Mic OPC Oli
#> ALL 2.391774e-18 -4.628917e-17 2.791554e-18
#> Ast 2.790367e-18 -6.422233e-17 3.420627e-18
#> End 3.070225e-18 -8.139281e-17 3.853016e-18
#> Exc 1.729675e-18 -9.654809e-18 1.703122e-18
#> Inh 3.021721e-18 -7.911124e-17 3.807192e-18
#> Mic 3.849673e-18 -8.050504e-17 3.793298e-18
#> OPC -8.050504e-17 1.671932e-02 -1.176857e-16
#> Oli 3.793298e-18 -1.176857e-16 5.965284e-18
#>
#> , , mash::mash::var7
#>
#> ALL Ast End Exc Inh
#> ALL 4.210216e-11 -1.900362e-10 4.091089e-11 2.061900e-11 3.853805e-11
#> Ast -1.900362e-10 1.874734e-03 -1.910806e-10 -9.599858e-11 -1.806020e-10
#> End 4.091089e-11 -1.910806e-10 4.578031e-11 2.123711e-11 3.953013e-11
#> Exc 2.061900e-11 -9.599858e-11 2.123711e-11 1.282913e-11 1.995468e-11
#> Inh 3.853805e-11 -1.806020e-10 3.953013e-11 1.995468e-11 4.029148e-11
#> Mic 4.238650e-11 -1.982127e-10 4.353519e-11 2.198199e-11 4.093253e-11
#> OPC 5.560787e-11 -2.617633e-10 5.687879e-11 2.868615e-11 5.353214e-11
#> Oli 5.468542e-11 -2.579022e-10 5.586970e-11 2.816944e-11 5.259757e-11
#> Mic OPC Oli
#> ALL 4.238650e-11 5.560787e-11 5.468542e-11
#> Ast -1.982127e-10 -2.617633e-10 -2.579022e-10
#> End 4.353519e-11 5.687879e-11 5.586970e-11
#> Exc 2.198199e-11 2.868615e-11 2.816944e-11
#> Inh 4.093253e-11 5.353214e-11 5.259757e-11
#> Mic 4.877789e-11 5.891194e-11 5.787281e-11
#> OPC 5.891194e-11 8.086453e-11 7.583870e-11
#> Oli 5.787281e-11 7.583870e-11 7.788301e-11
#>
#> , , mash::mash::var8
#>
#> ALL Ast End Exc Inh
#> ALL 9.813251e-09 3.431415e-08 1.747086e-08 7.475763e-09 1.438297e-08
#> Ast 3.431415e-08 1.427757e-03 7.678173e-08 3.296610e-08 6.355927e-08
#> End 1.747086e-08 7.678173e-08 4.375744e-08 1.673708e-08 3.219801e-08
#> Exc 7.475763e-09 3.296610e-08 1.673708e-08 8.864592e-09 1.379326e-08
#> Inh 1.438297e-08 6.355927e-08 3.219801e-08 1.379326e-08 2.947915e-08
#> Mic 8.155501e-09 3.544363e-08 1.827152e-08 7.796524e-09 1.497376e-08
#> OPC -1.214957e-07 -2.996184e-03 -2.717075e-07 -1.169823e-07 -2.259349e-07
#> Oli 1.749783e-08 7.702371e-08 3.917827e-08 1.676801e-08 3.226520e-08
#> Mic OPC Oli
#> ALL 8.155501e-09 -1.214957e-07 1.749783e-08
#> Ast 3.544363e-08 -2.996184e-03 7.702371e-08
#> End 1.827152e-08 -2.717075e-07 3.917827e-08
#> Exc 7.796524e-09 -1.169823e-07 1.676801e-08
#> Inh 1.497376e-08 -2.259349e-07 3.226520e-08
#> Mic 1.179513e-08 -1.242600e-07 1.827550e-08
#> OPC -1.242600e-07 9.280180e-03 -2.729245e-07
#> Oli 1.827550e-08 -2.729245e-07 4.356954e-08
#>
#> , , mash::mash::var9
#>
#> ALL Ast End Exc Inh
#> ALL 0.0003559253 0.0001295760 0.0001401115 0.0001153959 0.0001371165
#> Ast 0.0001295760 0.0004578580 0.0001695291 0.0001396257 0.0001659097
#> End 0.0001401115 0.0001695291 0.0005088667 0.0001509766 0.0001793957
#> Exc 0.0001153959 0.0001396257 0.0001509766 0.0003924719 0.0001477529
#> Inh 0.0001371165 0.0001659097 0.0001793957 0.0001477529 0.0004941634
#> Mic 0.0001208574 0.0001462205 0.0001581158 0.0001302199 0.0001547230
#> OPC 0.0001344312 0.0001626559 0.0001758803 0.0001448556 0.0001721223
#> Oli 0.0001425920 0.0001725310 0.0001865577 0.0001536500 0.0001825727
#> Mic OPC Oli
#> ALL 0.0001208574 0.0001344312 0.0001425920
#> Ast 0.0001462205 0.0001626559 0.0001725310
#> End 0.0001581158 0.0001758803 0.0001865577
#> Exc 0.0001302199 0.0001448556 0.0001536500
#> Inh 0.0001547230 0.0001721223 0.0001825727
#> Mic 0.0004172414 0.0001517066 0.0001609133
#> OPC 0.0001517066 0.0004811061 0.0001789942
#> Oli 0.0001609133 0.0001789942 0.0005211790
#>
#> , , mash::mash::var10
#>
#> ALL Ast End Exc Inh
#> ALL 0.0001638604 0.0001349313 0.0001351592 0.0001363130 0.0001347133
#> Ast 0.0001349313 0.0003074359 0.0001383156 0.0001343776 0.0001371735
#> End 0.0001351592 0.0001383156 0.0005291612 0.0001350802 0.0001371003
#> Exc 0.0001363130 0.0001343776 0.0001350802 0.0001702679 0.0001349952
#> Inh 0.0001347133 0.0001371735 0.0001371003 0.0001349952 0.0001672132
#> Mic 0.0001287492 0.0001383723 0.0001369741 0.0001217623 0.0001352446
#> OPC 0.0001591774 0.0001429146 0.0002016938 0.0001725927 0.0001521843
#> Oli 0.0001349813 0.0001384906 0.0001382440 0.0001342585 0.0001373453
#> Mic OPC Oli
#> ALL 0.0001287492 0.0001591774 0.0001349813
#> Ast 0.0001383723 0.0001429146 0.0001384906
#> End 0.0001369741 0.0002016938 0.0001382440
#> Exc 0.0001217623 0.0001725927 0.0001342585
#> Inh 0.0001352446 0.0001521843 0.0001373453
#> Mic 0.0001857111 0.0001144616 0.0001390134
#> OPC 0.0001144616 0.0128618223 0.0001481780
#> Oli 0.0001390134 0.0001481780 0.0001703273
#>
#> , , mash::mash::var11
#>
#> ALL Ast End Exc Inh
#> ALL 2.538008e-04 8.929927e-05 7.662591e-05 0.0001612659 0.0001184412
#> Ast 8.929927e-05 4.483214e-04 3.784194e-04 0.0001289392 0.0002688095
#> End 7.662591e-05 3.784194e-04 5.402028e-04 0.0001245756 0.0002924973
#> Exc 1.612659e-04 1.289392e-04 1.245756e-04 0.0002365730 0.0001419160
#> Inh 1.184412e-04 2.688095e-04 2.924973e-04 0.0001419160 0.0003356282
#> Mic 9.002180e-05 3.442321e-04 3.837985e-04 0.0001302757 0.0002725189
#> OPC 1.023871e-04 3.183596e-04 3.521095e-04 0.0001364263 0.0002582129
#> Oli 7.865035e-05 3.687384e-04 4.137661e-04 0.0001247310 0.0002861939
#> Mic OPC Oli
#> ALL 0.0000900218 0.0001023871 7.865035e-05
#> Ast 0.0003442321 0.0003183596 3.687384e-04
#> End 0.0003837985 0.0003521095 4.137661e-04
#> Exc 0.0001302757 0.0001364263 1.247310e-04
#> Inh 0.0002725189 0.0002582129 2.861939e-04
#> Mic 0.0004623247 0.0003228892 3.739318e-04
#> OPC 0.0003228892 0.0004167111 3.435047e-04
#> Oli 0.0003739318 0.0003435047 5.143134e-04
#>
#> , , mash::mash::var12
#>
#> ALL Ast End Exc Inh
#> ALL 1.813077e-04 4.693738e-05 5.190729e-05 2.024360e-05 3.377307e-05
#> Ast 4.693738e-05 3.899029e-04 9.934606e-05 3.874401e-05 6.463906e-05
#> End 5.190729e-05 9.934606e-05 4.417057e-04 4.284638e-05 7.148313e-05
#> Exc 2.024360e-05 3.874401e-05 4.284638e-05 1.461251e-04 2.787765e-05
#> Inh 3.377307e-05 6.463906e-05 7.148313e-05 2.787765e-05 2.624205e-04
#> Mic 2.897725e-05 5.546060e-05 6.133279e-05 2.391895e-05 3.990607e-05
#> OPC 5.048651e-05 9.662699e-05 1.068580e-04 4.167360e-05 6.952670e-05
#> Oli 5.013940e-05 9.596254e-05 1.061232e-04 4.138709e-05 6.904858e-05
#> Mic OPC Oli
#> ALL 2.897725e-05 5.048651e-05 5.013940e-05
#> Ast 5.546060e-05 9.662699e-05 9.596254e-05
#> End 6.133279e-05 1.068580e-04 1.061232e-04
#> Exc 2.391895e-05 4.167360e-05 4.138709e-05
#> Inh 3.990607e-05 6.952670e-05 6.904858e-05
#> Mic 2.194913e-04 5.965422e-05 5.924396e-05
#> OPC 5.965422e-05 4.266915e-04 1.032187e-04
#> Oli 5.924396e-05 1.032187e-04 4.230477e-04
#>
#> , , mash::mash::var13
#>
#> ALL Ast End Exc Inh
#> ALL 3.030740e-04 0.0001175361 8.750391e-05 8.881418e-05 9.690526e-05
#> Ast 1.175361e-04 0.0004522954 1.638871e-04 1.149303e-04 1.364276e-04
#> End 8.750391e-05 0.0001638871 5.488517e-04 1.061700e-04 1.471428e-04
#> Exc 8.881418e-05 0.0001149303 1.061700e-04 2.837803e-04 9.684780e-05
#> Inh 9.690526e-05 0.0001364276 1.471428e-04 9.684780e-05 3.583380e-04
#> Mic 7.323611e-05 0.0001662746 2.990968e-04 1.025349e-04 1.534208e-04
#> OPC 9.796187e-05 0.0001662735 2.317464e-04 1.107067e-04 1.467608e-04
#> Oli 8.988823e-05 0.0001602113 2.359681e-04 1.052942e-04 1.427509e-04
#> Mic OPC Oli
#> ALL 7.323611e-05 9.796187e-05 8.988823e-05
#> Ast 1.662746e-04 1.662735e-04 1.602113e-04
#> End 2.990968e-04 2.317464e-04 2.359681e-04
#> Exc 1.025349e-04 1.107067e-04 1.052942e-04
#> Inh 1.534208e-04 1.467608e-04 1.427509e-04
#> Mic 6.668493e-04 2.656678e-04 2.749941e-04
#> OPC 2.656678e-04 5.173515e-04 2.181491e-04
#> Oli 2.749941e-04 2.181491e-04 5.093209e-04
#>
#> , , mash::mash::var14
#>
#> ALL Ast End Exc Inh
#> ALL 1.584345e-05 1.162212e-05 -1.052131e-04 1.669171e-05 1.413864e-05
#> Ast 1.162212e-05 1.120499e-05 -8.155115e-05 1.273342e-05 1.135537e-05
#> End -1.052131e-04 -8.155115e-05 5.278346e-03 -1.199981e-04 -9.989065e-05
#> Exc 1.669171e-05 1.273342e-05 -1.199981e-04 2.029309e-05 1.591832e-05
#> Inh 1.413864e-05 1.135537e-05 -9.989065e-05 1.591832e-05 1.512406e-05
#> Mic 1.345136e-05 1.102103e-05 -9.392407e-05 1.501712e-05 1.306923e-05
#> OPC 1.253892e-05 1.047502e-05 -8.753495e-05 1.387123e-05 1.222362e-05
#> Oli 1.311570e-05 1.082683e-05 -9.147741e-05 1.459229e-05 1.276269e-05
#> Mic OPC Oli
#> ALL 1.345136e-05 1.253892e-05 1.311570e-05
#> Ast 1.102103e-05 1.047502e-05 1.082683e-05
#> End -9.392407e-05 -8.753495e-05 -9.147741e-05
#> Exc 1.501712e-05 1.387123e-05 1.459229e-05
#> Inh 1.306923e-05 1.222362e-05 1.276269e-05
#> Mic 1.436400e-05 1.182663e-05 1.231964e-05
#> OPC 1.182663e-05 1.278706e-05 1.159346e-05
#> Oli 1.231964e-05 1.159346e-05 1.381103e-05
#>
#> , , mash::mash::var15
#>
#> ALL Ast End Exc Inh
#> ALL 1.820217e-04 4.383462e-05 5.610726e-05 3.359970e-05 3.942686e-05
#> Ast 4.383462e-05 3.587041e-04 1.250699e-04 2.397509e-05 1.197085e-04
#> End 5.610726e-05 1.250699e-04 4.616632e-04 3.755636e-05 1.332618e-04
#> Exc 3.359970e-05 2.397509e-05 3.755636e-05 1.865156e-04 1.168093e-05
#> Inh 3.942686e-05 1.197085e-04 1.332618e-04 1.168093e-05 3.753818e-04
#> Mic 5.446991e-05 1.206004e-04 1.414434e-04 3.679894e-05 1.283077e-04
#> OPC 5.603076e-05 1.025911e-04 1.247316e-04 4.758806e-05 1.029973e-04
#> Oli 5.136942e-05 9.489293e-05 1.149965e-04 4.315983e-05 9.567130e-05
#> Mic OPC Oli
#> ALL 5.446991e-05 5.603076e-05 5.136942e-05
#> Ast 1.206004e-04 1.025911e-04 9.489293e-05
#> End 1.414434e-04 1.247316e-04 1.149965e-04
#> Exc 3.679894e-05 4.758806e-05 4.315983e-05
#> Inh 1.283077e-04 1.029973e-04 9.567130e-05
#> Mic 4.411154e-04 1.205894e-04 1.111771e-04
#> OPC 1.205894e-04 4.111389e-04 1.039607e-04
#> Oli 1.111771e-04 1.039607e-04 3.674209e-04
#>
#> , , mash::mash::var16
#>
#> ALL Ast End Exc Inh
#> ALL 2.881235e-07 1.880996e-07 2.148979e-07 2.063863e-07 2.051274e-07
#> Ast 1.880996e-07 4.421215e-07 2.055632e-07 1.870155e-07 2.028269e-07
#> End 2.148979e-07 2.055632e-07 1.467210e-06 2.168373e-07 2.197602e-07
#> Exc 2.063863e-07 1.870155e-07 2.168373e-07 2.932512e-07 2.063472e-07
#> Inh 2.051274e-07 2.028269e-07 2.197602e-07 2.063472e-07 3.101076e-07
#> Mic 2.055678e-07 2.098797e-07 2.220247e-07 2.062554e-07 2.148326e-07
#> OPC 1.972002e-06 7.639814e-08 4.553652e-06 2.097212e-06 1.534630e-06
#> Oli 2.148862e-07 1.952005e-07 2.271265e-07 2.176679e-07 2.156437e-07
#> Mic OPC Oli
#> ALL 2.055678e-07 1.972002e-06 2.148862e-07
#> Ast 2.098797e-07 7.639814e-08 1.952005e-07
#> End 2.220247e-07 4.553652e-06 2.271265e-07
#> Exc 2.062554e-07 2.097212e-06 2.176679e-07
#> Inh 2.148326e-07 1.534630e-06 2.156437e-07
#> Mic 3.244456e-07 1.322582e-06 2.159275e-07
#> OPC 1.322582e-06 9.312459e-03 2.161324e-06
#> Oli 2.159275e-07 2.161324e-06 3.237405e-07
#>
#> , , mash::mash::var17
#>
#> ALL Ast End Exc Inh
#> ALL 0.0005469253 0.0003879874 0.0004009435 0.0003974205 0.0003969694
#> Ast 0.0003879874 0.0011172850 0.0007967213 0.0007361418 0.0008056293
#> End 0.0004009435 0.0007967213 0.0009419880 0.0006729015 0.0007245252
#> Exc 0.0003974205 0.0007361418 0.0006729015 0.0008522451 0.0006766460
#> Inh 0.0003969694 0.0008056293 0.0007245252 0.0006766460 0.0009490760
#> Mic 0.0003976362 0.0007422121 0.0006775022 0.0006370248 0.0006814018
#> OPC 0.0003927180 0.0007503668 0.0006818251 0.0006398519 0.0006862030
#> Oli 0.0003867530 0.0002242706 0.0002876315 0.0003057678 0.0002779092
#> Mic OPC Oli
#> ALL 0.0003976362 0.0003927180 0.0003867530
#> Ast 0.0007422121 0.0007503668 0.0002242706
#> End 0.0006775022 0.0006818251 0.0002876315
#> Exc 0.0006370248 0.0006398519 0.0003057678
#> Inh 0.0006814018 0.0006862030 0.0002779092
#> Mic 0.0008603845 0.0006440029 0.0003037347
#> OPC 0.0006440029 0.0008621509 0.0002928240
#> Oli 0.0003037347 0.0002928240 0.0006966323
#>
#>