Generate beta matrix giving genetic effect per SNP on taxa abundances.
compute_beta_matrix_cluster.Rd
Generate beta matrix giving genetic effect per SNP on taxa abundances.
Usage
compute_beta_matrix_cluster(
n_b,
n_g,
n_clust,
n_qtl_o,
n_otus,
effect_size = 1,
correlation = 1,
beta_info = NULL
)
Arguments
- n_b
Number of taxa
- n_g
Number of SNPs
- n_clust
A vector with a length matching the total number of taxa with values from 0 to the number of clusters. Typically the
assign_taxa()
output.- n_qtl_o
Number of causative QTL on taxa abundances (per taxon)
- n_otus
Number of taxa under genetic control.
- effect_size
Vector giving the size of genetic effect to try
- correlation
A numeric value between 0 and 1 representing the target correlation between taxa within the same cluster.
- beta_info
Informations from Beta matrix of genetic effects
Value
A list
of two objects linked to the beta matrix.
The beta matrix itself : a matrix
of taxa (in rows) across SNPs (in columns) giving the multiplicative effect of genotype on taxa abundances.
A small data.frame
containing the parameters used during the beta matrix simulation, including cluster
, id_otu
and id_qtl_o
.