matplotlib histogram of nearest neighbor energies

[dave-doty]

I often run code like this in a Jupyter notebook to get a sense of the distribution of nearest neighbor energies of DNA strands of a certain length, to get a calibrated sense of good values to pick for low and high energy values to pass to NearestNeighborEnergyFilter:

import nuad.np as nn
import matplotlib.pyplot as plt

temperature = 37
length = 7
seqs = nn.DNASeqList(length=length)
energies = seqs.energies(temperature=temperature)
# help(plt.hist)
print(f'number of sequences of length {length}: {len(seqs):,}')
plt.figure(figsize=(15,6))
_=plt.hist(energies, bins=30)

number of sequences of length 7: 16,384

Add a helper function to the np submodule to automate this.