我的数组是这样的:
Slide Index A B C DoseGroup
482 778 l 0 0 2 13Gy_p_75_42wk
483 778 r 0 0 2 13Gy_p_75_42wk
484 779 l 0 0 2 13Gy_p_75_42wk
485 779 r 0 0 2 13Gy_p_75_42wk
486 4700 l 2 2 2 14.25Gy_C_50pl_42wk
487 4700 r 0 0 1 14.25Gy_C_50pl_42wk
488 4701 l 0 0 1 14.25Gy_C_50pl_42wk
我想使用DoseGroup
列的条目,以便能够在其他列中选择相应的条目。我希望能够告诉R,例如,"使用列C
在13Gy_p_75_42wk
和14.25Gy_C_50pl_42wk
数据集之间执行wilcox.test
。"
我怎么用R做这个?是否有某种方法可以选择所有具有14.25Gy_C_50pl_42wk
条目的列?
我修改了你的数据,在DoseGroup中添加了第三个级别,使其更真实。
txt <- "Slide Index A B C DoseGroup
778 l 0 0 2 13Gy_p_75_42wk
778 r 0 0 2 13Gy_p_75_42wk
779 l 0 0 2 13Gy_p_75_42wk
779 r 0 0 2 13Gy_p_75_42wk
4700 l 2 2 2 14.25Gy_C_50pl_42wk
4700 r 0 0 1 14.25Gy_C_50pl_42wk
4701 l 0 0 1 14.25Gy_C_50pl_42wk
4702 l 0 0 10 15Gy_C_50pl_42wk"
dat <- read.table(text = txt, header = TRUE)
wilcox.test(C ~ DoseGroup, data = dat,
subset = DoseGroup %in% c("13Gy_p_75_42wk", "14.25Gy_C_50pl_42wk"))
## Wilcoxon rank sum test with continuity correction
## data: C by DoseGroup
## W = 10, p-value = 0.1175
## alternative hypothesis: true location shift is not equal to 0
要选择数据,可以使用以下两个命令之一:
dat[dat$DoseGroup == "14.25Gy_C_50pl_42wk", ]
subset(dat, DoseGroup == "14.25Gy_C_50pl_42wk")
这些命令是R
的基础,如果你读过R
的任何介绍,你也能做到。
我强烈建议你们这么做,如果你们想真正享受R