造型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包装