> -Original Message-
> From: Deepayan Sarkar
> Sent: Friday, August 06, 2004 3:31 PM
> To: [EMAIL PROTECTED]
> Subject: Re: [R] Bug in qnorm or pnorm?
>
> On Friday 06 August 2004 08:13, Trenkler, Dietrich wrote:
>
> > Given that pnorm(8.30) d
On Friday 06 August 2004 08:13, Trenkler, Dietrich wrote:
> Given that pnorm(8.30) delivers 1 shouldn't we get Inf
> for x<-8.30;x-qnorm(pnorm(x)) ?
Why?
> pnorm(8.30)
[1] 1
> qnorm(pnorm(8.30)) ## same as qnorm(1)
[1] Inf
> 8.30 - qnorm(pnorm(8.30)) ## same as 8.30 - Inf
[1] -Inf
This seems
> -Original Message-
> From: Trenkler, Dietrich
> Sent: Friday, August 06, 2004 3:13 PM
> To: 'r-help'
> Subject: [R] Bug in qnorm or pnorm?
>
> I found the following strange behavior using qnorm() and pnorm():
>
> > x<-8.21;x-qnorm
I found the following strange behavior using qnorm() and pnorm():
> x<-8.21;x-qnorm(pnorm(x))
[1] 0.0004638484
> x<-8.22;x-qnorm(pnorm(x))
[1] 0.01046385
> x<-8.23;x-qnorm(pnorm(x))
[1] 0.02046385
> x<-8.24;x-qnorm(pnorm(x))
[1] 0.03046385
> x<-8.25;x-qnorm(pnorm(x))
[1] 0.04046385
> x<-8.26;x-qn