Số nhị phân âm trong GEE

Đối với các gói R triển khai GEE như gee, geepack, có vẻ như là họ nhị thức âm không được bao gồm. Tôi có hai câu hỏi:Số nhị phân âm trong GEE

Có gói R nào khác cho GEE mà tôi không biết?
Nếu không, là có một bước đơn giản để cho phép việc tạo ra một gia đình, tức là cung cấp các chức năng liên kết (log mu) và chức năng phương sai (mu + mu^2/theta), giả sử theta được quy định (nếu không thì NB không phải là một GLM) và sau đó để các mã gee hoặc geepack hoạt động kinh doanh theo cách tương tự với glm?

Nguồn

2012-12-19 alittleboy

Tìm kiếm nhanh, sản xuất các gói sau: GWAF, MuMIn, multgee, repolr, CorrBin, sếp, cưa, ape –

https://stat.ethz.ch/pipermail/r-help/2011-June/282116. html –

'thư viện (" sos "); findFn ("{generalized estimating equation}") 'rất hữu ích - tìm thấy một danh sách tương tự với @ RicardoSaporta ở trên - nhưng lưu ý rằng không phải tất cả những điều này thực sự phù hợp với GEE, một số chỉ tận dụng' geepack' hoặc những người khác như một phần của các thủ tục khác –

Bạn sẽ có thể sử dụng gia đình negative.binomial định nghĩa trong gói MASS để làm điều này (thiết lập một gia đình NB với một giá trị theta chỉ định). Có vẻ như geepack::geese (ít nhất) sẽ chấp nhận thông số kỹ thuật của gia đình trong biểu mẫu này. Để ước tính theta bạn có thể thử nhúng GEE phù hợp với một theta cố định vào một vòng lặp, hoặc thực hiện một chức năng geefit_NB(theta) và tối ưu hóa trên theta.

Nếu negative.binomial đã chưa tồn tại trong MASS, bạn có thể xác định gia đình riêng của mình (đây là thừa nhận một chút tiên tiến - Tôi sẽ bắt đầu bằng cách tải về mã nguồn của gói MASS và nhìn vào các tập tin R/neg.bin.R).

Nguồn

2012-12-21 19:40:29

Cảm ơn! Tôi định thêm tùy chọn 'negbin' vào hàm' geese.fit' thông qua 'VARIANCE', chỉ chứa gaussian, binomial, poisson và Gamma bây giờ. Tôi đoán phần quan trọng là để tìm một cách hiệu quả để cho 'geese.fit' lấy kết quả từ một sự phù hợp nhị thức âm từ gói MASS. Điều này, như bạn đã nói, có thể là một chút tiên tiến và tôi đã không tìm thấy bất kỳ nguồn cho một ví dụ đồ chơi ... – alittleboy

Bạn đang phải đó là khó hơn tôi nghĩ - Tôi nghĩ rằng sự tồn tại của "gia đình" có nghĩa là nó có thể tự động sử dụng phương thức do người dùng định nghĩa (hoặc ít nhất là không được định nghĩa), nhưng tôi thấy rằng các kiểu phương sai được mã hóa cứng trong các cuộc gọi đến mã C ... –

Nhìn lại, tôi không nghĩ câu trả lời này thực sự trả lời câu hỏi, bởi vì (như đã nêu trong bình luận ở trên), hàm sai số được mã hóa cứng hơn tôi nghĩ. Tồi tệ hơn, tôi không thấy một vị trí rõ ràng trong 'gee' hoặc' geepack' để chuyển một tham số vào hàm phương sai (tham số có thể điều chỉnh duy nhất liên quan đến phương sai là tham số tỷ lệ, là hằng số tỷ lệ tổng thể. ..) –

Số nhị phân âm trong GEE

Trả lời

Các vấn đề liên quan