cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 
Choose Language Hide Translation Bar
Marco_
Level III

Using the Prediction Profiler to optimize for a certain range

Hi everyone,

I'm trying to use the prediction profiler in JMP in combination with the desirability function in order to optimize several models simulatanously. 
There are several optimisation options in JMP implemented like maximizing, minimizing or matching a target. However, I need an option to match a range.

Since there are only 3 adjustable points in the desitability function and I need a desirability function where all values outside the target ranges are equal to 0 and within equal to 1, I need your help how I can achieve this ( the function should look similar to a rectangle)

Please note, I'm familiar with the option to set specification limits in order to perform simulations or to calculate defect rates. These static plots are unfortunately not a big help for me. What I need is the interactivity of the desirability plot.

 

I'm looking forward to your suggestions!

 

Thanks!

11 REPLIES 11

Re: Using the Prediction Profiler to optimize for a certain range

@Marco_ if you could go to the JMP Wish List https://community.jmp.com/t5/JMP-Wish-List/idb-p/jmp-wish-list and provide your feature request, that would be helpful for use to keep track of this and to be able to better evaluate it.  You can also reach out to me directly via PM on the Community and we could discuss this further.  I am interested to find out more specifics about what you want to do.  

Marco_
Level III

Re: Using the Prediction Profiler to optimize for a certain range

@SamGardner 
Instead of writing you a PM I thought it also might be more helpful for other JMP users if we keep the discussion public. I also google very frequently questions I have and I'm always happy if I find a post in the community

 

I spent now some time with prediction profiler and you guys have done a really great job there. This feature is really great to tackle QbD or PAR  requests from authorities such as the FDA or EMA. However, at the moment is see two big problems why I cannot use the feature.

 

1. As I already mentioned in my last post, I'm afraid that the model uncertainty in the simulation is underestimated. @John_Sall explained at the end of this video how to incorporate the model uncertainty within PredFormula. I could reproduce his example. But, the used RMSE of the model does not consider the correlation structure within the model. You anticipate that the model uncertainty across the whole investigated design space (lets call the design matrix X) is constant. This however is in general not true, since most of the models are not completely orthogonal regarding their input-parameters. See below a short equation how the uncertainty of an estimated coefficient is dependent on the design correlation structure.

Marco__0-1657873615966.png

Since JMP allows easily to save the PredSE formula of each model, where the correlation structure is considered. I wonder why the formula is not used for the simulation within the DesignSpace profiler? Instead of using a constant SD, simply connect the PredSE formula here and the issue should be solved.

 

2. The other point is that by using the RMSE of the model as the SD for the simulation within the DesignSpace Profiler, the uncertainty in the model mean prediction is neglected. From my understanding the simulation anticipates that the mean prediction of the model has no error and only considers the uncertainty around this model prediction. To make my point a bit clearer, see the following formula for the calculation of a prediction interval (just copy/pasted from a paper):

Marco__1-1657875685519.png

At the moment you are only considering the first part of the equation but not the second part (i.e. the standard error).

 

 

To sum up my points. I think at the present status, the DesignSpace profiler underestimates the true variability and thus allow to wide ranges for the input-parameters. I think both points could easily be fixed and I would much appreciate your feedback (also from @Phil_Kay if he has some time during his train trip )