Torah puzzle

Results are suspiciously close

(page 3)

Criticism

How to compute P2 ?

A formula is given in [WRR94, p.436] representing P2 by a finite sum, but the sum is long (n terms). If you have "Maple V" or another program for mathematical computation, just ask it to compute the incomplete gamma function divided by the (complete) gamma function:

     GAMMA (n,s)   GAMMA (163, 251.33)              -8
P  = ----------- = ------------------- = 0.1096 * 10   ;
 2    GAMMA (n)      GAMMA (163)

here s = 251.33 is the sum of logarithms of (inverted) corrected distances, and n = 163 is the number of summands, for the second experiment of [WRR94]. (In fact, the default precision of 10 digits is not enough for my "Maple V" to get the correct answer for so large numbers; I did it with increased precision.) Note that (n-1)! = GAMMA (n) . If you have no such program, or want to get the result more quickly, I recommend you the following asymptotic formula:

          n-1  -s
         s    e             1
P    =   -------- * ---------------------  .
 2a       (n-1)!        n-1         1
                    1 - --- ( 1 - ----- )
                         s        s-n+3

Its relative error is less than 0.002 whenever the result is less than 10-6. For examlpe, when s = 245 and n = 155 , the formula gives 0.2896*10-9 while the exact value is 0.2898*10-9. More examples (the first number is the approximate value, P2a, the second one is the exact value, P2):

145150155160165
230 .7513*10-9
.7518*10-9
.7458*10-8
.7464*10-8
.6284*10-7
.6291*10-7
.4519*10-6
.4525*10-6
.2789*10-5
.2795*10-5
235 .1084*10-9
.1084*10-9
.1195*10-8
.1195*10-8
.1117*10-7
.1118*10-7
.8913*10-7
.8923*10-7
.6101*10-6
.6110*10-6
240 .1468*10-10
.1469*10-10
.1793*10-9
.1794*10-9
.1857*10-8
.1859*10-8
.1641*10-7
.1643*10-7
.1244*10-6
.1245*10-6
245 .1872*10-11
.1873*10-11
.2529*10-10
.2530*10-10
.2896*10-9
.2898*10-9
.2829*10-8
.2831*10-8
.2369*10-7
.2371*10-7
250 .2252*10-12
.2253*10-12
.3358*10-11
.3359*10-11
.4245*10-10
.4247*10-10
.4574*10-9
.4577*10-9
.4225*10-8
.4229*10-8
255 .2562*10-13
.2563*10-13
.4209*10-12
.4211*10-12
.5861*10-11
.5864*10-11
.6957*10-10
.6961*10-10
.7077*10-9
.7081*10-9
260 .2762*10-14
.2763*10-14
.4991*10-13
.4992*10-13
.7643*10-12
.7646*10-12
.9975*10-11
.9979*10-11
.1115*10-9
.1116*10-9

A simple approximation for P2

ln P2b = - 20.55 + L - 0.0115 L2 where L = 17.19 - 0.365 ( s - 1.252 n ) .
The approximation is fitted to be applicable in a neighborhood of the values used in [WRR94]. Examples follow (the first number is the simple approximation, - ln P2b, the second one is the more exact approximation, - ln P2a ):

145150155160165
230 21.05
21.01
18.80
18.71
16.67
16.58
14.66
14.61
12.77
12.79
235 22.93
22.95
20.59
20.55
18.36
18.31
16.25
16.23
14.27
14.31
240 24.90
24.94
22.45
22.44
20.13
20.10
17.93
17.93
15.84
15.90
245 26.93
27.00
24.39
24.40
21.98
21.96
19.68
19.68
17.50
17.56
250 29.05
29.12
26.41
26.42
23.90
23.88
21.50
21.51
19.23
19.28
255 31.24
31.30
28.51
28.50
25.90
25.86
23.41
23.39
21.04
21.07
260 33.51
33.52
30.68
30.63
27.97
27.90
25.39
25.33
22.92
22.92

The simple approximation makes it easy, to estimate the contribution of any personality (or appellation, or word pair) into the overall proximity measure P2. For example, the personality considered on page 2 (Rabbi Avraham Av-Beit-Din of Narbonne) gives 12 word pairs, the corresponding sum of 12 logarithms being 4.828 + ... + 1.363 = 28.172 (as was said on page 2). Its contribution into L is determined by the value

s - 1.252 n = 28.172 - 1.252 * 12 = 13.15 .

Do you recognize the number? Look once again at page 1, the second list of outcomes. Yes, 13.15 appears there, it is just the last (maximal) element of the list. You may guess that each element of the list is s - 1.252 n for a personality, and you are right!

You may also guess that Rabbi Avraham Av-Beit-Din of Narbonne contributes more than any other personality to the overall proximity measure, and you are right, once again. Consider it in detail. The sum of the list (already calculated on page 2) is

s - 1.252 n = (-2.14) + (-2.00) + ... + 12.30 + 13.15 = 47.3 ,

which gives L = 17.19 - 0.365 * 47.3 = - 0.1 and ln P2 = - 20.55 + L - 0.0115 L2 = - 20.7 , that is, P2 = exp (-20.7) = 0.11 * 10-8 . It is just the overall proximity measure for the second experiment (calculated before as GAMMA(n,s)/GAMMA(n)).

What happens if we exclude the personality (Rabbi Avraham...)? Then s-1.252n decreases from 47.3 to 47.3 - 13.15 = 34.2, which gives L = 17.19 - 0.365 * 34.2 = 4.7 and ln P2 = - 20.55 + L - 0.0115 L2 = - 16.1 , that is, P2 = exp (-16.1) = 0.10 * 10-6 . The person contributes into P2 the factor 0.011 . The factor can be roughly estimated via the value s-1.252n of the personality: exp ( - 0.365 * 13.15 ) = 0.0082 . Why not 0.011 ? Because of the nonlinear correction 0.0115 L2 . The correction contributes the factor exp ( 0.0115 * ( 4.72 - 0.12 ) ) = 1.29 , and so, 0.0082 * 1.29 = 0.011 , as it should be.

After all, what is the relevance of the numbers on page 1 ?

Each number corresponds to a personality used by WRR. (Personalities that do not contribute are excluded.) The number is

s - 1.252 n ,

where (-s) is the logarithm of the product of all corrected distances for the personality, and n is the number of the distances. The relevant property of the linear form ( s - 1.252 n ) follows.

The close proximity of the sums (47.1 and 47.3) of the two lists of numbers (observed on page 1) is an equivalent reformulation of the close proximity (noted by Gil Kalai) of the two overall proximity measures.

back to page 2 back to Criticism back to Torah puzzle