help-gsl
[Top][All Lists]

## [Help-gsl] gsl_cdf_poisson_Q: strange values

 From: Stephan Lorenzen Subject: [Help-gsl] gsl_cdf_poisson_Q: strange values Date: Thu, 7 Mar 2019 05:58:35 +0100 User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:60.0) Gecko/20100101 Thunderbird/60.4.0

```Dear list,

I have N=1850923 drawings of characters which can be either A,T,C or G.
```
The probability of C is p=0.256903, so I would expect lambda=N*p=1850923*0.256903=475508 C's.
```Indeed, I got k=560073.
```
According to my understanding, the probability to get 560073 or more C's in 1850923 drawings should be gsl_cdf_poisson_Q(k-1, lambda)=gsl_cdf_poisson_Q(560072, 475508) I get slightly but not dramatically more than expected, so I would expect a p-value of, whatwever, 0.2 or so.
```However, I got 0.
```
I would expect the p-value in the same range than when getting 5 instead of 4 drawings, so I did a series:
```
gsl_cdf_poisson_Q(     5,      4) = 0.21487
gsl_cdf_poisson_Q(    56,     47) = 0.0859409
gsl_cdf_poisson_Q(   560,    475) = 6.65119e-05
gsl_cdf_poisson_Q(  5600,   4755) = 4.39096e-33
gsl_cdf_poisson_Q( 56007,  47550) = 1.46886e-311
gsl_cdf_poisson_Q(560072, 475507) = 0

The results starting from the second line are certainly not what I expect.
Probably I have a fundamental misunderstanding of gsl_cdf_poisson_Q?
```
What would be the proper way to use it, or which function should I use instead?
```
```