search for: dwilcox

...mproved [dpq]wilcox functions by implementing Harding's algorithm: Harding, E.F. (1984): An Efficient, Minimal-storage Procedure for Calculating the Mann-Whitney U, Generalized U and Similar Distributions, App. Statist., 33, 1-6 Results on my computer show (against R-2.9.1): > system.time( dwilcox( 800, 800, 80) ) user system elapsed 0.240 0.180 0.443 > system.time( dwilcox( (1:2800), 80, 80) ) user system elapsed 0.290 0.020 0.30 > system.time( dwilcox( (1:2800), 160, 160) ) user system elapsed 0.810 0.060 0.868 > system.time( dwilcox( (1:28000), 210,...

Full_Name: Mark J. Lamias Version: 1.7.0 OS: Windows 2000 Pro Submission from: (NULL) (65.222.84.72) I am running the qwilcox procedure and it is producing incorrect results. For example, dwilcox(.025, 3, 5) should equal 6, but it is equal to 1. Similarly, dwilcox(.025, 3, 6) should equal 7, but it equals 2. The critical values are not set being returned with the correct values. I've verified this with a program that performs direct enumeration to determine the appropriate critical v...

...39;t work for R-1.8.0 (bin wang) > 2. Re: var() doesn't work for R-1.8.0 (Uwe Ligges) > 3. strptime (PR#4204) (pliu3@ncsu.edu) > 4. possible bug in diag() (Marsland, John) > 5. Re: possible bug in diag() (Simon Urbanek) > 6. Re: R on BeOS (Simon Urbanek) > 7. dwilcox (PR#4212) (mark.lamias@grizzard.com) > 8. dwilcox (PR#4213) (mark.lamias@grizzard.com) > 9. Re: dwilcox (PR#4212) (Torsten Hothorn) > 10. Re: couldn't find function "setClass" (Kurt Hornik) > 11. Sweave: option family=... (Wolfram Fischer) > > > --------...

I?ve been looking at this for a couple hours--it ended up being trickier than I expected to implement well. I?ve attached a new patch here. This version scales significantly better than the existing method for both `pwilcox` and `dwilcox`. Examples are included below. I can?t think of any other ways to improve runtime at this point. That?s not to say there aren?t any, though?I?m hopeful inspection by someone else could identify more ways to save some time. This version uses Andreas?s algorithm for computing the value of `cwilcox`...

...ith_inexact -g -I../include -I../../src/include -c pnt.c -o pnt.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c ptukey.c -o ptukey.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qtukey.c -o qtukey.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c dwilcox.c -o dwilcox.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c pwilcox.c -o pwilcox.o gcc -ieee_with_inexact -g -I../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...

...src/include -c pnbeta.c cc -ieee_with_inexact -g -I../../src/include -c pnf.c cc -ieee_with_inexact -g -I../../src/include -c pnt.c cc -ieee_with_inexact -g -I../../src/include -c ptukey.c cc -ieee_with_inexact -g -I../../src/include -c qtukey.c cc -ieee_with_inexact -g -I../../src/include -c dwilcox.c cc -ieee_with_inexact -g -I../../src/include -c pwilcox.c cc -ieee_with_inexact -g -I../../src/include -c qwilcox.c cc -ieee_with_inexact -g -I../../src/include -c rwilcox.c ar cr ../lib/libmath.a mlutils.o d1mach.o i1mach.o fmax2.o fmin2.o fmod.o fprec.o fround.o ftrunc.o sign.o fsign.o fs...

Full_Name: Murray H Smith Version: 0.63.3 OS: Windows NT Submission from: (NULL) (130.216.5.57) The ppois function is displaced by -0.5. Try: > ppois(-0.5,1) [1] 0.3678794 > ppois(-0.51,1) [1] 0 > ppois(0,1) [1] 0.3678794 and > par(mfrow=c(2,1)) > x<-seq(-1,5,0.01) > plot(x,ppois(x,1),type="s",ylab="F(x)",main="Poisson CDF?") >

...I../include -I../../src/include -c pnf.c -o pnf.o gcc -g -O2 -I../include -I../../src/include -c pnt.c -o pnt.o gcc -g -O2 -I../include -I../../src/include -c ptukey.c -o ptukey.o gcc -g -O2 -I../include -I../../src/include -c qtukey.c -o qtukey.o gcc -g -O2 -I../include -I../../src/include -c dwilcox.c -o dwilcox.o gcc -g -O2 -I../include -I../../src/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 dsignran...

...I../include -I../../src/include -c pnf.c -o pnf.o gcc -g -O2 -I../include -I../../src/include -c pnt.c -o pnt.o gcc -g -O2 -I../include -I../../src/include -c ptukey.c -o ptukey.o gcc -g -O2 -I../include -I../../src/include -c qtukey.c -o qtukey.o gcc -g -O2 -I../include -I../../src/include -c dwilcox.c -o dwilcox.o gcc -g -O2 -I../include -I../../src/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 dsignran...

...-I../include -c pnchisq.c gcc -g -O2 -I../include -c qnchisq.c gcc -g -O2 -I../include -c dnbeta.c gcc -g -O2 -I../include -c pnbeta.c gcc -g -O2 -I../include -c pnf.c gcc -g -O2 -I../include -c pnt.c gcc -g -O2 -I../include -c ptukey.c gcc -g -O2 -I../include -c qtukey.c gcc -g -O2 -I../include -c dwilcox.c gcc -g -O2 -I../include -c pwilcox.c gcc -g -O2 -I../include -c qwilcox.c gcc -g -O2 -I../include -c rwilcox.c ar cr ../lib/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.o imin2.o chebyshev.o logrelerr.o lgammaco...

...t -g -I../include -I../../src/include -c pnt.c -o p= nt.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c ptukey.c -= o = ptukey.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qtukey.c -= o = qtukey.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c dwilcox.c = -o = dwilcox.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c pwilcox.c = -o = pwilcox.o gcc -ieee_with_inexact -g -I../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 -i...

I am a newbie to C (although I did some programming in Perl and Pascal) so forgive me for the language that follows. I am writing because I have a question concerning the *implementation of *(the internal function)* cwilcox* in https://svn.r-project.org/R/!svn/bc/81274/branches/R-DL/src/nmath/wilcox.c. This function is used to determine the test statistics of the Wilcoxon-Mann-Whitney U-test.

R-1.7.0 built on NetBSD 1.6, but the validation test suite failed: Machinetype: Intel Pentium III (600 MHz); NetBSD 1.6 (GENERIC) Remote gcc version: gcc (GCC) 3.2.2 Remote g++ version: g++ (GCC) 3.2.2 Configure environment: CC=gcc CXX=g++ LDFLAGS=-Wl,-rpath,/usr/local/lib make[5]: Entering directory `/local/build/R-1.7.0/src/library' >>> Building/Updating

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