Чтобы объяснить мою проблему, у меня есть эти смоделированные данные с использованием R.
require(splines)
x=rnorm(20 ,0,1)
y=rep(c(0,1),times=10)
Сначала я установил регулярную (линейные эффекты) модель логистической регрессии.
fit1=glm(y~x ,family = "binomial")
Затем, чтобы проверитьнелинейные эффекты, я приспособил эту естественную сплайн-модель.
fit2=glm(y~ns(x,df=2) ,family = "binomial")
Основываясь на моих моделях мышления, я считаю, что эти две модели не являются вложенными моделями.
Далее я хотел проверить, неЛинейная модель (fit2) имеет какие-либо существенные эффекты по сравнению с обычной логистической моделью (fit1).
Можно ли сравнить эти две модели?Я считаю, что не могу использовать функцию lrtest в пакете lmtest, потому что эти две модели не являются вложенными моделями.
Любое предложение будет высоко оценено
Спасибо.