Commit c086bdd6 authored by boussau's avatar boussau
Browse files

Replacing a deprecated function call in

parent 30d7f3d3
......@@ -36,7 +36,7 @@ import pandas as pd
from ete3 import Tree, NodeStyle, TreeStyle, TextFace
from Bio import AlignIO
from scipy.stats import chisqprob
from scipy.stats import chi2
from scipy.stats import multinomial
......@@ -284,7 +284,7 @@ def multinomial_lrt (vector1, vector2, vectorAll):
lk2 = rv2.logpmf(vector2)
lkAll = rvAll.logpmf([sum(x) for x in zip(vector1, vector2)] )
lr = 2*(-lkAll + lk1 + lk2)
lrt = chisqprob(lr, 1)
lrt = chi2.sf(lr, 1)
return lr, lrt
def applyMultinomialLRT(x, orderedAA):
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment