以下代码与我的数据类似:
data cast_comm_concat;
set cast_comm;
by id; /*make sure it's sorted*/
length comments $ 1000 ;
retain comments;
if (first.id) then do;
comments="";
end;
comments=catx(",",comments,comment);
if last.id then do;
output;
end;
run;
我需要使用quasipoisson分布运行混合模型:
uno<-rnorm(100, 0, 1)
due<-rpois(100, 2)
tre<-rnorm(100,0,1)
randomm<-as.factor(seq(from=1, to=length(uno), by=1))
data <- cbind.data.frame(uno=uno, due=due, randomm=randomm, tre=tre)
我收到以下错误:
mod <- glmmadmb(due~uno+tre+(1|randomm), data = data, family="nbinom1")
这究竟是什么意思?据我所知,它必须对我的数据类做一些事情。如何解决这个问题?