	model
	{
		smoking ~ dcat(p.smoking[1:2])
		tuberculosis ~ dcat(p.tuberculosis[asia,1:2])
		lung.cancer ~ dcat(p.lung.cancer[smoking,1:2])
		bronchitis ~ dcat(p.bronchitis[smoking,1:2])
		either <- max(tuberculosis,lung.cancer)
		xray ~ dcat(p.xray[either,1:2])
		dyspnoea ~ dcat(p.dyspnoea[either,bronchitis,1:2])
	}
