こんにちは
3 つのグループ間で変数の平均を比較するために、一元配置分散分析を実施しました。論文の査読者は、効果の大きさを測定するために Cohen の d 検定を求めました。JMP でこの測定値を取得するにはどうすればよいでしょうか? JMP pro V17 を使用しています。
CohenのdはJMP 17では直接利用できません。ただし、JMP 18では
分析=>YをXでフィット=>平均/分散分析/プールt
JMP 17では、JSLを使用して計算する必要があります。これは、すべてのグループのペアに対してCohenのdを生成するサンプルスクリプトです。
または
Names Default To Here( 1 );
dt = open("$SAMPLE_DATA/big class.jmp");
var = "height";
grp = "sex";
Eval(
Parse(
"Summarize(
dt,
bygp = by( :" || grp || " ),
stds = Std Dev( :" || var
|| " ),
ns = Count( :" || var || " ),
means = Mean( :" || var || " )
);"
)
);
nw = New Window( "Cohen's d", vlb = V List Box() );
For( i = 1, i <= N Items( bygp ) - 1, i++,
For( k = i + 1, k <= N Items( bygp ), k++,
pooled = Sqrt(
(stds[i] ^ 2 * (ns[i] - 1) + stds[k] ^ 2 * (ns[k] - 1)) / (ns[i] + ns[k] - 2)
);
dif = Abs( means[i] - means[k] );
cohen = (dif / pooled);
stats = Matrix( dif ) |/ Matrix( pooled ) |/ Matrix( cohen );
compare = {};
insert into( compare, bygp[i] || " vs. " || bygp[k]);
ob = Outline Box( grp || "=" || bygp[i] || " vs. " || grp || "=" || bygp[k],
Table Box(
String Col Box( "", compare ),
String Col Box( "", {"Difference", "Pooled Std Dev", "Cohen's d"} ),
Number Col Box( "", stats )
)
);
vlb << append( ob );
)
);
ありがとう
役に立ちました。JMP 18で2グループではなく3グループで取得する方法があるかどうか疑問に思っています。プールされたtは2グループの平均を比較する場合にのみ使用できますが、私の場合は3グループの平均を比較するために一元配置分散分析を実行しました。
こんにちはサミラ 、
2 つのグループ以外では、このオプションは JMP 18 ではネイティブでは利用できません。ウィッシュ リストをご覧ください。このニーズはすでに提起され、認識されています: 「Fit Y by X」の Cohen の d
将来の 18.X バージョンまたは JMP 19 でのリリースを期待しましょう。
こんにちはサミラ : ちなみに、ANOVA の 2 つ以上のグループの場合、RMSE (二乗平均平方根誤差) はプールされた SD です。
ありがとう
これは、Cohen の d を手動で計算するのに非常に役立ちました。さらに明確にするために、適合表の概要で RMSE を見つけました。
