similar to: Problem with pbirthday (PR#7837)

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)

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 > 0) 1 - exp(sum(log1p(-(0:(n-1))/c))) else 0. Simplifying this a bit further one

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

Please note that I am not on the mailing list. Please respond to Ken.Lynch@luxfer.net. I did not want to join the list just to ask one question. I need to know what UID's are being given to a domain user when they logon. Is there any way to view the winbindd_imap.tdb or winbindd_cache.tdb. Is there any utilities that would allow me to see the mappings. Thanks for any ideas, Ken Lynch

Hi, I am doing the following: 1. I have a list of files.. Files1=list.files("some directory",pattern="some pattern") 2. I define a list as res=vector("list", length(files1)) 3. I read all the files into this list: res=lapply(files1, read.csv) I now want to rowbind all the items in the list into one big mass (all files have same number of columns). I tried

Suppose I have: Book Value A 10 B 11 C 9 D 8 A 12 C 4 D 5 B 7 I want to summarize above not by Book but by groupings of Books as in (below) I have a list ... basic_map <- list(c("A",B"),c("C,D")) Big_names <- c("A1", "A2") Names(basic_map) <- big_names So I want to get : A1 40 A2 26 How do I use tapply AND the list to get my

As requested by Mark Williamson in BZ #681 Signed-off-by: Gawain Lynch <gawain.lynch@gmail.com> --- xendomains.orig 2006-06-17 14:35:48.000000000 +1000 +++ xendomains 2006-06-17 14:45:29.000000000 +1000 @@ -207,10 +207,11 @@ if [ "$XENDOMAINS_RESTORE" = "true" ] && contains_something "$XENDOMAINS_SAVE" then + XENDOMAINS_SAVED=`/bin/ls

Hi all: I have a directory of files as in -------------> bunchafiles <- list.files(path="/data/2.3/2006", pattern="returns", full.names=T,recursive=T) .... Each file is a bunch of returns for a particular date (unique). There are like 252 files or so. With a custom function myread (below), I define a vector ----- res <-

I have a dataframe as follows: Date time value 20110620 11:18:00 7 20110620 11:39:00 9 20110621 11:41:00 8 20110621 11:40:00 6 20110622 14:05:00 8 20110622 14:06:00 6 For every date, I want to extract the row that has the greatest time. Therefore, ending up like: 20110620 11:39:00 9 20110621 11:41:00 8 20110622 14:07:00 6 I am using for loops (for every date, find largest time value) to do

Re: Heteroscedasticity and mgcv. (Collin Lynch) The GAMLSS package can model heterogeneity in the scale parameter (e.g. standard deviastion) [and also heterogeity in skewness and kurtosis parameters].of the response variable distribution. For parametric models a generalized likelihood ratio test can be used to test whether the heterogeity is needed. Alternatively a generalized Akaike

I am using Version 1.3.1 (2001-08-31) on a PC running Win98 and have 384MB RAM The error that I get is Reached total allocation of 256Mb: see help(memory.size) According to the help I should issue a Command-line flag `--max-mem-size' Where do I do this and what is the statement ? Thanks in advance Steve Lynch School of Bioresources Engineering and Environmental Hydrology Univ of Natal

Hi, I'm looking for a command to extract lines from a large text file, a password file. A typical user has a username made from a letter followed by their id-number. m9718508:pw:9301:840: Lynch :/home/pgstud/m9718508:/bin/bash So for instance if I need to extract lines where; the 1st field, the username begins with an m and the 4th field, the group contains exactly 850 cat passwdfile |

I've rolled up R-2.1.1.tar.gz a short while ago. This is a maintenance release containing mainly bugfixes. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.1.1.tar.gz (give it some time to arrive there) or wait for it to be mirrored at a CRAN site nearer to you. If you're *really* impatient, www.biostat.ku.dk/~pd/R-release should work

I've rolled up R-2.1.1.tar.gz a short while ago. This is a maintenance release containing mainly bugfixes. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.1.1.tar.gz (give it some time to arrive there) or wait for it to be mirrored at a CRAN site nearer to you. If you're *really* impatient, www.biostat.ku.dk/~pd/R-release should work

I just upgraded to kernel 2.4.14 and applied the ext3 patch ext3-2.4-0.9.15-2414.gz. It appears (to me) that the patch breaks the loop block device module (loop.o). The problem is that the call to deactivate_page is unresolved. The routine used to be defined in swap.c. It looks like it no longer exists in 2.4.14. I'm not a kernel programmer so if I'm mistaken, sorry. I do get an

I am trying to install phpBB forum on version 5.2. It requires that I have imagemagick and gd graphics support. What rpm packages provide these. I have tried a yum install ImageMagicK with yum responding nothing todo. Is imageMagicK not available with 5.2? Any help would be greatly appreciated. -- Bo Lynch

My apologies for the statistical naivete of my question but... Is there an established method or calulating the statistical power of a correlation test? And if so is there a method in R for it? Thank you, Collin Lynch.

This is probably a really bad question to ask but here goes. Does asterisk work with any of the T1/E1 cards from SBE? I'm sure SBE is a competitor to Digium, but I have access to SBE cards and the linux driver. Just curious more than anything. Thanks. -- Michael J. Lynch What if the hokey pokey IS what it's all about -- author unknown

Just recently bought a Dell E6410 with a Intel? HD Graphics with ExpressCard with I7-640m processor. Having some issues getting a graphical env to work. By just the looks of lspci X is wanting to use the VGA Compatable controller: Intel Corp Core Processor Intergrated graphics controller. Everything else is using Intel Corp Series 5/3400 Series Chipset. Has anyone else had any luck getting

Hello Everyone, I have what I think is a complex text parsing task. I've provided some sample data below. There's a relatively simple version of the coding that needs to be done and a more complex version. If someone could help me out with either version, I'd greatly appreciate it. Here are my sample data. haveData <- structure(list(profile_key = structure(c(1L, 1L, 2L, 2L, 2L,