[R] Bug in qnorm or pnorm?

Trenkler, Dietrich dtrenkler at nts6.oec.uni-osnabrueck.de
Fri Aug 6 15:51:29 CEST 2004



> -----Original Message-----
> From:	Deepayan Sarkar 
> Sent:	Friday, August 06, 2004 3:31 PM
> To:	r-help at stat.math.ethz.ch
> Subject:	Re: [R] Bug in qnorm or pnorm?
> 
> 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 perfectly acceptable to me for all reasonable definitions of 
> Inf.
> 
> Deepayan
> 
> 
	[Dietrich Trenkler]  Yes of course, you're right. I meant 
	qnorm(pnorm(8.3))  should deliver Inf -- as it does.  
	Thank you.

	Dietrich.




More information about the R-help mailing list