Model with 2 states for amino acids ? about ccmgen HOT 1 OPEN

Hi @odannis,

sorry for the long wait - I have been busy lately. Generally speaking the two-state MRF is a special case of the standard 20 state mrf, where all the singleton potentials corresponding to the 18 "forbidden" states are set to -infinity.

By default training should make the singleton potentials small, given that there's no prior on the singletons. You need to be careful with pseudocounts, though! I would not add any, just to be sure.

If the "forbidden" amino acids still pop up during simulation, you could also modify the parameters manually and set the singleton potentials corresponding to the unwanted states to a large negative number.

from ccmgen.