This function simulates MC step for an entire matrix. Should not need to be used by user directly. NOTE: currently written in parallel
1 |
Y |
logratio matrix |
W |
corresponding count matrix |
eY |
current expected value of logratio matrix |
N |
number of samples, or nrow of Y |
Q |
number of OTUs minus the base, or ncol of Y |
base |
OTU index used for base |
sigma |
current estimate of sigma |
MCiter |
number of MC samples to generate |
stepsize |
variance used for MH samples, defaults to 1. Tweak to adjust acceptance ratio |
poorman |
boolean of whether to just use simple diagonal inverse to calculate sigma inverse. Defaults to FALSE |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.