作为一个在 gamlss 上相对较新的人，我可能会遗漏一些东西。我正在尝试对一些增长图表进行建模并评估不同的选项并发现（通过绘图、GAIC 评估和超参数优化）以下模型提供了最佳拟合：
mod_fp_poly2_gg <- gamlss(value~fp(age,npoly = 2),sigma.formula = ~fp(age,npoly = 2), family = GG, data = data)
现在鉴于这是一个增长图表，我希望能够为给定的任何新数据点计算 z 分数，所以我试试这个：
newdata <- data.frame(age = 1000, value = 600000, group = "new")
predictAll(mod_fp_gg2, newdata = newdata)
# as well as
centiles.pred(mod_fp_bct, xname = "age", xvalues = 1000, yval = 600000, type="z-scores")
这两个都给出了相同的错误，这似乎与这不是一个“简单”的 lm 拟合有关
lm.wfit (X[onlydata, , drop = FALSE], y, w) 中的错误：尺寸不兼容
有没有办法从给定模型的新数据中提取 z 分数（或百分位数）？