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.
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.
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]]),
study = "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 625 likelihood matrix.
#> - Likelihood calculations took 0.01 seconds.
#> - Fitting model with 625 mixture components.
#> - Model fitting took 0.02 seconds.
model <- mashModelFit(ssl, alpha = 0, priorCovariances = prior$U,
vhat = vhat)
#> - Computing 68 x 573 likelihood matrix.
#> - Likelihood calculations took 0.02 seconds.
#> - Fitting model with 573 mixture components.
#> - Model fitting took 1.20 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.723269e-01 1.723269e-01 1.723269e-01 1.723269e-01
#> mash::mash::var2 1.476307e-02 1.371301e-02 1.384520e-02 1.396853e-02
#> mash::mash::var3 1.961604e-01 1.961290e-01 1.961331e-01 1.961621e-01
#> mash::mash::var4 -1.061086e-01 -1.059492e-01 -1.059706e-01 -1.061259e-01
#> mash::mash::var5 -2.287148e-02 -2.450765e-02 -3.248295e-02 -2.316511e-02
#> mash::mash::var6 -2.011793e-14 -6.066732e-15 1.617767e-14 -5.756151e-14
#> mash::mash::var7 -2.092369e-07 -4.050663e-01 6.184644e-08 3.453836e-08
#> mash::mash::var8 -3.525636e-03 -1.618265e-01 -4.182271e-02 -6.332220e-03
#> mash::mash::var9 7.823844e-02 -1.172710e-02 3.728663e-02 2.341446e-02
#> mash::mash::var10 -3.964469e-03 1.505138e-04 -9.605846e-03 -7.621832e-03
#> mash::mash::var11 1.743615e-01 1.743612e-01 1.743612e-01 1.743615e-01
#> mash::mash::var12 -1.222330e-01 -2.866656e-03 -5.567241e-02 -1.177342e-01
#> mash::mash::var13 -1.470159e-01 -1.186523e-01 3.509555e-02 -1.301384e-01
#> mash::mash::var14 -3.370937e-02 5.214562e-03 -2.719534e-01 -5.322239e-02
#> mash::mash::var15 5.667875e-02 -4.342221e-03 -5.841683e-02 7.695519e-02
#> mash::mash::var16 9.302826e-06 -6.954558e-06 2.894470e-05 1.048894e-05
#> mash::mash::var17 1.880009e-01 2.800876e-02 2.478557e-03 6.808907e-02
#> Inh Mic OPC Oli
#> mash::mash::var1 1.723269e-01 1.723269e-01 1.723269e-01 1.723269e-01
#> mash::mash::var2 1.452357e-02 1.446077e-02 1.395165e-02 1.529599e-02
#> mash::mash::var3 1.961486e-01 1.961336e-01 1.961363e-01 1.961370e-01
#> mash::mash::var4 -1.060119e-01 -1.059598e-01 -1.058402e-01 -1.058739e-01
#> mash::mash::var5 -2.193268e-02 -2.234293e-02 -2.210340e-02 -2.243901e-02
#> mash::mash::var6 6.532494e-15 8.364501e-15 1.122940e+00 4.225834e-14
#> mash::mash::var7 -6.961329e-08 -2.454014e-08 -4.121774e-07 -3.817868e-07
#> mash::mash::var8 -2.390608e-02 1.819771e-02 -2.240378e-01 -2.591986e-02
#> mash::mash::var9 -9.609314e-02 2.854119e-01 3.698791e-02 -6.222324e-03
#> mash::mash::var10 -1.568961e-03 1.136253e-02 3.204517e-02 1.609956e-04
#> mash::mash::var11 1.743613e-01 1.743612e-01 1.743612e-01 1.743612e-01
#> mash::mash::var12 -1.182798e-02 2.219329e-02 2.481517e-02 -1.002487e-02
#> mash::mash::var13 -1.125961e-01 9.462162e-02 -1.490967e-02 -1.323115e-02
#> mash::mash::var14 -2.065683e-02 9.283119e-03 1.234156e-02 1.444185e-02
#> mash::mash::var15 -2.679104e-02 -1.302579e-02 3.325049e-02 3.469952e-02
#> mash::mash::var16 5.552969e-06 2.632456e-06 -4.194299e-01 1.294375e-05
#> mash::mash::var17 4.256589e-02 6.534125e-02 7.148197e-02 2.445175e-01
#>
#> $PosteriorSD
#> ALL Ast End Exc
#> mash::mash::var1 1.423632e-02 1.423632e-02 1.423632e-02 1.423632e-02
#> mash::mash::var2 2.628676e-02 3.010545e-02 3.075975e-02 2.648087e-02
#> mash::mash::var3 1.658302e-02 1.677592e-02 1.674652e-02 1.658427e-02
#> mash::mash::var4 1.559749e-02 1.601497e-02 1.613970e-02 1.560001e-02
#> mash::mash::var5 3.487428e-02 4.484312e-02 5.223474e-02 3.481717e-02
#> mash::mash::var6 3.887005e-08 3.100527e-08 6.777725e-08 8.264447e-08
#> mash::mash::var7 1.042198e-04 4.330002e-02 1.248490e-04 4.191389e-05
#> mash::mash::var8 1.430060e-02 4.840545e-02 5.895676e-02 1.273612e-02
#> mash::mash::var9 2.440548e-02 3.191805e-02 5.354348e-02 2.682470e-02
#> mash::mash::var10 1.954965e-02 2.261637e-02 4.336101e-02 2.369460e-02
#> mash::mash::var11 1.335228e-02 1.335333e-02 1.335427e-02 1.335227e-02
#> mash::mash::var12 1.578784e-02 3.167685e-02 5.859241e-02 1.365696e-02
#> mash::mash::var13 2.182130e-02 3.558855e-02 5.434310e-02 2.117686e-02
#> mash::mash::var14 2.054213e-02 2.416242e-02 6.609110e-02 2.022394e-02
#> mash::mash::var15 1.580368e-02 2.602597e-02 5.893582e-02 1.574040e-02
#> mash::mash::var16 7.646721e-04 9.693487e-04 1.983569e-03 7.890469e-04
#> mash::mash::var17 2.855566e-02 9.887204e-02 1.339184e-01 7.979914e-02
#> Inh Mic OPC Oli
#> mash::mash::var1 1.423632e-02 1.423632e-02 0.0142363204 1.423632e-02
#> mash::mash::var2 2.641911e-02 2.647316e-02 0.0326170301 2.670324e-02
#> mash::mash::var3 1.660691e-02 1.672227e-02 0.0166957176 1.668881e-02
#> mash::mash::var4 1.578191e-02 1.596081e-02 0.0166743961 1.642412e-02
#> mash::mash::var5 3.543364e-02 3.522021e-02 0.0409894636 3.516139e-02
#> mash::mash::var6 2.970933e-08 2.797418e-08 0.1293032270 6.432506e-08
#> mash::mash::var7 6.784124e-05 8.097093e-05 0.0002083411 1.898192e-04
#> mash::mash::var8 2.064958e-02 2.200773e-02 0.0797189797 3.335747e-02
#> mash::mash::var9 3.537902e-02 2.861635e-02 0.0340676717 3.839359e-02
#> mash::mash::var10 1.930900e-02 3.276899e-02 0.1085508186 2.003082e-02
#> mash::mash::var11 1.335270e-02 1.335366e-02 0.0133534775 1.335387e-02
#> mash::mash::var12 2.086878e-02 1.807112e-02 0.0362416918 3.572884e-02
#> mash::mash::var13 2.647006e-02 4.132784e-02 0.0389168201 3.594188e-02
#> mash::mash::var14 2.807698e-02 3.905170e-02 0.0325106131 3.804792e-02
#> mash::mash::var15 2.594235e-02 3.474575e-02 0.0329555709 2.851760e-02
#> mash::mash::var16 8.313165e-04 9.040518e-04 0.0964911607 9.307571e-04
#> mash::mash::var17 9.257107e-02 8.120049e-02 0.0772461546 4.536175e-02
#>
#> $lfdr
#> ALL Ast End Exc
#> mash::mash::var1 1.258875e-31 1.217716e-31 1.258875e-31 1.258875e-31
#> mash::mash::var2 5.951665e-01 5.879721e-01 5.951665e-01 5.951665e-01
#> mash::mash::var3 2.799192e-30 2.781516e-30 2.799192e-30 2.799192e-30
#> mash::mash::var4 5.910862e-10 5.815139e-10 5.910862e-10 5.910862e-10
#> mash::mash::var5 4.876709e-01 4.766487e-01 4.876709e-01 4.876709e-01
#> mash::mash::var6 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00
#> mash::mash::var7 9.999954e-01 4.899977e-16 9.999954e-01 9.999954e-01
#> mash::mash::var8 6.729062e-02 6.551635e-02 6.729062e-02 6.729062e-02
#> mash::mash::var9 1.475116e-23 1.469591e-23 1.475116e-23 1.475116e-23
#> mash::mash::var10 6.765962e-01 6.703465e-01 6.765962e-01 6.765962e-01
#> mash::mash::var11 8.788862e-37 8.246820e-37 8.788862e-37 8.788862e-37
#> mash::mash::var12 8.661927e-27 8.633273e-27 8.661927e-27 8.661927e-27
#> mash::mash::var13 5.535108e-20 2.132311e-23 5.535108e-20 5.535108e-20
#> mash::mash::var14 1.783882e-04 1.779844e-04 1.783882e-04 1.783882e-04
#> mash::mash::var15 2.516309e-07 2.510081e-07 2.516309e-07 2.516309e-07
#> mash::mash::var16 9.989070e-01 9.988822e-01 9.989070e-01 9.989070e-01
#> mash::mash::var17 8.271099e-18 8.191955e-18 8.271099e-18 8.271099e-18
#> Inh Mic OPC Oli
#> mash::mash::var1 1.258875e-31 1.258875e-31 1.221043e-31 1.258875e-31
#> mash::mash::var2 5.951665e-01 5.951665e-01 5.863796e-01 5.951665e-01
#> mash::mash::var3 2.799192e-30 2.799192e-30 2.704973e-30 2.799192e-30
#> mash::mash::var4 5.910862e-10 5.910862e-10 5.861839e-10 5.910862e-10
#> mash::mash::var5 4.876709e-01 4.876709e-01 4.795097e-01 4.876709e-01
#> mash::mash::var6 1.000000e+00 1.000000e+00 4.304383e-15 1.000000e+00
#> mash::mash::var7 9.999954e-01 9.999954e-01 9.999954e-01 9.999954e-01
#> mash::mash::var8 6.729062e-02 6.729062e-02 1.774268e-03 6.729062e-02
#> mash::mash::var9 1.475116e-23 1.475116e-23 1.455691e-23 1.475116e-23
#> mash::mash::var10 6.765962e-01 6.765962e-01 5.972514e-01 6.765962e-01
#> mash::mash::var11 8.788862e-37 8.788862e-37 8.263720e-37 8.788862e-37
#> mash::mash::var12 8.661927e-27 8.661927e-27 8.593514e-27 8.661927e-27
#> mash::mash::var13 5.535108e-20 5.535108e-20 5.535095e-20 5.535108e-20
#> mash::mash::var14 1.783882e-04 1.783882e-04 1.776436e-04 1.783882e-04
#> mash::mash::var15 2.516309e-07 2.516309e-07 2.482068e-07 2.516309e-07
#> mash::mash::var16 9.989070e-01 9.989070e-01 4.005620e-03 9.989070e-01
#> mash::mash::var17 8.271099e-18 8.271099e-18 8.071203e-18 8.271099e-18
#>
#> $NegativeProb
#> ALL Ast End Exc
#> mash::mash::var1 3.199621e-22 6.291960e-10 3.873869e-09 7.509221e-20
#> mash::mash::var2 4.573514e-02 5.831356e-02 5.417999e-02 5.445646e-02
#> mash::mash::var3 4.159652e-17 9.238328e-05 6.108686e-05 2.774288e-19
#> mash::mash::var4 1.000000e+00 9.995874e-01 9.994273e-01 1.000000e+00
#> mash::mash::var5 4.529558e-01 4.580060e-01 4.771944e-01 4.556996e-01
#> mash::mash::var6 4.680770e-13 3.421648e-13 2.080848e-13 5.603027e-13
#> mash::mash::var7 4.580838e-06 1.000000e+00 1.867783e-06 1.557175e-06
#> mash::mash::var8 5.605228e-01 9.344837e-01 7.208347e-01 6.512927e-01
#> mash::mash::var9 6.735009e-04 6.433445e-01 2.430953e-01 1.913672e-01
#> mash::mash::var10 2.057628e-01 1.663261e-01 2.069564e-01 2.349652e-01
#> mash::mash::var11 1.215720e-24 9.834967e-08 7.435085e-07 4.449304e-21
#> mash::mash::var12 1.000000e+00 5.360473e-01 8.290140e-01 1.000000e+00
#> mash::mash::var13 1.000000e+00 9.995607e-01 2.534150e-01 1.000000e+00
#> mash::mash::var14 9.499102e-01 4.138591e-01 9.998138e-01 9.958788e-01
#> mash::mash::var15 1.415221e-04 5.724803e-01 8.380580e-01 1.314542e-07
#> mash::mash::var16 3.835517e-04 6.021304e-04 3.440546e-04 3.683213e-04
#> mash::mash::var17 2.788582e-10 5.556565e-01 5.574849e-01 1.899452e-01
#> Inh Mic OPC Oli
#> mash::mash::var1 2.668123e-11 1.488648e-10 1.263674e-09 2.386847e-10
#> mash::mash::var2 4.857059e-02 4.929407e-02 5.789243e-02 4.189063e-02
#> mash::mash::var3 8.679296e-09 4.607035e-05 2.624577e-05 2.026810e-05
#> mash::mash::var4 9.998968e-01 9.996804e-01 9.985635e-01 9.988856e-01
#> mash::mash::var5 4.436108e-01 4.475963e-01 4.487960e-01 4.485377e-01
#> mash::mash::var6 2.144023e-13 1.875710e-13 1.384487e-15 1.983544e-14
#> mash::mash::var7 3.251676e-06 2.560317e-06 4.559971e-06 4.582854e-06
#> mash::mash::var8 8.360747e-01 1.771586e-01 9.982257e-01 7.411856e-01
#> mash::mash::var9 9.966973e-01 1.367742e-23 1.388017e-01 5.643733e-01
#> mash::mash::var10 1.795819e-01 8.728935e-02 1.315836e-01 1.621528e-01
#> mash::mash::var11 2.398536e-09 2.844480e-07 1.985082e-07 4.163929e-07
#> mash::mash::var12 7.145664e-01 1.096829e-01 2.467593e-01 6.104611e-01
#> mash::mash::var13 9.999891e-01 1.141784e-02 6.480826e-01 6.424704e-01
#> mash::mash::var14 7.690508e-01 4.066413e-01 3.521535e-01 3.526381e-01
#> mash::mash::var15 8.640381e-01 6.487946e-01 1.553625e-01 1.108950e-01
#> mash::mash::var16 4.461323e-04 4.817317e-04 9.956001e-01 3.599389e-04
#> mash::mash::var17 4.136149e-01 2.113395e-01 1.564198e-01 1.462485e-12
#>
#> $lfsr
#> ALL Ast End Exc
#> mash::mash::var1 3.199621e-22 6.291960e-10 3.873869e-09 7.509221e-20
#> mash::mash::var2 6.409016e-01 6.462856e-01 6.493465e-01 6.496229e-01
#> mash::mash::var3 4.159652e-17 9.238328e-05 6.108686e-05 2.774288e-19
#> mash::mash::var4 4.381088e-08 4.126006e-04 5.726751e-04 7.720361e-09
#> mash::mash::var5 5.470442e-01 5.419940e-01 5.228056e-01 5.443004e-01
#> mash::mash::var6 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00
#> mash::mash::var7 9.999954e-01 4.440892e-16 9.999973e-01 9.999970e-01
#> mash::mash::var8 4.394772e-01 6.551635e-02 2.791653e-01 3.487073e-01
#> mash::mash::var9 6.735009e-04 3.566555e-01 2.430953e-01 1.913672e-01
#> mash::mash::var10 7.942372e-01 8.336739e-01 7.930436e-01 7.650348e-01
#> mash::mash::var11 1.215720e-24 9.834967e-08 7.435085e-07 4.449304e-21
#> mash::mash::var12 4.773959e-15 4.639527e-01 1.709860e-01 0.000000e+00
#> mash::mash::var13 7.269851e-12 4.392677e-04 2.534150e-01 4.175823e-10
#> mash::mash::var14 5.008977e-02 4.140371e-01 1.862322e-04 4.121195e-03
#> mash::mash::var15 1.417737e-04 4.275197e-01 1.619420e-01 3.830852e-07
#> mash::mash::var16 9.992905e-01 9.993979e-01 9.992510e-01 9.992753e-01
#> mash::mash::var17 2.788583e-10 4.443435e-01 4.425151e-01 1.899452e-01
#> Inh Mic OPC Oli
#> mash::mash::var1 2.668123e-11 1.488648e-10 1.263674e-09 2.386847e-10
#> mash::mash::var2 6.437371e-01 6.444606e-01 6.442720e-01 6.370571e-01
#> mash::mash::var3 8.679296e-09 4.607035e-05 2.624577e-05 2.026810e-05
#> mash::mash::var4 1.032166e-04 3.195730e-04 1.436454e-03 1.114368e-03
#> mash::mash::var5 5.563892e-01 5.524037e-01 5.512040e-01 5.514623e-01
#> mash::mash::var6 1.000000e+00 1.000000e+00 5.688870e-15 1.000000e+00
#> mash::mash::var7 9.999967e-01 9.999974e-01 9.999954e-01 9.999954e-01
#> mash::mash::var8 1.639253e-01 2.444492e-01 1.774271e-03 2.588144e-01
#> mash::mash::var9 3.302743e-03 2.842858e-23 1.388017e-01 4.356267e-01
#> mash::mash::var10 8.204181e-01 7.638855e-01 7.288351e-01 8.378472e-01
#> mash::mash::var11 2.398536e-09 2.844480e-07 1.985082e-07 4.163929e-07
#> mash::mash::var12 2.854336e-01 1.096829e-01 2.467593e-01 3.895389e-01
#> mash::mash::var13 1.090280e-05 1.141784e-02 3.519174e-01 3.575296e-01
#> mash::mash::var14 2.309492e-01 4.068197e-01 3.523311e-01 3.528165e-01
#> mash::mash::var15 1.359619e-01 3.512054e-01 1.553628e-01 1.108952e-01
#> mash::mash::var16 9.993531e-01 9.993887e-01 4.399939e-03 9.992669e-01
#> mash::mash::var17 4.136149e-01 2.113395e-01 1.564198e-01 1.462493e-12
#>
#> $PosteriorCov
#> , , mash::mash::var1
#>
#> ALL Ast End Exc Inh
#> ALL 0.0002026727 0.0002026727 0.0002026727 0.0002026727 0.0002026727
#> Ast 0.0002026727 0.0002026728 0.0002026728 0.0002026727 0.0002026727
#> End 0.0002026727 0.0002026728 0.0002026729 0.0002026727 0.0002026728
#> Exc 0.0002026727 0.0002026727 0.0002026727 0.0002026727 0.0002026727
#> Inh 0.0002026727 0.0002026727 0.0002026728 0.0002026727 0.0002026727
#> Mic 0.0002026727 0.0002026728 0.0002026728 0.0002026727 0.0002026727
#> OPC 0.0002026727 0.0002026728 0.0002026728 0.0002026727 0.0002026728
#> Oli 0.0002026727 0.0002026728 0.0002026728 0.0002026727 0.0002026727
#> Mic OPC Oli
#> ALL 0.0002026727 0.0002026727 0.0002026727
#> Ast 0.0002026728 0.0002026728 0.0002026728
#> End 0.0002026728 0.0002026728 0.0002026728
#> Exc 0.0002026727 0.0002026727 0.0002026727
#> Inh 0.0002026727 0.0002026728 0.0002026727
#> Mic 0.0002026728 0.0002026728 0.0002026727
#> OPC 0.0002026728 0.0002026728 0.0002026728
#> Oli 0.0002026727 0.0002026728 0.0002026728
#>
#> , , mash::mash::var2
#>
#> ALL Ast End Exc Inh
#> ALL 0.0006909938 0.0006458895 0.0006409679 0.0006415301 0.0006439243
#> Ast 0.0006458895 0.0009063380 0.0006608438 0.0006588277 0.0006497759
#> End 0.0006409679 0.0006608438 0.0009461624 0.0006575217 0.0006459222
#> Exc 0.0006415301 0.0006588277 0.0006575217 0.0007012365 0.0006456850
#> Inh 0.0006439243 0.0006497759 0.0006459222 0.0006456850 0.0006979693
#> Mic 0.0006436743 0.0006507968 0.0006472260 0.0006467765 0.0006446651
#> OPC 0.0006437204 0.0006587955 0.0006570388 0.0006552477 0.0006471921
#> Oli 0.0006473391 0.0006371660 0.0006297452 0.0006321984 0.0006429540
#> Mic OPC Oli
#> ALL 0.0006436743 0.0006437204 0.0006473391
#> Ast 0.0006507968 0.0006587955 0.0006371660
#> End 0.0006472260 0.0006570388 0.0006297452
#> Exc 0.0006467765 0.0006552477 0.0006321984
#> Inh 0.0006446651 0.0006471921 0.0006429540
#> Mic 0.0007008280 0.0006481028 0.0006417810
#> OPC 0.0006481028 0.0010638707 0.0006359785
#> Oli 0.0006417810 0.0006359785 0.0007130628
#>
#> , , mash::mash::var3
#>
#> ALL Ast End Exc Inh
#> ALL 0.0002749965 0.0002745978 0.0002746345 0.0002748951 0.0002747741
#> Ast 0.0002745978 0.0002814315 0.0002803586 0.0002742504 0.0002770865
#> End 0.0002746345 0.0002803586 0.0002804458 0.0002743340 0.0002767864
#> Exc 0.0002748951 0.0002742504 0.0002743340 0.0002750379 0.0002746527
#> Inh 0.0002747741 0.0002770865 0.0002767864 0.0002746527 0.0002757893
#> Mic 0.0002746394 0.0002802435 0.0002795163 0.0002743452 0.0002767462
#> OPC 0.0002746636 0.0002796766 0.0002790260 0.0002744004 0.0002765482
#> Oli 0.0002746695 0.0002795388 0.0002789069 0.0002744139 0.0002765000
#> Mic OPC Oli
#> ALL 0.0002746394 0.0002746636 0.0002746695
#> Ast 0.0002802435 0.0002796766 0.0002795388
#> End 0.0002795163 0.0002790260 0.0002789069
#> Exc 0.0002743452 0.0002744004 0.0002744139
#> Inh 0.0002767462 0.0002765482 0.0002765000
#> Mic 0.0002796343 0.0002789389 0.0002788223
#> OPC 0.0002789389 0.0002787470 0.0002784052
#> Oli 0.0002788223 0.0002784052 0.0002785163
#>
#> , , mash::mash::var4
#>
#> ALL Ast End Exc Inh
#> ALL 0.0002432816 0.0002417954 0.0002417974 0.0002418336 0.0002418153
#> Ast 0.0002417954 0.0002564791 0.0002517695 0.0002405471 0.0002487880
#> End 0.0002417974 0.0002517695 0.0002604898 0.0002407170 0.0002478493
#> Exc 0.0002418336 0.0002405471 0.0002407170 0.0002433603 0.0002410597
#> Inh 0.0002418153 0.0002487880 0.0002478493 0.0002410597 0.0002490688
#> Mic 0.0002417981 0.0002525536 0.0002511082 0.0002406328 0.0002483258
#> OPC 0.0002417593 0.0002611878 0.0002585803 0.0002396547 0.0002535424
#> Oli 0.0002417708 0.0002587611 0.0002564802 0.0002399302 0.0002520766
#> Mic OPC Oli
#> ALL 0.0002417981 0.0002417593 0.0002417708
#> Ast 0.0002525536 0.0002611878 0.0002587611
#> End 0.0002511082 0.0002585803 0.0002564802
#> Exc 0.0002406328 0.0002396547 0.0002399302
#> Inh 0.0002483258 0.0002535424 0.0002520766
#> Mic 0.0002547474 0.0002598989 0.0002576340
#> OPC 0.0002598989 0.0002780355 0.0002703759
#> Oli 0.0002576340 0.0002703759 0.0002697516
#>
#> , , mash::mash::var5
#>
#> ALL Ast End Exc Inh
#> ALL 0.0012162151 0.0010584148 0.0009978922 0.001097171 0.0011102730
#> Ast 0.0010584148 0.0020109051 0.0008905208 0.001053293 0.0010748045
#> End 0.0009978922 0.0008905208 0.0027284684 0.001032203 0.0008884603
#> Exc 0.0010971707 0.0010532928 0.0010322027 0.001212236 0.0011034956
#> Inh 0.0011102730 0.0010748045 0.0008884603 0.001103496 0.0012555427
#> Mic 0.0011059307 0.0010676553 0.0009365401 0.001100761 0.0011224696
#> OPC 0.0011076040 0.0010713326 0.0008979656 0.001101200 0.0011280800
#> Oli 0.0011049116 0.0010659797 0.0009477745 0.001100117 0.0011202519
#> Mic OPC Oli
#> ALL 0.0011059307 0.0011076040 0.0011049116
#> Ast 0.0010676553 0.0010713326 0.0010659797
#> End 0.0009365401 0.0008979656 0.0009477745
#> Exc 0.0011007607 0.0011011999 0.0011001173
#> Inh 0.0011224696 0.0011280800 0.0011202519
#> Mic 0.0012404630 0.0011191326 0.0011135595
#> OPC 0.0011191326 0.0016801361 0.0011170372
#> Oli 0.0011135595 0.0011170372 0.0012363236
#>
#> , , mash::mash::var6
#>
#> ALL Ast End Exc Inh
#> ALL 1.510881e-15 2.199846e-16 -5.789633e-16 2.070331e-15 -2.325893e-16
#> Ast 2.199846e-16 9.613269e-16 -1.743369e-16 6.270680e-16 -6.956340e-17
#> End -5.789633e-16 -1.743369e-16 4.593755e-15 -1.657957e-15 1.885407e-16
#> Exc 2.070331e-15 6.270680e-16 -1.657957e-15 6.830108e-15 -6.673305e-16
#> Inh -2.325893e-16 -6.956340e-17 1.885407e-16 -6.673305e-16 8.826444e-16
#> Mic -2.982280e-16 -8.941461e-17 2.411722e-16 -8.550172e-16 9.784729e-17
#> OPC 2.001059e-14 6.032698e-15 -1.609534e-14 5.725472e-14 -6.500396e-15
#> Oli -1.516551e-15 -4.579683e-16 1.218027e-15 -4.339707e-15 4.914371e-16
#> Mic OPC Oli
#> ALL -2.982280e-16 2.001059e-14 -1.516551e-15
#> Ast -8.941461e-17 6.032698e-15 -4.579683e-16
#> End 2.411722e-16 -1.609534e-14 1.218027e-15
#> Exc -8.550172e-16 5.725472e-14 -4.339707e-15
#> Inh 9.784729e-17 -6.500396e-15 4.914371e-16
#> Mic 7.825547e-16 -8.323507e-15 6.290898e-16
#> OPC -8.323507e-15 1.671932e-02 -4.203599e-14
#> Oli 6.290898e-16 -4.203599e-14 4.137713e-15
#>
#> , , mash::mash::var7
#>
#> ALL Ast End Exc Inh
#> ALL 1.086176e-08 -3.924292e-08 -2.814017e-09 -1.571490e-09 3.167434e-09
#> Ast -3.924292e-08 1.874891e-03 1.159923e-08 6.477631e-09 -1.305607e-08
#> End -2.814017e-09 1.159923e-08 1.558728e-08 4.645615e-10 -9.361999e-10
#> Exc -1.571490e-09 6.477631e-09 4.645615e-10 1.756774e-09 -5.228084e-10
#> Inh 3.167434e-09 -1.305607e-08 -9.361999e-10 -5.228084e-10 4.602434e-09
#> Mic 1.116597e-09 -4.602556e-09 -3.300007e-10 -1.842727e-10 3.715214e-10
#> OPC 1.875426e-08 -7.730311e-08 -5.543474e-09 -3.095758e-09 6.239623e-09
#> Oli 1.737143e-08 -7.160393e-08 -5.134717e-09 -2.867491e-09 5.779536e-09
#> Mic OPC Oli
#> ALL 1.116597e-09 1.875426e-08 1.737143e-08
#> Ast -4.602556e-09 -7.730311e-08 -7.160393e-08
#> End -3.300007e-10 -5.543474e-09 -5.134717e-09
#> Exc -1.842727e-10 -3.095758e-09 -2.867491e-09
#> Inh 3.715214e-10 6.239623e-09 5.779536e-09
#> Mic 6.556292e-09 2.199610e-09 2.037415e-09
#> OPC 2.199610e-09 4.340601e-08 3.422076e-08
#> Oli 2.037415e-09 3.422076e-08 3.603132e-08
#>
#> , , mash::mash::var8
#>
#> ALL Ast End Exc Inh
#> ALL 2.045072e-04 3.962181e-05 1.063809e-05 1.610659e-06 6.080761e-06
#> Ast 3.962181e-05 2.343087e-03 4.700206e-04 7.116258e-05 2.686636e-04
#> End 1.063809e-05 4.700206e-04 3.475900e-03 1.910649e-05 7.213380e-05
#> Exc 1.610659e-06 7.116258e-05 1.910649e-05 1.622088e-04 1.092133e-05
#> Inh 6.080761e-06 2.686636e-04 7.213380e-05 1.092133e-05 4.264051e-04
#> Mic -4.628785e-06 -2.045119e-04 -5.490967e-05 -8.313512e-06 -3.138638e-05
#> OPC -6.403711e-05 -3.023150e-03 -7.596182e-04 -1.150142e-04 -4.342075e-04
#> Oli 6.593058e-06 2.913020e-04 7.821228e-05 1.184142e-05 4.470582e-05
#> Mic OPC Oli
#> ALL -4.628785e-06 -6.403711e-05 6.593058e-06
#> Ast -2.045119e-04 -3.023150e-03 2.913020e-04
#> End -5.490967e-05 -7.596182e-04 7.821228e-05
#> Exc -8.313512e-06 -1.150142e-04 1.184142e-05
#> Inh -3.138638e-05 -4.342075e-04 4.470582e-05
#> Mic 4.843401e-04 3.305246e-04 -3.403101e-05
#> OPC 3.305246e-04 6.355116e-03 -4.707689e-04
#> Oli -3.403101e-05 -4.707689e-04 1.112721e-03
#>
#> , , mash::mash::var9
#>
#> ALL Ast End Exc Inh
#> ALL 5.956276e-04 -1.983508e-11 3.913313e-11 3.091547e-11 -1.859103e-10
#> Ast -1.983508e-11 1.018762e-03 -8.681732e-12 -6.858629e-12 4.124448e-11
#> End 3.913313e-11 -8.681732e-12 2.866904e-03 1.353159e-11 -8.137222e-11
#> Exc 3.091547e-11 -6.858629e-12 1.353159e-11 7.195647e-04 -6.428465e-11
#> Inh -1.859103e-10 4.124448e-11 -8.137222e-11 -6.428465e-11 1.251675e-03
#> Mic 4.143896e-10 -9.193294e-11 1.813768e-10 1.432890e-10 -8.616692e-10
#> OPC 6.819339e-11 -1.512879e-11 2.984799e-11 2.358014e-11 -1.417992e-10
#> Oli -1.330085e-11 2.950825e-12 -5.821724e-12 -4.599202e-12 2.765740e-11
#> Mic OPC Oli
#> ALL 4.143896e-10 6.819339e-11 -1.330085e-11
#> Ast -9.193294e-11 -1.512879e-11 2.950825e-12
#> End 1.813768e-10 2.984799e-11 -5.821724e-12
#> Exc 1.432890e-10 2.358014e-11 -4.599202e-12
#> Inh -8.616692e-10 -1.417992e-10 2.765740e-11
#> Mic 8.188957e-04 3.160671e-10 -6.164766e-11
#> OPC 3.160671e-10 1.160606e-03 -1.014494e-11
#> Oli -6.164766e-11 -1.014494e-11 1.474068e-03
#>
#> , , mash::mash::var10
#>
#> ALL Ast End Exc Inh
#> ALL 3.821886e-04 9.201332e-05 2.544001e-04 2.197672e-04 1.181715e-04
#> Ast 9.201332e-05 5.115001e-04 8.040490e-05 8.459545e-05 9.686519e-05
#> End 2.544001e-04 8.040490e-05 1.880177e-03 4.153373e-04 1.491398e-04
#> Exc 2.197672e-04 8.459545e-05 4.153373e-04 5.614340e-04 1.379234e-04
#> Inh 1.181715e-04 9.686519e-05 1.491398e-04 1.379234e-04 3.728374e-04
#> Mic -1.023382e-04 1.235168e-04 -4.294190e-04 -3.110269e-04 3.414050e-05
#> OPC 1.420161e-04 1.033217e-04 1.890650e-04 1.737007e-04 1.211899e-04
#> Oli 8.872148e-05 1.004243e-04 7.188187e-05 7.796072e-05 9.575987e-05
#> Mic OPC Oli
#> ALL -1.023382e-04 1.420161e-04 8.872148e-05
#> Ast 1.235168e-04 1.033217e-04 1.004243e-04
#> End -4.294190e-04 1.890650e-04 7.188187e-05
#> Exc -3.110269e-04 1.737007e-04 7.796072e-05
#> Inh 3.414050e-05 1.211899e-04 9.575987e-05
#> Mic 1.073807e-03 1.380259e-05 1.344216e-04
#> OPC 1.380259e-05 1.178328e-02 1.067552e-04
#> Oli 1.344216e-04 1.067552e-04 4.012337e-04
#>
#> , , mash::mash::var11
#>
#> ALL Ast End Exc Inh
#> ALL 0.0001782835 0.0001782782 0.0001782769 0.0001782821 0.0001782797
#> Ast 0.0001782782 0.0001783114 0.0001783181 0.0001782812 0.0001782988
#> End 0.0001782769 0.0001783181 0.0001783364 0.0001782809 0.0001783043
#> Exc 0.0001782821 0.0001782812 0.0001782809 0.0001782830 0.0001782815
#> Inh 0.0001782797 0.0001782988 0.0001783043 0.0001782815 0.0001782945
#> Mic 0.0001782776 0.0001783131 0.0001783233 0.0001782811 0.0001783012
#> OPC 0.0001782780 0.0001783108 0.0001783202 0.0001782812 0.0001782997
#> Oli 0.0001782773 0.0001783155 0.0001783265 0.0001782810 0.0001783027
#> Mic OPC Oli
#> ALL 0.0001782776 0.0001782780 0.0001782773
#> Ast 0.0001783131 0.0001783108 0.0001783155
#> End 0.0001783233 0.0001783202 0.0001783265
#> Exc 0.0001782811 0.0001782812 0.0001782810
#> Inh 0.0001783012 0.0001782997 0.0001783027
#> Mic 0.0001783202 0.0001783149 0.0001783203
#> OPC 0.0001783149 0.0001783154 0.0001783175
#> Oli 0.0001783203 0.0001783175 0.0001783258
#>
#> , , mash::mash::var12
#>
#> ALL Ast End Exc Inh
#> ALL 2.492557e-04 1.408806e-07 1.971550e-06 1.441377e-06 3.073965e-07
#> Ast 1.408806e-07 1.003423e-03 1.428308e-07 1.044210e-07 2.227050e-08
#> End 1.971550e-06 1.428308e-07 3.433071e-03 1.461319e-06 3.116503e-07
#> Exc 1.441377e-06 1.044210e-07 1.461319e-06 1.865126e-04 2.278431e-07
#> Inh 3.073965e-07 2.227050e-08 3.116503e-07 2.278431e-07 4.355061e-04
#> Mic -4.513040e-07 -3.269363e-08 -4.575468e-07 -3.345080e-07 -7.133792e-08
#> OPC -1.453028e-06 -1.052643e-07 -1.473130e-06 -1.076989e-06 -2.296840e-07
#> Oli 5.768817e-07 4.179360e-08 5.848641e-07 4.275864e-07 9.119069e-08
#> Mic OPC Oli
#> ALL -4.513040e-07 -1.453028e-06 5.768817e-07
#> Ast -3.269363e-08 -1.052643e-07 4.179360e-08
#> End -4.575468e-07 -1.473130e-06 5.848641e-07
#> Exc -3.345080e-07 -1.076989e-06 4.275864e-07
#> Inh -7.133792e-08 -2.296840e-07 9.119069e-08
#> Mic 3.265655e-04 3.372128e-07 -1.338786e-07
#> OPC 3.372128e-07 1.313460e-03 -4.310418e-07
#> Oli -1.338786e-07 -4.310418e-07 1.276550e-03
#>
#> , , mash::mash::var13
#>
#> ALL Ast End Exc Inh
#> ALL 4.761690e-04 1.955052e-07 -1.964941e-05 1.670598e-06 -5.904133e-07
#> Ast 1.955052e-07 1.266545e-03 5.613971e-06 6.971937e-07 1.223253e-06
#> End -1.964941e-05 5.613971e-06 2.953172e-03 -4.615568e-06 1.100979e-05
#> Exc 1.670598e-06 6.971937e-07 -4.615568e-06 4.484592e-04 4.843519e-07
#> Inh -5.904133e-07 1.223253e-06 1.100979e-05 4.843519e-07 7.006638e-04
#> Mic -2.733612e-05 7.375147e-06 1.955936e-04 -6.675391e-06 1.479524e-05
#> OPC -1.319613e-05 4.127887e-06 9.802544e-05 -2.888580e-06 7.825638e-06
#> Oli -1.341297e-05 4.177405e-06 9.952067e-05 -2.946736e-06 7.932290e-06
#> Mic OPC Oli
#> ALL -2.733612e-05 -1.319613e-05 -1.341297e-05
#> Ast 7.375147e-06 4.127887e-06 4.177405e-06
#> End 1.955936e-04 9.802544e-05 9.952067e-05
#> Exc -6.675391e-06 -2.888580e-06 -2.946736e-06
#> Inh 1.479524e-05 7.825638e-06 7.932290e-06
#> Mic 1.707990e-03 1.343787e-04 1.364338e-04
#> OPC 1.343787e-04 1.514519e-03 6.851087e-05
#> Oli 1.364338e-04 6.851087e-05 1.291819e-03
#>
#> , , mash::mash::var14
#>
#> ALL Ast End Exc Inh
#> ALL 4.219792e-04 1.135024e-07 4.735343e-05 6.656436e-06 4.902050e-06
#> Ast 1.135024e-07 5.838225e-04 -4.379314e-05 -3.158948e-06 1.989563e-06
#> End 4.735343e-05 -4.379314e-05 4.368033e-03 8.874998e-05 3.763359e-05
#> Exc 6.656436e-06 -3.158948e-06 8.874998e-05 4.090077e-04 5.995103e-06
#> Inh 4.902050e-06 1.989563e-06 3.763359e-05 5.995103e-06 7.883170e-04
#> Mic -2.024314e-06 7.835428e-06 -7.324767e-05 -6.699345e-06 -1.713435e-08
#> OPC -2.102918e-06 8.300953e-06 -7.701500e-05 -7.041911e-06 4.393699e-08
#> Oli -3.476550e-06 8.896871e-06 -9.546455e-05 -9.279262e-06 -1.269071e-06
#> Mic OPC Oli
#> ALL -2.024314e-06 -2.102918e-06 -3.476550e-06
#> Ast 7.835428e-06 8.300953e-06 8.896871e-06
#> End -7.324767e-05 -7.701500e-05 -9.546455e-05
#> Exc -6.699345e-06 -7.041911e-06 -9.279262e-06
#> Inh -1.713435e-08 4.393699e-08 -1.269071e-06
#> Mic 1.525035e-03 1.067127e-05 1.204646e-05
#> OPC 1.067127e-05 1.056940e-03 1.270271e-05
#> Oli 1.204646e-05 1.270271e-05 1.447644e-03
#>
#> , , mash::mash::var15
#>
#> ALL Ast End Exc Inh
#> ALL 2.497564e-04 -8.334163e-06 -3.448466e-05 1.587193e-05 -1.862465e-05
#> Ast -8.334163e-06 6.773510e-04 1.065095e-04 -2.792673e-05 7.418229e-05
#> End -3.448466e-05 1.065095e-04 3.473431e-03 -7.832001e-05 1.664662e-04
#> Exc 1.587193e-05 -2.792673e-05 -7.832001e-05 2.477603e-04 -4.813811e-05
#> Inh -1.862465e-05 7.418229e-05 1.664662e-04 -4.813811e-05 6.730053e-04
#> Mic -1.561690e-05 6.098681e-05 1.378687e-04 -3.988391e-05 9.235293e-05
#> OPC 1.620102e-05 1.174983e-05 -1.195485e-05 1.507269e-05 3.118979e-06
#> Oli 1.351432e-05 1.103544e-05 -7.537432e-06 1.217714e-05 4.327816e-06
#> Mic OPC Oli
#> ALL -1.561690e-05 1.620102e-05 1.351432e-05
#> Ast 6.098681e-05 1.174983e-05 1.103544e-05
#> End 1.378687e-04 -1.195485e-05 -7.537432e-06
#> Exc -3.988391e-05 1.507269e-05 1.217714e-05
#> Inh 9.235293e-05 3.118979e-06 4.327816e-06
#> Mic 1.207267e-03 1.348600e-06 2.587658e-06
#> OPC 1.348600e-06 1.086070e-03 3.521218e-05
#> Oli 2.587658e-06 3.521218e-05 8.132535e-04
#>
#> , , mash::mash::var16
#>
#> ALL Ast End Exc Inh
#> ALL 5.847234e-07 5.668066e-08 5.453314e-07 2.847402e-07 2.161922e-07
#> Ast 5.668066e-08 9.396369e-07 -3.065605e-07 3.502090e-08 1.242021e-07
#> End 5.453314e-07 -3.065605e-07 3.934545e-06 6.120703e-07 3.369853e-07
#> Exc 2.847402e-07 3.502090e-08 6.120703e-07 6.225951e-07 2.233917e-07
#> Inh 2.161922e-07 1.242021e-07 3.369853e-07 2.233917e-07 6.910872e-07
#> Mic 1.748719e-07 1.784132e-07 1.702209e-07 1.745955e-07 1.757303e-07
#> OPC 3.621336e-06 -2.139007e-06 1.046817e-05 4.036148e-06 2.305237e-06
#> Oli 3.204191e-07 -1.234227e-08 7.572199e-07 3.464740e-07 2.391654e-07
#> Mic OPC Oli
#> ALL 1.748719e-07 3.621336e-06 3.204191e-07
#> Ast 1.784132e-07 -2.139007e-06 -1.234227e-08
#> End 1.702209e-07 1.046817e-05 7.572199e-07
#> Exc 1.745955e-07 4.036148e-06 3.464740e-07
#> Inh 1.757303e-07 2.305237e-06 2.391654e-07
#> Mic 8.173097e-07 1.288184e-06 1.739824e-07
#> OPC 1.288184e-06 9.310544e-03 4.882203e-06
#> Oli 1.739824e-07 4.882203e-06 8.663089e-07
#>
#> , , mash::mash::var17
#>
#> ALL Ast End Exc Inh
#> ALL 0.0008154257 0.0003229833 0.000355803 0.0002714591 0.0003042697
#> Ast 0.0003229833 0.0097756802 0.009961036 0.0065525108 0.0078784629
#> End 0.0003558030 0.0099610360 0.017934131 0.0075547883 0.0090870891
#> Exc 0.0002714591 0.0065525108 0.007554788 0.0063679027 0.0059810196
#> Inh 0.0003042697 0.0078784629 0.009087089 0.0059810196 0.0085694021
#> Mic 0.0002749915 0.0066952628 0.007719756 0.0050868934 0.0061111046
#> OPC 0.0002670974 0.0063762465 0.007351093 0.0048458177 0.0058203962
#> Oli 0.0000446557 -0.0026131053 -0.003037209 -0.0019472981 -0.0023712848
#> Mic OPC Oli
#> ALL 0.0002749915 0.0002670974 0.0000446557
#> Ast 0.0066952628 0.0063762465 -0.0026131053
#> End 0.0077197559 0.0073510933 -0.0030372086
#> Exc 0.0050868934 0.0048458177 -0.0019472981
#> Inh 0.0061111046 0.0058203962 -0.0023712848
#> Mic 0.0065935200 0.0049507409 -0.0019929446
#> OPC 0.0049507409 0.0059669684 -0.0018909359
#> Oli -0.0019929446 -0.0018909359 0.0020576882
#>
#>