global shapeParameterP = 1;
global shapeParameterQ = 1;
shapeParameterP:>0.05;shapeParameterP:<85;
shapeParameterQ:>0.05;shapeParameterQ:<85;
category categoryPlaceholder = (rateClassCount, EQUAL, MEAN, 
							 _x_^(shapeParameterP-1)*(1-_x_)^(shapeParameterQ-1)/Beta(shapeParameterP,shapeParameterQ), /* density */
							 IBeta(_x_,shapeParameterP,shapeParameterQ), /*CDF*/
							 0, 				   /*left bound*/
							 1, 			   /*right bound*/
							 IBeta(_x_,shapeParameterP+1,shapeParameterQ)*shapeParameterP/(shapeParameterP+shapeParameterQ)
		  	 				);
