「統計学関連なんでもあり」の過去ログ---036

★ t分布の両側確率について ★

9015.　t分布の両側確率について　龍　2006/01/13 (金) 14:02
└9017.　Re: t分布の両側確率について　青木繁伸　2006/01/13 (金) 14:35
　└9022.　Re^2: t分布の両側確率について　龍　2006/01/13 (金) 15:32
　　└9028.　Re^3: t分布の両側確率について　青木繁伸　2006/01/13 (金) 15:51
　　　└9032.　Re^4: t分布の両側確率について　龍　2006/01/13 (金) 16:25

9015.　t分布の両側確率について　龍　　2006/01/13 (金) 14:02

はじめまして。こちらでお伺いしてよいものかの判断も付かない程，統計学に無知なものですが，どうぞよろしくお願いします。

以下のプログラムの意味が分からないので，教えていただきませんか？
# t分布の両側確率　<a name="txp">txp.awk</a>

function txp(t, df, i, t1, t2, p, w, m2pi)
{
m2pi = 0.636619772367581343076
t1 = abs(t)/sqrt(df)
t2 = 1.0/(1.0+t1*t1)
if ((df%2) == 0) {
w = t1*sqrt(t2)
p = 1.0-w
for (i = 2; i <= df-2; i += 2) {
p -= w *= t2*(i-1)/i
}
}
else {
p = 1.0-m2pi*atan2(t1, 1)
if (df >= 3) {
w = m2pi*t1*t2
p -= w
for (i = 3; i <= df-2; i += 2) {
p -= w *= t2*(i-1)/i
}
}
}
return (p < 0.0 && abs(p) < 1e-10) ? 0.0 : p
}

　　　　　[このページのトップへ]

9017.　Re: t分布の両側確率について　青木繁伸　　2006/01/13 (金) 14:35

> 以下のプログラムの意味が分からないので，教えていただきませんか？

何をどのようにお教えすればよいかわかりません。

せめて，何がしたいのかを質問し，あなたが利用できる環境が書いてあれば，代替法をお教えすることはできますが。

　　　　　[このページのトップへ]

9022.　Re^2: t分布の両側確率について　龍　　2006/01/13 (金) 15:32

> 何をどのようにお教えすればよいかわかりません。

ご返答ありがとうございます。

プログラムの使い方は分かっているのですが，どのような公式からプログラムが作成されたかがわかりません。

もしよろしければ，その公式を教えていただけませんか？

　　　　　[このページのトップへ]

9028.　Re^3: t分布の両側確率について　青木繁伸　　2006/01/13 (金) 15:51

近似公式がどのようなものであるかは，プログラムを読めばわかると思いますが，その出所がどこであったかはちょと忘れました。はるか25年ほども前の話なので。

　　　　　[このページのトップへ]

9032.　Re^4: t分布の両側確率について　龍　　2006/01/13 (金) 16:25

ご返答ありがとうございます。

　　　　　[このページのトップへ]