我的数据帧很长,如下所示:
data = data.frame(
id = rep(c(1, 2, 3, 4, 5, 6, 7), each = 5),
pt = rep(c("A", "A", "A", "B", "B", "B", "B"), each = 5),
value = rnorm(35, 10, 1)
)
id
是样本ID,而pt
是患者ID。不需要的信息,但可能有助于了解目的。一个pt
具有多个id
,因为它们是重复项。
我想对重复样本进行成对散点图。这意味着对于单个pt
,我想要id
的所有散点图组合。
所需的输出将是带有A_1 vs A_2, A_2 vs A_3, A_1 vs A_3, B_4 vs B_5, B_4 vs B_6, B_4 vs B_7, B_5 vs B_6, B_5 vs B_7, B_6 vs B_7
的散点图的构面网格,格式为pt_id
。
这怎么办?如果我将其从长转换为宽,然后运行对,我将获得所有组合,但是我只对来自同一pt
的组合感兴趣,并且想要一个facet_wrap对象。
谢谢, 杰克