Can I define/constrain number of runs for each level of a categorical factor in custom design?

modelcitizen89

Dear JMP community,

 

I am in the processing of exploring the JMP and trying to make it work out a custom design for a bioprocess involved with 3 continuous factors and two categorical factors and some factor constrains. I would like to group same number of runs into each level of one of two categorical factors (6 levels in total, it is a "vessel"-like factor to contain different runs essentially).  Is there a possible way I can define or constrain the categorical factor to make it work.  It seems to me a more logical way is to define this categorical factor as a block effect instead. I have tried that option but then I would not allow me to keep the factor constrains that I had set up. Any guidance and explanation would be truly appreciated! 

Thanks in advance!

SDF1

Hi @modelcitizen89 ,

 

  Thanks for sharing your script, that helps to see things a little more clearly.  I can confirm your original error message: JMP does not allow disallowed combinations when explicitly including a Blocking variable in the model. You can include disallowed combinations if you leave Container as a categorical variable, then you can have the constraints present, but you get a different number of runs for each Container, which you want to keep fixed.

 

  I don't know if this is an appropriate way to get around that, but you can do what you are looking for by removing Container as a variable (it's not in your model, except as a main effect (by default); you are not looking at any interactions between the container and other factors). Instead of having Container as a categorical variable or blocking factor, you can set the DOE to generate it such that you have random block sizes (which you would choose as the number of vessels, 6). By default, the DOE wants to do 28, but this doesn't allow for an equal number of runs within each block size, so you need to change this to 36. This will generate a DOE for you with the constraints and the appropriate "containers". Again, I'm not sure this is the appropriate use, but it might get you to a place that allows you to set up your experiments.

 

  Lastly, I'm not positive, but I think there might have be an error in the disallowed script. Below is the JSL output from the Custom DOE that I ran with your constraints and using the runs per random block as the containers. I think this is what you want.

DOE(
	Custom Design,
	{Add Response( Maximize, "Y", ., ., . ),
	Add Factor( Continuous, 0.005, 0.03, "A", 0 ),
	Add Factor( Continuous, 0, 0.03, "B/C", 0 ),
	Add Factor( Continuous, -1, 1, "D", 0 ), Add Factor( Continuous, -1, 1, "E", 0 ),
	Add Factor( Categorical, {"B", "C"}, "Ingredient", 0 ),
	Set Random Seed( 531473239 ), Number of Starts( 11 ), Add Term( {1, 0} ),
	Add Term( {1, 1} ), Add Term( {2, 1} ), Add Term( {3, 1} ), Add Term( {4, 1} ),
	Add Term( {5, 1} ), Add Term( {1, 1}, {2, 1} ), Add Term( {1, 1}, {3, 1} ),
	Add Term( {1, 1}, {4, 1} ), Add Term( {1, 1}, {5, 1} ),
	Add Term( {2, 1}, {3, 1} ), Add Term( {2, 1}, {4, 1} ),
	Add Term( {2, 1}, {5, 1} ), Add Term( {3, 1}, {4, 1} ),
	Add Term( {3, 1}, {5, 1} ), Add Term( {4, 1}, {5, 1} ), Add Term( {1, 2} ),
	Add Term( {2, 2} ), Add Term( {3, 2} ), Add Term( {4, 2} ),
	Set Runs Per Random Block( 6 ), Set Sample Size( 36 ),
	Disallowed Combinations(
		Ingredient == "B" & ("B/C"n >= 0.02 & "B/C"n <= 0.03) | "B/C"n == 0 & (A >=
		0 & A <= 0.015)
	), Simulate Responses( 0 ), Save X Matrix( 0 ), Make Design,
	Set Run Order( Keep the Same ), Make Table}
);

  Hope this helps,

DS