#r
#r
Вопрос:
Я выполняю анализ логистической регрессии / выживаемости, в котором я группирую стандартные ошибки по переменной в наборе данных. Я использую R.
Поскольку это не так просто, как в STATA, я использую решение, которое я нашел в прошлом : https://www.rdocumentation.org/packages/miceadds/versions/3.0-16/topics/lm.cluster
В качестве наглядного примера того, о чем я говорю:
model <- miceadds::glm.cluster(data = data, formula = outcome ~ a b c years years^2 years^3, cluster = "cluster.id", family = "binomial")
Это хорошо работает для получения важных значений, при этом создаются коэффициенты, std. ошибки (кластеризованные) и z-значения. Мне потребовалась вечность, чтобы просто найти это решение; и даже сейчас оно не идеально (например, невозможность вывода в Stargazer). Я изучил множество других распространенных предложений по этому вопросу, таких как решение экономической теории (https://economictheoryblog.com/2016/12/13/clustered-standard-errors-in-r /); однако это для lm(), и я не могу заставить его работать для логистической регрессии.
Я не ограничиваюсь простым запуском двух моделей, одной с glm () и одной с glm.cluster() и заменой стандартных ошибок в stargazer вручную.
Меня беспокоит то, что я не знаю, как бы я построил график вышеупомянутой функции, скажем, если бы я вместо этого сделал следующее:
model <- miceadds::glm.cluster(data = data, formula = outcome ~ a*b c years years^2 years^3, cluster = "cluster.id", family = "binomial")
В этом случае я хочу построить график прогнозируемой вероятности, чтобы посмотреть на взаимодействие между a * b в моем результате; однако я не могу этого сделать с объектом glm.cluster(). Я должен сделать это с помощью модели glm (), но тогда мои доверительные интервалы переполнены.
Я рассматривал множество вариантов кластеризации стандартных ошибок для логистической регрессии здесь, но нахожусь в полной растерянности.
Кто-нибудь нашел какие-либо недавние разработки о том, как это сделать в r? Существуют ли какие-либо пакеты, которые позволяют кластеризировать SE по переменной в наборе данных и отображать объекты? (Бонусные баллы за взаимодействия)
Любая информация была бы оценена. Спасибо!