search for: pbirthday

Hi, pbirthday(, coincident = 2) starts to issue warnings (see (*) below) for larger number of classes (R 4.0.0, R-devel ./src/library/stats/R/birthday.R:47). The default coincident = 2 is computed as 1 - prod((c:(c - n + 1))/rep(c, n)) where c = classes. Using exp(log(...)), one can derive the return value if(n...

Full_Name: Andy Lynch Version: 1.9.1 OS: Windows Submission from: (NULL) (131.111.86.211) As I understand it, pbirthday(n,c,k) gives the approximate probability that we see a class with k coicident people in it when n people are sorted into c classes. so the command > pbirthday(4,classes=3,coincident=4) should give the approximate probability that when four people fall into three classes, all four end up i...

Recent news articles concerning an article from The Lancet with fabricated data indicate that in the sample containing some 900 or so patients, more than 200 had the same birthday. I was curious and tried out the p and q birthday functions but pbirthday could not handle 250 coincidences with n = 1000. The calculation of upper prior to using uniroot produces NaN, upper<-min(n^k/(c^(k-1)),1) I was able to get it to work by using logs, however, as in the following version function(n, classes = 365, coincident = 2){ k <- coincident c...

Ralf I don't doubt this is expected with the current implementation, I doubt the implementation is desirable. Suggesting to turn this to pbirthday(1e6, classes = 2^53) ## [1] 5.550956e-05 (which is still non-zero, but much less likely to cause confusion.) Best regards Kirill On 26.02.19 10:18, Ralf Stubner wrote: > Kirill, > > I think some level of collision is actually expected! R uses a 32bit MT > that can produce 2^32 diff...

Gabe As mentioned on Twitter, I think the following behavior should be fixed as part of the upcoming changes: R.version.string ## [1] "R Under development (unstable) (2019-02-25 r76160)" .Machine$double.digits ## [1] 53 set.seed(123) RNGkind() ## [1] "Mersenne-Twister" "Inversion"??????? "Rejection" length(table(runif(1e6))) ## [1] 999863 I don't

...ault value of cex.axis was used. (Related to PR#14550.) ? dataframe[[row,col]] now dispatches on [[ methods for the selected column (spotted by Bill Dunlap). ? sort.int() would strip the class of an object, but leave its object bit set. (Reported by Bill Dunlap.) ? pbirthday() and qbirthday() did not implement the algorithm exactly as given in their reference and so were unnecessarily inaccurate. pbirthday() now solves the approximate formula analytically rather than using uniroot() on a discontinuous function. The description of the pro...

...ault value of cex.axis was used. (Related to PR#14550.) ? dataframe[[row,col]] now dispatches on [[ methods for the selected column (spotted by Bill Dunlap). ? sort.int() would strip the class of an object, but leave its object bit set. (Reported by Bill Dunlap.) ? pbirthday() and qbirthday() did not implement the algorithm exactly as given in their reference and so were unnecessarily inaccurate. pbirthday() now solves the approximate formula analytically rather than using uniroot() on a discontinuous function. The description of the pro...

Hi, I tried to calculate the formula for the birthday problem (the probability that at least two people out of a group of n people share the same birthday) But the factorial-function allows me only to calculate factorials up to 170. So is there a way to push that limit? to solve this formula: (factorial(365) / factorial((365-23))) / (365^23) (n=23)

Kirill, I think some level of collision is actually expected! R uses a 32bit MT that can produce 2^32 different doubles. The probability for a collision within a million draws is > pbirthday(1e6, classes = 2^32) [1] 1 Greetings Ralf On 26.02.19 07:06, Kirill M?ller wrote: > Gabe > > > As mentioned on Twitter, I think the following behavior should be fixed > as part of the upcoming changes: > > R.version.string > ## [1] "R Under development (unstable)...

...work correctly on objects fitted from a data frame with row names. o The coding standards recommendations had nuke-trailing-whitespace where newer versions of ESS need ess-nuke-trailing-whitespace. (PR#7888) o package.skeleton() missed the first newline in the DESCRIPTION file. o pbirthday() reported p = 1 too often when coincident > 2. o plot(1:3, exp(1:3), log = "y", ylim = c(30,1)) {reversed log-scale axis} now works, based on Uwe Ligges' suggestions. (PR#7894) o install.packages() was aborting when a package in a bundle was chosen from a menu. It fa...

...work correctly on objects fitted from a data frame with row names. o The coding standards recommendations had nuke-trailing-whitespace where newer versions of ESS need ess-nuke-trailing-whitespace. (PR#7888) o package.skeleton() missed the first newline in the DESCRIPTION file. o pbirthday() reported p = 1 too often when coincident > 2. o plot(1:3, exp(1:3), log = "y", ylim = c(30,1)) {reversed log-scale axis} now works, based on Uwe Ligges' suggestions. (PR#7894) o install.packages() was aborting when a package in a bundle was chosen from a menu. It fa...

...PR#14452.) o gl() runs faster by avoiding calling factor(). o The print() method for object.size() accepts B as well as b as an abbreviation for 'bytes'. o unlink() gains a force argument to work like rm -f and if possible override restrictive permissions. o pbirthday() and qbirthday() now use exact calculations for coincident = 2. o unzip() and unz() connections have been updated with support for more recent Zip64 features (including large file sizes and bzip2 compression, but not UTF-8 file names). unzip() has a new option to rest...

...PR#14452.) o gl() runs faster by avoiding calling factor(). o The print() method for object.size() accepts B as well as b as an abbreviation for 'bytes'. o unlink() gains a force argument to work like rm -f and if possible override restrictive permissions. o pbirthday() and qbirthday() now use exact calculations for coincident = 2. o unzip() and unz() connections have been updated with support for more recent Zip64 features (including large file sizes and bzip2 compression, but not UTF-8 file names). unzip() has a new option to rest...

...Page never being called on the device, which in turn resulted in incorrect output for postcript() and pdf() devices. (Reported by Marc Schwartz in discussion of the non-bug PR#7820.) o terms.formula needed to add parentheses to formulae with terms containing '|'. (PR#8462) o pbirthday() and qbirthday() now also work for very improbable events {those you are typically *not* interested in}. o Only source help files starting with an upper- or lower-case letter or digit and extension .Rd or .rd are documented to be processed. This is more liberal in that starting with a dig...

...Page never being called on the device, which in turn resulted in incorrect output for postcript() and pdf() devices. (Reported by Marc Schwartz in discussion of the non-bug PR#7820.) o terms.formula needed to add parentheses to formulae with terms containing '|'. (PR#8462) o pbirthday() and qbirthday() now also work for very improbable events {those you are typically *not* interested in}. o Only source help files starting with an upper- or lower-case letter or digit and extension .Rd or .rd are documented to be processed. This is more liberal in that starting with a dig...