model
{
	Tau.noninformative <- 1.0E-3
	Hyper.gamma <- 1.0E-2
	Tau.dummy <- 2.0E+4

	# observation
	for (i in 1:N.sample) {
		O[i] ~ dpois(nobs[i])
		nobs[i] <- exp(n[i])
	}

	# time change of density
	g[1] ~ dnorm(0.0, Tau.noninformative) # initial g
	n1 ~ dnorm(0.0, Tau.noninformative) # initial n
	n[1] <- n1
	for (j in 2:N.year) {
		g[j] ~ dnorm(g[j - 1], tau)
		mean.n[j] <- n[j - 1] + g[j]
		n[j] ~ dnorm(mean.n[j], Tau.dummy)
	}

	# priors and hyper-priors
	tau ~ dgamma(Hyper.gamma, Hyper.gamma)
}