问题标签 [vgam]

这个

VGAM加载包时发生故障：

两者都实现s()了功能，但这不应该发生对吗？这是错误mgcv还是VGAM包装错误？

为什么添加第二个解释变量后拟合会变差？

我正在尝试将 tobit 模型从 Stata 迁移到 R。

稳健的 Stata 命令将只是添加,vce(robust)到模型中。而对于聚类它会是,vce(cluster idvar).

可重现的状态示例：

可重现的 R 示例：

我的理解是，这coeftest(m, vcov = sandwich)应该给我强大的 se。

但我得到以下信息：Error: $ operator not defined for this S4 class.

有人可以建议一种方法来估计来自 vglm 模型的稳健 se 以及使用 vglm 聚类 se 吗？

VGAM 版本 0.93

问题在这里：

这里log1p(exp(1000))变成Inf

因此，与正常工作的基础相比，它使用的数值方法不能处理大数字plogis。

是否值得提交错误报告，我在哪里可以提交？

我有一个矩阵Expr，其中行代表变量和列样本。我有一个分类向量groups（包含“A”、“B”或“C”）我想测试哪些变量 'Expr' 可以通过样本属于 a 的事实来解释group。

我的策略是使用广义加法模型（具有负二项分布）对问题进行建模。然后我想以变量明智的方式使用似然比测试来获得每个变量的 ap 值。我愿意：

接着：

最后一步是错误的，因为它测试的是两个模型的整体似然比，而不是变量。我想获得每个变量的 p 值向量，而不是单个 p 值。

我该怎么做？（我对 R 很陌生，所以请原谅我的愚蠢）

我正在使用序数显式拟合累积 logit 模型，得到的系数如下：

我想知道函数放在变量名称之后的 L 和 Q 代表什么。

先感谢您

我无法通过常规方法获得它们。

这里它指的是AIC.vglm，但它甚至不存在。

我找到的最接近的函数是AICrrvglmor AICvgam，两者都不起作用。

我正在使用包 VGAM 的帮助页面中的以下示例

我得到了一个合适的图表：

如何避免从图中绘制点？此外，我需要在 20、30、40...80 岁时打印出这些百分位数的值（单独作为表格）。如何才能做到这一点？我可以在这里使用 ggplot() 格式命令而不是 qtplot() 命令吗？谢谢你的帮助。

我正在从 VGAM 的帮助文件中运行以下代码：

我怎样才能防止在绘图上放置点，以便图表仅显示百分位曲线？qtplot 中是否有一些选项可以抑制该图上的点，这样我就不需要像本页那样使用长 ggplot 路线：来自 VGAM 的百分位数？在我之前的问题中，还有其他问题，所以这一点被忽略了。谢谢你的帮助。

TL;博士

我正在使用 R 中的 VGAM 包运行 Tobit 回归——这是一个玩具数据集，它一直给我一个我无法诊断的错误：

尝试诊断

带回溯：

我查看了源代码lm.wfit并找到了错误的来源：

我在源代码中找到了以下内容vglm：

该方法默认为vglm.fit.

我仍然无法找到use.i11创建排除标准的位置、它在做什么以及为什么它会导致权重、回归量和回归量之间的维度冲突。

我观察到将minSand舍入maxS到十个或更少的位置会导致成功运行，但这是因为maxS增加了所以第一次观察不再被右删失，并且minS增加了所以第二、第四和第五次观察不再被左删失。两者都改变​​了最大似然函数中的观察处理，所以我怀疑我会用错误的结果污染回归。

有人可以帮助诊断为什么会发生这种类型的错误吗？