コード例 #1
0
ファイル: neg_binom_test.py プロジェクト: 16NWallace/phylowgs 
def testNegBinomWSd(k, r, mu, sd):
	truth = mean([nbinom.logpmf(k, i, mu) for i in range(int(r-0.5*sd),int(r+0.5*sd)+1)])
	val = u.log_neg_binom_likelihood(k,r,mu,sd)
	print truth, val
	return truth==val
コード例 #2
0
ファイル: neg_binom_test.py プロジェクト: 16NWallace/phylowgs 
def testNegBinomNoSd(k, r, mu):
	truth = nbinom.logpmf(k,r,mu)
	val = u.log_neg_binom_likelihood(k,r,mu)
	return truth==val
コード例 #3
0
ファイル: data.py プロジェクト: 16NWallace/phylowgs 
	def calcLogLikelihoodNegBinom(self, mu, tp, sd_d=0):
		return u.log_neg_binom_likelihood(self.a[tp], self.d[tp]-self.a[tp], mu, sd_d)