这是一个难以表述的问题。以下是我想要做的事情示例,以及我现在所拥有的示例:
set.seed(0)
dt <- data.table(dr1.d=rnorm(5), dr1.p=abs(rnorm(5, sd=0.08)),
dr2.d=rnorm(5), dr2.p=abs(rnorm(5, sd=0.08)),
dr3.d=rnorm(5), dr3.p=abs(rnorm(5, sd=0.08)),
dr4.d=rnorm(5), dr4.p=abs(rnorm(5, sd=0.08)),
sym = paste("sym", c(1,1,1,2,2)))
dt
dr1.d dr1.p dr2.d dr2.p dr3.d dr3.p dr4.d dr4.p sym
1: 1.2629543 0.1231960034 0.7635935 0.03292087 -0.22426789 0.040288638 -0.2357066 0.09215294 sym 1
2: -0.3262334 0.0742853628 -0.7990092 0.02017788 0.37739565 0.086861549 -0.5428883 0.07937283 sym 1
3: 1.3297993 0.0235776357 -1.1476570 0.07135369 0.13333636 0.055276307 -0.4333103 0.03436105 sym 1
4: 1.2724293 0.0004613738 -0.2894616 0.03485466 0.80418951 0.102767948 -0.6494716 0.09906433 sym 2
5: 0.4146414 0.1923722711 -0.2992151 0.09900307 -0.05710677 0.003738094 0.7267507 0.02234770 sym 2
对于所有共享药物的列对(例如“dr1”),我想按“sym”分组行,然后在每个组中选择具有最小p值(以“.p”结尾)的行。以上数据表的最终结果如下:
dr1.d dr1.p dr2.d dr2.p dr3.d dr3.p dr4.d dr4.p sym
1: 1.3297993 0.0235776357 -0.7990092 0.02017788 -0.22426789 0.040288638 -0.4333103 0.03436105 sym 1
2: 1.2724293 0.0004613738 -0.2894616 0.03485466 -0.05710677 0.003738094 0.7267507 0.02234770 sym 2
我尝试使用.SD和lapply来完成这个任务,但我无法理解它。谢谢!