<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />我很好奇是否有人知道并且可以向我解释JMP如何为自定义设计计算其非中心性参数。 我知道计算它的一般形式，但我不能完全按照JMP使用的数学计算。 原因是df 在JMP 10中，默认情况下，用于计算分割图设计功率的误差是错误的。 虽然我可以调整df 手动，如果我可以编写一个脚本来从一开始就正确地进行计算，但是设计评估平台在JMP 10中不能编写脚本（不是我无论如何都要弄清楚 - 甚至尝试使用<< Get Scriptable Object'返回一个空集）所以我必须能够在JSL中“手动”进行数学运算。
提前致谢！
琼斯博士，
谢谢您的回复。 您是否介意在裂缝结构存在的情况下解释它是如何工作的？我使用自定义设计器创建了一个简单的5因子（2个HTC，3个ETC因子），并在每个WP中获得了6个WP和8个重复的设计，以实现完美复制的设计。 我使用REML拟合模型，然后保存编码表以获得X矩阵。 但是，我不能用6个整体图的编码解决X'X，因为它们（故意）与HTC因素混淆。 如果我按照你为ETC因素指定的那样做，我会得到与JMP相同的ncp。 如果我将HTC因子降低到只有6个独立设置来对待它们就好像只有6次运行（实际上是正在发生的事情），我得到了一个很好的甚至6为ncp。 但是，JMP给了我5.333和4.741（确切的值并不重要）。 所以，还有其他事情发生，我不确定它是什么或如何解释它。
再次，我感谢您花时间提供上述代码来获取ncp。
JMP9和JMP10之间的功率计算有所不同。在JMP9中，计算信噪比的功率，其中回归系数是信号，误差标准偏差是噪声。在JMP 10中，2级因子的信号是从因子的低值到因子的高值时响应的变化。所以现在信号是1自由度效应的回归系数的2倍。
以下是您需要在JMP 9中计算1自由度效应的功率计算的JSL
设X为模型的X矩阵。
然后，系数的方差是
varBetaHat = vecdiag（inverse（X` * X））;
alpha = 0.05; //默认显着性水平
denominatorDF = 11; // df用于4个因子中的16次正交设计
Fcrit = F分位数（1-α，1，分母DF）; // F分布的临界值
SNRatio = 1; //默认SNRatio
noncentrality =（SNRatio ^ 2）/ varBetaHat [2]; //第一个主要效果的非中心性假设第一个效果是截距
power = 1 - F Distribution（Fcrit，1，denominatorDF，noncentrality）;
下面是JSL，您需要使用相同的X矩阵计算JMP 10中的1自由度效应的功率计算。
varContrast = 4 * vecdiag（inverse（X` * X））; //由于对比度是系数的两倍，因此其方差大4倍。
alpha = 0.05; //默认显着性水平
denominatorDF = 11; // df用于4个因子中的16次正交设计
Fcrit = F分位数（1-α，1，分母DF）; // F分布的临界值
SNRatio = 2; //默认SNRatio是容纳的两倍大，并为电源提供相同的默认值。
noncentrality =（SNRatio ^ 2）/ varContrast [2]; //第一个主要效果的非中心性假设第一个效果是截距
power = 1 - F Distribution（Fcrit，1，denominatorDF，noncentrality）;
