Тест Вальда для сравнения двух моделей в Stata

Я запускаю логистическую регрессию с результатом (ser) и рядом переменных-предикторов, со стандартными ошибками, кластеризованными (by provider_id). Одним из предикторов (years_cat) является категориальная переменная, которую я создал из непрерывной переменной (срок владения); другими словами, years_cat - это категории владения недвижимостью. Коэффициенты шансов для категорий years_cat монотонно растут, и я хотел бы проверить, является ли переменная владения фактически линейной (хотя я буду сообщать о категориальной переменной в моей статье).

Обычно я запускаю модель с категориальной переменной (cat1), сохраняю оценки, а затем запускаю модель с линейной переменной (lin1), сохраняю оценки и затем использую тест отношения правдоподобия для сравнения двух (lrtest cat1 lin1). Поскольку у меня есть кластеризованные стандартные ошибки, тест отношения правдоподобия не подходит, и я хотел бы использовать тест Вальда. Однако, когда я запускаю тест Уолда, я получаю ошибки (подробности см. Ниже), что Stata не может найти модели, которые я сохранил (несмотря на подтверждение того, что они на самом деле хранятся). Как я могу запустить тест Вальда, чтобы сравнить равенство этих двух хранимых моделей?

*Run the model with categorical variable (years_cat), store.
logistic ser i.age_ra i.sex ib5.indication ib2.gender_doc ib1.specialty1 ///
    i.vol ib4.years_cat, cluster(provider_id)
estimates store cat1

*Run the model with continuous variable (tenure), store.
logistic ser i.age_ra i.sex ib5.indication ib2.gender_doc ib1.specialty1 ///        
    i.vol tenure, cluster(provider_id)
estimates store lin1

*Confirm that stata has stored the estimates    
estimates table cat1 lin1

*Run wald test    
test lin1 cat1
*returns error "r(111) lin1 not found"

*Run wald test with alternative syntax
test [lin1 cat1]
*returns error "r(303) equation lin1 not found"

/*The same errors are reutrned if I put an = between lin1 and cat1. If I put         
cat1 before lin1,the error tells me instead taht cat1 is not found.*/

*I've also tried using the suest command, but I get the same errors
logistic ser i.age_ra i.sex ib5.indication ib2.gender_doc ///
ib1.specialty1 i.vol ib4.years_cat
estimates store cat2

logistic ser i.age_ra i.sex ib5.indication ib2.gender_doc ib1.specialty1 ///     
tenure i.vol ib4.years_cat
estimates store lin2

suest cat2 lin2, cluster(provider_id)

test lin2 cat2
*Error: r(111) lin2 not found

0 ответов

Другие вопросы по тегам