В настоящий момент у меня проблемы с созданием кубического сплайна.
Кто-нибудь знает, как создать кубический сплайн для скорректированной модели?Моя нескорректированная модель выглядит следующим образом:
M1 <- lm(menSpmet127ttscore1~ bs(G_alcre, knots = knots), data = MR).
У меня нет проблем с запуском модели.
Мой код для скорректированной модели:
M1 <- lm(menSpmet127ttscore1~ bs(G_alcre+Pa_Index, knots = knots), data = MR)
Это делаетне работает, и я получаю следующую ошибку:
Ошибка в splineDesign (Aknots, x, ord): длина «производных» больше длины «x»
Было бы здорово, если бы кто-то мог мне помочь.Заранее спасибо.