Kruskal Wallis检验的p值不正确;可能与NA值有关?

时间:2013-08-09 21:33:56

标签: r na p-value

我试图在数据框中的两列之间执行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

谁能告诉我我做错了什么?

0 个答案:

没有答案