model {
	A1 ~ dcat(p_A1[input_A1, ])
		input_A1 <- 2*(B2 - 1) + B1
	B1 ~ dcat(p_B1[input_B1, ])
		input_B1 <- 2*(C2 - 1) + C1
	B2 ~ dcat(p_B2[input_B2, ])
		input_B2 <- 2*(C4 - 1) + C3
	C1 ~ dcat(p_C1[input_C1, ])
		input_C1 <- 2*(D2 - 1) + D1
	C2 ~ dcat(p_C2[input_C2, ])
		input_C2 <- 2*(D4 - 1) + D3
	C3 ~ dcat(p_C3[input_C3, ])
		input_C3 <- 2*(D6 - 1) + D5
	C4 ~ dcat(p_C4[input_C4, ])
		input_C4 <- 2*(D8 - 1) + D7
	D1 ~ dcat(p_D1[])
	D2 ~ dcat(p_D2[])
	D3 ~ dcat(p_D3[])
	D4 ~ dcat(p_D4[])
	D5 ~ dcat(p_D5[])
	D6 ~ dcat(p_D6[])
	D7 ~ dcat(p_D7[])
	D8 ~ dcat(p_D8[])
}
