[R] pchisq error

[Bernardo Rangel Tura]

On Mon, 2009-01-19 at 09:54 +0100, Jeremy Silver wrote:
> Dear R experts,
(...)
> > pchisq(5.464342,1,lower.tail = FALSE)
> 
> [1] 0.01940836
> 
> > reproduceError(5.464342)
> 
> stat = 5.464342, p = 0.019408
> 
> > pchisq(5.464342,1,lower.tail = FALSE)
> 
> [1] NaN
> 
> Warning messages:
> 
> 1: In pchisq(5.464342, 1, lower.tail = FALSE) :
> 
>    ** NON-convergence in pgamma()'s pd_lower_cf() f= nan.
> 
> 2: In pchisq(q, df, lower.tail, log.p) : NaNs produced
> > 
(...)
> 
> > sessionInfo()
> 
> R version 2.8.0 (2008-10-20) 
> 
> i686-pc-linux-gnu 
> 
> locale:
> 
> LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C
> 
> attached base packages:
> 
> [1] stats     graphics  grDevices utils     datasets  methods   base     

Hi Jeremy,

In my computer your error is not occur. Look This:

> pchisq(5.464342,1,lower.tail = FALSE)
[1] 0.01940836

> pchisq(5.464342,1,lower.tail = FALSE)
[1] 0.01940836

> sessionInfo()
R version 2.8.1 Patched (2009-01-16 r47630) 
x86_64-unknown-linux-gnu 

locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods
base     


Well, do you already try update your R?


-- 
Bernardo Rangel Tura, M.D,MPH,Ph.D
National Institute of Cardiology
Brazil