r - 二项式 glmer() 预测的分类精度

Question

我一直在用我的(非 r-savvy)大脑来让 R 产生二项式 glmer 模型的正确预测的百分比。我知道这不是统计上的 super 信息，但经常被报道；所以我也想举报。

数据:

因变量:Tipo，它有 2 个值:'s' 或 'p'。
一堆因子预测变量，而不是单个连续变量。
2 随机截取:测试对象和他/她回应的名词

用于模型的代码:

model <- glmer(Tipo ~ agency + tense + 
               co2pr + pr2pr + socialclass + 
               (1|muestra) + (1|nouns), 
               data=datafile, family="binomial",
               control=glmerControl(optimizer="bobyqa"), 
               contrasts=c("sum", "poly"))

我知道有一个函数 predict()它需要一个模型对象并根据该模型制定预测，但我似乎无法让它为我工作。如果您愿意分享代码，我将不胜感激。

提前致谢。

Answer 1

为了进行预测，您需要一个阈值(关于此主题有大量文献 [搜索“ROC 曲线”或“AUC”]...)天真地选择 0.5 截止值(如果您不这样做，这是一个合理的默认值) t 知道或想要假设任何关于假阳性与假阴性的相对成本，或等效的敏感性与特异性的值(value))，然后

p <- as.numeric(predict(model, type="response")>0.5)

应该给出预测的概率并将它们分别转换为 0 或 1。然后

mean(p==datafile$Tipo)

应该给你正确的比例。

table(p,datafile$Tipo)

应该给你一个预测与观察表。