Displaying 4 results from an estimated 4 matches for "lfastchoose".
2000 Mar 24
3
quantiles of the hypergeometric distribution (PR#502)
Hello!
I use R-version 1.0.0
To get the 0.95 quantile of the hypergeometric
distribution with the parameters m=45000,n=5000 and
k=600 I use the R-command
> qhyper(0.95,45000,5000,600).
The value obtained is 600. However, the true value
is 552. The latter can be obtained for example by
calling the corresponding distribution function
with the R commands
> x<-540:580
>
1998 Feb 23
0
R-beta: Hypergeometric Probabilities
...abilities as the
parameters of the hypergeometric distribution get large. For example, when
N1=1345, N2=1055, and n=1330, phyper returns either 0 or 1, and nothing in
between.
Looking at phyper.c, it's clear what's happening. First a term (called
"term") is computed:
term = exp(lfastchoose(NR, xr) + lfastchoose(NB, xb)
- lfastchoose(N, n));
Then in a while loop, this is updated by multiplication, and each
subsequent "term" is added:
while(xr <= x) {
sum += term;
xr++;
NB++;...
2001 Feb 08
2
dnbinom(,size<1,)=0 (PR#842)
This came up on r-help but indicates a bug.
dnbinom(x,n,p) calls dbinom_raw(n-1,...)
which returns 0 for n<1.
-thomas
---------- Forwarded message ----------
Date: Thu, 08 Feb 2001 17:10:23 +0000
From: Yudi Pawitan <yudi@stat.ucc.ie>
To: Mark Myatt <mark@myatt.demon.co.uk>
Cc: R-Help <r-help@stat.math.ethz.ch>
Subject: Re: [R] Goodness of fit to Poisson / NegBinomial
2006 Mar 28
2
R 2.3.0 (alpha) on FreeBSD 6.1 fails make check-all
Hi Developers,
The alpha, compiles successfully, but it is failing make check-all (on
two seperate machines, both FreeBSD 6.1).
Here is the version string:
platform i386-unknown-freebsd6.1
arch i386
os freebsd6.1
system i386, freebsd6.1
status alpha
major 2
minor 3.0
year 2006
month 03
day 27
svn rev