以下数据代表四种药物对某症状的影响评分，做kruska-Walis检验，得到P<0.05，请问应如何进一步作两两比较

x <- c(2, 3,3, 3, 3,0)
y <- c(1, 3, 1, 2, 2,3)     
z <- c(3, 1, 2, 1,3, 3)
k<-c(1,0,0,0,0,-1)
kruskal.test(list(x, y, z,k))

另外一个问题是如果希望看到kruskal.test的程序编写过程，有何办法？我试着在提示符下输入kruskal.test，只得到以下回复

function (x, ...)
UseMethod("kruskal.test")
<environment: namespace:stats>

谢谢！