Đầu ra của mô hình logistic trong R

Tôi đang cố gắng diễn giải loại mô hình logistic sau:

mdl <- glm(c(suc,fail) ~ fac1 + fac2, data=df, family=binomial)

Là đầu ra của predict(mdl)tỷ lệ thành công dự kiến ​​cho mỗi điểm dữ liệu? Có một cách đơn giản để lập bảng tỷ lệ cược cho từng cấp độ yếu tố của mô hình, thay vì tất cả các điểm dữ liệu?

Các trang trợ giúp cho

predict.glm

trạng thái: "Do đó, đối với mô hình nhị thức mặc định, các dự đoán mặc định là tỷ lệ cược log (xác suất theo thang đo logit) và 'type =" reply "' đưa ra xác suất dự đoán". Vì vậy, predict(mdl)trả về nhật ký (tỷ lệ cược) và sử dụng "type =" reply sẽ trả về xác suất dự đoán. Bạn có thể tìm thấy ví dụ đồ chơi này mang tính hướng dẫn:

> y <- c(0,0,0,1,1,1,1,1,1,1)
> prop.table(table(y))
y
  0   1 
0.3 0.7 
> glm.y <- glm(y~1, family = "binomial")
> ## predicted log(odds)
> predict(glm.y)
        1         2         3         4         5         6         7         8 
0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 0.8472979 
        9        10 
0.8472979 0.8472979 
> ## predicted probabilities (p = odds/(1+odds))
> exp(predict(glm.y))/(1+exp(predict(glm.y)))
  1   2   3   4   5   6   7   8   9  10 
0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 
> predict(glm.y, type = "response")
  1   2   3   4   5   6   7   8   9  10 
0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 

Về câu hỏi thứ hai của bạn, bạn có thể muốn xem gói hiệu ứng http://socserv.socsci.mcmaster.ca/jfox/Misc/effects/index.html của John Fox; xem thêm bài viết JSS của ông "Hiệu ứng hiển thị trong R cho các mô hình tuyến tính tổng quát" (trang 8-10).