def setUp(self):
res = random_multinomial_model(
num_samples=50,
num_features=5,
reps=1,
low=-1,
high=1,
beta_mean=0,
beta_scale=1,
mu=1000, # sequencing depth
sigma=0.5,
seed=0)
self.table, self.md, self.beta = res
tf.set_random_seed(0)
md = self.md
md.name = 'sampleid'
md = qiime2.Metadata(md)
self.ref_beta, self.ref_stats = multinomial(table=self.table,
metadata=md,
summary_interval=1,
formula="X",
epoch=50000)
self.base_beta, self.base_stats = multinomial(table=self.table,
metadata=md,
summary_interval=1,
formula="1",
epoch=50000)
self.results = "results"
if not os.path.exists(self.results):
os.mkdir(self.results)
def setUp(self):
res = random_multinomial_model(
num_samples=50, num_features=5,
reps=1,
low=-1, high=1,
beta_mean=0,
beta_scale=1,
mu=1000, # sequencing depth
sigma=0.5,
seed=0)
self.table, self.md, self.beta = res