Участок для прогнозируемых вероятностей после команды blogit

Question

Я ищу способ построения прогнозируемых вероятностей двоичного результата после оценки модели логита для сгруппированных данных с использованием команды blogit в Stata. Сюжет должен показывать прогнозируемую вероятность на оси y и непрерывную переменную (дни) по оси x.

Я пробовал команду xi3 и postgr3 из UC LA, но они, похоже, не работают с командой blogit.

Справка действительно оценена! Пожалуйста, дайте мне знать, если вам нужна дополнительная информация.

Спасибо!

Answer 1

Если я вас правильно понимаю, один простое решение было бы сделать Контрафактно набор данных, в котором pop_var является 1, а затем просто использовать команду predict, чтобы получить вероятности, как следующее:

webuse xmpl2 
blogit deaths pop agecat exposed 
replace pop = 1 
predict probability 

Затем вы можете нарисовать график, используя различные команды twoway, однако вы хотите. Хотя мой пример не имеет сплошную ось х, это может быть что-то вроде:

twoway scatter probability agecat, by(exposed) 

Answer 2

Вы также можете соответствовать этому типу модели в рамках GLM, как это:

webuse xmpl2, clear 
list, clean noobs 
blogit deaths pop agecat exposed 
glm deaths agecat exposed, link(logit) family(binomial pop) nolog 
replace pop = 1 
margins, at(agecat = (0 1)) 
marginsplot