Будучи относительно новичком в игре, я, вероятно, что-то упустил. Я пытаюсь смоделировать некоторые диаграммы роста и оценил различные варианты и нахожу (путем построения графиков, оценки GAI C и гиперпараметрической оптимизации), что следующая модель обеспечивает наилучшее соответствие:
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
Error in lm.wfit(X[onlydata, , drop = FALSE], y, w) :
incompatible dimensions
Есть ли способ извлечь z-показатели (или, в этом отношении, процентили) из новых данных с данной моделью?