similar to: R-alpha: planned update of ctest

This was filed as a bug report on the Debian r-base package. It is more properly a bug report on the ctest package in R. The default method for wilcox.test manipulates x and y without checking the class or data.class of these objects. Possible solutions are - create wilcox.test.factor (if appropriate) - check the class and/or data.class of x and y in wilcox.test.default and produce error

Sorry for leaving this one in my mail box for so long, but - well, I suppose you know what I mean. (I'm shifting it over to r-devel, so I'll include all your original text) Kurt Hornik <hornik@ci.tuwien.ac.at> writes: > Well, ctest is not making progress as quickly as I wanted it ... > Anyway, here are a few questions/remarks. > > * I am still a bit confused about

I have the following problem. Consider a `classical' test which works for k .ge. 2 samples. Possible interfaces are e.g. xxx.test(x, g) x ... all data, g ... corresponding groups xxx.test(x1, ..., xk) xxx.test(list(x1, ..., xk)) etc etc. Clearly, the first and the second one are nice, but cannot be combined without making `g' (i.e., `group') a named argument. Hence, in

I frequently want to test for differences between animal size frequency distributions. The obvious test (I think) to use is the Kolmogorov-Smirnov two sample test (provided in R as the function ks.test in package ctest). The KS test is for continuous variables and this obviously includes length, weight etc. However, limitations in measuring (e.g length to the nearest cm/mm, weight to the nearest

Appology in advance if this has been asked and answered. I am getting different answers using Fisher's exact test in the package ctest. For example: x_cbind(c(8, 22), c(31, 29)) fisher.test(x) Fisher's Exact Test for Count Data data: x p-value = 0.04024 alternative hypothesis: two.sided However, in SAS or other packages on the net the answer I get is p-value= 0.02664. Is this

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

A simple question - using the following fishers test it appears that the P value is significant, but the CI includes 1. Is this result correct? > data.50p10min <- matrix(c(16,15, 8, 24),nrow=2) > fisher.test(data.50p10min) Fisher's Exact Test for Count Data data: data.50p10min p-value = 0.03941 alternative hypothesis: true odds ratio is not equal to 1 95

o <- (structure(list(sand.silt = structure(list(statistic = structure(185, .Names = "W"), parameter = NULL, p.value = 0.0478835773838087, null.value = structure(0, .Names = "location shift"), alternative = "two.sided", method = "Wilcoxon rank sum test with continuity correction", data.name = ".column by site"), .Names =

We've seen the integer overflow problem in ks.test before, easily solved. The help page says x and y must be numeric, so this is user error. I've added tests to the code. Why do people file bug reports without reading the help/man page? On Tue, 14 Jan 2003 bates@stat.wisc.edu wrote: > This was filed as a bug report on the Debian r-base package. It is > more properly a bug

R 0630 for windows > library(ctest) > binom.test(7,10,p=0.3, alternative="two.sided") returns a p-value of =< 2.2e-016 and a warning In Splus 3.4 > binom.test(7,10,p=0.3, alternative="two.sided") returns a p-value of 0.0106 I think it is the max(v[v<=(1+eps)*PVAL]) causing the problem... max() of an empty vector....... Mai Z

Hi, I have a data set: > Dataset X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 X11 X12 X13 X14 X15 X16 X17 1 user1 m 22 19 28 24 12 18 9 7 4 5 4 7 5 7 9 2 user2 f 25 19 23 18 18 15 6 8 6 6 7 10 7 7 7 3 user3 f 28 21 24 18 15 12 10 6 7 9 5 10 5 9 5 4 user4 f 26 19 26 21 12 18 6 6 5 1 3 8 6 5 6 5 user5 m 21 22 26 18 9 6 4 6 1

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

Hi all, I have a distribution, and take a sample of it. Then I compare that sample with the mean of the population like here in "Wilcoxon signed rank test with continuity correction": > wilcox.test(Sample,mu=mean(All), alt="two.sided") Wilcoxon signed rank test with continuity correction data: AlphaNoteOnsetDists V = 63855, p-value = 0.0002093 alternative hypothesis:

Hi! I have some data that looks like this up down percentaje uew_21 20 14 58.82 uew_20_5 27 40 40.29 uew_20 8 13 38.09 uew_19_5 17 42 28.81 So I have 4 experimental conditions and I am counting number of animals in the up and down compartment and the calculating the percentage, I want to know which one of the conditions is different from each other. If the data wouldn't be percentage

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, I' am doing a stats project using R to work out the size of a t-test and wilcoxon test depending on the distribution and sample size. I just can't get it to work - I want to put my results from the function size() into an array.At the moment I keep getting the error message:Error in res[distribution, test, samplesize] <- results : subscript out of boundsCan anyone tell me where

I've put all ``current'' add-on packages into CRAN's source/contrib tree and created an INDEX file (attached below). As you can see, currently we have acepack bootstrap ctest date e1071 fracdiff gee jpn snns splines survival4 (Yes, e1071 and jpn are new ... more on the latter in a later mail.) In the near future, I am hoping for the following: oz (Bill

fisher.test(matrix(c(1,20,246,6873),2),hybrid=F) Error in if (f(lower, ...) * f(upper, ...) >= 0) stop("f() values at end points not of opposite sign") : missing value where logical needed Thomas Lumley Assistant Professor, Biostatistics University of Washington, Seattle -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list --

Por si alguno pudiera ayudarme. Al realizar el t.test para una muestra, junto con el valor de t y el p-valor, la función proporciona la estimación de la media y su INTERVALO DE CONFIANZA. Desde el punto de vista de la estadística de rangos esto se puede hacer mediante: > iris$MEDIANA <- with(iris, 2.95) > median(iris$Sepal.Width - iris$MEDIANA, na.rm=TRUE) # median difference [1]

1. There is, I believe, some redundant code in the calculation of the test statistic in ks.test in the package ctest. Lines 34-37 of the code read x <- y(sort(x), ...) - (0:(n - 1))/n STATISTIC <- switch(alternative, two.sided = max(abs(c(x, x - 1/n))), greater = max(c(x, x - 1/n)), less = -min(c(x, x - 1/n))) Lines 35-37 could read