This function estimates the LNM model fit from Xia et al.
1 2 | LNM.EM.nocov(W, base, EMiter = 10, EMburn = 5, MCiter = 1000,
MCburn = 500, stepsize = 0.01, p = 0.05, poorman = FALSE)
|
W |
count matrix, with OTUs as columns |
base |
OTU index to be used for base |
EMiter |
number of EM iterations, defaults to 10 |
EMburn |
number of EM iterations to burn, defaults to 5 |
MCiter |
number MC iterations, defaults to 1000 |
MCburn |
number of MC iterations to burn, defaults to 500 |
stepsize |
variance used for MH samples, defaults to 0.01. Tweak to adjust acceptance ratio |
p |
size of purturbation used for logratios, defaults to 0.05 |
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.