使用时变泊松伽马模型非齐次 lambda进行注册预测

如何解决使用时变泊松伽马模型非齐次 lambda进行注册预测

我正在从事招生预测硬件项目，我正在尝试使用 Poisson-gamma 模型绘制招生曲线。我指的主要是 Anisimov 关于使用固定 lambda 进行不同持续时间的临床试验预测的论文。在接下来的步骤中，我将尝试实现随时间变化的非齐次 lambda。公式如下所示：enter image description here

其中 t_in,s 是启动/开始时间。我不知道如何在 r 中编写此函数并将其集成以查找入组患者的平均数量。以前，对于固定的 lambda，我使用了以下内容：enter image description here

注册遵循泊松，率 lambda 遵循均值“m”的伽马分布。 (t-ai) 是持续时间。任何帮助表示赞赏。我正在尝试在 R studio 中实现这一点。谢谢。