Я регулярно запускаю свой скрипт со статистикой, каждый раз подгоняю модель со случайными эффектами и проверяю значимость предиктора. Каждый раз модель показывает, что ее нет. Я подгонял уже раз 150 и 150 раз значимости не было. Смысл в том, что скрипт и данные всегда одни и те же. А я каждый раз с неизменной надеждой подгоняю снова:
- Эффект есть?
- Эффекта нет.

- Значим ли предиктор?
- Предиктор не значим.

- Изменяется ли шанс исхода в зависимости от регрессора?
- Нет, не изменяется.

- Есть ли статистическая значимость?
- p-value > 0.05.

И ведь этот скрипт, зараза, знает, что я подгоню модель и знает, что она рассчитает. Но ещё ни разу ни один из нас ни предупреждением, ни ошибкой не показал, что каждый из нас знает сценарий. Бывает идешь ставить ПЦР, а комп стоит без дела, тогда я запускаю расчет, что бы просто так не простаивал, он равнодушно подгоняет модель, рассчитывает p-value и печатает результаты:
- Generalized linear mixed-model: fitting.
- Значимо?
- OR 1.01, p-value 0.837.
- Жаль.

Иногда я читаю статьи, тогда я запускаю расчет параллельно и терпеливо жду, когда он закончится:
- glmm(outcomeVar ~ ., data = mydata, family = binomial)
- OR 1.0, p-value 0.936.

Это очень суровое, по-настоящему мужское противостояние, исход которого не ясен. Очевидно, что каждая сторона рассчитывает на победу. Впрочем, я уже теряю надежду.