<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />JMPがカスタムデザインの非中心性パラメータを計算する方法を誰かが知っていて説明できるのであれば、私は興味があります。 私はそれを計算するための一般的な形式を知っています、しかし私はJMPが使用しているものに合算するように私の数学を得ることはできません。 その理由は、dfが JMP 10では、分割プロット設計の検出力の計算に使用される誤差がデフォルトで間違っています。 私はdfを調整することができますが 手動でget-goから計算を正しく行うためのスクリプトを書くことができれば設計評価プラットフォームはJMP 10ではスクリプト化できません（とにかく理解できない - << Get Scriptable Objectを試しても）空のセットを返したので、私はJSLで「手作業で」数学を実行できなければなりません。
前もって感謝します！
ジョーンズ博士、
お返事ありがとうございます。 スプリットプロット構造の存在下でこれがどのように機能するかを説明していただけませんか。カスタムデザイナーを使用して単純な5つの要素（2つのHTC、3つのETC要素）を作成し、6つのWPと各WP内に8つの複製を含むデザインを取得し、完全に複製されたデザインを作成しました。 REMLを使用してモデルを近似した後、X行列を取得するためにコーディングテーブルを保存しました。 ただし、6つのプロット全体に使用されているコーディングではX'Xを解くことはできません。これらはHTC係数と（意図的に）混同されているためです。 あなたがETCファクタのために指定したように私がするならば、私はJMPがするのと同じncpを得ます。 HTC係数を6つの独立した設定にまとめて6つの実行しかないかのように扱うと（つまり、実際には何が進行しているのか）、ncpには6という素晴らしい値が得られます。 しかし、JMPは私に5.333と4.741を与えます（正確な値は重要ではありません）。 それで、他に何かが起こっています、そして、それが何であるか、またはそれをどう説明するかについてはわかりません。
繰り返しになりますが、ncpを取得するために上記のコードを提供していただいてありがとうございます。
電力計算に関してはJMP9とJMP10の間に違いがあります。JMP9では、回帰係数が信号で、誤差標準偏差がノイズであるSN比に対して、電力が計算されました。JMP 10では、2レベルの因子のシグナルは、因子の低い値から高い値に変化したときの応答の変化です。そのため、信号は1自由度の効果の回帰係数の2倍になります。
これは、自由度1の効果を得るためにJMP 9でべき乗計算を計算するために必要なJSLです。
XをモデルのX行列とします。
次に、係数の分散は
varBetaHat = vecdiag（inverse（X` * X））;
アルファ＝ ０．０５。//デフォルトの有意水準
分母ＤＦ ＝ １１。// 4因子の16ラン直交計画の場合はdf
Ｆｃｒｉｔ ＝ Ｆ四分位数（１ −α、１、分母ＤＦ）； // F分布の臨界値
ＳＮ比＝ １。//デフォルトのSN比
非中心度=（SN比^ 2）/ varBetaHat [2]; // 1番目の効果が切片であると仮定して、1番目の主効果の非心度
度数= 1 - F分布（Fcrit、1、分母DF、非心度）。
これは、同じX行列を使用して自由度1の効果を得るためにJMP 10でべき乗計算を計算するために必要なJSLです。
varContrast = 4 * vecdiag（inverse（X` * X））; //コントラストは係数の2倍なので、その分散は4倍大きくなります。
アルファ＝ ０．０５。//デフォルトの有意水準
分母ＤＦ ＝ １１。// 4因子の16ラン直交計画の場合はdf
Ｆｃｒｉｔ ＝ Ｆ四分位数（１ −α、１、分母ＤＦ）； // F分布の臨界値
ＳＮ比＝ ２。// default SNRatioは、収容するために2倍の大きさで、電力に同じデフォルト値を与えます。
非心度=（SN比^ 2）/ varContrast [2]; // 1番目の効果が切片であると仮定して、1番目の主効果の非心度
度数= 1 - F分布（Fcrit、1、分母DF、非心度）。
