read.lefse.Rd
This reads an input file in a LEfSe input format
read.lefse(filename = "data.txt", rows.meta = 1, row.samples = 2)
name of the input file in a LEfSe input format
specifies in which rows medata variables are stored
specifies in which row sample names are stored
a list with two elements:
feat
a features matrix
meta
a metadate matrix
fn.in.lefse<- system.file("extdata",
"LEfSe_crc_zeller_msb_mocat_specI.tsv",package = "SIAMCAT")
meta.and.features <- read.lefse(fn.in.lefse, rows.meta = 1:6,
row.samples = 7)
meta <- meta.and.features$meta
feat <- meta.and.features$feat
label <- create.label(meta=meta, label="label", case = "cancer")
#> Label used as case:
#> cancer
#> Label used as control:
#> healthy
#> + finished create.label.from.metadata in 0 s
siamcat <- siamcat(feat=feat, label=label, meta=meta)
#> + starting validate.data
#> +++ checking overlap between labels and features
#> + Keeping labels of 141 sample(s).
#> +++ checking sample number per class
#> +++ checking overlap between samples and metadata
#> + finished validate.data in 0.283 s