Bu hesaplayıcı ne işe yarar?
Bu araç, a ve b şekil parametreleri için bir x aralığında birbiriyle yakından ilişkili iki özel fonksiyonu tablo haline getirir. Alt tamamlanmamış beta fonksiyonu \(B_x(a,b)\), beta çekirdeğinin 0'dan x'e kadar olan integralidir; düzenlenmiş tamamlanmamış beta fonksiyonu \(I_x(a,b)\) ise bu integralin tam beta fonksiyonu \(B(a,b)\)'ye bölünmüş halidir. \(I_x(a,b)\) aynı zamanda Beta(a,b) dağılımının kümülatif dağılım fonksiyonudur (CDF) ve binom, Student t, F ile diğer pek çok dağılımın CDF'lerinin temelini oluşturur.
Nasıl kullanılır?
a ve b şekil parametrelerini girin (her ikisi de pozitif olmalıdır). Bir başlangıç x değeri (x'in başlangıç değeri), bir Artış adımı ve Tekrar sayısı (satır sayısı) seçin. i. satır şu formülü kullanır: \(x = \text{başlangıçX} + i \times \text{adım}\). Değerler [0, 1] aralığında tutulur; x değeri 1'e ulaştığında tablo durur. Üstteki kutu tam beta \(B(a,b)\) ile son satırı gösterirken, kaydırılabilir tablo her x değerini \(B_x(a,b)\) ve \(I_x(a,b)\) ile birlikte listeler.
Formülün açıklaması
Tam beta, log-gama fonksiyonu için Lanczos yaklaşımı kullanılarak kararlı biçimde şöyle hesaplanır:
$$B(a,b) = \exp(\text{lgamma}(a) + \text{lgamma}(b) - \text{lgamma}(a+b))$$Düzenlenmiş değer \(I_x(a,b)\), klasik Numerical Recipes sürekli kesir yöntemiyle (Lentz algoritması) değerlendirilir: \(bt = \exp(\text{lgamma}(a+b) - \text{lgamma}(a) - \text{lgamma}(b) + a\cdot\ln x + b\cdot\ln(1-x))\) olmak üzere, \(x < (a+1)/(a+b+2)\) iken \(I_x = bt\cdot\text{betacf}(a,b,x)/a\), aksi halde \(I_x = 1 - bt\cdot\text{betacf}(b,a,1-x)/b\) kullanılır. Son olarak \(B_x(a,b) = I_x(a,b) \times B(a,b)\). \(\log(0)\) durumundan kaçınmak için uç noktalar \(I_0 = 0\) ve \(I_1 = 1\) tam olarak ele alınır.
Çözümlü örnek
a = 1, b = 3 için:
$$B(1,3) = \frac{\Gamma(1)\Gamma(3)}{\Gamma(4)} = \frac{1\cdot 2}{6} = \frac{1}{3} \approx 0{,}333333$$x = 0,5 noktasında kapalı form
$$B_x = \frac{1 - (1-x)^3}{3} = \frac{1 - 0{,}125}{3} = 0{,}291667$$verir; dolayısıyla
$$I_x = \frac{0{,}291667}{0{,}333333} = 0{,}875$$olur. Simetri kontrolü \(I_x(1,3) = 1 - (1-x)^3 = 1 - 0{,}125 = 0{,}875\) sonucu doğrular.
Sıkça sorulan sorular
Bx ile Ix arasındaki fark nedir? \(B_x(a,b)\) ham integraldir; \(I_x(a,b) = B_x(a,b)/B(a,b)\) ise 0 ile 1 arasında olacak şekilde normalleştirilmiştir.
a ve b neden pozitif olmalı? İntegral ve gama fonksiyonları yalnızca \(a > 0\) ve \(b > 0\) için yakınsar.
Adımım x'i 1'in ötesine taşırsa ne olur? Her x değeri 1'e sabitlenir ve tablo x = 1 noktasında durur; burada \(B_x(a,b) = B(a,b)\) ve \(I_x(a,b) = 1\) olur.