Тест Вальда для сравнения двух моделей в 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