SimDirFactorBlock: Generate synthetic OTU abundance table from DirFactor model...

Description Usage Arguments Value Examples

Description

Generate synthetic OTU abundance table from DirFactor model with block diagonal between-sample Gram matrix.

Usage

1
SimDirFactorBlock(vcounts, n, p, m, hyper, K = 1, a = 2)

Arguments

vcounts

Required. A vector of total counts per biological sample. Example: vcounts=c(1e5,1e6).

n

Required. Number of biological samples in the simualted table.

p

Required. Number of species in the simulated dataset.

m

Required. True number of factors in the simulated dataset.

hyper

Required. A list of hyper-parameters in the priors. See Details in DirFactor for the fields in this list.

K

Number of diagonal blocks in the Gram matrix of biological samples. Default is 1 (not block diagonal).

a

Random weights of species are generated by sigma*Q^a*(Q>0). Default is 2.

Value

A list contains the following fields:

Examples

1
2
3
my.hyper = list( nv = 3, a.er = 1, b.er = 0.3, a1 = 3, 
                 a2 = 4, m = 10, alpha = 10, beta = 0 )
SimDirFactorBlock( vcounts=c(1e5,1e6), n = 22, p = 68, m = 3, hyper = my.hyper )

boyuren158/DirFactor documentation built on May 13, 2019, 1:38 a.m.