问题标签 [gamlss]

我正在使用 gamlss 包中的示例代码来绘制百分位曲线：

它正在绘制自己的一组百分位曲线。如何选择仅绘制第 10、第 50 和第 90 个百分位曲线？另外我想避免绘制点，以便只绘制曲线。谢谢你的帮助。

这个问题与：Selecting Percentile curve using gamlss::lms in R

我可以从以下数据和代码中得到百分位曲线：

我现在如何获得曲线上每个 x 值 (5:15) 的 yvar，这将代表平滑后数据的第 90 个百分位数？

我试图阅读帮助页面并找到 fit(h) 和 fv(h) 来获取整个数据的拟合值。但是如何在第 90 个百分位曲线水平上获得每个年龄水平的值？谢谢你的帮助。

编辑：下图显示我需要什么：

我尝试了以下但它是正确的，因为值不正确：

这些值与直接来自数据的第 90 个分位数有很大不同：

我想从 gamlss R 函数预测二进制类概率/类标签，如何使用预测函数来获取它们？

我有以下示例代码

predict.gamlss(fit, newdata = data.frame(X1, X2)[testing, ], 中的错误：使用选项 data 定义原始数据

任何想法？

上下文：tidyverse和dplyr环境/工作流程。

我希望能深入了解如何解决以下问题，这是我在尝试使用回归结果集合时遇到的。

这个最小的可重复性显示了问题

下面的示例说明了更广泛的上下文并按预期工作（生成显示的图像）。如果我所做的只是更改~lm(...)为~glm(...) or ，它也可以工作~gam(...)：

但是，如果我尝试~gamlss(...)如下使用：

我观察到以下错误：

是否有人建议需要更改哪些内容才能使此示例按预期工作？

我会很感激任何关于出了什么问题的见解。为什么它不起作用。如何诊断此类问题。

我正在努力将gamlss结果收集到数据框中。这将继续此处的示例

工作示例使用lm

破碎的例子使用gamlss

感谢任何提示或提示。

我正在尝试获得留一的预测值。请帮我解决这个“找不到对象”的问题。我已经搜索了类似的问题，但没有设法弄清楚。这是在 Windows 10 上。

提前致谢

我想估计候选人的名字对她获得的负面偏好票的比例的影响。我的因变量“负面偏好投票”衡量候选人从政党选票中划掉的百分比。变量在 [0,1] 范围内。我已经使用 GAMLSS 包来适应 Beta-zero-one 膨胀分布 (BEINF)。我的观察还集中在市政当局和候选人中——但数据结构是交叉嵌套的。我想在模型中添加交叉分类的随机效应，但我还没有找到任何关于如何在 GAMLSS 中建模的文档。

据我所知，有两个函数可以在 GAMLSS 模型中拟合随机效应，random() 和 re()。我使用了 re() 函数并估计了具有两个随机效应的以下模型，但我怀疑这是否指定了交叉分类的随机效应。

有没有办法在 GAMLSS 中指定交叉分类的随机效应？

我一直在处理这些数据

year rango_edad Sexo zona_2016 conteo siniestros expuestos upc_millon valor_millon freq 1 2010 01. < 1 F Alejada 180 87 75 121 111 0.48 2 2010 01. < 1 F Ciudades 103453 76219 40228 60755 84981 0.74 3 2010 01. < 1 F Especial 5129 3194 2078 3289 3013 0.62 4 2010 01. < 1 F Normal 27393 18436 10735 15656 16692 0.67 5 2010 01. < 1 M Alejada 185 98 73 116 110 0.53 6 2010 01. < 1 M Ciudades 106915 80731 41719 62991 105135 0.76 costo.medio ratio 1 1.27 0.92 2 1.11 1.40 3 0.94 0.92 4 0.91 1.07 5 1.12 0.94 6 1.30 1.67 我正在尝试用 gamlss 对频率进行建模

gamlss(freq~Sexo+zona_2016+rango_edad,family=PO(mu.link = "log"),data=na.omit(subset(datos,is.na(freq)==FALSE ))) gamlss(freq~Sexo+zona_2016+rango_edad,family=NBI(mu.link = "log"),data=na.omit(subset(datos,is.na(freq)==FALSE )))

但我收到了这个错误信息

Error in while (abs(G.dev.old - G.dev) > c.crit && iter < n.cyc) { : missing value where TRUE/FALSE needed

我该如何解决？

我正在尝试执行 gamlss 的零膨胀 beta 回归模型，如下所示>它适用于变量 V255，但通过 simialr 变量 V256 的错误

有什么建议可以解决这个问题吗？代码和模型。

我已经建立了一个 GAMLSS 模型来基于现有数据集调查鱼类的种群驱动因素，但是在尝试预测新数据集的响应值时出现了问题。

要创建一个虚拟数据集：

使用 GAMLSS 模型：

当我在predict()现有数据集上运行该函数时，我能够毫无问题地获得所有模型参数。但是，当我尝试使用新数据集 ( MC.df.newdata <- data.frame(Datej=c(1:100), Year=2010, Turbidity=1.959954, Site="Retallacks Hole")) 获取预测响应值时，我收到以下错误消息：

我使用了上面所有相同的代码，+random(Site)一切正常，所以我的问题是如何让predict()函数在模型中使用随机效应？