[Rd] choose fails a fundamental property of binomial coefficients (PR#11035)
maechler at stat.math.ethz.ch
maechler at stat.math.ethz.ch
Wed Mar 26 11:02:37 CET 2008
>>>>> "JL" == Jerry Lewis <Jerry.Lewis at biogenidec.com>
>>>>> on Wed, 26 Mar 2008 05:10:04 +0100 (CET) writes:
JL> Full_Name: Jerry W. Lewis Version: 2.7.0 (2008-03-23
JL> r44847) OS: Windows XP Professional Submission from:
JL> (NULL) (71.184.230.48)
JL> choose(n,k) = choose(n,n-k) is not satisfied if either
JL> 1. n is a negative integer with k a positive integer
JL> (due to automatically returning 0 for n-k<0)
JL> 2. n is not an integer with k a positive integer (due to
JL> rounding n-k to an integer, compounded by automatically
JL> returning 0 if n<0 which implies n-k<0)
I think you have not really understood from help(choose)
that we have
choose(n,k) := n(n-1)...(n-k+1) / k!
defined for non-negative integer k and arbitrary real n
This definition has the very important property that it is used
in the 'binomial theorem' and everything depending on that
theorem,
(1 + x) ^ \alpha = \sum_{k=0}^\infty choose(x,k) x^k
which holds for all |x| < 1 for all \alpha \in \R
-------
Your "fundamental" property does only hold when n is
non-negative integer.
Definitely no bug in R here!
Martin Maechler, ETH Zurich
More information about the R-devel
mailing list