search for: exactranktests

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070629/17f7d932/attachment.pl

Hi, In the exactRankTest package, I've become aware that you can get anomalous p-values (i.e., above 1) from the wilcox.exact method, as in: > wilcox.exact(c(-0.6,0.8,-0.5)) Exact Wilcoxon signed rank test data: c(-0.6, 0.8, -0.5) V = 3, p-value = 1.25 alternative hypothesis: true mu is not equal to 0 This is disturbing. Has anyone encountered this before, and if so

Full_Name: Michael W. Grant Version: 1.7.1(2003-06-16) OS: Windows 98 Second Edition 4.10.2222 A Submission from: (NULL) (24.55.165.101) I just installed version 1.7.1 tonight. In the process I installed via CRAN from the menu several packages. One installation, exactRankTests, had a distinct anomaly. A DLL attempt to change FPU control word as seem in the last line of the following listing: START LISTING: > install.packages(choose.files('',filters=Filters[c('zip','All'),]), .libPaths()[1], CRAN = NULL) > local({a <- CRAN.packages() +...

..."cannot allocate vector of length ...". On http://rapidshare.com/files/298621893/wilcox.exact_trouble.Rdata.html you can download the data that cause the trouble. (Attention: 10 Downloads only, so you may have to send me an e-mail so I can mail you the Rdata file directly.) > library(exactRankTests) Package ?exactRankTests? is no longer under development. Please consider using package ?coin? instead. > wilcox.exact(group_1,group_2,exact=T) Error in .Call("cpermdist2", ma = as.integer(m), mb = as.integer(col), : negative length vectors are not allowed > wilcox.exact(gro...

Given the following data, and hypothesized median M.0 I've found a method to implement the Wilcoxon signed-rank test. Data: (with one zero difference and tied ranks) x <- c(136, 103, 91, 122, 96, 145, 140, 138, 126, 120, 99, 125, 91,142, 119, 137) M.0 <- 119 > library(exactRankTests) Package ?exactRankTests? is no longer under development. Please consider using package ?coin? instead. > wilcox.exact(x, mu=M.0) Exact Wilcoxon signed rank test data: x V = 65.5, p-value = 0.771 alternative hypothesis: true mu is not equal to 119 I've been unable to implement thi...

...how to take the ties into account? (Or a chance that the correction is taken into account for the stats package?) Stefan Grosse Take the following example from Bortz/Lienert/Boehnke: > x1<-c(9,14,8,11,14,10,8,14,12,14,13,9,15,12,9) > x2<-c(13,15,9,12,16,10,8,13,12,16,9,10,16,12,9) # exactRankTests package: > wilcox.exact(x1,x2,paired=TRUE) Exact Wilcoxon signed rank test data: x1 and x2 V = 13, p-value = 0.1367 alternative hypothesis: true mu is not equal to 0 # wilcox.test by stats package: > wilcox.test(x1,x2,paired=TRUE,exact=TRUE) Wilcoxon signed rank test wi...

Hi, I'm trying to run a permutation test on paired samples. First I tried the package "exactRankTests": require("exactRankTests") x <- c(1.83,0.50,1.62,2.48,1.68,1.88,1.55,3.06,1.30) y <- c(0.878,0.647,0.598,2.05,1.06,1.29,1.06,3.14,1.29) wilcox.test(x,y,paired = TRUE,alternative = "greater") perm.test(y,x,paired = TRUE,exact = TRUE,alternative = "greater")...

Dear collegues, Using maxstat I am getting the following: > blood <- maxstat.test(Surv(SUPER, FV)~ZAP,data=zap70, smethod="LogRank") Error in maxstat(y = structure(c(24.4301369863014, 26.4164383561644, 18.7835616438356, : couldn't find function "cscores" I do not know the meaning of this problem. Could you please help me on dat? Thank you in advance for your

Apropos the question about permutation tests in multiple regression: We do have perm.test in package ExactRankTests, but it does one- and two-sample tests, as in t.test, wilcox.test, etc. There doesn't seem to be an exact version of the permutation test for correlations, i.e., the one that could be estimated using replicate(10000, cor(x,sample(y))) # or other values of 10000 or, of course, computed exactly...

Hi, I am running the following example which can be found on page 12 of the pdf file of COIN package wt<-wilcox_test(pd~age,data=water_transfer,distribution="exact", conf.int=TRUE) "wt" actually contains the estimate of difference in location and the confidence interval of it. I am just wondering how can I extract these values? From the examples, I understand that the

Dear R friends, I want to compare two datasets and I get the message Error in .Call("cpermdist2", ma = as.integer(m), mb = as.integer(col), : negative length vectors are not allowed after specifying the exact test. I'm using the exactRankTests package. Do you suggest me using the coin library, or is there anything "wrong" with my data? Kind regards, David -- Jetzt kostenlos herunterladen: Internet Explorer 8 und Mozilla Firefox 3.5 - sicherer, schneller und einfacher! http://portal.gmx.net/de/go/atbrowser

Using the package ?exactRankTests? one can execute a one-sample permutation test for a hypothesized location parameter of 0 like: perm.test(rnorm(30,0)) The package ?exactRankTests? seems now to be deprecated in favor of the ?coin? package which as I understand is a superset of ?exactRankTests? in terms of functionality. The ?co...

Dear All, I have two simple data samples (no groups or factors, etc.) and would just like to compute the two-sample Wilcoxon Rank Sum test using the wilcox_test function contained in the coin package, which is reportedly better than the regular wilcox.test function because it performs some adjustment for ties. Would anyone know how to craft a script to perform this task? Much appreciated. Janh

...KernSmooth Matrix NISTnls Oarray PHYLOGR PTAk RArcInfo RMySQL RODBC RPgSQL RandomFields RmSQL Rstreams Rwave SASmixed SuppDists VR XML acepack adapt akima ash bindata blighty boot bootstrap bqtl car cclust cfa chron cluster cmprsk coda conf.design cramer date diamonds dse e1071 ellipse event.chart exactRankTests fastICA fdim fields foreign fracdiff gafit gee geoR gld gregmisc gss ineq leaps lgtdl lmtest locfit logspline lokern lpridge maptree maxstat mclust mda meanscore mgcv mlbench muhaz multiv mvnmle mvtnorm netCDF nlme nlrq norm odesolve oz panel pcurve permax pinktoe pixmap polymars polynom princurve...

...ant to compute the Hodges Lehmann exact confidence interval. So far, I know that "pairwiseCI" has the function "HL.diff". The description is as follows : HL.diff calculates the Hodges-Lehmann confidence interval for the difference of locations by calling wilcox.exact in package exactRankTests ; But when I check wilcox.exact, it says " an exact p-value is computed if the samples contain less than 50 finite values and there are no ties " However, there are many ties in the dataset, so is there any other way I can compute the exact CI for a small dataset even with ties? Any he...

...this function (Myles & Hollander 1973, as well as 2nd ed. 1999 by same authors) and it is clear that W is the sum of the ranks of the data set with the smaller number of members. Here is sample input/output: > x<-c(1,2,3,4,5,6,7,8,9,10) > y<-c(11,12,13,14,15) > library(exactRankTests) > wil<-wilcox.exact(x,y) > wil Exact Wilcoxon rank sum test data: x and y W = 0, p-value = 0.000666 alternative hypothesis: true mu is not equal to 0 In this case W should be the sum of the ranks of y (coincidentally also the values of the y data)= 11+12+13+14+15 = 65. W of...

...nts). so far, I've figured out how to: -create a marked ppp object: ms.ppp -calculate my test statistic: teststat <- mean(abs(markstat(ms.ppp, diff, N=2))) -randomly allocate marks to a point pattern: Y <- rlabel (ms.ppp, labels=ms.ppp$marks, permute=TRUE) I have looked at perm.test{exactRankTests} and perused the R help archive but haven't been able to find how to work the permutation test with a marked ppp object. I thank you in advance for any help, J-F Savard Doctoral Candidate University of Maryland

Hello, can someone please help me with coding a function for a bootstrap WMW test (package boot, R under Windows, version 1.6.2)?

Help needed! Has anyone access to StatXact? I just hacked exact two-sided p-values for rank tests (for package exactDistr, which will move to CRAN/contrib as exactRankTests soon ;-) and would like to compare the results of my implementation to that of StatXact. Could someone please calculate the exact one-sided (both greater and less) and two-sided p-values? # Data from the StatXact-4 manual, page 221 treat <- c(94, 108, 110, 90) contr <- c(80, 94, 85, 90, 90,...

...not calculating the statistics correctly. The following text shows the commands I use in R to calculate the signed-rank test using these different R commands: Thanks in advance for any assistance. --Mike ################################################################################# library(exactRankTests) #this loads the package for calculating the modified signed-rank test library(coin) #this adds additional routines for the wilcoxon signed-rank test and the Pratt signed-rank test # # Data from Page 658 of Pratt x <- c(1, 1, 1, 1, 1, 7, 10, 12, 13, 16, 17) y <- c(1, 1, 3, 4, 6, 1, 1, 1,...