poisson相关内容

确定飞机到达之间的等待时间

我正在写一个程序,它需要200架分布在泊松分布的飞机在12个小时内,这些飞机需要降落在一个只有一条跑道的机场。利用指数分布的反向CDF方法确定到达间隔时间。然而,我似乎无法计算出在空中等待的时间。 例如,一架飞机100秒到达,75秒着陆,175秒完成。 飞机2在150秒到达,必须等待175-150=25秒。考虑到等待时间也可以是0秒,我如何编程我的函数可以输出的变量? import m ..
发布时间:2022-04-06 22:47:25 Python

手工泊松回归

我想手工进行泊松回归,并定义一个函数,该函数可用于估计任意数量的系数。我有两个问题: 首先:如何才能获得Beta矩阵,而不必编写每个Beta显式。我想这样编写lambda=exp(t(X)%*%beta)。我以为我可以做一个for循环,为x、a、beta中的每一列创建一个for,然后在矩阵中将它们相加ab,但我不知道如何编码。 第二个: 因为我不知道怎么写,所以我试着写了估算6个贝塔的 ..
发布时间:2022-02-26 18:53:26 其他开发

Python:endog 已评估为具有多列的数组

我正在尝试运行泊松模型,如下所示: poisson_model_xg = smf.glm(公式=“xG~主场+球队+对手", data=xg_model_data,family=sm.families.Poisson()).fit() 我收到以下错误: ValueError:endog 已评估为包含多个列的数组具有形状 (760, 9).当变量转换为 endg 时会发生这种情况是非数字的 ..
发布时间:2021-09-11 19:44:21 Python

计算泊松概率百分比

当您在 Excel(或 OpenOffice Calc)中使用 POISSON 函数时,它需要两个参数: 一个整数 “平均"数字 并返回一个浮点数. 在 Python 中(我尝试过 RandomArray 和 NumPy),它返回一个随机泊松数数组.我真正想要的是这个事件发生的百分比(它是一个常数,数组每次都有不同的数字 - 那么它是平均值吗?). 例如: 打印泊松(2 ..
发布时间:2021-08-30 18:43:08 Python

具有上限的 Scipy 泊松分布

我正在使用 scipy stats 生成一个随机数.我使用了泊松分布.下面是一个例子: 将 scipy.stats 导入为 sctA = 2.5Pos = sct.poisson.rvs(A,size = 20) 当我打印 Pos 时,我得到以下数字: array([1, 3, 2, 3, 1, 2, 1, 2, 2, 3, 6, 0, 0, 4, 0, 1, 1, 3, 1, 5]) 从 ..
发布时间:2021-07-16 21:07:36 Python

Python SciPy chisquare 测试从 Excel 和 LibreOffice 返回不同的 p 值

阅读最近的博客文章 关于泊松分布的应用,我尝试使用 Python 的“scipy.stats"模块以及 Excel/LibreOffice 的“POISSON"和“CHITEST"重现其发现功能. 对于文章中显示的期望值,我简单地使用了: 导入scipy.stats对于范围内的 i (8):打印(scipy.stats.poisson.pmf(i,2)* 31) 这重现了博客文章中显示的 ..
发布时间:2021-07-16 20:41:17 其他开发

如何用seaborn拟合泊松分布?

我尝试将我的数据拟合为泊松分布: 导入 seaborn 为 sns导入 scipy.stats 作为统计信息sns.distplot(x, kde = False, fit = stats.poisson) 但我收到此错误: AttributeError: 'poisson_gen' 对象没有属性 'fit' 其他分布(伽马等)运行良好. 解决方案 泊松分布 (在 scip ..
发布时间:2021-07-16 20:33:03 其他开发

如何在 Python 中将数据帧的一列拟合为泊松分布

我一直在尝试找到一种方法,使我的一些列(包含用户点击数据)适合 python 中的泊松分布.这些列(例如,click_website_1、click_website_2)可能包含从 1 到数千的值.我正在尝试这样做,因为某些 资源: 我们建议不要分析计数数据对它进行对数转换,而是基于泊松和负数的模型应该使用二项分布. 我在 scipy 和 numpy 中找到了一些方法,但是这些方法似乎 ..
发布时间:2021-07-16 20:13:31 Python

zeroinfl“系统在计算上是奇异的"而在预测变量中没有相关性

我正在尝试对一年中工人缺勤天数的计数数据进行建模(因变量).我有一组预测变量,包括关于工人的信息、关于他们的工作等......,其中大部分是分类变量.因此,有大量系数需要估计(83),但由于我有超过 600 000 行,我认为这应该没有问题.此外,我的数据集中没有缺失值. 我的因变量包含很多零值,所以我想估计一个零膨胀模型(泊松或负二项式),使用 pscl 的函数 zeroinfl包,代码: ..
发布时间:2021-07-07 18:51:45 其他开发

如何在python中使用random.expovariate()生成整数到达时间

在 python random 模块中,expovariate() 函数生成浮点数,可用于对泊松过程的到达间隔时间进行建模.我如何利用它来生成到达之间的整数时间而不是浮点数? 解决方案 jonrsharpe 已经有点提到了,你可以让函数生成浮点数,并使用 int() 这个 >>>随机导入>>>[random.expovariate(0.2) for i in range(10)][ ..
发布时间:2021-07-02 19:50:49 Python

Statsmodels Poisson glm与R不同

我正在尝试根据R中提供的某些代码来拟合某些模型(空间交互模型).我已经能够在python框架中使用statsmodels使某些代码工作.完全没有.我相信我为R和Python使用的代码应能得到相同的结果.有没有人看到任何差异?还是有一些根本的差异可能会导致事情失败?R代码是与教程中给出的数字相匹配的原始代码(位于此处:http://www.bartlett.ucl.ac.uk/casa/pdf/pa ..
发布时间:2021-05-11 18:58:22 Python

使用带有矩形域的FFTW的泊松方程

我试图通过使用具有矩形域(-4 <= x <= 4和-2 <= y <= 2)的FFTW库来求解泊松方程.如果域是正方形,我有正确的结果,如果域是矩形,那是错误的.请给我一些建议.非常感谢.这是我的代码. #include#include#include#include#include#includ ..
发布时间:2021-04-19 19:50:21 C/C++开发

向图像添加泊松噪声

我有一些图像需要添加更多的泊松噪声,以便对其进行更彻底的分析.我知道您可以在MATLAB中执行此操作,但是如何在Python中执行此操作呢?到目前为止,搜索没有任何结果. 解决方案 如果您可以使用numpy/scipy,则以下内容应会有所帮助.我建议您将数组转换为浮点数以进行中间计算,然后将其转换回uint8以进行输出/显示.由于泊松噪声都> = 0,因此您需要决定在转换回uint8时如何 ..
发布时间:2020-11-26 22:38:54 Python

检查R中的GLM摘要的解释

只想检查我在做什么都是正确的! 我在几个分类为两个栖息地的地点(农田和湿地)都有鸟类计数.我只是想看看哪个栖息地的数量更高. 我正在使用具有Poisson函数的GLM(因为它们是计数数据): > mod summary(mod) Call: glm(formula = count ~ ..
发布时间:2020-11-16 00:23:30 其他开发

R中Rate变量的回归

我的任务是开发一个回归模型,研究不同课程的学生入学情况.这是一个非常好的,干净的数据集,其注册计数很好地遵循了Poisson分布.我在R中拟合了模型(同时使用GLM和零膨胀泊松.)产生的残差似乎是合理的. 但是,然后我被指示将学生人数更改为“费率",该费率是按学生/school_population计算的(每所学校都有自己的人口.)现在,这不再是一个计数变量,而是一个介于0和1.这被视为程序 ..
发布时间:2020-11-16 00:12:49 其他开发

从R中的零膨胀分布中随机选择值

您好,谢谢您的帮助! 前一段时间,我问了 我在上一个问题中创建的向量是二进制的,现在我想生成一个加权向量(即与有界整数)。我从具有长尾巴的零膨胀或拟泊松分布中进行采样,因此选择零的可能性比另一个值高得多,但是选择大值(例如63)的可能性有限。 / p> 我可以使用rpois从泊松分布中选择值并创建给定长度的向量。这与我想做的类似,因此我将以它为例。 e = seq(0: ..
发布时间:2020-10-22 19:35:20 其他开发