当我看到这一点时,我试图在数据框中的两列之间执行 Kruskal-Wallis 测试,如果我将一列与其自身进行比较,我不会得到 1 的 p 值:
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
谁能告诉我我做错了什么?