similar to: Errors in wilcox family functions

While using the Hodges Lehmann Mean in DescTools (DescTools::HodgesLehmann), I found that it generated incorrect answers (see <https://github.com/AndriSignorell/DescTools/issues/97> https://github.com/AndriSignorell/DescTools/issues/97). The error is driven by the existence of tied values forcing wilcox.test in Base R to switch to an approximate algorithm that returns incorrect results - see

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

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 everyone, My apologies if this isn?t the right place to submit this?I?m new to the R-devel community and still figuring out what is where. If people want to skip my writeup and just look at the code, I?ve made a repository for it here: https://github.com/ahl27/new_dendrapply/tree/master. I?m not quite sure how to integrate it into a fork of R-devel; the package structure is different from

It looks like a homework assignment. It also looks like you didn't read the documentation carefully enough. The 'len.out' argument in seq is solely for specifying the length of a sequence. The 'quantile' function omputes the empirical quantile of raw data in the vector 'x' at cumulative probabilit(y)(ies) given in the weights' argument, with interpolation I'm

When I try to install the R RPM under Mandrake Linux 7.1 (intel), I get the following error: [root at edgeworth programs]# rpm -ivh R-base-1.1.1-1.i386.rpm error: failed dependencies: libreadline.so.3 is needed by R-base-1.1.1-1 [root at edgeworth programs]# I do have realine installed: [root at edgeworth programs]# rpm -qa | grep readline readline-4.1-4mdk and

Greetings, I am having a problem with DescTools::Quantile (a function computing quantiles from weighted samples): # these sum to one probWeights = c( 0.0043, 0.0062, 0.0087, 0.0119, 0.0157, 0.0204, 0.0257, 0.0315, 0.0378, 0.0441, 0.0501, 0.0556, 0.06, 0.0632, 0.0648, 0.0648, 0.0632, 0.06, 0.0556, 0.0501, 0.0441, 0.0378, 0.0315, 0.0257, 0.0204, 0.0157, 0.0119, 0.0087,

Hello Aidan, Sorry for dropping this for a while. ? Thu, 2 Mar 2023 21:03:59 +0000 "Lakshman, Aidan H" <AHL27 at pitt.edu> ?????: > //after > curnode = eval(lang3(R_Bracket2Symbol, parent->node, DEND_IND), env); lang3() always constructs a new language object. If you do end up using eval(), it may make sense to move lang3() out of the loop and reuse the existing object

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

Dear R users, I have a statistical question regarding on the Edgeworth expansion of a random vector. Assume x is a d-dimensional random vector, and whats the expression of its Edgeworth expansion up to order 4? I only get the formulation for a scalar random variable. So anyone can help me out on this problem? Thanks for your comments or advice. Fred [[alternative HTML version deleted]]

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

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

I have just upgraded from Red Hat Linux 5.2 to version 6.0 (with GNOME) on an intel box. I have also upgraded to the latest R version (0.64.1-2 rpm) and have also installed the corresponding GNOME rpm. However, I am not able to start R: [cribari at edgeworth install]$ R /usr/lib/R/bin/R.binary: error in loading shared libraries: /usr/lib/R/bin/R.binary: undefined symbol: __setfpucw [cribari at

Hi All, I am trying to estimate a functional form in terms of cumulants / moments from the results of a MC simulation. My first idea was to use an Edgeworth expansion, but this looks pretty difficult to implement. Before I go ahead and write some code to do this, could anyone tell me if this has already been done, and which package I could find it in if it has? If it hasn't, does anyone know

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

Hello, I have three numerical variables and I would like to test if their correlation is significantly different. I have seen that there is a package that "Test the difference between two (paired or unpaired) correlations". [https://www.personality-project.org/r/html/paired.r.html] However, there is the need to convert the correlations to "z scores using the Fisher r-z

Full_Name: David Simcha Version: 2.10 OS: Windows XP Home Submission from: (NULL) (173.3.208.5) > a <- c(1:10) > b <- c(1:10) > cor.test(a, b, method = "spearman", alternative = "greater", exact = TRUE) Spearman's rank correlation rho data: a and b S = 0, p-value < 2.2e-16 alternative hypothesis: true rho is greater than 0 sample estimates:

Hello! I want to calculate simultaneous confidence intervals for a nominal variable with three categories: "yes", "no", "partially" and I expect that far more than 5 samples fall into each category. I have read that Glaz & Sison's method is only appropriate for variables with 7 or more categories. Therefore, the Goodman method seems like a good idea. I have

I am trying to install R under Red Hat Linux 5.2, but the following error happens: [root at edgeworth cribari]# rpm -ivh R-base-0.63.2-1.i386.rpm failed dependencies: libncurses.so.3.0 is needed by R-base-0.63.2-1 I must be missing something very simple. I would appreciate any guidance. Thanks. FC. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help

how do I check the rank of a matrix ? say A= 1 0 0 0 1 0 then rank(A)=2 what is this function? thanks I did try help.search("rank"), but all the returned help information seem irrelevant to what I want. I would like to know how people search for help information like this. rank(base) Sample Ranks SignRank(stats) Distribution of the