metadata_groups: Easily extract groups of individuals from metadata variables.

In microSTASIS: Microbiota STability ASsessment via Iterative cluStering

Description

Easily extract groups of individuals from metadata variables.

Usage

metadata_groups(metadata, samples, individuals, col_number)

Arguments

metadata	data frame with one column of samples matching with the rownames of the original input matrix to paired_times().
samples	vector from metadata corresponding to the samples ID.
individuals	vector of individuals located in the first column of the st_previz() output.
col_number	number of the column with the variable used for grouping individuals.

Value

A vector with the same length as the number of rows in the st_previz() output.

Examples

t1_t2 <- paired_times(data = clr[,1:25], first = "_1",
                      second = "_25", common = "_0_")
t2_t3 <- paired_times(data = clr[,1:25], first = "_25",
                      second = "_26", common = "_0_")
klist_t1_t2 <- iterative_clustering(data = t1_t2, parallel = FALSE)
klist_t2_t3 <- iterative_clustering(data = t2_t3, parallel = FALSE)
result_t1_t2 <- stabilitas(klist_t1_t2, common = "_0_")
result_t2_t3 <- stabilitas(klist_t2_t3, common = "_0_")
results <- st_previz(results = list(result_t1_t2, result_t2_t3),
                     times = c("t1_t2", "t2_t3"))
metadata <- data.frame(Sample = rownames(clr),
                       age = c(rep("youth", 65), rep("old", 131-65)))
group <- metadata_groups(metadata = metadata, samples = metadata$Sample,
                         individuals = results$individual, col_number = 2)

microSTASIS documentation built on Sept. 14, 2021, 9:07 a.m.