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

★ Re: カイ自乗分布の確率密度関数について ★

　166　Re: カイ自乗分布の確率密度関数について　　名無しさん　　2001/07/09 (月) 11:23
　　167　Re^2: カイ自乗分布の確率密度関数について　　名無しさん　　2001/07/09 (月) 11:27
　　　169　Re^4: カイ自乗分布の確率密度関数について　　名無しさん　　2001/07/09 (月) 13:37

166.　Re: カイ自乗分布の確率密度関数について　　名無しさん　　2001/07/09 (月) 11:23

> ところが，こいつが見事にカイ自乗分布になりません!!

どこがどういう風に変なのか，書いてください。

> Public Function CHIDIST(X As Long, DF As Integer) As Double

第一引数は，double だと思うが?
それと，CHIDIST というのは Excel の関数と名前が競合しているので，Excel の関数として使うと，まずいことが起こるのでは?
別の名前を付けてみればいいと思いますが。
また，GammaLn は別の所で定義済みなんでしょうか?もし，Excel の関数を援用しようとするなら，Application.GammaLn としないといけないのではないかな。

Public Function ccc(X As Double, DF As Integer) As Double
ccc = ((X ^ (DF / 2 - 1)) * (Exp(-X / 2))) / _
((2 ^ (DF / 2)) * (Exp(Application.GammaLn(DF / 2))))
End Function

で動くと思いますが。

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

167.　Re^2: カイ自乗分布の確率密度関数について　　名無しさん　　2001/07/09 (月) 11:27

> ccc = ((X ^ (DF / 2 - 1)) * (Exp(-X / 2))) / _
> ((2 ^ (DF / 2)) * (Exp(Application.GammaLn(DF / 2))))

不要なカッコが多くて，わかりづらいですね。
LISP みたいです(^_^;)

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

169.　Re^4: カイ自乗分布の確率密度関数について　　名無しさん　　2001/07/09 (月) 13:37

> ところで，上記のようにコードを修正したのですが，自由度(上記関数の引数
> でいうところの「DF」)を「1」にしたところ，「引数が不正です」とのエラー
> がでてしまいました。

df が 1 で，x が 0 のときの話でしょう?
df = 1 でも，x > 0 ならエラーにならないと思いますが?

0 ^(1/2-1) の計算結果がどうなるかを考えると，上記エラーの原因がわかります。

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