R中泊松过程的回归

机器算法验证 r 回归 软件 泊松过程
2022-04-03 17:41:51

我有一系列不同长度的样本,以及在这些时间样本中创建的错误数量。阅读文献,这通常被建模为泊松过程。如果你这样写:

P(k,t)=eλt(λt)kk!

我知道kt并想回归寻找λ.

是否有一个 R 包可以做到这一点?该标准glm适用于泊松回归,但我找不到任何允许样本具有非单位时间的东西。

1个回答

您仍然可以glm在 R 中使用,但您包含的日志t作为考虑它的“偏移量”,例如:

fit <- glm( k ~ 1 + offset(log(t)), data=mydata, family=poisson)

这将适合一个截距,该截距将是λ,但如果需要,您也可以包含协变量。