Module: Distribution::F::Ruby_

Defined in:
lib/distribution/f/ruby.rb

Class Method Summary collapse

Class Method Details

.c_pdf(f, df) ⇒ Object 



6
7
8
 
# File 'lib/distribution/f/ruby.rb', line 6

def c_pdf(f,df)
  Distribution::ChiSquare.pdf(f,df)
end

.cdf(f, n1, n2) ⇒ Object

F-distribution interface



107
108
109
 
# File 'lib/distribution/f/ruby.rb', line 107

def cdf(f,n1, n2)
  1.0 - q_f(n1, n2, f)
end

.p_value(y, n1, n2) ⇒ Object 



110
111
112
 
# File 'lib/distribution/f/ruby.rb', line 110

def p_value(y, n1, n2)
  pf(1.0 - y, n1, n2)
end

.pdf(x, d1, d2) ⇒ Object 



9
10
11
 
# File 'lib/distribution/f/ruby.rb', line 9

def pdf(x,d1,d2)
  Math.sqrt(((d1*x)**d1*(d2**d2)).quo((d1*x+d2)**(d1+d2))).quo( x*Math.beta(d1/2.0, d2/2.0))
end

.pf(q, n1, n2) ⇒ Object

Inverse CDF [x, infty)



58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
 
# File 'lib/distribution/f/ruby.rb', line 58

def pf(q, n1, n2)
  if(q < 0.0 || q > 1.0 || n1 < 1 || n2 < 1)
  $stderr.printf("Error : Illegal parameter in pf()!\n")
  return 0.0
end

if n1 <= 240 || n2 <= 240
eps = 1.0e-5
if(n2 == 1) then eps = 1.0e-4 end
fw = 0.0
s = 1000.0
loop do
fw += s
if s <= eps  then return fw end
if (qe = q_f(n1, n2, fw) - q) == 0.0 then return fw end
if qe < 0.0
  fw -= s
  s /= 10.0 #/
end
end
end

eps = 1.0e-6
qn = q
if q < 0.5 then qn = 1.0 - q
u = pnorm(qn)
w1 = 2.0 / n1 / 9.0
w2 = 2.0 / n2 / 9.0
w3 = 1.0 - w1
w4 = 1.0 - w2
u2 = u * u
a = w4 * w4 - u2 * w2
b = -2. * w3 * w4
c = w3 * w3 - u2 * w1
d = b * b - 4 * a * c
if(d < 0.0)
fw = pfsub(a, b, 0.0)
else
if(a.abs > eps)
  fw = pfsub(a, b, d)
else
  if(b.abs > eps) then return -c / b end
  fw = pfsub(a, b, 0.0)
end
end
fw * fw * fw
end
end

.pfsub(x, y, z) ⇒ Object

inverse of F-distribution ([2])



52
53
54
 
# File 'lib/distribution/f/ruby.rb', line 52

def pfsub(x, y, z)
(Math.sqrt(z) - y) / x / 2.0
end

.q_f(df1, df2, f) ⇒ Object

F-distribution ([1]) Integral over [x, infty)



14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
 
# File 'lib/distribution/f/ruby.rb', line 14

def q_f(df1, df2, f)
if (f <= 0.0) then return 1.0; end
if (df1 % 2 != 0 && df2 % 2 == 0)
return 1.0 - q_f(df2, df1, 1.0 / f)
end
cos2 = 1.0 / (1.0 + df1.to_f * f / df2.to_f)
sin2 = 1.0 - cos2

if (df1 % 2 == 0)
prob = cos2 ** (df2.to_f / 2.0)
temp = prob
i = 2
while i < df1
  temp *= (df2.to_f + i - 2) * sin2 / i
  prob += temp
  i += 2
end
return prob
end
prob = Math.atan(Math.sqrt(df2.to_f / (df1.to_f * f)))
temp = Math.sqrt(sin2 * cos2)
i = 3
while i <= df1
prob += temp
temp *= (i - 1).to_f * sin2 / i.to_f;
i += 2.0
end
temp *= df1.to_f
i = 3
while i <= df2
prob -= temp
temp *= (df1.to_f + i - 2) * cos2 / i.to_f
i += 2
end
prob * 2.0 / Math::PI
end