Evaluating coevolution

A method on the alignment provides an interface to the simpler (and yet robust and fast) methods for estimating coevolution. The default measure is normalised mutual information (NMI).

Display coevolution as a heatmap

from cogent3 import load_aligned_seqs

aln = load_aligned_seqs("../data/brca1.fasta", moltype='dna')
aln = aln.no_degenerates(motif_length=3)
aln = aln.get_translation()
aln = aln[:100]  # for compute speed in testing the documentation
coevo = aln.coevolution(show_progress=False, drawable='heatmap')
coevo.show()

Display coevolution scores as a Violin plot

coevo = aln.coevolution(show_progress=False, drawable='violin')
coevo.show(width=300)

Display coevolution scores as a Boxplot

coevo = aln.coevolution(show_progress=False, drawable='box')
coevo.show(width=300)