Как использовать функцию Бокса-Тидвелла с логистической регрессией в R - PullRequest
Новая
       21

Как использовать функцию Бокса-Тидвелла с логистической регрессией в R

1 голос
/ 29 мая 2019

Я получаю ошибки, используя функцию boxTidwell из пакета 'car' для модели логистической регрессии.

Я хочу смоделировать

fatalCancer ~ globy1, где fatalCancer - это фактор с двумя уровнямии globy1 является числовым (все положительные).Я проверяю это, чтобы проверить предположение о линейности globy1 с логитом результата.

Глядя на сообщения об ошибках (ниже) и код функции boxTidwell, кажется, что, возможно, существует проблема с то, что fatalCancer являетсяфактор.В документации boxTidwell ничего не указано о том, что это логистическая модель.В примере, приведенном в разделе 6.4 «Fox Companion к прикладной регрессии» Fox (стр.312), пример логистической регрессии не требует какой-либо спецификации.

Есть ли способ исправить синтаксис функции boxTidwell ниже?

> library(car)
Loading required package: carData
> 
> load("m2dat.RData")
> m2dat <- na.omit(m2dat)
> dim(m2dat) 
[1] 116   3    
> head(m2dat)
   dog globy1 fatalCancer
1 101A    3.1          No
2 102A    2.9          No
3 103A    4.9          No
4 104A    3.1         Yes     
5 105A    2.8         Yes
6 106A    3.5          No
> boxTidwell(fatalCancer ~ globy1, data=m2dat)
 MLE of lambda Score Statistic (z) Pr(>|z|)    
        6.5694                  NA       NA

iterations =  21      
There were 48 warnings (use warnings() to see them)
> warnings()    
Warning messages:
1: In model.response(mf, "numeric") :
  using type = "numeric" with a factor response will be ignored
2: In model.response(mf, "numeric") :
  using type = "numeric" with a factor response will be ignored
3: In Ops.factor(y, z$residuals) : ‘-’ not meaningful for factors
4: In model.response(mf, "numeric") :
  using type = "numeric" with a factor response will be ignored
5: In Ops.factor(y, z$residuals) : ‘-’ not meaningful for factors
6: In Ops.factor(r, 2) : ‘^’ not meaningful for factors
7: In Ops.factor(y, z$residuals) : ‘-’ not meaningful for factors
8: In Ops.factor(y, z$residuals) : ‘-’ not meaningful for factors    
9: In Ops.factor(y, z$residuals) : ‘-’ not meaningful for factors
10: In Ops.factor(y, z$residuals) : ‘-’ not meaningful for factors
...
46: In Ops.factor(y, z$residuals) : ‘-’ not meaningful for factors 
47: In Ops.factor(y, z$residuals) : ‘-’ not meaningful for factors
48: In Ops.factor(y, z$residuals) : ‘-’ not meaningful for factors

Статистика результатов заканчивается NA, и я хотел бы успешно выполнить тест.

...