配方实验设计我必须为 2 个连续因子和 1 个定性因子定义 DOE,并在不同浓度下测试 2 个水平(A 和 B)。我的限制如下:
- 如果 A 浓度 >0,则 B=0
- 如果 B 浓度 >0,则 A=0
目的是优化配方。
我尝试进行自定义设计,但找不到定义约束的方法。
最佳值应选择定性因子(A 或 B)的最佳水平以及 -1 到 1 范围内的浓度水平。
我将感谢您对这个主题的帮助。
提前致谢。
回复:DOE 配方如果范围非常不同,可以使用的另一种方法是使用从 -1 到 1 的编码范围。当您遇到 (A,-1) 处的处理时,这意味着您在 10 处运行 A。如果是(B,-1),您以 800 运行 B。您需要记住将结果转换回适当的比例。
回复:DOE 配方感谢您提供此附加信息。所以,如果我很好理解的话,当射程确实不同时,并不真正建议使用这个母鹿尖端。回复:DOE 配方定制设计对我来说没有收敛,所以我无法检查最终的设计。我相信你对 A 和 B 水平不平衡的观察。
我删除了约束并进行了设计,启用了模拟响应,制作了数据表,手动将 A 的 505 更改为 12.5,将 1000 更改为 15,将 B 的 10 更改为 800,将 505 更改为 900,然后模拟响应。
我怀疑由于存在奇点,这些数据实际上并不支持估计模型参数。事实证明,情况并没有那么糟糕,只是非常糟糕。 A 和 B 的所需浓度水平在 X3、浓度以及与 X3 的任何相互作用的估计值之间产生非常强的相关性。这种相关性反过来又导致估计不稳定。我在模拟响应中加入了 Conc 的强烈非线性效应和较大的 X3*Conc 相互作用,但根本没有检测到它们。
我附上这个数据表供您检查。
回复:DOE 配方非常感谢您的回复,这是解决我的问题的一个有趣的方法(我有点绝望)。只是一条评论,也许是您的一条建议,当 A 浓度范围 (10-15) 与 B 浓度范围 (800-1000) 确实不同时,您如何处理。 A 和 B 之间的运行次数确实不平衡(A 运行 3 次,B 运行 17 次)。我觉得很奇怪。下面是代码示例。
你怎么认为?
再次感谢您的及时回复。
美国能源部(
定制设计,
{添加响应(最大化, "Y", ., ., . ),
添加因子(连续, -1, 1, "X1", 0 ),
添加因子(连续, -1, 1, "X2", 0 ),
添加因子( 分类, {"A", "B"}, "X3", 0 ),
添加因子( Continuous, 10, 1000, "Conc", 0 ), 设置随机种子( 2111591509 ),
开始次数( 1453 ),添加术语( {1, 0} ),添加术语( {1, 1} ),
添加术语( {2, 1} ), 添加术语( {3, 1} ), 添加术语( {4, 1} ), 添加术语( {1, 2} ),
添加术语( {1, 1}, {2, 1} ), 添加术语( {2, 2} ), 添加术语( {1, 1}, {3, 1} ),
添加术语( {2, 1}, {3, 1} ), 添加术语( {1, 1}, {4, 1} ),
添加术语( {2, 1}, {4, 1} ), 添加术语( {3, 1}, {4, 1} ), 添加术语( {4, 2} ),
添加别名项( {1, 1}, {2, 1} ), 添加别名项( {1, 1}, {3, 1} ),
添加别名项( {1, 1}, {4, 1} ), 添加别名项( {2, 1}, {3, 1} ),
添加别名项( {2, 1}, {4, 1} ), 添加别名项( {3, 1}, {4, 1} ),
设置样本大小( 20 ),不允许的组合(
X3 == 1 & 浓 > 15 | X3 == 2 & 浓度 < 800
)、最优性准则( 2 )、进行设计、设置运行顺序(随机化)、
制作表格}
)回复:DOE 配方请参阅此讨论。将浓度定义为跨越 A 和 B 的较宽范围。设置一个约束,以便当因子设置为 A 时,排除浓度范围的一部分。对设置为 B 的因子执行相同操作。
我认为它应该回答您的问题,但如果您仍然需要帮助,请告诉我们。
回复:DOE 配方我忘了补充一点,两种成分的浓度范围是不同的。
事实上我的限制如下:
如果 CA>0 则 CB=0
或者如果 CB>0 则 CA=0
其中 CA(分别 CB)是 A(分别 B)成分的浓度。
当我尝试实现此约束时,我收到以下消息:
“模型大纲中选定的项与之前的项线性相关。请按“删除项”按钮,然后重试。”
我希望它有帮助。
在此先感谢您的帮助。
回复:DOE 配方需要约束吗?似乎您可以有两个因素:成分(A 或 B)的分类因素和浓度的连续因素。这样,A 和 B 就不会出现在同一个配方中。