model {
	A ~ dcat(p_A[input_A, ])
		input_A <- B
	B ~ dcat(p_B[])
}
