Я создал матрицу дизайна X, используя
X <- model.matrix( ~ age + BMI + sex, data=df)
, и результат выглядит так, как я ожидаю.Я хочу проверить мультиколлинеарность, поэтому я использую
library(mctest)
mctest(X, y)
, что дает
Error in mctest(X, y) : X must be a numeric matrix
Я что-то здесь упускаю?
Редактировать: Вот выводул ().Я немного сократил его, моя фактическая модель намного больше, чем та, что показана выше.
num [1:210, 1:41] 1 1 1 1 1 1 1 1 1 1 ...
- attr(*, "assign")= int [1:41] 0 1 2 3 4 4 4 5 5 5 ...
- attr(*, "contrasts")=List of 27
..$ Sex : chr "contr.treatment"
[...]
attr(*, "dimnames")=List of 2
..$ : chr [1:210] "1" "2" "3" "4" ...
..$ : chr [1:41] "(Intercept)" "age" "BMI" "Sex1" ...
- attr(*, "scaled:center")= Named num [1:41] 1 65.0048 25.9557 0.2857 0.0381 ...
..- attr(*, "names")= chr [1:41] "(Intercept)" "age" "BMI" "Sex1" ...
- attr(*, "scaled:scale")= Named num [1:41] 0 11.317 5.129 0.453 0.192 ...
..- attr(*, "names")= chr [1:41] "(Intercept)" "age" "BMI" "Sex1" ...