similar to: hybrid in fisher.test broken?

Hi! I am Marta Colombo, student in Mathematical Engineering at "Politecnico di Milano". For my master degree thesis I have to analyze some categorical data. My dataset is composed by 327 individuals and 16 variables. I am using Fisher exact test to test independence on IxJ contingency tables, but I have a problem with one variable. R gives me this error message: FEXACT error 7.

Dear all, I’m having an awkward problem in R. When I write the command Fisher.test(school.data,workspace=2e+07), where school.data is the matrix corresponding to the data set, I get the error message: FEXACT error 7. LDSTP is too small for this problem. Try increasing the size of the workspace. Increasing the workspace: Fisher.test(school.data,workspace=1e+10), I

(CC'ed to R-bugs ``for the record'') >>>>> "BDR" == Prof Brian D Ripley <ripley@stats.ox.ac.uk> writes: BDR> On Thu, 13 Jun 2002, Martin Maechler wrote: >> >>>>> "MM" == Martin Maechler >> <maechler@stat.math.ethz.ch> writes: >> >> >>>>> "BDR" ==

This is a bad bug as reported by Robin Hankin, it is still in "R-patched" ... ##- From: Robin Hankin <r.hankin@auckland.ac.nz> ##- To: r-help@stat.math.ethz.ch ##- Subject: [R] possum sleeping: thanks and fisher.test() FEXACT error ##- Date: Thu, 13 Jun 2002 16:46:26 +1200 ## ..... ## Example slighlty modified (MM) d4 <- matrix(c(0, 0, 0, 0, 0, 0, 3, 0, 1, 0, 0, 0, 0, 0,

Hello, I have an error message that doesn't seem to make sense in that having read the R documentation I was under the impression that R was able to grab as much memory as it needed, and has been able to do so for some time, so the advice given below about increasing the size of the workspace is redundant. If I'm right, does anyone have a solution to the problem of the size of

Sorry, but I lost this thread, so I sending this as a new message. This is really a follow-up to a post from a couple days ago saying that fisher.test from the ctest library crashed on the following data set: > T [,1] [,2] [1,] 2 1 [2,] 2 1 [3,] 4 0 [4,] 8 0 [5,] 6 0 [6,] 0 0 [7,] 1 0 [8,] 1 1 [9,] 7 1 [10,] 8 2 [11,]

Hi folks, Forgive me if this question is a trivial issue. I was doing a series of Fishers' exact test using the fisher.test function in stats package. Since the counts I have were quite large (c(64, 3070, 2868, 4961135)), R suggested me to use *other algorithms* for the test which can be specified through the 'control' argument of the fisher.test function as I understood. But where

I posted this last night but I think I figured out the problem. I checked on the underlying equation for the Fisher Exact Test for tables greater than 2X2. It looks like I have an insane amount of factorials to be multiplied. Is that problem? Is it an overflow issue? Jeff Last night’s post: I just a ran a fisher.test on a 9x5 table and received the following message

I want to ascertain the basis of the table ranking, i.e. the meaning of "extreme", in Fisher's Exact Test as implemented in 'fisher.test', when applied to RxC tables which are larger than 2x2. One can summarise a strategy for the test as 1) For each table compatible with the margins of the observed table, compute the probability of this table conditional on the

Hi, I want to analyse a contigency table (3 x 12) with a fisher.test beacause there are cells that are less than 5. ?mmen Anken Baf Belchen H?chi Hof Porti R?m Schmutz Sch?n Sissa Tann class14 7 26 150 2 46 68 126 66 3 31 7 61 class24 7 6 55 5 49 71 93 90 1 18 16 79 class34 1 1 4 3 19 8 29 61

Full_Name: Anja von Heydebreck Version: 1.9.1 OS: Linux Submission from: (NULL) (155.250.128.25) I obtained a negative p-value from the fisher.test() function: > fisher.test(matrix(c(14576,3023,89,68),2), alternative="g")$p.value [1] -8.426593e-13 With R 1.8.1, I got a p-value of -6.239231e-12 for this example. Anja von Heydebreck

I am running R 2.2.1 on a University-supported linux installation based on Redhat EL3. I am sorry that it did not occur to me to mention this before; I updated R very recently, with the most recent version available for EL3 at http://cran.cnr.berkeley.edu/bin/linux/redhat/el3/. Looking at the gridBase documentation, I find that 2.2.1 is not in fact the most recent version. I have now spent

I have two matrices, generated by R functions that I don't understand. I want to confirm that they're the same, but I know that they have different attributes. If I want to compare the dimnames, I can say > identical(attr(tm, "dimnames"), attr(tmm, "dimnames")) [1] FALSE or even: > identical(dimnames(tm), dimnames(tmm)) [1] FALSE But I can't find any

Hi folks, I'm trying to call an R function (fisher.test) in my program for like a billion times! Though my program is in Python and I feel that using rpy2 to interface R to python doesn't give me satisfactory performance. So I looked into R code and found out that fisher.test is actually a wrapper around another function called fexact which is implemented in C. Using Cython I managed to

Hi Ian/Stefano: I changed over to the PV clock for hybrid liked we talked at the hackathon. I still have the hang in update_wall_time() after dom0 switches to xen as clocksource. The source of hang seems to be in xen stime_local_stamp in cpu_time that suddenly jumps to a large 64bit value. I''ve been chasing to figure where that happens, and why for the hybrid and not PV. It appears the

Hi, I have hybrid smp running with autoxlate. However, without autoxlate, I am running into issues realted to TLB flush. The guest in this case makes multicalls as part of which cache is flushed (__do_update_va_mapping, etc..). However, the guest is using VPIDs and it is getting complicated. I can just xen not do any TLB management and let the guest just do it after return from the hypercall.

Under Details, the documentation says "[if Rowv or Colv] is 'NULL', _no reordering_ will be done for the corresponding side." In fact, as explained elsewhere in the documentation, 'NA' is required, not 'NULL'. An anonymous expert suggested that I should attach an svn diff file to this report, so I will. JD -------------- next part -------------- Index:

The last line gets duplicated when a file is read like this: a <- scan(file=filename, what="", sep="\n", strip.white=c(TRUE), quiet=TRUE) (This error does not occur on Linux, the only other platform I tested.) Version: platform = hppa2.0-hp-hpux10.20 arch = hppa2.0 os = hpux10.20 system = hppa2.0, hpux10.20 Actually, all binaries are

Hi guys, At a loss trying to figure why if (has_eflag(X86_EFLAGS_ID)) returns false in my HVM domU. Standard function has_eflag() in cpucheck.c running in real mode. Works fine on PV dom0, but fails when guest is booting on my hybrid dom0. LMK if any ideas. I''ll keep digging in the manuals, but nothing so far. thanks, Mukesh

Dear helplist Many many thanks to everyone who helped me. The trick was to use tabulate() or, better, tab <- rep(0,50) names(tab) <- 1:50 tab[names(table(sleeps))] <- table(c) My original dataset was a list of 50 trees and a length 12 vector recording which tree a certain possum slept in on 12 nights. As Professor Ripley points out, a Monte-Carlo simulation is easy to set up, and it