Calculate tree mapping quality call — tree_quality

This function returns the "Core", "I1", "I2", "I3", "PoorQ" calls used to help assess Patch-seq quality. These scores are based on a combination of KL divergence, correlation to the best matching type, and tree mapping to the top two types.

tree_quality_call(AIT.anndata, query.mapping)

Arguments

AIT.anndata: A reference taxonomy object on which buildPatchseqTaxonomy has already been run
query.mapping: The output from taxonomy_mapping with both corr.map and tree.map set to TRUE, and with the same mode as currently set in AIT.anndata

Value

A matrix with "Tree_call" and some other useful columns, one per cell in query.mapping