Re: [Rd] proposed pbirthday fix

2006-01-23 Thread Martin Maechler
> "MM" == Martin Maechler <[EMAIL PROTECTED]> > on Mon, 23 Jan 2006 11:52:55 +0100 writes: > "ken" == ken knoblauch <[EMAIL PROTECTED]> > on Mon, 23 Jan 2006 09:43:28 +0100 writes: ken> Actually, since NaN's are also detected in na.action ken> operations, a simpler

Re: [Rd] proposed pbirthday fix

2006-01-23 Thread Martin Maechler
> "ken" == ken knoblauch <[EMAIL PROTECTED]> > on Mon, 23 Jan 2006 09:43:28 +0100 writes: ken> Actually, since NaN's are also detected in na.action ken> operations, a simpler fix might just be to use the ken> na.rm = TRUE option of min ken> upper <- min(n^k/(c^(k - 1))

[Rd] proposed pbirthday fix

2006-01-23 Thread ken knoblauch
Actually, since NaN's are also detected in na.action operations, a simpler fix might just be to use the na.rm = TRUE option of min upper <- min(n^k/(c^(k - 1)), 1, na.rm = TRUE) > Recent news articles concerning an article from The Lancet with > fabricated > data indicate > that in the sample

[Rd] proposed pbirthday fix

2006-01-22 Thread Ken Knoblauch
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 = 1