similar to: wilcox.test p-value = 0

i'm somehow embarrassed to even ask this, but is there any built-in method for doing this: my_list <- list() my_list[[1]] <- matrix(1:20, ncol = 5) my_list[[2]] <- matrix(20:1, ncol = 5) now, knowing that these matrices are identical in dimension, i'd like to unfold the list to a 2x4x5 (or some other permutation of the dim sizes) array. i know i can initialize the array, then

hi all - i've seen versions of this question before, but none seem to get directly at my solving my (probably very simple) issue: i simply want to annotate the tick marks on an axis with (superscripted) 10^x notation, and tried this: axis(1, at = axTicks(1), as.expression(substitute(10^foo, list(foo = axTicks(1)))) thinking the as.expression/substitute would create the appropriate

Hi folks, sorry for this beginner question but what means a p-value = NA on a menn-whitney test? v1 <- c(0.022, 0.022, 0.022, 0.022, 0.022, 0.022) v2 <- c(0.022, 0.022, 0.022, 0.022, 0.022) wilcox.test(v1, v2, PAIRED=FALSE) W = 15, p-value = NA I know that there's no statistical difference between v1 and v2, so why my p-value is not ONE ?? Can I consider p-value= NA as ONE ?? what i

i've written a function to coerce a matrix (e.g. from numeric to logical), but i'd like to know if someone has a more elegant method for this: > m <- matrix(c(0, 1, 1, 0), ncol = 2) > m <- as.logical(m) > m [1] FALSE TRUE TRUE FALSE i'd like 'm' to still be a matrix with the original dimensions. in my function to do this, i coerce 'm' to a logical,

hi all - i'm having some difficulty figuring out how to convert between "user units" (which i can't find a definition for in the plotrix package) and either (a) device units (e.g. inches with PDFs) or (b) user coordinates along any particular axis. as an example, suppose i set up a PDF device with inches, the device has both outer and inner magins, and the plot region has

I have a stupid question on how to get the real p-values for wilcox.test and correlation. the minmun can be reached is 2.2E-16 using the R version 2.6.2. I do not think it is the R version causing this but other issues. Any help is highly appreciated. Charles -- View this message in context: http://www.nabble.com/minmun-p-value-for-wilcox.test-and-correlation-tp23637568p23637568.html Sent

does anybody know why there are the two warnings in the example above? Regards Knut > day_4 [1] 540 1 1 1 1 1 1 300 720 480 > day_1 [1] 438 343 1 475 1 562 500 435 1045 890 > is.vector (day_1) [1] TRUE > is.vector (day_4) [1] TRUE > wilcox.test(day_4 ,day_1,paired=TRUE,alternative="two.sided",exact=TRUE,conf.int=TRUE) Wilcoxon

hi all - i'm trying to 'R CMD build' a package, but i have what appears to be a bootstrapping problem: i've included a vignette in my package, with R code interwoven (and built using Sweave), but in this documentation i have a code line: > library(MyPackage) now, when trying to build a .tar.gz install-able version of my package in a clean state, i remove the original

hi all - i'm just wondering what sort of code people write to essentially performa an aggregate call, but with different functions being applied to the various columns. for example, if i have a data frame x and would like to marginalize by a factor f for the rows, but apply mean() to col1 and median() to col2. if i wanted to apply mean() to both columns, i would call: aggregate(x, list(f),

hi, i'm trying to prepend some plain (i.e. unevalutated) text to a typographically evaluated R expression in a legend(...) call. here's a working legend(...) call that is close to what i'd like (where x and y are returned from an lm(...) call): legend("topleft", legend = c(bquote(R^2 == .(summary(x)$r.squared)), bquote(R^2 == .(summary(y)$r.squared))), fill =

hi all - i'm working on an R package that makes use of my own shared library written in C. but i also am making use of another C-written library. (my package is for facilitating biological namespace translations via online (i.e. up-to-date) biological databases.) problem is, the library i'm using is not a standard library (i.e. i doubt it will be installed on most users' machines). i

Hi all, I have a question about finalizers... I have a package that manages state for a few connections, and I'd like to ensure that these connections are 'cleanly' closed upon either (i) R quitting or (ii) an unloading of the package. So, in a pared-down example package with a single R file, it looks something like: ##### BEGIN PACKAGE CODE ##### .CONNS <- new.env(parent =

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 the summed up ranks of either the positive or negative differences,

hi all - i have a source package i'm writing that i'd like to be able to install with a special library that my R src files rely on. to be more precise, i have a normal R package directory structure (i.e. src/ R/ man/ DESCRIPTION NAMESPACE etc.). i also have another directory here called depPkg, and it has it's own configure file for the canonical './configure && make

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), c(0,9,8,7), paired=TRUE)

hi, i am looking for a way to correct for ties in the wilcoxon signed rank test -> e.g. wilcox.test(x,mu=.5) one way i have heard of is to look up the p value in a table that has been produced by Buck (1975). obviously i need to know the T-value to do that -> how do i get the T-value from the wilcox.test() function. is there any other (already implemented) way to correct for ties in

The point estimates produced by wilcox.test are perverse (not wrong, just brain damaged). The Hodges-Lehmann estimator that goes with the signed rank test is the median of the Walsh averages. The Hodges-Lehmann estimator that goes with the rank sum test is the median of the pairwise differences. wilcox.test agrees except that it uses the following very peculiar definition of "sample

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, alternative = "g")? or wilcox.test(test,

Good day all, My first posting to this list. It looked like the best place to post this question. When running the wilcox.test(), I noticed that the output values change if you change the ordering of the levels (example below which includes a t.test for comparison). I think this has something to do with the change in ranking order, but this doesn't make much sense as I would expect the

Hi all --- this is less a specific question and more general regarding S3 classes. I've noticed that quite a few very common default implementations of generic functions (e.g. `unique`, `[`, `as.data.frame`) strip away class information. In some cases, it appears conditionals have been created to re-assign the class, but only for a few special types. For example, in `unique.default`, if the