Re: [Rd] Are r2dtable and C_r2dtable behaving correctly?

2017-08-25 Thread Martin Maechler
> Gustavo Fernandez Bayon > on Thu, 24 Aug 2017 16:42:36 +0200 writes: > Hello, > While doing some enrichment tests using chisq.test() with simulated > p-values, I noticed some strange behaviour. The computed p-value was > extremely small, so I decided to dig a little

Re: [Rd] Are r2dtable and C_r2dtable behaving correctly?

2017-08-25 Thread Jari Oksanen
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 14 There are three possible matrices, and these come out in propo

Re: [Rd] Are r2dtable and C_r2dtable behaving correctly?

2017-08-25 Thread Peter Dalgaard
> On 25 Aug 2017, at 10:30 , Martin Maechler wrote: > [...] > https://stackoverflow.com/questions/37309276/r-r2dtable-contingency-tables-are-too-concentrated > > >> set.seed(1); system.time(tabs <- r2dtable(1e6, c(100, 100), c(100, 100))); >> A11 <- vapply(tabs, function(x) x[1, 1], numeric(1

Re: [Rd] Are r2dtable and C_r2dtable behaving correctly?

2017-08-25 Thread Peter Dalgaard
> On 25 Aug 2017, at 11:23 , Jari Oksanen wrote: > > 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 16

Re: [Rd] Are r2dtable and C_r2dtable behaving correctly?

2017-08-25 Thread Martin Maechler
> Peter Dalgaard > on Fri, 25 Aug 2017 11:43:40 +0200 writes: >> On 25 Aug 2017, at 10:30 , Martin Maechler wrote: >> > [...] >> https://stackoverflow.com/questions/37309276/r-r2dtable-contingency-tables-are-too-concentrated >> >> >>> set.seed(1); sys

Re: [Rd] Are r2dtable and C_r2dtable behaving correctly?

2017-08-25 Thread peter dalgaard
> On 25 Aug 2017, at 12:04 , Peter Dalgaard wrote: > >> 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. > > ... and > >> dhyper(0:2,2,2,2) > [1] 0.167 0.667 0.167 >> dhyper(0:2,2,2,2) *6 >

[Rd] A problem with headless R on WIndows

2017-08-25 Thread Eric Hopper
I need a version of R that has no dependency on GDI32.dll. Ideally, this would be a version of R that had no dependency on rgraphapp.dll. Is there a way to configure the build system to build a version of R like this? I've looked through the ./configure --help for seemingly appropriate flags, and