Я помогаю коллеге подобрать обобщенно-пуассоновскую обобщенную линейную смешанную модель в R, используя функцию cpglmm из пакета cplm ( ссылка ).Модель предполагает трехстороннее взаимодействие, и я хотел бы вычислить некоторые интерпретируемые величины.До сих пор я пытался вычислить некоторые отношения шансов, но я не уверен, что это правильный способ сделать это?
# Fit model with three-way interaction in fixed effects #
m <- cpglmm(ncs ~ diversity_index*diversity_speciality*n_authors + selfcit +
n_refs + (1|region), data = diversity)
# Calculate Odds-ratio #
se <- sqrt(diag(vcov(m)))
tab <- cbind(Est = m$fixef,
S.E. = se,
LL = m$fixef - 1.96 * se,
UL = m$fixef + 1.96 * se)
print(exp(tab), digits=3)
Я также хочу вычислить некоторые прогнозируемые значения, например, прогнозируемые вероятности или тому подобное, но я не могу заставить predict() работать на cpglmm.Могу ли я использовать какие-либо функции?