Я довольно новичок в R (хочу переключиться с Stata), и пока мне нравится!
Однако я столкнулся с проблемой, для которой не могу найти решение.
Я пытаюсь построить прогнозируемые значения зависимой переменной на основе конкретных значений одной независимой переменной после запуска модели Heckit (двухэтапная Хекмана) из пакета sampleSelection.
Общая модель выглядит так:
Mod1 <- Heckit(selection = binaryindicator ~ x1 + z , outcome = y ~ x1, df, method)
где структура данных выглядит следующим образом:
z x1 y binaryindicator
0.5 4 300 1
0.1 8 400 1
0.2 10 500 1
0.2 18 NA 0
0.4 20 50 1
0.3 30 1000 1
После заливки документации для sampleSelection я могу только извлечь прогнозы из функции предиката.выбор и получить соответствующие остатки. Однако этот метод не позволяет мне определять прогнозы для одной переменной (скажем, x1), и полученные прогнозы взяты из всей модели.
Поэтому я ищу способ (постоянная функция или руководство) построить прогноз по y из x1 с доверительными интервалами. Все попытки сделать это (используя plot, cplot, gpredict, effect_plot) потерпели неудачу. Все вышеперечисленные функции работают с линейными моделями, но не с моделью выбора.
В идеале я хотел бы также указать, из какого уравнения (исход или выбор) я хочу получить прогноз, но это также ускользнуло от меня, поскольку я не могу получить никаких прогнозов для начала. Есть ли способ добиться этого? В Stata это можно сделать с помощью команды margins, но эквивалентный пакет полей в R не работает для моделей sampleSelection (для меня пока). Любая помощь будет высоко ценится!