rnorm() with sd=0

[Ben Bolker]

  Is there a compelling reason why

rnorm(1,mean=m,sd=0)

can't return m instead of NaN?

  The following (trivial) patch makes the change.

*** rnorm.c.orig	Sat Jan  9 15:39:41 1999
--- rnorm.c	Sat Jan  9 15:39:50 1999
***************
*** 35,43 ****
  #ifdef IEEE_754
          !finite(mu) || !finite(sigma) ||
  #endif
! 	sigma <= 0.0) {
  	ML_ERROR(ME_DOMAIN);
  	return ML_NAN; 
!     }
!     return mu + sigma * snorm();
  }
--- 35,46 ----
  #ifdef IEEE_754
          !finite(mu) || !finite(sigma) ||
  #endif
! 	sigma < 0.0) {
  	ML_ERROR(ME_DOMAIN);
  	return ML_NAN; 
!     } else
!     if (sigma==0.0)
!        return mu;
!     else 
!        return mu + sigma * snorm();
  }

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._