问题标签 [gamlss]

我正在尝试使用emmeans库计算gamlss对象的估计边际均值。我的gamlss对象来自零膨胀的 beta 回归。假设我的模型称为 m1 并且我的变量之一是internationaltreaty，所以我调用：

emmeans(m1,"internationaltreaty",type="response")

我收到以下错误消息：

如果我使用不同的模型对象（例如 glm），emmeans 可以使用此代码。对我来说，emmeans 似乎无法识别我的类型参数。有没有人经历过类似的事情？

谢谢。

我正在通过随机效应在 R 建模中运行 GAMLSS。但是，会出现此错误：

solve.default(pdMatrix(a, factor = TRUE)) 中的错误：系统在计算上是奇异的：倒数条件数 = 0

我已经尝试过没有随机效应的模型并且它有效。

我是 r 中 GAMLSS 函数的新手，当我运行我的代码时，我总是收到这个错误：Response Variable out of range

查看数据框后，我意识到问题是响应变量之一是 0.0000。

我想知道是否有人可以向我解释为什么 0 超出范围以及可能的解决方案（例如替换值）？

如何对随机效应 beta 分布进行事后检验？

为了构建具有 beta 分布的混合模型，我使用了库 gamlss（我发现在 R 中没有其他方法可以做到这一点）。在这个模型中，有一个具有三个水平的因子。假设 level1、level2 和 level3。摘要显示了 level1 x level2 和 level1 x level3 之间的比较，因为 level1 是基线。

library(gamlss) mymodel <- gamlss(response~ph+dose, family = BE, random=~1|id, data = mydata) 摘要(mymodel)

它显示了 level1 x level2 和 level1 x level3 之间的比较

对于所有三个级别，摘要显示 level1 x level2 和 level1 x level3 之间的比较，因为 level1 是基线。

我想测试level2 x level3之间是否有区别。我怎样才能做到这一点？也许可能有另一个图书馆。谁能帮帮我吗？

我希望为一组线性模型制作一个汇总表。模型是 lmer()、glmer() 或 gamlss()。我正在尝试将其中 6 个模型的结果打印到一个表中。但是，当我尝试使用 sjPlot::tab_model 执行此操作时，我收到以下错误消息，“错误：$ 运算符对原子向量无效”。

我删除了两个 gamlss() 模型以查看表格是否会打印并且确实如此。在我看来，问题在于 sjPlot::tab_model 不能与 gamlss() 模型一起使用。

数据库

分析

表代码

有没有办法使用 sjPlot::tab_model 来获取带有 gamlss 模型对象的 html 表输出，或者是否有其他软件包可以推荐用于为 gamlss 以及 lmer 和 glmer 的线性模型对象制作出版质量表？谢谢！

我的y变量 ( n=30,000) 分布有非常重的尾部（正面和负面），fitDistGAMLSS 函数为此选择 ST4 系列。

我尝试使用解释变量x（pb 平滑）来评估基于 GAMLSS 的回归，但尾部y太重以至于在 50 个周期后收敛不会达到，即使在 refit ( time consuming+++) 之后也是如此。

因此，我y使用orderNorm转换（bestNormalize包）进行归一化，它可以轻松快速地达到收敛，然后从 GAMLSS 对象预测拟合值。

然而，这些拟合的“ orderNormalized”值是一个 GAMLSS 对象，因此不能使用预测函数从bestNormalize（因为后者似乎无法识别 GAMLSS 对象）来反转。

我的问题：无论采用何种方式，是否有可能对orderNorm来自 GAMLSS 对象的拟合值应用逆变换？

我收到ERROR : NA's in the working vector or weights for parameter mu gamlss 函数的以下错误：fit <- gamlss::gamlss(x ~ y, sigma.formula = ~ y, data = d, family = GA)

其中 d 是包含 x 和 y 的数据框

这个错误的原因是什么？

我目前正在尝试为心血管研究项目生成列线图和 z 分数。

但是，我对gamlss包有疑问，尤其是centiles.pred()函数。

例如，如果我运行以下代码，我会陷入错误：

当我运行最后一行时，我收到以下消息：

我该如何克服这个错误？

作为一个在 gamlss 上相对较新的人，我可能会遗漏一些东西。我正在尝试对一些增长图表进行建模并评估不同的选项并发现（通过绘图、GAIC 评估和超参数优化）以下模型提供了最佳拟合：

现在鉴于这是一个增长图表，我希望能够为给定的任何新数据点计算 z 分数，所以我试试这个：

这两个都给出了相同的错误，这似乎与这不是一个“简单”的 lm 拟合有关

lm.wfit (X[onlydata, , drop = FALSE], y, w) 中的错误：尺寸不兼容

有没有办法从给定模型的新数据中提取 z 分数（或百分位数）？

我gamlss在 R 中遇到了一些令人困惑的行为。文档列出了 predict.gamlss 并?predict.gamlss返回 RStudio 中的函数文档。但是，该函数在输入时不会自动完成predict.gamlss，并且尝试运行它会返回Error: 'predict.gamlss' is not an exported object from 'namespace:gamlss'。这是怎么发生的？该功能是否以某种方式停用？有一个单独的功能predictAll可以工作。

该文档确实说明了

此功能正在开发中

我正在尝试访问该功能，因为我在使用 和 时遇到了一些令人困惑的predict结果predictAll。

R 版本是 4.0.0。gamlss版本是 5.1.6。