Compute carriers-per-trajectory matrix \(k_{jT}\) for many phylotypes

Compute carriers-per-trajectory matrix \(k_{jT}\) for many phylotypes

Usage

compute.carriers.per.trajectory.mat(Z, trajectories, eps = 0)

Arguments

Z

Numeric matrix (samples x phylotypes) of relative abundance.

trajectories

List of integer vectors; each vector are sample indices (columns of Z) belonging to a trajectory.

eps

Numeric scalar; treat z > eps as "detected". Default 0.

Value

Integer matrix k.mat of dimension (n.phylotypes x n.trajectories), where k.mat[j, t] is the number of samples on trajectory t with Z[, j] > eps.