私達がJMPかProでこのテストを実行できるかどうか誰でも知っています... 顧客から依頼されました。ホワイトテスト - Wikipedia
私はSPSSのWebサイトで同様のことをするための手順に関する記事を見つけました。手順
チームの意見は?
M
私がAmgenにいたとき、私はBreusch-Paganのために短い脚本を書きました。それがどんな用途でもあるならば、ここにあなたは行きます:
tbl = Current Data Table();
cDlog = Column Dialog(
rCol = Col List("Residual", Max Col(1)),
pCol = Col List("Predicted Value", Max Col(1))
);
residCol = Column(cDlog["rCol"]);
predCol = Column(cDlog["pCol"]);
tbl << Fit Y by X(
Y(residCol),
X(predCol),
Fit Line({Save Residuals})
);
Column(tbl,N Col(tbl)) << Set Name("BP Resids");
bpResidCol = Column(tbl,N Col(tbl));
n = N Row(tbl);
zVec = predCol << Get as matrix;
eVec = (:BP Resids << Get as matrix);
e2Vec = eVec^2;
iVec = J(n,1,1);
zVec = iVec || zVec;
uProd = (e2Vec-Mean(e2Vec)*iVec);
eSqr = eVec`*eVec/n;
outVal = Summation(i=1,N Row(),
(residCol[i]^2 - eSqr)^2
);
e2Var = outVal/n;
bp = (uProd`*zVec*Inverse(zVec`*zVec)*zVec`*uProd)/e2Var;
pVal = 1-Chi Square Distribution(bp,2);
Show(bp,pVal)これは、直線性とバイアスを調べるときに含めると便利なテストです。
試験方法の線形範囲にわたる反復測定試験では、変動はしばしば入力変数の大きさに比例します。このテストで私は答えることができたと思います、「範囲にわたる異分散性は重要ですか?」
それは知っておくとよいでしょう。
残念ながら、JMP Proはそれを持っていないか、Whiteの異分散性の一貫した標準エラーを持っていません。JMPはWhiteテストとMacKinnon and White Estimatorの両方を実装することを強くお勧めします。