Сравнение не вложенных моделей с использованием R
Чтобы объяснить мою проблему, у меня есть эти смоделированные данные с использованием 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")
Основываясь на моих моделях мышления, я считаю, что эти 2 модели не являются вложенными моделями.
Далее я хотел проверить, работает ли нелинейная модель (fit2) имеет какие-либо существенные эффекты по сравнению с обычной логистической моделью (fit1).
Есть ли способ сравнить эти две модели? Я считаю, что я не могу использовать lrtest функция в lmtest пакет, потому что эти две модели не являются вложенными моделями.
Любое предложение будет высоко оценено
Спасибо.