小编tor*_*gos的帖子

使用scipy.stats计算条件期望

假设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的误解,将不胜感激。

statistics scipy

3
推荐指数
1
解决办法
345
查看次数

标签 统计

scipy ×1

statistics ×1