Displaying 1 result from an estimated 1 matches for "presnew".
Did you mean:
resnew
2012 Nov 28
2
Error message R2Jags
...I3[i]
log(mu[i]) <- max(-20, min(20, eta[i]))
PRes[i] <- (Y[i] - mu[i]) / sqrt(mu[i]) #estimate residuals (=
observed - expected / sqrt(variance))
#Discrepancy measures (used for checking overdispersion)
YNew[i] ~ dpois(mu[i]) #New data
PResNew[i] <- (YNew[i] - mu[i]) / sqrt(mu[i])
D[i] <- pow(PRes[i], 2)
DNew[i] <- pow(PResNew[i], 2)
}
# Add up discrepancy measures
fit <- sum(D[1:N])
fit.new <- sum(DNew[1:N])
for (i in 1:6) { b[i] ~ dnorm(0.0, 0.01)
# b[i] ~ dun...