[R] test wilcoxon sur R help!

arun smartpink111 at yahoo.com
Fri Oct 25 06:38:19 CEST 2013


It looks much better than mine.


with p value adjustment:
p.adjust(fun1(dat1), method = "holm", n = 153)
#########

dat1$id <- 1:10
library(reshape2)
dat2 <- melt(dat1,id.var="id")
with(dat2,pairwise.wilcox.test(value,variable))
 with(dat2,pairwise.wilcox.test(value,variable,p.adj="none")) 


A.K.




On Friday, October 25, 2013 12:05 AM, vikram ranga <babuawara at gmail.com> wrote:
Hi,
Check out this function:-
pairwise.wilcox.test {package=stats}.

example(pairwise.wilcox.test)


On Fri, Oct 25, 2013 at 2:15 AM, Rui Barradas <ruipbarradas at sapo.pt> wrote:
> Hello,
>
> There's a bug in your function, it should be 'dat', not 'dat1'. In the line
> marked, below.
>
> fun1 <- function(dat){
>         mat1 <- combn(colnames(dat),2)  # Here, 'dat' not 'dat1'
>
>         res <- sapply(seq_len(ncol(mat1)),function(i) {x1<- dat[,mat1[,i]];
> wilcox.test(x1[,1],x1[,2])$p.value})
>         names(res) <- apply(mat1,2,paste,collapse="_")
>         res
> }
>
>
> Hope this helps,
>
> Rui Barradas
>
> Em 24-10-2013 20:16, arun escreveu:
>
>> Hi,
>> Try:
>> fun1 <- function(dat){
>> mat1 <- combn(colnames(dat1),2)
>>   res <- sapply(seq_len(ncol(mat1)),function(i) {x1<- dat[,mat1[,i]];
>> wilcox.test(x1[,1],x1[,2])$p.value})
>> names(res) <- apply(mat1,2,paste,collapse="_")
>> res
>> }
>>
>> set.seed(432)
>> dat1 <- as.data.frame(matrix(sample(18*10,18*10,replace=FALSE),ncol=18))
>>
>>    fun1(dat1) #gives the p-value for each pair of columns
>>
>>
>>
>>
>> Hi,
>>
>> I want to make a wilcoxon test, i have 18 columns each column
>> corresponds to a different sample and i want to compare one to each
>> other with a wilcoxon test in one step this is possible ? or do i
>> compare two by tow?
>>
>> Does it exist a code for automation this test? like this i dont have to
>> type the code for each couple.
>>
>> thanks!
>> denisse
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.

>>
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.




More information about the R-help mailing list