Displaying 1 result from an estimated 1 matches for "totaln1".
Did you mean:
totaln
2010 Nov 15
2
Zero truncated Poisson distribution & R2WinBUGS
...l model for true abundance
for (i in 1:R) { # Loops over R sites
N1[i] ~ dpois(lambda1[i])
lambda1[i] <- lam1
}
# Observation model for replicated counts
for (i in 1:n) { # Loops over all n observations
C1[i] ~ dbin(p1[i], N1[site.p[i]])
p1[i] <-p0
}
# Derived quantities
totalN1 <- sum(N1[]) # Estimate total population size across all sites
}
",fill=TRUE)
sink()
# Package data for WinBUGS
R = dim(y)[1] # number of sites: 200
n = dim(y)[1] * dim(y)[2]#number of observations (sites*surveys)
win.data <- list(R = R, n = n, C1 = C1, site.p = site.p)
# Inits
Nst <...