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