Skip to content

Instantly share code, notes, and snippets.

@kyleabeauchamp
Created November 1, 2012 04:25
Show Gist options
  • Save kyleabeauchamp/3991712 to your computer and use it in GitHub Desktop.
Save kyleabeauchamp/3991712 to your computer and use it in GitHub Desktop.
Test driver for new mle_code
from msmbuilder import io
from msmbuilder import MSMLib
ass = io.loadh("~/dat/redesign/NTL9/Ward-2500/Assignments.Fixed.h5")["Data"]
c = MSMLib.get_count_matrix_from_assignments(ass)
m = MSMLib.reversible_MLE_estimator(c)
u = m.optimize()
c0 = c + c.transpose()
c.data[:] = 1.
c3 = MSMLib.mle_reversible_count_matrix(c,initial_guess=c0)
c2 = m.log_vector_to_matrix(u)
c2 /= c2.sum()
c3 /= c3.sum()
d = c2 - c3
u2 = m.matrix_to_log_vector(c3)
m.log_likelihood(u2) - m.log_likelihood(u)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment