test hooks
Bugfix in likelihood.py The problem was that scipy.sparse matrices use np.matrices as opposed to np.arrays, so when I was using the '*' operator it was doing matrix multiplication when I really wanted element-wise multiplication. The solution was to cast everything to arrays so that it is obvious what is going on.