similar to: Bug in wilcox.test

Full_Name: Charles J. Geyer Version: 1.3.1 OS: SuSE 7.2 Submission from: (NULL) (134.84.86.22) in line 90 of signrank.c and lines 114 and 273 of wilcox.c calloc is used and the return is not checked to see if it is null (violating numbers 2 and 6 of the ten commandments for C programmers). The following code thus makes R segfault (for suitably large values of 1000) for (n in seq(100, 1000,

The note to 1004 says "fixed for 1.3.1" Uh. No. It ain't. The problem was more serious than guessed as even the simplest testing would show. For example, Example 5.4 in Hollander and Wolfe (Nonparametric Statistical, Methods, 2nd ed., Wiley, 1999, pp. 180-181) R Version 1.3.1 (SuSE Linux 7.1) > X <-

The point estimates produced by wilcox.test are perverse (not wrong, just brain damaged). The Hodges-Lehmann estimator that goes with the signed rank test is the median of the Walsh averages. The Hodges-Lehmann estimator that goes with the rank sum test is the median of the pairwise differences. wilcox.test agrees except that it uses the following very peculiar definition of "sample

I took some time and liberty and tried to improve existing implementation of SignRank functions in R. (dsignrank, ...) As I have seen they've been based on csignrank. So I modified csignrank and, I believe, improved calculation time and memory efficiency. The idea is basically the same. I use the same recursion as original author used with one slight modification. I am generating Wilcoxon

Hi Everyone, Following the previous discussion on optimizing *wilcox functions, Andreas Loeffler brought to my attention a few other bugs in `wilcox` family functions. It seems like these issues have been discussed online in the past few months, but I haven?t seen discussion on R-devel...unless I missed an email, it seems like discussion never made it to the mailing list. I haven?t seen any bug

I am puzzled at the use of regression. I have a categorical variable ClassePop33000 which factors a Population variable into 3 levels. I want to investigate whether that categorical variable has some relation with my dependent variable, so I go : lm(Cout.ton ~ ClassePop33000, data=ech2) Call: lm(formula = Cout.ton ~ ClassePop33000, data = ech2) Residuals: Min 1Q Median 3Q

The R_registerRoutines C function allows the number of parameters to a .Call function to be registered. For example, the tools package function md5sum() calls "Rmd5", which has been registered to require just one parameter. But if it is called with the wrong number of parameters, only the first error gets caught: > library(tools) >

This is my first post to the mailing list and I guess it's a pretty stupid question but I can't figure it out. I hope this is the right forum for these kind of questions. Before I started using R I was using STATA to run a Wilcoxon signed-rank test on two variables. See data below:

Hi, I noticed this problem on my home desktop running FC4 and again on my laptop running FC5. Both have previously compiled and passed make check-all on 2.3.1 svn revisions from 10 days ago or so. On both these machines, make check-all is consistently failing (4 out of 4 attempts on the FC 4 desktop and 3 out of 3 on the FC 5 laptop) in the p-r-random-tests tests. This is with both default

Full_Name: Patrick Hodgson Version: 2.0 OS: solaris 2.9 Submission from: (NULL) (65.94.128.161) The value reported for the parameter W in the function wilcox.exact appears to be incorrect. I have checked the reference in the help file for this function (Myles & Hollander 1973, as well as 2nd ed. 1999 by same authors) and it is clear that W is the sum of the ranks of the data set with the

Full_Name: Charles J. Geyer Version: 1.8.0 OS: i686-pc-linux-gnu (Suse 8.2) Submission from: (NULL) (134.84.86.22) Two bugs (perhaps related, perhaps independent) revealed by the same Poisson regression with offset mydata <- read.table(url("http://www.stat.umn.edu/geyer/5931/mle/seeds.txt")) out.fubar <- glm(seedlings ~ burn01 + vegtype * burn02 + offset(log(totalseeds)),

Dear R friends, here I write again about the wilcox.exact() problem. I want to compare two sets of categorical data, and in one case it says "negative length vectors not allowed", and in the other one I get the error "cannot allocate vector of length ...". On http://rapidshare.com/files/298621893/wilcox.exact_trouble.Rdata.html you can download the data that cause the

Here's a simpler example that should reproduce that error for you: ans <- utils::select.list(c("hello", "world", "again"), graphics=TRUE) Does it? FYI, I installed R 3.6.3 from source on CentOS 7 a few hours ago, and for me the above works just fine. For your immediate needs of selecting a CRAN mirror, you can set: options(menu.graphics = FALSE) as a

On a recent FreeBSD 8.0-CURRENT (i386) building R (any version) breaks with the following messages: ---------------------------------------------------------------------- [...snip...] gcc -std=gnu99 -I. -I../../src/include -I../../src/include -I/usr/local/include -DHAVE_CONFIG_H -g -O2 -c wilcox.c -o wilcox.o gcc -std=gnu99 -I. -I../../src/include -I../../src/include -I/usr/local/include

Dear Developers Team, I have two items: 1. wilcox.test with the paired=T option appears to delete zeros before ranking absolute differences. Would it be possible to add the feature of removing zeros after ranking, which is given in Lehmann's Nonparametrics as the preferred choice. See also Pratt (1959), JASA 54, 655-667. It is given in wilcoxsign_test of the coin package as an option

OK. I give up. I'll ask a stupid question. How do I get the $!#@*$ signif stars line printed by summaries to not look extremely bizarre in the latex produced by Sweave? For example, see p. 7 of http://www.stat.umn.edu/geyer/aster/library/aster/doc/tutor.pdf I can see what the problem is. R emits non-ascii characters (as it is supposed to do), Sweave puts them in the tex file, and

The subject line says it all. How can I find what Error: cannot set length of non-vector means? RTFS is no help. I can find out of course that it comes from "lengthgets", but who called that? Not me! The situation is as follows. I am trying to get a package ready for CRAN. Every time I run "R CMD check" I get this error when building the package vignette. But (!)

Can anyone tell me if the following C code (which proved very useful when using the snow package -- use it to nice slaves) compiles and dyn.loads under Windoze or Mac? It is (apparently) POSIX, so I suppose it is fine in OS X, but does Windoze's advertized POSIX compliance mean anything here? If it doesn't work on some platform, how do I deal with that? I want to submit as contributed

I knew I could work around. But this shouldn't happen. And yes. Same problem with your example. blurfle$ R --vanilla R version 3.6.3 (2020-02-29) -- "Holding the Windsock" Copyright (C) 2020 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain

This is just to alert the list about a new technical report http://www.stat.umn.edu/geyer/gdor/ that shows how to determine, using exact infinite-precision rational arithmetic when the MLE exists in a GLM (or for that matter in any exponential family, such as aster models -- contributed package aster) and what to do when it does not. This suggests many changes to glm, glm.fit, predict.glm,