model {
	A1 ~ dcat(p_A1[])
	A2 ~ dcat(p_A2[input_A2, ])
		input_A2 <- A1
	A3 ~ dcat(p_A3[input_A3, ])
		input_A3 <- 2*(A1 - 1) + A2
	B1 ~ dcat(p_B1[input_B1, ])
		input_B1 <- 4*(A1 - 1) + 2*(A2 - 1) + A3
	B2 ~ dcat(p_B2[input_B2, ])
		input_B2 <- 4*(A3 - 1) + 2*(A2 - 1) + A1
	B3 ~ dcat(p_B3[input_B3, ])
		input_B3 <- 4*(A3 - 1) + 2*(A1 - 1) + A2
	C1 ~ dcat(p_C1[input_C1, ])
		input_C1 <- 4*(A2 - 1) + 2*(B3 - 1) + B2
	C2 ~ dcat(p_C2[input_C2, ])
		input_C2 <- 4*(A1 - 1) + 2*(B3 - 1) + B1
	C3 ~ dcat(p_C3[input_C3, ])
		input_C3 <- 4*(B1 - 1) + 2*(A2 - 1) + B2
	D1 ~ dcat(p_D1[input_D1, ])
		input_D1 <- 4*(B2 - 1) + 2*(C3 - 1) + C2
	D2 ~ dcat(p_D2[input_D2, ])
		input_D2 <- 4*(A1 - 1) + 2*(C3 - 1) + C1
	D3 ~ dcat(p_D3[input_D3, ])
		input_D3 <- 4*(C1 - 1) + 2*(A2 - 1) + C2
	E1 ~ dcat(p_E1[input_E1, ])
		input_E1 <- 4*(A3 - 1) + 2*(D3 - 1) + D2
	E2 ~ dcat(p_E2[input_E2, ])
		input_E2 <- 4*(D3 - 1) + 2*(D1 - 1) + A1
	E3 ~ dcat(p_E3[input_E3, ])
		input_E3 <- 4*(D2 - 1) + 2*(A2 - 1) + D1
}
