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; } } }

Fitting scheffe model in mixture + non mixture design

akrishna39

Hi,

I have a custom design with 4 mixture variables and 13 non mixture variables. I am performing a screening experiment with 24 runs. When I use the "Fit Model" functionality with only main effects and choose "No intercept" to fit a model,

1) Do i get a scheffe model?

2) I am not including any interaction terms in the model as I am performing screening and so am interested in only the main effects. So, is there any benefit of developing a scheffe model in this case?

Accepted answers

Dan_Obermiller

Anytime you are building a model with mixture factors that add up to 1, you should be fitting a mixture model, such as the Scheffe mixture model. If you don't, you will get a singular matrix message. This is true even with just a main effects model. You really don't have a choice.

Remember that a Scheffe mixture model is not the same as a no-intercept model. When fitting a Scheffe mixture model in JMP, the Analysis of Variance report will state: Tested against reduced model Y=mean. If the model is a no-intercept model, the message would be: Tested against reduced model Y=0.

All comments

Dan_Obermiller

akrishna39

Thanks Dan!

But when I click on "Fit model" under "Analyze" I do not get an Analysis of Variance report. I get the Least square fit parameter estimates. I observe that the mixture variabels are scaled (mean = 0 and range = 1) and their parameter estimates are also present. So, are these estimates the co-efficients of the Scheffe's model?

Dan_Obermiller

With version 13 of JMP the Analysis of Variance report is closed by default. You need to open it by going to the red pop-up menu at the top of your report next to "Response Y" with Y being the name of your response variable. From there choose Regression Reports > Analysis of Variance. Remember that you sometimes will need to request the items that you need for your analysis from those red popup menus.

As for your parameter estimates, if JMP created the design, then it will automatically turn on pseudocomponent coding for your mixture factors as well as the "regular" coding for the non-mixture factors. This assumes that when you created the design you specified the natural units for your factors. This coding does not indicate if a mixture model was fit to the data or not.

akrishna39

Thanks Dan!

I understand that for scheffe model the calculation of p-values for mixture variables is not straightforward as for a regular model because the t-test do not assess the statistical significance of the component effects. Is there any documentation which shows how the p-values are being calculated for mixture variables in case of a scheffe model?