Какой пакет R может взаимодействовать между сглаживанием и контролем родословного эффекта?
Я застрял на несколько недель. Я пытаюсь смоделировать набор данных, в который я хочу включить взаимодействия между сглаживаниями [как в пакете GAM: s(a, b, by = c) или ti(a, b, by = c)] и одновременно контролировать влияние родословной животных из набора данных (т.е. мать и отец каждого животного на протяжении многих поколений). Используя пакет GAM, я могу приспособить взаимодействия между сглаживаниями, но не могу найти способ включить структуру корреляции родословной каждого животного. Может быть, это можно сделать с помощью GAM, но я не нашел, как это сделать в документации GAM или где-то еще. MCMCglmm может работать с родословной, а также со сглаживанием с помощью функции spl () (сплайн с тонкими пластинами) и альтернативной версии spl2 () ( https://stat.ethz.ch/pipermail/r-sig-mixed-models/2014q1/021810.html), однако я полагаю, что эти две функции не рассчитаны на то, чтобы соответствовать взаимодействиям между гладкими (как те, что я описал выше), или, по крайней мере, подходящим образом для нестатиков вроде меня. Кто-нибудь знает, как решить эту проблему, используя GAM, MCMCglmm или любой другой пакет R? Хорошего дня