8512_zh-CN

liliconcom

JMP 上的平衡不完全区组设计<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />

如何在具有 13 个属性的 JMP 软件上执行平衡不完全区组设计 (BIBD)，形成 13 个区组或选择集，每个属性在所有选择集中出现 4 次，每个选择集包含 4 个属性，每个属性在每个属性中出现一次其他： (b, r, k, lambda) (13, 4,4,1) 其中 b 是块数，r 是每个级别的重复次数，k 是每个选择集中的属性数（块大小）和 lambda是配对频率？

任何帮助将不胜感激，谢谢。

布莱恩.

回复：JMP 上的平衡不完全区组设计

感谢您的及时回复。

回复：JMP 上的平衡不完全区组设计

沃尔特·帕茨科夫斯基

数据分析公司总裁

使用 DOE 平台创建 BIBD——我该怎么做？

你好，

我有一个关于使用 DOE 平台创建平衡不完整区组设计 (BIBD) 的问题。 我有 13 个处理，我想要一个包含 13 个块或游程的设计，每个块 4 个处理，每个处理重复 4 次，块中没有重复，每对重复一次。 我知道存在这样的设计（请参阅 Cochrane 和 Cox 的书或使用 SAS 宏），但我无法让 JMP 来做到这一点。 如果您是设计专家并且了解 DOE 平台，您能否告诉我这在 JMP 中是否可行以及如何实现？

谢谢，

沃尔特·帕茨科夫斯基

2年

• 肖恩·舒伯特可能有一种方法可以使用自定义设计功能进行设置，将运行分组为块大小 4，难以更改和易于更改的因素，但我可能会尝试这些选项，但不会得到满足所有要求的设计限制，所以我最终会手动构建它。
  
  我会创建一个 13 次处理 x 4 次完全随机重复的设计。 然后，我从第 1 行开始将处理分配给块，并询问该处理是否遵循规则。 如果不是，我会跳过它并将其添加到下一个块。 最终所有的治疗都会被分配。
  
  不是最好的答案，但我至少知道我会很快得到满足要求的东西。

  2年
• 
  

  沃尔特·帕茨科夫斯基早上好，
  
  我几乎决定 JMP 无法创建 BIBD —— 似乎很少
  包可以。 您关于“蛮力”方法的建议可能会奏效。 它
  编写脚本需要付出一些努力，但这样做可能很有趣。
  
  谢谢你的建议，
  
  沃尔特
  
  
  _____________________________________________________
  
  JMP Professional Network 于 2013 年 9 月 21 日上午 9:39 写道：

  2年
• 
  

  肖恩·舒伯特沃尔特，
  
  这是一个您将严格按照指定运行多次的设计吗？ 我喜欢编程挑战，但这可能不是一个容易实现的算法。

  2年

• 沃尔特·帕茨科夫斯基你好，
  
  它会运行很多次，这就是为什么我不打算尝试它
  我自己有一段时间了。 我需要 BIBD 来进行不同的市场研究
  项目。 JMP 的最佳方法是使用 SAS 接口
  SAS 宏 %mktBIBD 生成一个（如果存在）。
  
  再次感谢，
  
  沃尔特
  
  
  _____________________________________________________
  
  
  JMP Professional Network 于 2013 年 9 月 22 日中午 12:30 写道：

  
  

  2年
• 
  

  布雷迪·布雷迪嗨沃尔特，
  
  此时，JMP 的算法肯定不会生成 BIBD。 不过，我将分享一种在某些情况下有效的替代方法。 （例如，在您上面描述的情况下，它反复对我和其他人有效。）
  
  JMP的DOE开发人员一直在思考如何在JMP的未来版本中增强DOE平台。 客户的要求在很大程度上影响了他们的决策；如果您愿意，请与他们分享您对 BIBD 支持的渴望： support@jmp.com
  
  与此同时，也许这种方法会对您有所帮助。 再说一次，这不是保证，但绝对值得一试：
  
  * 在 DOE 菜单中调用自定义设计器。
  
  * 创建 13 个连续因子，低水平 = 0，高水平 = 1。 （这对于下一步的约束很重要。）
  
  * 添加2个约束：Sum(factors) >=4；总和（因子）<=4。 为此，请为所有系数输入“1”，选择适当的不等式，然后在不等式的右侧输入 4。 （每次要添加约束时，请按“添加约束”按钮，该按钮在展开“定义因子约束”轮廓框时会显示。
  
  * 单击截距的可估计性“必要”并将其更改为“如果可能”，从而使最小游程为 13。 选择此最小运行次数。 （对于块数多于因子数的设计，无需执行此步骤。）
  
  * 创建设计并将设计点舍入为整数（某些值不完全是 0 或 1。）
  
  
  干杯，感谢您使用 JMP！
  
  布雷迪

  2年

• 布雷迪·布雷迪上述内容的后续：
  
  在步骤 4 中，只需继续并完全删除截距项（无论是否确实需要，它似乎提高了成功的可能性；通过使用此调整，在另一个示例中，我能够找到我失败的 BIBD前。）

  
  

  2年
• 
  

  沃尔特·帕茨科夫斯基你好，
  
  感谢您对此的帮助。 我一定会尝试一下。 我想
  请参阅 JMP 中添加的 BIBD，因为它是市场研究中的常见设计。
  
  谢谢，
  
  沃尔特
  
  
  _____________________________________________________
  
  
  JMP Professional Network 于 2013 年 9 月 23 日下午 4:12 写道：

  2年

• 比阿特丽斯·布鲁姆嗨沃尔特，
  
  我的经验是：如果 SAS 可以完成 DoE 的工作，那么 JMP 也可以完成。
  以下内容应该可以满足您的所有需求。
  
  美国能源部 - 定制设计
  添加 1 个分类因素“治疗”，共 13 个级别
  添加 1 个阻塞因子“Block”，其中“每块运行 4 次”
  （重要提示：不要使用“组随机进入 4 块”。 这将生成与 Brady 的方法相同的设计）
  运行次数：用户指定：52（每块 13 trt * 4 = JMP 总共 52 次运行）
  单击继续
  制作表格
  
  使用制表检查结果表。 它是平衡的。 据我所知，每一对都是重复的（我没有检查全部），而布雷迪的方法没有提供这一点。 但我仍然不确定我在这里错过了什么，因为布雷迪的答案往往是优秀和终极的。
  
  希望这可以帮助，
  比亚

  2年

• 布雷迪·布雷迪大家好，
  
  Beatrice提到的方法是我最初尝试的方法。 这种设计是平衡的，因为每个处理出现在相同数量的块中，但 BIBD 还要求每对处理出现在相同数量的块中。
  
  由于创建定制设计的算法固有的随机性，实际上没有什么可以阻止 13 级分类方法创建 BIBD，但发生这种情况的可能性似乎比连续方法低得多。 我（和其他几个一直在研究这个问题的人）还没有使用分类方法获得 BIBD（无论是添加阻塞因子，还是简单地单击单选按钮将运行分组为块），但已经创建了 BIBD （对于本例）每次都使用连续方法。
  
  下面，我展示了 3 个有代表性的 x`*x 矩阵以进行比较。 （对于分类示例，您首先必须拆分设计表以形成 13x13 矩阵。）
  
  沿着主对角线，所有条目均为 4。 这意味着无论使用哪种方法，13 个因子中的每一个都会出现在 4 个块中。
  
  非对角线条目记录了给定因子对出现的块数。 请注意，在前两个示例中，主对角线之外有 0、1 和 2。 这意味着一些处理对出现在零个块中，一些出现在一个块中，一些出现在两个块中。 仅当所有非对角线条目都相同时，该设计才是 BIBD。 所以，只有在第 3 种情况下我们才有 BIBD。
  
  示例 1：来自自定义设计的 x`*x 矩阵，具有 52 次运行的 13 级分类因子和 4 次运行的每块分块因子：
  [4 1 2 1 1 2 0 1 0 0 2 1 1,
  1. 4 1 1 0 0 2 2 1 1 1 2 0,
  
  1. 1 4 0 0 1 1 2 0 1 1 1 2,
  
  1. 1 0 4 3 1 2 0 1 1 0 1 1,
  
  1. 0 0 3 4 2 2 0 1 1 1 0 1,
  
  1. 0 1 1 2 4 0 1 1 1 1 1 1,
  
  0 2 1 2 2 0 4 1 1 1 1 1 0,
  1. 2 2 0 0 1 1 4 2 1 1 0 1,
  
  0 1 0 1 1 1 1 2 4 2 1 1 1,
  0 1 1 1 1 1 1 1 2 4 0 2 1,
  1. 1 1 0 1 1 1 1 1 0 4 1 2,
  
  1. 2 1 1 0 1 1 0 1 2 1 4 1,
  
  1.0 2 1 1 1 0 1 1 1 2 1 4]
  
  示例 2：来自自定义设计的 x`*x 矩阵，具有 13 级分类因子的 52 次运行，使用对话框按钮“将运行分组为大小为 4 的随机块”：
  [4 2 1 2 1 1 0 1 2 1 0 1 0,
  1. 4 1 1 1 1 1 2 1 1 0 1 0,
  
  1. 1 4 1 1 1 2 0 1 1 1 1 1,
  
  1. 1 1 4 0 0 1 2 1 1 1 1 1,
  
  1. 1 1 0 4 2 1 2 1 0 1 1 1,
  
  1. 1 1 0 2 4 1 1 1 1 1 1 1,
  
  0 1 2 1 1 1 4 0 1 0 2 1 2,
  1. 2 0 2 2 1 0 4 0 1 1 1 1,
  
  1. 1 1 1 1 1 1 0 4 2 0 0 2,
  
  1. 1 1 1 0 1 0 1 2 4 2 1 1,
  
  0 0 1 1 1 1 2 1 0 2 4 2 1,
  1. 1 1 1 1 1 1 1 0 1 2 4 1,
  
  0 0 1 1 1 1 2 1 2 1 1 1 4]
  
  
  示例 3：来自自定义设计的 x`*x 矩阵，具有 13 个连续因子的 13 次约束游程（如本线程前面所述）：
  [4 1 1 1 1 1 1 1 1 1 1 1 1,
  1. 4 1 1 1 1 1 1 1 1 1 1 1,
  
  1. 1 4 1 1 1 1 1 1 1 1 1 1,
  
  1. 1 1 4 1 1 1 1 1 1 1 1 1,
  
  1. 1 1 1 4 1 1 1 1 1 1 1 1,
  
  1. 1 1 1 1 4 1 1 1 1 1 1 1,
  
  1. 1 1 1 1 1 4 1 1 1 1 1 1,
  
  1. 1 1 1 1 1 1 4 1 1 1 1 1,
  
  1. 1 1 1 1 1 1 1 4 1 1 1 1,
  
  1. 1 1 1 1 1 1 1 1 4 1 1 1,
  
  1. 1 1 1 1 1 1 1 1 1 4 1 1,
  
  1. 1 1 1 1 1 1 1 1 1 1 4 1,
  
  1. 1 1 1 1 1 1 1 1 1 1 1 4]
  
  我希望这有帮助！
  
  干杯，
  布雷迪

  2年

• 肖恩·舒伯特我可以设想将标准设计表从 JMP 复制到 Excel 中，然后在 VBA 中对约束进行编程以获得阻塞的运行顺序。
  
  一旦设置正确，每次我按 F9 时，我都会得到一个新的运行命令。

回复：JMP 上的平衡不完全区组设计

LinkedIn 页面不再可用。是否有其他资源定义如何在 JMP 中开发 BIBD？

回复：JMP 上的平衡不完全区组设计

嗨布莱恩，

这个问题不久前在 LinkedIn 上出现并得到了关注相当不错的一组回应。

-杰夫