造型beta-二項式迴歸

Question

假設這個簡單的例子：

treatment <- factor(rep(c(1, 2), c(43, 41)), levels = c(1, 2),labels = c("placebo", "treated")) 
improved <- factor(rep(c(1, 2, 3, 1, 2, 3), c(29, 7, 7, 13, 7, 21)),levels = c(1, 2, 3),labels = >c("none", "some", "marked")) 
numberofdrugs<-rpois(84, 50)+1 
healthvalue<-rpois(84,5) 
y<-data.frame(healthvalue,numberofdrugs, treatment, improved) 
test<-lm(healthvalue~numberofdrugs+treatment+improved, y) 

我是什麼supossed辦時，我想估計有Rβ-二項式迴歸？有人熟悉它嗎？任何想法是讚賞！

Answer 1

我不明白這個例子是如何與beta-binomial迴歸有關的（例如，你已經產生了計數數據，而不是數量（總數可能））。爲了模擬的β-二項式數據，參見rbetabinom在任的emdbook或rmutil包...

library(sos); findFn("beta-binomial")發現了許多有用的出發點，包括

• AOD（overdispersed數據的分析）， betabin功能
• betabinomial家庭VGAM
• hglm包
• emdbook包裝（對於dbetabinom）加mle2包裝