问题标签 [pscl]

我真的很喜欢政治和选举，我刚刚开始学习 R，我想用来自当地县的新数据重新创建这个博客中列出的过程。在进行区域分析之前，我已经能够使用修改后的代码可靠地完成博客中的大部分过程。

作者分析了特定的住宅区，而我更愿意分析整个县。我修改了代码以使用 US House 作为我的目标区域，因为它包含整个县。

我想知道是否有人对我为什么无法从该县数据中获取辖区级别摘要提出建议。我收到这样的错误：

当我在数据中有“NA”时，我只会收到这个错误。当我用“0”代替空白时， District.analyze 起作用，但是“0”会抛出所有方程。

我可以重现的最少代码量是：

我的数据是我从 .csv 读入 R 的大型数据集：

我们正在尝试使用零膨胀泊松（在 pscl 包中实现）对具有过多零的计数变量进行建模。这是一个（简化的）输出，显示了分类和连续解释变量：

summary 给了我们一些很好的答案，但我们正在寻找一个类似 ANOVA 的表格。所以，问题是：可以使用 car::Anova 来获取这样的表吗？

它似乎工作正常，但我不确定是否是一种有效的方法，因为缺少文档（似乎只考虑“计数模型”部分？）。您建议遵循这种方法还是有更好的方法？

我正在使用函数ideal和rollcallfrom psclpackage 但相同的选民根据数据集的顺序（在列表中）接收交替的负号值或正号值。当我计算几个时期时，我需要它的行为每年都保持一致。有没有办法控制这方面？

我正在尝试针对 3 种不同的效果大小运行偏差校正的百分位引导程序，并从前任中运行一些代码。虽然此代码有效，但运行时间太长了。完成 3 种不同的效果大小中的每一种大约需要 3 天。

它还显示了我不需要的 1000 次迭代中的每一次。理想情况下，我想最大限度地优化运行时间。我还希望控制台上的所有 3 个输出或 Word/Excel 文档中的所有输出，而不必一次运行它们。

我知道代码有点笨拙，但我们将不胜感激。

我对 R 相当陌生，并且正在寻找有关如何适当进行的一些建议。代码如下：

在我看来，每个效果单独运行的问题来自于将每个循环的结果命名为“结果”。我不知道在不干扰循环的情况下打印所有结果（对于每个效果大小）的最佳方法。

长度显然来自于大量的迭代和缺乏足够快的内存来处理它们。这是可以补救的吗？当程序运行时，时间并不像没有得到所有效果大小的结果那样令人担忧。

我正在使用 R 版本 3.6.0，并且我正在遵循本指南以便对霍乱数据运行零膨胀泊松回归：https ://stats.idre.ucla.edu/r/dae/zip/

我已经按照说明安装了 pscl 包，但是当我运行我的代码时，会显示一条错误消息。

我尝试使用不同版本的 R 并且它有同样的问题。

错误消息：zeroinfl 中的错误（cases ~ wateraccess + atpc，data = cholera2）：找不到函数“zeroinfl”

我的数据是零膨胀的，所以我正在使用 glmmamdb 运行一个零膨胀模型：

Model3z <- glmmadmb(Count3 ~ Light3 + (1|Site3), zeroInflation = T, family= "poisson", data = dframe3)

但是，当我尝试对此模型中的不同灯光类型进行成对比较时pwcs3 <- lsmeans(Model3z, "Light")我收到错误消息：

ref_grid(object, ...) 中的错误：无法处理“glmmadmb”类的对象 使用 help("models", package = "emmeans") 获取有关受支持模型的信息。

当我访问 emmeans 软件包网站时，它说 glmmadmb 不再受支持。

我已切换到 pscl 和 zeroinfl 函数，但不确定如何重组我的代码以适应 pscl 格式。键入P <- zeroinfl(Count3 ~ Light3 + (1|Site3), family = poisson, data = dframe3)得到错误消息：

contrasts<-( , value = contr.funs *tmp*[1 + isOF[nn]]) 中的错误：对比只能应用于具有 2 个或更多级别的因子此外：警告消息：在 Ops.factor(1, Site3) 中：'|' 对因子没有意义

是否有另一种方法可以将 glmmadmb 与 lsmeans 一起使用？如果没有，有谁知道 pscl 中的零膨胀模型代码应该是什么样子？谢谢。

我正在尝试对一些数据运行零膨胀负二项式计数模型，这些数据包含按县划分的政治家的竞选访问次数。（对数似然检验表明负二项式是正确的，Vuong 检验表明零膨胀，尽管这可能会因为我的零膨胀模型显然没有收敛的事实而被抛弃。）我在 R 中使用 pscl 包。问题是当我跑步时

显然，所有这些 NA 对我的帮助都不大。任何建议将不胜感激！我是 R、StackOverflow 和 Statistics 的新手，但我正在努力学习。我正在尝试提供最小可重现示例所需的一切，但我看不到任何地方可以分享我的实际数据......所以如果你需要这些东西来回答这个问题，请告诉我我可以把它放在哪里！

我正在使用“pscl”包来构建障碍模型，以便将鸟类的丰度建模为环境变量的函数。我有兴趣预测非采样区域的发生概率和丰度，以便为该物种生成种群估计。我能够生成模型和预测，但我很难为这些估计生成可信区间。

我正在使用包中的hurdle()函数pscl在 R 中运行负二项式障碍模型。

我想只从模型输出中提取 z 值列表并将这些值分配给一个向量。但我不知道该怎么做？

这是我的模型的基本轮廓：

例如，我尝试过：

coef(summary(real.mod))[,3]

但这不起作用。

任何有关如何执行此操作的指导将不胜感激。