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

	# observation
	for (i in 1:N.sample) {
		Landing[i] ~ dpois(landing[i])
		landing[i] <- exp(log.landing[i])
		log.landing[i] ~ dnorm(m[i], tau[1])
		m[i] <- x[LV[i]] + q.year[Year[i]]
	}

	# time change of density
	# LV[i]: pointer to lataent variable
	for (i in 1:N.sample) {
		x[LV[i]] ~ dnorm(mean.x[i], tau[2])
		mean.x[i] <- x[LV[i] - 1] + g[LV[i]]
		g[LV[i]] ~ dnorm(g[LV[i] - 1], tau[3])
	}
	for (k in 1:N.LV.init) {
		g[LV.init[k]] ~ dnorm(0.0, Tau.noninformative)
		x[LV.init[k]] ~ dnorm(0.0, Tau.noninformative)
	}
	for (y in 1:N.year) {
		q.year[y] ~ dnorm(0.0, tau[4])
	}

	# priors and hyper-priors
	for (k in 1:N.tau) {
		tau[k] ~ dgamma(Hyper.gamma, Hyper.gamma)
	}
}