Tried to fix memory allocation issue related to fixed CI_MaxSites / SITESINC constant.   Added methods to compute likelihood of a model for a given sample HCN and to find the optimal model within a parametrized set.
Added basic test for parameter estimation.