nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Custom design - Help to split-plot/blocking

Bertelsen92

Hi JMP community,

I am currently working on a custom DOE design. Overall I have 3 factors with varying levels.

Stabilizer solution (I have 7 different stabilizer solution I would like to screen) + potentially adding a reference sample (no stabilizer added).

Dosage of stabilizer - 2 levels (high and low) => (0,01% and 0,02%)

Process (2 different heat treatments/autoclavations).

The pilot plant is able to mix 4 unique codes pr. day. (we have 4 days in the pilot plant booked). This means we can maximum mix 16 unique codes that can be run through both heat treatments (giving a total of 32 codes).

However the fact that we can only mix 4 codes pr. day and these need to be processed with both heat treatments (giving 8 codes pr. day) gives me some troubles. The ideal order will off course be so that:

Day 1:

Stabilizer 1: 0,01% - heat treatment 1

Stabilizer 4: 0,02% - heat treatment 1

Stabilizer 5: 0,01% - heat treatment 1

Stabilizer 6: 0,02% - heat treatment 1

Stabilizer 1: 0,01% - heat treatment 2

Stabilizer 4: 0,02% - heat treatment 2

Stabilizer 5: 0,01% - heat treatment 2

Stabilizer 6: 0,02% - heat treatment 2

etc...

However I can not get JMP to do this (have tried using blocking, and editing the "changes" from easy to hard etc.).

This best answer I got myself is attached.

I am not familiar with scripting (JSL), so those solutions will be hard for me. I am using JMP 17.

Thank you in advance.

Best regards

Kristian Bertelsen

Design 1 - run order.PNG

Custom design - 1.PNG

Accepted answers

All comments

Victor_G

Hi @Bertelsen92,

Which effects would you like to see in the models ?

For the screening of main effects, I'm able to create a 32-runs in 4 whole plots (so 8 runs per day like you mentioned), where you only have one process per day.
Here is the script :

DOE(
	Custom Design,
	{Add Response( Maximize, "Y", ., ., . ),
	Add Factor( Continuous, -1, 1, "X1", 0 ),
	Add Factor( Continuous, -1, 1, "X2", 0 ),
	Add Factor( Continuous, -1, 1, "X3", 0 ),
	Add Factor( Continuous, -1, 1, "X4", 0 ),
	Add Factor( Continuous, -1, 1, "X5", 0 ),
	Add Factor( Continuous, -1, 1, "X6", 0 ),
	Add Factor( Continuous, -1, 1, "X7", 0 ),
	Add Factor( Continuous, -1, 1, "X8", 0 ),
	Add Factor( Categorical, {"L1", "L2"}, "Process", 1 ),
	Set Random Seed( 235758137 ), Number of Starts( 1574 ), 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( {6, 1} ), Add Term( {7, 1} ), Add Term( {8, 1} ),
	Add Term( {9, 1} ), Add Alias Term( {1, 1}, {2, 1} ),
	Add Alias Term( {1, 1}, {3, 1} ), Add Alias Term( {1, 1}, {4, 1} ),
	Add Alias Term( {1, 1}, {5, 1} ), Add Alias Term( {1, 1}, {6, 1} ),
	Add Alias Term( {1, 1}, {7, 1} ), Add Alias Term( {1, 1}, {8, 1} ),
	Add Alias Term( {1, 1}, {9, 1} ), Add Alias Term( {2, 1}, {3, 1} ),
	Add Alias Term( {2, 1}, {4, 1} ), Add Alias Term( {2, 1}, {5, 1} ),
	Add Alias Term( {2, 1}, {6, 1} ), Add Alias Term( {2, 1}, {7, 1} ),
	Add Alias Term( {2, 1}, {8, 1} ), Add Alias Term( {2, 1}, {9, 1} ),
	Add Alias Term( {3, 1}, {4, 1} ), Add Alias Term( {3, 1}, {5, 1} ),
	Add Alias Term( {3, 1}, {6, 1} ), Add Alias Term( {3, 1}, {7, 1} ),
	Add Alias Term( {3, 1}, {8, 1} ), Add Alias Term( {3, 1}, {9, 1} ),
	Add Alias Term( {4, 1}, {5, 1} ), Add Alias Term( {4, 1}, {6, 1} ),
	Add Alias Term( {4, 1}, {7, 1} ), Add Alias Term( {4, 1}, {8, 1} ),
	Add Alias Term( {4, 1}, {9, 1} ), Add Alias Term( {5, 1}, {6, 1} ),
	Add Alias Term( {5, 1}, {7, 1} ), Add Alias Term( {5, 1}, {8, 1} ),
	Add Alias Term( {5, 1}, {9, 1} ), Add Alias Term( {6, 1}, {7, 1} ),
	Add Alias Term( {6, 1}, {8, 1} ), Add Alias Term( {6, 1}, {9, 1} ),
	Add Alias Term( {7, 1}, {8, 1} ), Add Alias Term( {7, 1}, {9, 1} ),
	Add Alias Term( {8, 1}, {9, 1} ), Set N Whole Plots( 4 ), Set Sample Size( 32 ),
	Simulate Responses( 0 ), Save X Matrix( 0 ), Make Design}
);

And attached you can find the datatable.

The trick is to specify your 2-level categorical factor "Process" as a hard-to-change factor, and specify 4 as the number of whole plot (4 days) and 32 as the number of runs (8 runs per day). Note that for screening main effects, you could reduce the number of runs if needed (but it should be a multiple of the number of experimentation day, so if you keep 4 days, you can reduce to 28 or 24 runs for example) :

If I misunderstand your objective and/or design specifications, don't hesitate to add more details and explanations, so we can find a better solution.
I hope this will help you,

Edit: This design doesn't seem to respect your constraint of having only 4 different mixtures/day, so increasing the number of whole plots to 8 (and doing 2 whole plots of 4 experiments each per day) with the same number of runs could be an alternative, I have to check later if it works.

Sorry for the misunderstanding,

Custom Design_Split-plot.jmp

Bertelsen92

Hi Victor,

Thanks a lot for helping out.

Keep in mind I want to test 7 different stabilizers + a reference. So shouldn't these be put in as a categorical 8 level factor?

Because we can only add one stabilizer solution to the blend at a time?

Yes the 4 mixtures pr day is the limiting factor. However the process (heat treatment - 2 different ones) can be run at the same day.

That's why it would be perfect if you could run 4 mixtures pr. day that will then be heat-treated by both treatments.

I have succeeded in getting this order (please see the attached pictures). However I am afraid that the power score is way to low for the stabilizers which is the main effect I want to investigate. What do you think?

Custom design - bad power for stabilizer (main effect).PNG

Custom design - bad power (correct order).PNG

Custom design - bad power.PNG

Victor_G

Hi Kristian,

Sorry for my misunderstanding, you're completely right about the settings of your factors.

Now that everything seems clear (thanks a lot for your snapshots that helped me define the factors and model needed), I was able to find a design able to fit your constraints : 4 mixtures per day, with the 2 processes involved.

In this case, I didn't set a Factor as "Hard to change", but created a blocking factor "Day", with 8 runs per day, and "naturally" 4 mixtures are done each day.

I hope this would be a better solution that fits your needs ?

Please find the datatable attached, you can then compare this design with the one you already have prepared.

Custom Design_Stabilizers.jmp

Bertelsen92

Hi Victor,

No worries. Great, however I am open to other suggestions if it could be done smarter.

However I am afraid your suggestions does not solve my problem. Because L6 = 0,01 and L6 = 0,02 are unique mixes and just as different as L7=0,01 and L3=0,02 would be. So within Day 1 there would need to be mixed 8 codes (4 different stabilizers but a 2 levels), which is not possible.

Maybe it cannot be done in JMP or in general. I might need to rethink the design completely.