0

我的数组如下所示：

     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，例如，“wilcox.test在 the13Gy_p_75_42wk和14.25Gy_C_50pl_42wk使用 column 的数据集之间做一个C”。

我怎么能用 R 做到这一点？是否有某种方法可以选择所有具有条目的列14.25Gy_C_50pl_42wk？

4

1 回答 1

1

我修改了您的数据以在 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。

于 2013-04-25T13:18:33.330 回答