count_references: Count alignment hits to a reference database.

In Gibbons-Lab/mbtools: Provides helper functions for microbiome preprocessing and analysis

Count alignment hits to a reference database.

Description

This will correct for effective reference lengths as done by almost any good tool those days. So the returned counts are not correlated with feature lengths. By default an expectation maximization algorithm is used to resolve multiple mappings of one read to many references which pretty much always happens in metagenomics data sets. The optimized likelihood function is very similar to the one in kallisto (https://doi.org/10.1038/nbt.3519).

Usage

count_references(object, ...)

Arguments

object	An experiment data table as returned by any alignment method like align_short_reads or align_long_reads.
...	A configuration as generated by config_count.

Details

Note that for the EM method there will be a NA reference reported which corresponds to the approximate abundance of references not contained in the database.

Value

A list containing the used alignments and the transcript counts in 'counts'.

Gibbons-Lab/mbtools documentation built on Jan. 28, 2024, 11:08 a.m.