Как создать таблицу классификации логита в R при удалении значений ...?

Question

Я французский шпиль, и я пытаюсь сделать классификационную таблицу на модели Logit.

Проблема, R возвращает мне ошибку, когда я пытаюсь ее построить, потому что длина glm не совпадает с длиной моей зависимой переменной.

Что я пытался сделать ...:

1.

>reg4=glm(buyEV3~lvlEtu+typecar1+typecar2+typecar3+bonus+goodForEnv+cheaper+knowEV,data=bdd,family=binomial(link="logit")) 
> reg4 
[...] 
Degrees of Freedom: 230 Total (i.e. Null); 222 Residual 
    (85 observations deleted due to missingness) 
[...] 
>mysteps=step(reg4,buyEV3~lvlEtu+typecar1+typecar2+typecar3+bonus+goodForEnv+cheaper+knowEV,data=bdd,family=binomial(link="logit")) 
[...] 
> classDF=data.frame(response=bdd$buyEV,predicted=round(fitted(mysteps),0)) 
    Error in data.frame(response = bdd$buyEV, predicted = round(fitted(mysteps), : 
    arguments imply differing number of rows: 316, 231 

2.

> library(QuantPsyc) 
> ClassLog(reg4,bdd$buyEV,cut=.5) 
Error in table(predict(MOD, type = "response") > cut, resp) : 
    all arguments must have the same length 

Как исправить эту проблему?

спасибо :)

Answer 1

Отсносящихся Пытаться 1. Поскольку только полные наблюдения включены в логистической модели «(85 наблюдений удаленных из-за missingness)» длину вашего зависимого varable в исходном наборе данных и тому количество значений, заданных при запуске fitted(), будет отличаться на ту же сумму.

Экстрагирование зависимой переменной из наблюдений фактически включены в модель будет делать трюк

classDF=data.frame(response=reg4$model$buyEV,predicted=round(fitted(mysteps),0)) 

glmobject$model$variable Где это то, что важно в подобных ситуациях.