similar to: Ansari vs wilcoxon vs ks test

The help page for ansari.test() says (emphasis added): By default (if exact is not specified), an exact p-value is computed if both samples contain less than 50 finite values **and there are no ties**. Otherwise, a normal approximation is used. However, this does not appear to be the case in R 2.4.1 or R 2.5.0. In fact, even example(ansari.test) produces a warning Warning message:

A recent message on ansari.test() prompted me to play with the examples. This doesn't work for me in R version 2.4.1 R> ansari.test(rnorm(100), rnorm(100, 0, 2), conf.int = TRUE) Error in uniroot(ab, srange, tol = 1e-04, zq = qnorm(alpha/2, lower = FALSE)) : object "ab" not found It looks like there's a small typo in ccia() inside ansari.test.default() in which

Hello. I am probably wrong... I am wondering if it could have a mistake in the code of the ansari.test function. For me, it seems that the function do not recover the p value at the correct side of the normal law N(0, 1) when it use the normal approximation (presence of ties) in a one tailed test. Exemple : quanti<-c(197, 205, 228, 234, 237, 195, 233, 226, 244, 227, 259, 185, 198, 253, 207,

Full_Name: Gael Millot Version: 2.2.0. OS: XP Submission from: (NULL) (195.220.102.20) Hello. I sent an Email in r-help without answer for the moment. I am wondering if it could have a mistake in the code of the ansari.test function. For me, it seems that the function do not recover the p value at the correct side of the normal law N(0, 1) when it use the normal approximation (presence of

Full_Name: Wei Xu Version: 1.5.1 OS: WindowsME Submission from: (NULL) (63.215.238.92) The P-value for a two.sided test is not consistent with the confidence interval. For example, P-value=0.1372, but the 95% CI doesn't include the H0 value(1). > x [1] 0.80 0.83 1.89 1.04 1.45 1.38 1.91 1.64 0.73 1.46 > y [1] 1.15 0.88 0.90 0.74 1.21 >

Probably, it's an obvious info, but I have not found anything in R FAQ concerning this feature/bug. The results of ks.test and wilcoxon.test (in the Mann-Whitney version, paired = 'FALSE') don't coincide with the results from the other statistical packages, e.g. Statistica, Medcalc, and (as for MW test) from the numerous online MW tests. E.g. Statistica p-value=0.0435353 Medcalc

Dear members of the R help list, I want to do a hierarchical glm with binomial family but am unsure about how to write the syntax which involves nesting. I want to test whether the risk of being attacked by Herbivores for Meadowsweet plants is significantly dependent on the Distance to heterospecific source plants. Dependent variable = Herbivory (yes/no) Explanatory continuous variable =

Hi, is there a way to analyze subpopulations (e.g. women over 50, those who answered "yes" to a particular question) in a survey using Survey package? Other packages (e.g. Stata, SUDAAN) do this with a subpopulation option to identify the subpopulation for which the analysis shoud be done. I did not see this option in the Survey package. Is there another way to do this?

I am trying to get R 1.0 running on the Mac. The main target is MacX. Anyone else working on that? For the recent Macintosh system, I am trying to compile R using MachTen. The R core compiles and runs without any changes, using the Unix make files. The libraries give some problems, presumably due to handling of shared libraries. It seems I am missing something to link. Before I spend my time

Hi, I'm trying to compile R-1.7.1 from source (on a RedHat 8.0) instead of using the binary version, as it has often been advised. However I don't manage to find a solution to the following error which occurs during the make procedure : /!\-------------------------------------------/!\ building package 'ctest' mkdir -p -- ../../../library/ctest/R mkdir -p --

I would like to be able to use two functions; qansari and pansari which are found in the function ansari.test. How can I evaluate these functions independently? Thanks. --Dale For example, when I load the function ... qansari <- function(p, m, n) { .C(R_qansari, as.integer(length(p)), q = as.double(p), as.integer(m), as.integer(n))$q } and

Hi all, Trying to install R-1.7.1 on a RedHat 8.0 platform, I have a few problems. R gets installed without default packages (but base and ctest) : make script fails at the end of the procedure (configure is made successfully). This is the (translated) log I have : /gcc -I../../../../include /usr/include/mysql -D__NO_MATH_INLINES -mieee-fp -fPIC -g -O2 -c ansari.c -o ansari.o gcc: cannot

Hi! I need to run a wilcoxon (Mann-whitly, in fact) test with bonferroni correction, as I am running 10 consecutive wilcoxon test not independent, and I know that bonferroni will partially correct for this problem, but I have no idea how to do it with R, I have been looking in the archive but couldn't understand how to do it. The format I am using at the moment is r4_o <- [1]

Hello there, I would like to do a Wilcoxon matched pairs signed rank sum test in R and tried the function wilcoxon.test. Is it in the "base" library? If not, please let me know which library it is in. Thank you in advance for your help Lisa Wang Biostatistician Princess Margaret Hospital Toronto, Ca

I tried to build R-1.1.1 on my HP-UX 10.20 box. For this I installed gcc 2.95.2 and built g77 (version 2.95.2 19991024 (release)). Making R runs to the point, where shared libs are linked. There it reports: make[4]: Entering directory `/disk4/hoffmann/Work/rwork/src/library/ctest/src' gcc -I../../../../include -I/usr/local/include -fPIC -g -O2 -c

Hi, R allows me to run a one sample Wilcoxon test like this: wilcox.test(c(1,3.5,2.1,4,1.5,5), mu=2, exact=TRUE) The function 'wilcoxsign_test' from the package 'coin' should (I suppose) be able to calculate exact p values even if there are ties in the ranks. However, I couldn't find information on how to run a one sample test using 'wilcoxsign_test' like in the

Hi everyone! I need to perform a Wilcoxon rank sum test, but I have some ties and the groups have different size also. When I deal with ties I use the wilcox.exact function, how can I solve the different size problem using this function? thanks net -- View this message in context: http://r.789695.n4.nabble.com/Wilcoxon-rank-sum-in-unbalanced-design-tp3447400p3447400.html Sent from the R help

Dear R-ers, I am currently running some Wilcoxon tests in R-64. How do I find the degrees of freedom in the output I am receiving? > wilcox.test(good$TRUE, good$x4a, paired=FALSE) Wilcoxon rank sum test with continuity correction data: good$TRUE and good$x4a W = 2455, p-value < 2.2e-16 alternative hypothesis: true location shift is not equal to 0 Thank you, Stephen.

Hello All, I am writing to ask your opinion on how to interpret this case. I have two vectors "a" and "b" that I am trying to compare. The wilcoxon test is giving me a pvalue of 5.139217e-303 of a over b with the alternative "greater". Now if I make a summary on each of them I have the following > summary(a) Min. 1st Qu. Median Mean 3rd Qu.

Hallo, is there any extension of the pairwise Wilcoxon test to a dependent samples layout with replicates (or, in other terms, a one-way layout with blocking and replicates)? The Wilcoxon method with matched pairs works for the case of dependent samples with one observation per block, while the Mann-Whitney test works for independent samples, thus one single block and replicated observations. Is