我试图在数据框中的两列之间执行Kruskal-Wallis测试时,如果我将列与自身进行比较,我的p值不是1:
controlvcontrol<- kruskal.test(neurite.5RGD[,1], neurite.5RGD[,1], na.rm=TRUE)
controlvcontrol$p.value
[1] 7.967453e-99
如果我比较两个包含相同数据的不同列,则会发生同样的事情:
(1:length(neurite.5RGD[,1]))[which((neurite.5RGD[,1]-neurite.5RGD[,3])!=0)]
integer(0)
> controlvBDNF<- kruskal.test(neurite.5RGD[,1], neurite.5RGD[,3], na.rm=T)
> controlvBDNF$p.value
[1] 7.967453e-99
我认为这是由于列中的NA值所致,所以我尝试了一个简单的例子:
> x <- c(1, 2, 3, 4,5, NA, NA)
> y <- c(1, 2, 3, 4,5, NA, NA)
> krus <- kruskal.test(x,y, na.rm = TRUE)
> krus$p.value
[1] 0.4060058
> ttest <- t.test(x,y, na.rm=T)
> ttest$p.value
[1] 1
谁能告诉我我做错了什么?