Cảnh báo: #successes không phải số nguyên trong một nhị thức nhị thức! (gói khảo sát)

Question

Tôi đang sử dụng gói twang để tạo điểm số xu hướng, được sử dụng làm trọng số trong một nhị thức bằng cách sử dụng survey::svyglm. Mã này trông giống như sau:

pscore <- ps(ppci ~ var1+var2+.........., data=dt....) 

dt$w <- get.weights(pscore, stop.method="es.mean") 

design.ps <- svydesign(ids=~1, weights=~w, data=dt,) 

glm1 <- svyglm(m30 ~ ppci, design=design.ps,family=binomial) 

Điều này tạo ra những cảnh báo sau đây:

Warning message: 
    In eval(expr, envir, enclos) : non-integer #successes in a binomial glm! 

Có ai biết những gì tôi có thể làm sai?

Tôi không chắc liệu thông báo này có tốt hơn trên số liệu thống kê.SE, nhưng về số dư, tôi nghĩ tôi sẽ thử ở đây trước tiên.

Answer 1

Không có gì sai, glm chỉ là cầu kỳ khi nói đến việc chỉ định mô hình nhị thức (và Poisson). Nó cảnh báo nếu nó phát hiện ra không. các thử nghiệm hoặc thành công là không tách rời, nhưng nó vẫn đi trước và phù hợp với mô hình. Nếu bạn muốn chặn cảnh báo (và bạn chắc chắn đó không phải là vấn đề), hãy sử dụng family=quasibinomial để thay thế.

Answer 2

Không có gì sai, tính toán, nhưng thống kê bạn có thể không làm điều gì đó có ý nghĩa nhiều. Trong trường hợp này, tốt hơn là nên sử dụng phương pháp hồi quy mạnh mẽ, thường là ý tưởng tốt cho dữ liệu đáp ứng tỷ lệ nếu dữ liệu của bạn bao gồm đơn vị chính xác 1 hoặc chính xác 0.

Answer 3

Xin lỗi, nhưng nó mạnh mẽ hơn ý nghĩa rằng nếu cơ chế bên dưới là một mô hình nhị thức quá mức, thì nhị thức thừa sẽ tính đến nó khi ước lượng erorr chuẩn. Do đó, bạn sẽ nhận được bảo hiểm tốt hơn, mặc dù điểm ước tính là như nhau.