search for: rhyper

I reported bugs in "rhyper" and the TASKS.OLD file claims that they have been fixed. But R continues to get stuck (the prompt never comes back) if I use rhyper with N1, N2, n large, for example "rhyper(5,200,250,60)". E. S. Venkatraman (venkat at biosta.mskcc.org) =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-...

...I wrote a package, which in the NAMESPACE file exports functions like this: exportPattern("^\\rh") On R-2.8 (Linux, 64), upon loading the package I have the rh functions present. On R-2.10, Mac OS X, (32 bit), it builds, loads, but the functions are not loaded, i.e the only function is rhyper (which is not from my package). Is there something wrong with my package setup? Regards Saptarshi

Hi Developers, The alpha, compiles successfully, but it is failing make check-all (on two seperate machines, both FreeBSD 6.1). Here is the version string: platform i386-unknown-freebsd6.1 arch i386 os freebsd6.1 system i386, freebsd6.1 status alpha major 2 minor 3.0 year 2006 month 03 day 27 svn rev

There is a problem with `rpois'. It does seem to take care about the order of the arguments. This is an example: > rpois(n=1,lambda=2) [1] 3 > rpois(lambda=2,n=1) [1] 2 0 It obviously uses the first argument as the number of samples to be drawn, which is wrong. I used Version 0.49 Beta (April 23, 1997). Fredrik

There is a problem with `rpois'. It does seem to take care about the order of the arguments. This is an example: > rpois(n=1,lambda=2) [1] 3 > rpois(lambda=2,n=1) [1] 2 0 It obviously uses the first argument as the number of samples to be drawn, which is wrong. I used Version 0.49 Beta (April 23, 1997). Fredrik

On 1/20/20 4:26 AM, Martin Maechler wrote: > Coming late here -- after enjoying a proper weekend ;-) -- > I have been agreeing (with Spencer, IIUC) on this for a long > time (~ 3 yrs, or more?), namely that I've come to see it as a > "design bug" that rpois() {and similar} must return return typeof() "integer". > > More strongly, I'm actually pretty

...=20 should mean for non-integer m, n, k. =09-thomas > # OK > dhyper(1,10,90,10) > > # NaN > dhyper(1,10.5,90,10) > > # OK > phyper(1,10,90,10) > phyper(1,10.5,90,10) > > # OK > qhyper(0.71,10,90,10) > qhyper(0.71,10.5,90,10) > > # OK > rhyper(1,10,90,10) > rhyper(1,10.5,90,10) > > > > # this "naive" version of dhyper() works > > dhyperx<-function(x,m,n,k,log=3DFALSE){ > lp<-lchoose(m, x) + lchoose(n, k-x) - lchoose( m+n, k) > if (log) {lp} else {exp(lp)} > } > > # O...

...e with R were 32 - not 64 bit. I was able to install them and so now the update.packages() works. I still have a couple of issues: 1) a lot of packages are not available for R 2.14.0. 2) When compiling I occasionally get compilation errors - for example: "../inst/include/Rcpp/stats/random/rhyper.h", line 31: Warning (Anachronism): Formal argument gen of type double(*)(double,double,double) in call to Rcpp::Vector<14>::Vector<double, double, double>(const int&, double(*)(double,double,double), const double&, const double&, const double&) is being passed exte...

...Monitor the history of this process, assuming that one individual is infected at the beginning. The code is: * simulation <- function(number, prob){ cumulative.time <- 0 current.time <- 0 number.sick <- 1 while(number.sick<number){ current.time <- current.time + 1 meetings <- rhyper(nn=1, m=number.sick, n=number-number.sick, k=2) if(meetings==1){ one.sick <- rbinom(n=1, size=1, prob) if(one.sick==1){ cumulative.time <- c(cumulative.time, current.time) number.sick <- number.sick + 1 }}} cumulative.time } number <- 1000 prob <- .05 model <- simulate(number, p...

On 1/20/20 12:33 PM, Martin Maechler wrote: > > It's really something that should be discussed (possibly not > here, .. but then I've started it here ...). > > The NEWS for R 3.0.0 contain (in NEW FEATURES) : > > * Functions rbinom(), rgeom(), rhyper(), rpois(), rnbinom(), > rsignrank() and rwilcox() now return integer (not double) > vectors. This halves the storage requirements for large > simulations. > > and what I've been suggesting is to revert this change > (svn rev r60225-6) which was purposefu...

...t; > yes, of course there are several such distributions. > > It's really something that should be discussed (possibly not > here, .. but then I've started it here ...). > > The NEWS for R 3.0.0 contain (in NEW FEATURES) : > > * Functions rbinom(), rgeom(), rhyper(), rpois(), rnbinom(), > rsignrank() and rwilcox() now return integer (not double) > vectors. This halves the storage requirements for large > simulations. > > and what I've been suggesting is to revert this change > (svn rev r60225-6) which was purposefull...

....1666667 0.6666667 0.1666667 > dhyper(0:2,2,2,2) *6 [1] 1 4 1 so that is exactly what you would expect. However, > dhyper(10782,209410, 276167, 25000) [1] 0.005230889 so you wouldn't expect 10782 to recur. It is close to the mean of the hypergeometric distribution, but > sim <- rhyper(1e6,209410, 276167, 25000) > mean(sim) [1] 10781.53 > sd(sim) [1] 76.14209 (and incidentally, rhyper is plenty fast enough that you don't really need r2dtable for the 2x2 case) -pd > > Cheers, Jari O. > ________________________________________ > From: R-devel <r-devel-b...

...>>> It's really something that should be discussed (possibly not >>> here, .. but then I've started it here ...). >>> >>> The NEWS for R 3.0.0 contain (in NEW FEATURES) : >>> >>> * Functions rbinom(), rgeom(), rhyper(), rpois(), rnbinom(), >>> rsignrank() and rwilcox() now return integer (not double) >>> vectors. This halves the storage requirements for large >>> simulations. >>> >>> and what I've been suggesting is to revert this change...

...cf8475861b3690bc8bf1cefaa8f R-3/R-3.6.3.tar.gz This is the relevant part of the NEWS file CHANGES IN R 3.6.3: NEW FEATURES: * The included LAPACK has been updated to version 3.9.0 (for the included routines, just bug fixes). BUG FIXES: * Fixed a C level integer overflow in rhyper(); reported by Benjamin Tyner in PR#17694. * Uses of url(gzcon(.)) needing to extend buffer size have failed (with HTTP/2 servers), reported by G'abor Cs'ardi. * predict(loess(..), se=TRUE) now errors out (instead of seg.faulting etc) for large sample sizes, than...

...cf8475861b3690bc8bf1cefaa8f R-3/R-3.6.3.tar.gz This is the relevant part of the NEWS file CHANGES IN R 3.6.3: NEW FEATURES: * The included LAPACK has been updated to version 3.9.0 (for the included routines, just bug fixes). BUG FIXES: * Fixed a C level integer overflow in rhyper(); reported by Benjamin Tyner in PR#17694. * Uses of url(gzcon(.)) needing to extend buffer size have failed (with HTTP/2 servers), reported by G'abor Cs'ardi. * predict(loess(..), se=TRUE) now errors out (instead of seg.faulting etc) for large sample sizes, than...

...cf8475861b3690bc8bf1cefaa8f R-3/R-3.6.3.tar.gz This is the relevant part of the NEWS file CHANGES IN R 3.6.3: NEW FEATURES: * The included LAPACK has been updated to version 3.9.0 (for the included routines, just bug fixes). BUG FIXES: * Fixed a C level integer overflow in rhyper(); reported by Benjamin Tyner in PR#17694. * Uses of url(gzcon(.)) needing to extend buffer size have failed (with HTTP/2 servers), reported by G'abor Cs'ardi. * predict(loess(..), se=TRUE) now errors out (instead of seg.faulting etc) for large sample sizes, than...

...lnorm.o qlnorm.o rlnorm.o df.o pf.o qf.o rf.o dnf.o dt.o pt.o qt.o rt.o dnt.o dchisq.o pchisq.o qchisq.o rchisq.o rnchisq.o dbinom.o pbinom.o qbinom.o rbinom.o rmultinom.o dcauchy.o pcauchy.o qcauchy.o rcauchy.o dexp.o pexp.o qexp.o rexp.o dgeom.o pgeom.o qgeom.o rgeom.o dhyper.o phyper.o qhyper.o rhyper.o dnbinom.o pnbinom.o qnbinom.o rnbinom.o dpois.o ppois.o qpois.o rpois.o dweibull.o pweibull.o qweibull.o rweibull.o dlogis.o plogis.o qlogis.o rlogis.o dnchisq.o pnchisq.o qnchisq.o dnbeta.o pnbeta.o qnbeta.o pnf.o pnt.o qnf.o qnt.o ptukey.o qtukey.o toms708.o wilcox.o signrank.o ranlib libnma...

It is not about "really arge total number of observations", but: set.seed(4711);tabs <- r2dtable(1e6, c(2, 2), c(2, 2)); A11 <- vapply(tabs, function(x) x[1, 1], numeric(1));table(A11) A11 0 1 2 166483 666853 166664 There are three possible matrices, and these come out in proportions 1:4:1, the one with all cells filled with ones being most common. Cheers, Jari

...> In rgeom(10L, 1e-10) : NAs produced yes, of course there are several such distributions. It's really something that should be discussed (possibly not here, .. but then I've started it here ...). The NEWS for R 3.0.0 contain (in NEW FEATURES) : * Functions rbinom(), rgeom(), rhyper(), rpois(), rnbinom(), rsignrank() and rwilcox() now return integer (not double) vectors. This halves the storage requirements for large simulations. and what I've been suggesting is to revert this change (svn rev r60225-6) which was purposefully and diligently done by a fel...

...f many...] BUG FIXES o Remove "Alpha" revision from opening blurb. Didn't make sense with our development model o set SHELL=/bin/sh in Makefile and use $(MAKE). Avoid use of wildcards in $(MANSRC) for non-gnu "make" (notably Solaris). o rhyper gave wrong results in degenerate cases o Using get with pos=length(search()) or envir=NULL was actually actually starting the search in .GlobalEnv. o first=0 (or negative) argument to substr() caused trouble. o 1%o%1 would segfault o .Dyn.libs would linger in .RD...