0

我的数据有三列Time, Interval, Count。我有一个像这样的泊松混合物

mod_string = " model{

for(i in 2:length(Count)){

Count[i] ~ dpois(lambda.hacked[i]*z[i]+0.0001)
z[i] ~dbern(p)
lambda.hacked[i] <- mu[ clust[i] ]
Prob <- p^-(1:i) * (1-p) / p
mu <- (Time[1:i] - Interval[1:i])*lambda
clust[i] ~ dcat( Prob)

}

## Priors

lambda ~ dgamma(0.01,0.02)
p ~ dbeta(1,1)
}"

mu每次迭代都会改变大小。随着i增长,集群的数量也在增长。

我该如何适应这个?

4

0 回答 0