turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- JMP User Community
- :
- Discussions
- :
- Power and DoE with mixtures (JMP 12)

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Apr 30, 2015 1:17 AM
(1035 views)

Hi,

I'm trying to make a DoE with three mixture factors (x1, x2, x3) which sum up to some constant, for example x1+x2+x3=50.

I've noticed that, when I make an optimal mixture design (with proportions - i.e. I suppose that x1+x2+x3=1) that the power in design evaluation window, for all coefficients is equal to 0.073.

On the other hand, If I change "Mixture Sum" parameter to 50 (i.e. x1+x2+x3=50) and then make the same design, the power for all coefficients is 1.

This is puzzling to me, since model matrix isn't changed, that is the range of all factors is translated to [0,1] interval.

Can someone explain the reason for that?

From the technical details of power calculation (DESIGN OF EXPERIMENTS GUIDE • EVALUATE DESIGNS • TECHNICAL DETAILS) I don't see the reason for a change.

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Solution

The model matrix X has changed because the design matrix D has changed. Here is a simple example for illustration. I am using Custom Design with three mixture components X1-X3. For the simple model with only the main effects and intercept, I get this design in 6 runs with its associated power:

[Missing image]

If I click Back and change the Mixture Sum to 50 and then click Make Design I get a different design matrix, which in turn is expanded into a different model matrix:

[Missing image]

The power depends on the model matrix, which is determined by the design matrix and the terms in the model. Comparing these two designs, the one that spans a smaller factor space will have less power to estimate the same sized effect (coefficient).

The same thing happens if you leave the mixture sum set to the default value of 1, but restrict the range of any of the components.

Learn it once, use it forever!

4 REPLIES

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Apr 30, 2015 5:58 AM
(778 views)

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

May 4, 2015 1:21 AM
(778 views)

Hi LouV, thank you for your answer. I'm aware of that power of test for coefficients == 0 isn't helpful here.

Information about prediction variance and stuff like G-efficiency is in my focus also, but still I'm trying to understand how things are calculated because that is the only way I can remember what and why is important to look at

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

The model matrix X has changed because the design matrix D has changed. Here is a simple example for illustration. I am using Custom Design with three mixture components X1-X3. For the simple model with only the main effects and intercept, I get this design in 6 runs with its associated power:

[Missing image]

If I click Back and change the Mixture Sum to 50 and then click Make Design I get a different design matrix, which in turn is expanded into a different model matrix:

[Missing image]

The power depends on the model matrix, which is determined by the design matrix and the terms in the model. Comparing these two designs, the one that spans a smaller factor space will have less power to estimate the same sized effect (coefficient).

The same thing happens if you leave the mixture sum set to the default value of 1, but restrict the range of any of the components.

Learn it once, use it forever!

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

May 4, 2015 1:51 AM
(778 views)

Hi Mark, thx for your explanation.

I figured out that the reason must be in model matrix X, but I looked at coding table (i.e. wrong table) of both designs and didn't see the difference. Now when I look at the correct table, everything is clear now (although, at this moment I don't understand why model matrix doesn't have coded values (as documentation states) instead of 0-50 ?)

This widget could not be displayed.