similar to: documentation bug in (ctest) chisq.test (PR#2346)

Dear all, I have a question about the chisq.test command. As an option one can chose the computation of p-values by Monte-Carlo simulation (simulate.p.value=T). Is there any documentation available how this calculations are done and how this simulation based test behaves in small samples? Thanks Klaus Abberger University of Konstanz, Germany [[alternate HTML version deleted]]

When R prints warnings, they often go "out of the line", it would be better if they where wrapped with writeLines(strwrap ... I tried to do that , changing the function warnings, but it has only effect when called explicitely, not when R prints the warnings unasked. Anyhow, here is the changed warnings: warnings <- function (...) { if (!(n <- length(last.warning)))

(I tried to send this earlier, but it doesnt seem to have come through, due to problems on my system) Hola: Both cannot be correct: > m <- matrix(1:4, 2) > svd(m) $d [1] 5.4649857 0.3659662 $u [,1] [,2] [1,] -0.5760484 -0.8174156 [2,] -0.8174156 0.5760484 $v [,1] [,2] [1,] -0.4045536 0.9145143 [2,] -0.9145143 -0.4045536 > La.svd(m) $d [1]

....in example data <- read.spss("c:/data.sav",use.value.label=F,to.data.frame=T) try help(read.spss) ! regards,christian kjetilh at umsanet.edu.bo schrieb am 28.08.02 02:08:39: > Hola! > > Can the foreign package now read data from spss 10? > > Kjetil Halvorsen > -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- > r-help

The following code in clusplot.default (package cluster) is in error: x1 <- cmdscale(x, k = 2, eig = TRUE) var.dec <- sum(x1$eig)/sum(diag(x1$x)) if (var.dec < 0) var.dec <- 0 if (var.dec > 1) var.dec <- 1 x1 <- x1$points x1 has components with names "points" and "eig", not "x", so

Hola! look at the following: > test <- function() { + options(warn=-1) + on.exit( options(warn=0)) + log(-1) + } > test() [1] NaN but what should be better, doesnt work: > test1 <- function() { + oldwarn <- options("warn") + options(warn=-1) + on.exit( options(warn=oldwarn)) + log(-1) + } > test1() Error in options(...) : warn parameter

(I sent this earlier, but it seems not to have come through, due to problems witkh my system) The following command from appendix A, "a sample session", isnt correct: contour(x, y, fa, nint=15) when used R protests: Warning message: parameter "nint" couldn't be set in high-level plot() function it should probably be nlevels, as used a few lines before. This is

I am resubmitting this to r-bugs, since Thomas Lumley indicates that it might be an error: On Wed, 5 Jul 2000, Thomas Lumley wrote: > On Wed, 5 Jul 2000, halvorsen wrote: > > > Hola! > > > > I have the following simple function: > > > > > testcar > > function(pow){ > > ob <- glm(Pound~CG+Age+Vage,data=car,weights=No, > >

Hola! > rm(list=ls(all=TRUE)) > x <- 1:20 > y <- 1+x+rnorm(x) > xy.coords(y ~ x,NULL) ... expected output, correct, but when called from inside lowess: > lowess(y ~ x) Error in xy.coords(x, y) : x and y lengths differ > debug(xy.coords) > lowess(y ~ x) debugging in: xy.coords(x, y) ... long listing deleted if (is.language(x)) { if (inherits(x,

I reported this earlier, and have got reports that others have the same promlem on UNIX machines, so it is not only a windows problem S I file a bug report. The function (boot.stat) given at the end produces a postscript file, which cannot be included correctly in LaTeX. Specifically, the image in LaTeX (when translated by dvips to postscript) becomes very small, not using the bounding box, and

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

Dear useRs I am looking for a way to randomise the values within a matrix: the conditions are that the sums of the rows and the sums of the columns should remain the same as in the original matrix. Any help would be appreciated Cheers Yann

This is a message for whoever maintains "chisq.test": For an outcome more extreme than 2000 simulations, a Monte Carlo p-value of "< 2.2e-16" was printed. Ripley said the proper p-value for such cases should be 1/(B+1) = 1/2001. This can be easily fixed by adding "if(PVAL==0)PVAL <- 1/(B+1)" right after the following line in the code for chisq.test (in R

>>>>> Gustavo Fernandez Bayon <gbayon at gmail.com> >>>>> 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 deeper and debug >

> On 25 Aug 2017, at 11:23 , Jari Oksanen <jari.oksanen at oulu.fi> 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 166664 > > There are

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 deeper and debug chisq.test(). I noticed then that the simulated statistics returned by the following call tmp <- .Call(C_chisq_sim, sr, sc, B, E) were all the same, very small numbers. This, at first,

This morning I had 37 messages from the r-help list in my mailbox. I think its becoming excessive for an e-mail list. I wonder if whoever looks after this list could either move or gateway it to a usenet group? That would also eliminate the need for special purpose archiving and searching facilities since the site: http://groups.google.com would automatically provide that. That site also

I am in the process of setting up R1.3.1 on a new computer running windows 2000. I am having problems running the PFE text editor simultaneously within R for editing functions and outside R for editing ordinary text files. The PFE editor behaves as I expected, if it is opened in R AND but no other PFE window is open outside R. Similarly, it also works fine if I am editing a text file outside

> From: kjetil halvorsen [mailto:kjetilh at umsanet.edu.bo] > > Martin Hoyle wrote: > > > > Dear Greg, > > Thanks for the new release. The decomposition of the SSQ is > just what I need! > > > I thought decomposition of the SSQ can be obtained with the split= > argument to summary.aov? That argument seems tro be new --- it didnt > exist two years

Full_Name: Kjetil Kjernsmo Version: 1.1.1.1 (patched the rpois :-)) OS: osf4.0e Submission from: (NULL) (129.240.28.227) Dear all, I have just gotten through a debugging session of my code, being confused for a few days. I sent what I thought was a straightforward matrix to chisq.test() and once in a while got the error message: > chisq.test(t2) Error in chisq.test(t2) : missing value where