similar to: dwilcox (PR#4212)

Has there been a recent change in the behavior of binary operators? In SparseM it was, until quite recently ok to do scalar multiplication but now, > A*4 Error in A * 4 : non-numeric argument to binary operator > 4*A Error in 4 * A : non-numeric argument to binary operator > A%*%A An object of class "matrix.csr" Slot "ra": [1] 1 1 1 Slot "ja": [1] 1 2 3

Hi, I believe I have significantly improved [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

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

It seems benign, but has implications since checking time is actually not a cheap operation: adding jus ta time check alone incurs a penalty of ca. 700% compared with the time it takes to call R_CheckUserInterrupt(). Generally, it makes no sense to check interrupts at every iteration - you'll find code like if (++i % 10000 == 0) R_CheckUserInterrupt(); in loops to make sure it's not called

Does anyone have a piece of code or know how I can use R to generate a table of critical values for the Mann-Whitney (aka Wilcoxon Rank Sum) test. Ideally, I'd like a table that contains the critical values for any two samples of size 3 through 30. I could use Monte Carlo simulation or the normal approximation when n1 and n2 are greater than, 10, but I figured someone may know how to

tl;dr: R_CheckUserInterrupt() can be a performance bottleneck within GUIs. This also affects functions in the 'stats' package, which could be improved by changing the position of calls to R_CheckUserInterrupt(). Dear all, Recently I was puzzled because some code in a package under development, which consisted almost entirely of a .Call() to a function written in

A more generic solution would be for R to throttle calls to R_CheckUserInterrupt(), because it makes no sense to check 1000 times per second if a user has interrupted, but it is difficult for the caller to know when R_CheckUserInterrupt() has been last called, or do it regularly without over-doing it. Here is a simple patch: https://github.com/r-devel/r-svn/pull/125 See also:

This seems like a great idea. Would it help to escalate this to a post on R-bugzilla, so it is less likely to fall through the cracks? On 12/17/24 09:51, Jeroen Ooms wrote: > A more generic solution would be for R to throttle calls to > R_CheckUserInterrupt(), because it makes no sense to check 1000 times > per second if a user has interrupted, but it is difficult for the > caller to

On 12/18/24 01:19, Simon Urbanek wrote: > It seems benign, but has implications since checking time is actually not a cheap operation: adding jus ta time check alone incurs a penalty of ca. 700% compared with the time it takes to call R_CheckUserInterrupt(). Generally, it makes no sense to check interrupts at every iteration - you'll find code like if (++i % 10000 == 0)

The function 'wilcox.test' in R and S gives (almost) identical results (see below). 'qwilcox' however, does not: > qwilcox(p,5,5) p: 0.025 0.975 -------------------- R> 3 22 S> 18 37 I originally wanted to ask a questions, but then I found the answer. Given the confusion I run into, I wonder if this experience is worth reporting. The

Hi again ! Thanks for the info on updating the config.site file which I have done. I have also added -lm in the Makeconf manually because this is needed explicitly for DEC cc. However, there are still a few problems when linking some of the files as you can see from the enclosed log. Ciao, Andreas ------------------------------------------------------- R-0.63.3>make make[1]: Entering

Hi, I would like to do a MannWitney test. Can anyone help me with the propper command? Thanks, Margarida

Hi, I am trying to convert a colleague from using SPSS to R, but am having trouble generating a result that is similar enough to a bootstrapped cox regression analysis that was run in SPSS. I tried unsuccessfully with bootcens, but have had some success with the bootcov function in the rms package, which at least generates confidence intervals similar to what is observed in SPSS. However, the

Hi, Just downloaded the R-0.62.4 of R and tried to install it. With the standard procedure : ./configure make At the end of the compilation I got the following message : ld: ../lib/libunix.a(system.o): main: multiply defined fort: Severe: Failed while trying to link. *** Exit 1 Stop. *** Exit 1 Stop. *** Exit 1 Stop. *** Exit 1 Stop. I attach the printout after the ./configure and make

How does R compute the p-value in the Wilcoxon Rank Sum Test? If I have the test statistics, can I get the p-value using the function pwilcox? Thanks -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

I am trying to graph two planes on the same graph using persp(). I can only get one plane to plot at a time. Can someone explain how I can graph two planes on the same graph using persp? I've looked throught the documentation, but cannot find any references to appending a persp plot. Thanks. Sincerely yours, Mark J. Lamias Statistical Consultant

R-users, I have the following piece of code which I am trying to run on a dataframe (aga2) with about a half million records.  While the code works, it is extremely slow.  I've read some of the help archives indicating that I should allocate space to the p1 and ags1 vectors, which I have done, but this doesn't seem to improve speed much.  Would anyone be able to provide me with advice on

> > > Performance statistics are interesting. If we assume the two populations > have a total of `m` members, then this implementation runs slightly slower > for m < 20, and much slower for 50 < m < 100. However, this implementation > works significantly *faster* for m > 200. The breakpoint is precisely when > each population has a size of 50; `qwilcox(0.5,50,50)`

Hello, I've allways had some problems compiling R (R-62.x and R-63.x) on a AIX 4.2.1 System to the end: make produces following output: ... ... mkdir ../../../library/eda/R ld: 0711-327 WARNING: Entry point not found: __start ld: 0711-244 ERROR: No csects or exported symbols have been saved. gnumake[4]: *** [eda.so] Error 8 gnumake[3]: *** [all] Error 1 gnumake[2]: *** [R] Error 1 gnumake[1]:

Hello, I've allways had some problems compiling R (R-62.x and R-63.x) on a AIX 4.2.1 System to the end: make produces following output: ... ... mkdir ../../../library/eda/R ld: 0711-327 WARNING: Entry point not found: __start ld: 0711-244 ERROR: No csects or exported symbols have been saved. gnumake[4]: *** [eda.so] Error 8 gnumake[3]: *** [all] Error 1 gnumake[2]: *** [R] Error 1 gnumake[1]: