similar to: permutation test on paired samples

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 ?coin? package allows one to run a two-sample permutation test using

hello, sorry for posting what may be a simple question: i do have a matrix of coordinates (positional judgments, see below) and now want to calculate and plot the corresponding error ellipse. can anyone help me with the exact steps/syntax? > xyDat X Y 1 -0.49 -2.13 2 0.91 0.48 3 0.20 -2.80 4 -0.76 -3.23 5 -0.36 2.50 6 1.38 1.24 7

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?

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,

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 same dataset, I get different p-values. I'm hoping that

Dear expeRts, I would like to use a oneway_test (from package coin) to test whether two groups differ on various variables. The variables are encoded within a data frame. Unfortunately, I obtained an error, that I don't understand. Could you please help me ? Example: library(coin) y <- as.data.frame(matrix(rnorm(200), ncol=2)) group <- as.factor(unif(100)) lapply(y,

Dear R users, I am trying to perform a hypothesis test on a marked point pattern. I would like to calculate the mean of the absolute value of the difference of marks between nearest neigbours, randomize the marks among points, then calculate this mean again. Ideally, I would test whether random mean values smaller than the observed mean value occur less than 5% of the time. I suppose 1000

Dear fellow R users, Keywords: Kruskal-Wallis, Post-Hoc, pair-wise comparisons, Nemenyi-Damico-Wolfe-Dunn test, coin package, oneway_test I am using the "oneway_test" function in the R package "coin" and I am obtaining results which I cannot believe are accurate. I do not wish to waste anyone's time and so if the following problem is rather trivial, I apologize, however I

Hi, I have a little problem that puzzles me about contradictory results returned by a bootstraped t-test (as in MASS 3rd ed p. 146) versus a permutation t-test (as in MASS 3rd ed, p 147). Data are measurements done on 100 cells in 9 slides of fish blood. With one method, cells are randomly sampled, and with the other method, the operator selects cells arbitrarily (in a way it is done usually

Hi all Consider the classic data below from Darwin on the heights of 15 pairs of zea mays (corn) plants either cross-fertilized or self-fertilized, where the goal is to see if it makes a difference. > head(ZeaMays) pair pot cross self diff 1 1 1 23.500 17.375 6.125 2 2 1 12.000 20.375 -8.375 3 3 1 21.000 20.000 1.000 4 4 2 22.000 20.000 2.000 5 5 2 19.125

Dear R friends, here I write again about the wilcox.exact() problem. I want to compare two sets of categorical data, and in one case it says "negative length vectors not allowed", and in the other one I get the error "cannot allocate vector of length ...". On http://rapidshare.com/files/298621893/wilcox.exact_trouble.Rdata.html you can download the data that cause the

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:

Dear all, I've discovered the possibility to do the Nemenyi-Damico-Wolfe-Dunn test in the library(coin); oneway_test() With the given example I am unfortunately not able to reproduce the test. What does trafo and contrMat mean? I have a dataframe with 176 elements in 7 classes. It may be a problem, that my model isn't balanced? In Class1 I happen to have 4 elements, while there are 90 in

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 --

Dear List, after updating the exactRanksumTests package I receive a warning that the package is not developed any further and that one should consider the coin package. I don't find the signed rank test in the coin package, only the Wilcoxon Mann Whitney U-Test. I only found a signed rank test in the stats package (wilcox.test) which is able to calculate the exact pvalues but unfortunately

Dear r-helpers, I have a small dataset (n<50), and I want 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

Full_Name: Patrick Hodgson Version: 2.0 OS: solaris 2.9 Submission from: (NULL) (65.94.128.161) The value reported for the parameter W in the function wilcox.exact appears to be incorrect. I have checked the reference in the help file for 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

I am trying to do a Nonparametric Tukey-type multiple comparison post-hoc test to determine which groups are significantly different. I have read the dialogue on this topic from the R-help, and am still not clear why no statistical packages include this test as an option? Is it not an appropriate test to conduct on non-normally distributed data? Is the only option to calculate it by hand

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,

Hi R users! I have the following problem: how appropriate is my aov model under the violation of anova assumptions? Example: a<-c(1,1,1,1,1,1,1,1,1,1,2,2,2,3,3,3,3,3,3,3) b<-c(101,1010,200,300,400, 202, 121, 234, 55,555,66,76,88,34,239, 30, 40, 50,50,60) z<-data.frame(a, b) fligner.test(z$b, factor(z$a)) aov(z$b~factor(z$a))->ll TukeyHSD(ll) Now from the aov i found that my model