search for: psignrank

...ers 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, 100)) { mu <- n * (n + 1) / 4 sigma <- sqrt(n * (n + 1) * (2 * n + 1) / 24) z <- floor(qnorm(0.05, mu, sigma)) p <- psignrank(z, n) cat("z = ", z, ", n = ", n, ", p = ", p, "\n", sep="") } for (n in seq(100, 1000, 100)) { m <- n mu <- n * m / 2 sigma <- sqrt(m * n * (m + n + 1) / 12) z <- floor(qnorm(0.05, mu, sigma)) p <- pwilcox(z, m, n)...

There is a fairly new bug in wilcox.test in R-2.2.1 (stable). It wasn't there when I last taught nonparametrics in fall 2003. Line 86 of wilcox.test.R achieved.alpha<-2*psignrank(trunc(qu),n) It should be achieved.alpha<-2*psignrank(trunc(qu)-1,n) If you don't see why, decode the cookbook instructions p. 56 in Hollander and Wolfe (2nd ed.) or see http://www.stat.umn.edu/geyer/5601/examp/signrank.html#conf or just do a sanity check: does this to the right...

...ion formula. What is changed? There is no need for SINGRANK_MAX in src/nmath/nmath.h anymore. Only functions: static void w_free() void signrank_free() static void w_init_maybe(int n) static double csignrank(int k, int n) in src/nmath/signrank.c have been altered. There was no change to dsignrank, psignrank, ... I've tried to make as little changes as possible. So, to compile with new functions only src/nmath/signrank.c needs to be changed with the one given in attachment. I hope it really is an improvement. With respect, -- Ivo Ugrina ICQ: 47508335 | www.iugrina.com --------------------------...

Hi, I am using R to perform wilcox.test and wondering if you know how the p-value in wilcox.test is calculated? Thank you for your help:-) Regards, Latha

...../include -I../../src/include -c qwilcox.c -o qwilcox.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c rwilcox.c -o rwilcox.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c dsignrank.c -o dsignrank.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c psignrank.c -o psignrank.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qsignrank.c -o qsignrank.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c rsignrank.c -o rsignrank.o ar cr libmath.a mlutils.o d1mach.o i1mach.o fmax2.o fmin2.o fmod.o fprec.o fround.o ftrunc.o sign....

...I../../src/include -c qwilcox.c = -o = qwilcox.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c rwilcox.c = -o = rwilcox.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c dsignrank.= c -o = dsignrank.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c psignrank.= c -o = psignrank.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qsignrank.= c -o = qsignrank.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c rsignrank.= c -o = rsignrank.o ar cr libmath.a mlutils.o d1mach.o i1mach.o fmax2.o fmin2.o fmod.o fprec.= o = fround....

...include -c pwilcox.c -o pwilcox.o gcc -g -O2 -I../include -I../../src/include -c qwilcox.c -o qwilcox.o gcc -g -O2 -I../include -I../../src/include -c rwilcox.c -o rwilcox.o gcc -g -O2 -I../include -I../../src/include -c dsignrank.c -o dsignrank.o gcc -g -O2 -I../include -I../../src/include -c psignrank.c -o psignrank.o gcc -g -O2 -I../include -I../../src/include -c qsignrank.c -o qsignrank.o gcc -g -O2 -I../include -I../../src/include -c rsignrank.c -o rsignrank.o ar cr libmath.a mlutils.o d1mach.o i1mach.o fmax2.o fmin2.o fmod.o fprec.o fround.o ftrunc.o sign.o fsign.o fsquare.o fcube.o imax2....

...include -c pwilcox.c -o pwilcox.o gcc -g -O2 -I../include -I../../src/include -c qwilcox.c -o qwilcox.o gcc -g -O2 -I../include -I../../src/include -c rwilcox.c -o rwilcox.o gcc -g -O2 -I../include -I../../src/include -c dsignrank.c -o dsignrank.o gcc -g -O2 -I../include -I../../src/include -c psignrank.c -o psignrank.o gcc -g -O2 -I../include -I../../src/include -c qsignrank.c -o qsignrank.o gcc -g -O2 -I../include -I../../src/include -c rsignrank.c -o rsignrank.o ar cr libmath.a mlutils.o d1mach.o i1mach.o fmax2.o fmin2.o fmod.o fprec.o fround.o ftrunc.o sign.o fsign.o fsquare.o fcube.o imax2....

Problem compiling R. version: 0.64.2 machine: SGI O2 OS: IRIX6.5 CC cc FC f77 MAKE GNU make 3.75 Here is part of the the output for make. make[2]: Leaving directory `/usr/people/faculty/math/mgass/stage/R/R-0.64.2/src/include' make[2]: Entering directory `/usr/people/faculty/math/mgass/stage/R/R-0.64.2/src/appl' cc -g -OPT:IEEE_NaN_inf=ON -I../include