Я хотел бы использовать фиксированные эффекты ранее в blme :: bglmer () для одного из моих ковариат, показывающих полное разделение, следуя указаниям предыдущего поста: Биномиальный блеск с категориальной переменной с полнымуспехи .Однако в моем случае есть дополнительный ковариат, не показывающий полного разделения.
Ниже приведен обзор моего набора данных (RAMP_Subset).ID идентифицирует каждую наблюдаемую рыбу, выборка идентифицирует каждую отобранную выборку (на промысловом судне, с несколькими рыбами из одного и того же вылова), смерть - это моя переменная реакции с 2 уровнями (смертность на асимптоте, 0 для выживания и 1 для мертвых), R1 иR2 - мои ковариаты, каждый из которых имеет 2 уровня (рефлекс для каждой рыбы, 0 при рефлексе и 1 при отсутствии рефлекса):
ID haul death R1 R2
2704 9 1 1 1
2705 14 1 1 1
2715 13 1 1 1
2718 13 1 1 1
2719 8 1 1 1
2722 9 1 1 1
...
Скажем, R1 - ковариата с полным разделениемдля которого я хотел бы добавить фиксированные эффекты ранее.Поэтому моя модель будет выглядеть следующим образом:
bglmer(death ~ R1 + (1|haul) + (1|ID), data=RAMP_Subset, family=binomial, fixef.prior = normal(cov = diag(9,2)))
Теперь предположим, что мой второй ковариат R2 не показывает полного разделения.Справочное руководство для пакета R "blme" гласит, что априоры для фиксированных эффектов определены так же, как и для ковариационных априоров, поэтому я попробовал следующее:
bglmer(death ~ R1 + R2 + (1|haul) + (1|ID), data=RAMP_Subset, family=binomial, fixef.prior = R1 ~ normal(cov = diag(9,2)))
Тем не менее, моя попытка получила сообщение об ошибке:
Ошибка вvaluFixefPrior (fixefPrior, defnEnv, evalEnv): нераспознанное распределение fixef: '~'
Я не смог найти ни одногодругой пример онлайн, и был бы очень признателен за совет по правильному форматированию.
Спасибо, Эстер