Poisson lda works but is giving incorrect results, need to debug