This is a wrapper arround function SpiecEasi::spiec.easi
with argument method = "mb"
.
correlate_mb( data, pulsar.params = list(thresh = 0.05, subsample.ratio = 0.8, ncores = getOption("mc.cores"), rep.num = 20, seed = 1337), ... )
data | integer matrix of abundance count data. One sample per row and one taxon per column |
---|---|
pulsar.params | list of options passed to pulsar.params |
... | further options passed to spiec.easi |
Kurtz ZD, M?ller CL, Miraldi ER, Littman DR, Blaser MJ, Bonneau RA (2015). “Sparse and compositionally robust inference of microbial ecological networks.” PLoS Comput Biol, 11(5), e1004226.
Meinshausen N, Bühlmann P, others (2006). “High-dimensional graphs and variable selection with the lasso.” The annals of statistics, 34(3), 1436--1462.