model {
	A ~ dcat(p_A[])
	B ~ dcat(p_B[input_B, ])
		input_B <- 3*(C - 1) + A
	C ~ dcat(p_C[])
}
