search for: wilcox

Hello, Writing to share some things I've found about wilcox.test() that seem a a bit inconsistent. 1. Inf values are not removed if paired=TRUE # returns different results (Inf is removed): wilcox.test(c(1,2,3,4), c(0,9,8,7)) wilcox.test(c(1,2,3,4), c(0,9,8,Inf)) # returns the same result (Inf is left as value with highest rank): wilcox.test(c(1,2,3,4),...

...ack and guidelines about issues 3 and 4 as well. On 2019-12-07 21:59, Martin Maechler wrote: >>>>>> Karolis Koncevi?ius >>>>>> on Sat, 7 Dec 2019 20:55:36 +0200 writes: > > > Hello, > > Writing to share some things I've found about wilcox.test() that seem a > > a bit inconsistent. > > > 1. Inf values are not removed if paired=TRUE > > > # returns different results (Inf is removed): > > wilcox.test(c(1,2,3,4), c(0,9,8,7)) > > wilcox.test(c(1,2,3,4), c(0,9,8,Inf)) > > > #...

...finities, we are computationally looking for roots over > infinite interval which results in an error. I suspect this is the > reason Inf values were removed in the first place. Maybe. It's still wrong to be done "up front". I'm sure 98% (or so ;-) of all calls to wilcox.test() do *not* use conf.int = TRUE > Just a note, I found a few more errors/inconsistencies when requesting > confidence intervals with paired=TRUE (due to Infinities being left in). > Current error in Inf-Inf scenario: > wilcox.test(c(1,2,Inf), c(4,8,Inf), paired...

hi, when performing a wilcox.test or a wilcox.exact i get results that looks like this: wilcox.exact(x, mu=.5) Exact Wilcoxon signed rank test data: x V = 207, p-value = 0.0006905 alternative hypothesis: true mu is not equal to 0.5 the way i understand the wilcox.test (or wilcox.exact) the V-value represents th...

I'm having trouble with the wilcox.test command in R. To demonstrate the anomalous behavior of wilcox.test, consider > wilcox.test(c(1.5,5.5), c(1:10000), exact = F)$p.value [1] 0.01438390 > wilcox.test(c(1.5,5.5), c(1:10000), exact = T)$p.value [1] 6.39808e-07 (this calculation takes noticeably longer). > wilcox.test(c(...

I get these BIZARRE results from wilcox.test() When INCREASING the number of samples i get INCREASED p-values. When increasing the number of samples further, the p-values goes down again. This seems really bizarre! Can anyone explain why this is so?! Example: > w <- wilcox.test(c(1:40),(c(1:40)+100)) > w$p.value [1] 1.8603...

>>>>> Karolis Koncevi?ius >>>>> on Sat, 7 Dec 2019 20:55:36 +0200 writes: > Hello, > Writing to share some things I've found about wilcox.test() that seem a > a bit inconsistent. > 1. Inf values are not removed if paired=TRUE > # returns different results (Inf is removed): > wilcox.test(c(1,2,3,4), c(0,9,8,7)) > wilcox.test(c(1,2,3,4), c(0,9,8,Inf)) > # returns the same result (Inf is left...

...lt;-read.table(file="C:\\Documents .....txt") dat<-as.data.frame(dat) dimnames(dat)[[1]]<-paste("g",c(1:nrow(dat)),sep="") l<-c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1) rbind(l,dat) dat<-as.matrix(dat) x<-dat[1:100,] wilcox.test<-function(x,s1,s2) { x1<-x[s1] x2<-x[s2] x1<-as.numeric(x1) x2<-as.numeric(x2) wilcox.out<-wilcox.test(x,y=NULL,exact=FALSE,mu=0,correct=TRUE,alternative="two.sided",paired=FALSE) out<-as.numeric(wilcox.out$statistic) return(out) } wilcox.test.run<-apply(x,1...

...ptions(digits=20) > unique(c(4,3,2)-c(3,2,1)) [1] 1 > unique(c(0.4,0.3,0.2)-c(0.3,0.2,0.1)) [1] 0.100000000000000033307 0.099999999999999977796 0.100000000000000005551 On 2019-12-07 1:55 p.m., Karolis Koncevi?ius wrote: > Hello, > > Writing to share some things I've found about wilcox.test() that seem a > a bit inconsistent. > > 1. Inf values are not removed if paired=TRUE > > # returns different results (Inf is removed): > wilcox.test(c(1,2,3,4), c(0,9,8,7)) > wilcox.test(c(1,2,3,4), c(0,9,8,Inf)) > > # returns the same result (Inf is left as value...

Hi, I became a little bit confused when working with the Wilcoxon test in R. As far as I understood, there are mainly two versions: 1) wilcox.test{stats}, which is the default and an approximation, especially, when ties are involved 2) wilcox_test{coin}, which does calculate the distribution _exactly_ even, with ties. I have the following scenario: #---Begi...

...ionally looking for roots over infinite interval >> which results in an error. I suspect this is the reason >> Inf values were removed in the first place. > Maybe. It's still wrong to be done "up front". I'm sure > 98% (or so ;-) of all calls to wilcox.test() do *not* use > conf.int = TRUE >> Just a note, I found a few more errors/inconsistencies >> when requesting confidence intervals with paired=TRUE >> (due to Infinities being left in). >> Current error in Inf-Inf scenario: >> wilcox...

Full_Name: David Wooff Version: 1.7.0 OS: i686-pc-linux-gnu Submission from: (NULL) (129.234.4.10) wilcox.test exits with error message when confidence interval required, under some situations. I suspect this occurs when the data contain a zero and for some data lengths only: print(wilcox.test(c(2,1,4,3,6,-5,0),conf.int=T)) fails print(wilcox.test(c(2,1,4,3,6,-5,0,1),conf.int=T)) works print(wilco...

...finite interval which results in an error. I suspect this is the reason Inf values were removed in the first place. Just a note, I found a few more errors/inconsistencies when requesting confidence intervals with paired=TRUE (due to Infinities being left in). Current error in Inf-Inf scenario: wilcox.test(c(1,2,Inf), c(4,8,Inf), paired=TRUE, conf.int=TRUE) Error in if (ZEROES) x <- x[x != 0] : missing value where TRUE/FALSE needed NaN confidence intervals: wilcox.test(c(1:9,Inf), c(21:28,Inf,30), paired=TRUE, conf.int=TRUE) Wilcoxon signed rank test with continuity correction...

Hi, I am very confused with constructing the wilcox.test in R. I have two populations 'original' and 'test'. I want to know if the 'test' is generally 'lower' than original. I use alpha of 0.05. So do I write the function as wilcox.test(original, test, alternative="l")? or wlcox.test(original, test, alterna...

Dear All, I am using wilcox.test to test two samples, data_a and data_b, earch sample has 3 replicates, suppose data_a and data_b are 20*3 matrix. Then I used the following to test the null hypothesis (they are from same distribution.): wilcox.test(x=data_a, y=data_b, alternative="g") I got pvalue = 1.9080617086331...

...message("Group number for t-test is out of range.") data.norm.test = data.norm.zo #normalized data nrow <- nrow(data.norm.test) ncol <- ncol(data.norm.test) p.ttest = rep(NA, nrow) #parametric test -- t test p.wilcox = rep(NA, nrow) #non-parametric test -- Mann-Whitney test fold.mean <- rep(NA, nrow) mean.na <- rep(NA, nrow) mean.nb <- rep(NA, nrow) sig.peak.all = data.norm.test g_1 = sig_test[1] g_2 = sig_test[2] col_s1 = 1; col_s...

Dear R Core, pairwise.wilcox.test does not handle "paired = TRUE" correctly; e.g. set.seed(13) x <- rnorm(20) g <- c(rep(1, 10), rep(2, 10)) wilcox.test(x ~ g)$p.value # 0.075 pairwise.wilcox.test(x, g)$p.value # 0.075, o.k wilcox.test(x ~ g, paired = TRUE)$p.value # 0.105 pairwise.wilcox.test(x, g, paired = T...

I have a question concerning the Wilcoxon signed-rank test, and specifically, which R subroutine I should use for my particular dataset. There are three different commands in R (that I'm aware of) that calculate the Wilcoxon signed-rank test; wilcox.test, wilcox.exact, and wilcoxsign_test. When I run the three commands on the sam...

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

Does anyone know why wilcox.exact gives W-statistic 6 instead of 12 as indicated below. 12 is the rank sum of group 0 of x, which is the linear statistic computed by wilcox_test. y<-c(1,2,3,4,5) x<-c(1,1,0,0,0) (a) wilcox.exact wilcox.exact(y~x) Exact Wilcoxon rank sum test data: y by x W = 6, p-value = 0.2 alter...