假设x〜Poisson(2.5); 我想计算类似E(x | x> 2)的值。
我假设可以使用.dist.expect运算符完成此操作,即:
D = stats.poisson(2.5)
cond_expect = D.dist.expect(lambda x: x, D.args,lb=2)
Run Code Online (Sandbox Code Playgroud)
计算结果为cond_expect = 2.29478750344
但是,如果我只是根据该分布计算随机样本的平均值
D = stats.poisson(2.5)
test = D.rvs(size = 100000)
empirical_expectation = np.mean(test[test>=2])
Run Code Online (Sandbox Code Playgroud)
empirical_expectation的计算结果为3.20875563063。
如果有人可以澄清我对API的误解,将不胜感激。