Gia đình nhị thức R với liên kết danh tính

Tôi muốn điều chỉnh mô hình tuyến tính bằng R family=binomial(link="identity"), tuy nhiên, họ nhị thức không có liên kết danh tính. Tôi nên làm gì?

r generalized-linear-model binomial link-function

— yêu
nguồn

Tôi nghĩ rằng có một vấn đề thống kê cơ bản ở đây.

— Glen_b -Reinstate Monica

có, câu hỏi sau đây sẽ hỏi về Điều chỉnh lỗi tiêu chuẩn cho quá mức.

— David

Nhưng đối với người đầu tiên, tôi cần sử dụng liên kết danh tính trong họ nhị thức, nhưng R không cho phép.

— David

Lưu ý rằng nếu bạn cần cả hai xử lý quá mức và sử dụng một liên kết nhận dạng, bạn nên xem xét việc chuyển trực tiếp đến một mô hình gần đúng với chức năng phương sai nhị thức. Một mô hình nhị thức chỉ đánh chặn có thể được trang bị bằng tay .

— Glen_b -Reinstate Monica

(Và nếu bạn muốn tham số phân tán cố định tại một : family=binomial(link=make.link("identity")).

— Scortchi - Tái lập Monica

Xem Wikipedia về mô hình xác suất tuyến tính và các bài đăng CV tại đây & tại đây để biết nền tảng thống kê. Mặc dù không "sai", bạn muốn có một lý do chính đáng để sử dụng liên kết nhận dạng để mô hình xác suất Bernoulli.

Theo familyhướng dẫn

gia đình nhị thức [chấp nhận] các liên kết logit, probit, cauchit, (tương ứng với logistic, bình thường và Cauchy CDF tương ứng) log và cloglog(bổ sung log-log)

Nhưng

Các đối số liên kết và phương sai có ngữ nghĩa khá lúng túng cho khả năng tương thích ngược. Cách được đề xuất là cung cấp chúng dưới dạng các chuỗi ký tự được trích dẫn, nhưng chúng cũng có thể được cung cấp không trích dẫn (dưới dạng tên hoặc biểu thức). Ngoài ra, chúng cũng có thể được cung cấp dưới dạng một vectơ dài một ký tự cho biết tên của một trong các tùy chọn hoặc dưới dạng danh sách (cho link, của lớp "link-glm"). Các hạn chế chỉ áp dụng cho các liên kết được đặt dưới dạng tên: khi được cung cấp dưới dạng chuỗi ký tự, tất cả các liên kết được biết make.linklà được chấp nhận.

Vì vậy, family=binomial(link="identity")hoạt động nhưng family=binomial(link=identity)không. (Nếu bạn thấy khác đi thì có thể làm với phiên bản R.) Để cho phép phân tán quá mức, sau đó sử dụng family=quasi(link="identity", variance = "mu(1-mu)").

— Scortchi - Tái lập Monica
nguồn

liên kết = "danh tính" so với liên kết = sửa lỗi danh tính là một trợ giúp rất lớn. Đây là một bài tập trong sách giáo khoa CDA của Agresti. Mã anh ấy cung cấp là gần đúng (liên kết ...) mà bạn thảo luận, tuy nhiên tính đơn giản của việc thêm "" là sửa lỗi thanh lịch. Theo hiểu biết của tôi, cuộc gọi link = "nhận dạng" đại diện cho nhị thức như một mô hình tuyến tính.

— Justin Peterson