rnorm returning NA's (PR#8635)

classic Classic list List threaded Threaded
3 messages Options
Reply | Threaded
Open this post in threaded view
|

rnorm returning NA's (PR#8635)

jjmcnutt
Full_Name: Josh McNutt
Version: 2.2.1
OS: Win XP
Submission from: (NULL) (192.88.209.232)


> which(is.na(rnorm(20000000)))
[1] 15242377

> which(is.na(rnorm(10000000)))
[1] 3692029

> which(is.na(rnorm(40000000)))
[1]  5560337  5938719 33888822

> which(is.na(rnorm(50000000)))
[1] 25231754 42397181 45085564 45363557

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
Reply | Threaded
Open this post in threaded view
|

Re: rnorm returning NA's (PR#8635)

Brian Ripley
This is not reproducible, and you have not told us the seed you used which
is what is needed to make random results reproducible.

Please supply a reproducible example (if you can).  I ran the first line
100 times (about 15 minutes) without encountering any NAs.  If you found
these fairly easily it indicates a bug in your Windows installation, not
one in R.  Further, if it were at all common then other people would have
encountered it and surely mentioned it.

On Thu, 23 Feb 2006, [hidden email] wrote:

> Full_Name: Josh McNutt
> Version: 2.2.1
> OS: Win XP
> Submission from: (NULL) (192.88.209.232)
>
>
>> which(is.na(rnorm(20000000)))
> [1] 15242377
>
>> which(is.na(rnorm(10000000)))
> [1] 3692029
>
>> which(is.na(rnorm(40000000)))
> [1]  5560337  5938719 33888822
>
>> which(is.na(rnorm(50000000)))
> [1] 25231754 42397181 45085564 45363557
>
> ______________________________________________
> [hidden email] mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>
>

--
Brian D. Ripley,                  [hidden email]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
Reply | Threaded
Open this post in threaded view
|

Re: rnorm returning NA's (PR#8635)

Duncan Murdoch
In reply to this post by jjmcnutt
On 2/23/2006 1:04 PM, [hidden email] wrote:

> Full_Name: Josh McNutt
> Version: 2.2.1
> OS: Win XP
> Submission from: (NULL) (192.88.209.232)
>
>
>> which(is.na(rnorm(20000000)))
> [1] 15242377
>
>> which(is.na(rnorm(10000000)))
> [1] 3692029
>
>> which(is.na(rnorm(40000000)))
> [1]  5560337  5938719 33888822
>
>> which(is.na(rnorm(50000000)))
> [1] 25231754 42397181 45085564 45363557
>
> ______________________________________________
> [hidden email] mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel

I don't see these.  Could you use set.seed() to put the RNG in a
reproducible state, and show us the results then?

Duncan Murdoch

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel