如何用R中的样本ID创建样本类型值差异的数据框架？

Question

我想要为每个样本ID创建两种类型的样本之间的值差异的数据框架，问题是有些样本ID不具有这两种样本类型。我想我正在寻找一个for循环/条件语句，但是我没有成功地创建一个

我在下面创建了一些示例数据：

ID <- c(1,3,3,4,6,6,7,7,8,8)
Type <- c('soil', 'soil', 'water','soil','soil', 'water','soil', 'water','soil', 'water')
data1 <- as.data.frame(cbind(ID, Type))
data1$value <- runif(length(data1$ID), min=0, max=1)
head(data1)

结束结果将为每个具有两种样本类型的样本ID提供差异(水-土壤)的数据框架。只有一种数据类型的示例ID将被记录为NA。

Answer 1

您可以使用if/else在dplyr中尝试条件检查：

library(dplyr)

data1 %>%
  group_by(ID) %>%
  summarise(diff = if(n_distinct(Type) == 1) NA else 
                      value[Type == 'water'] - value[Type == 'soil'])

如果NA中有一个唯一值，则返回Type，否则返回两者之间的差异。

不使用if/else的诀窍是将差异子集，因此当没有值时，它会自动返回NA。

data1 %>%
  group_by(ID) %>%
  summarise(diff = (value[Type == 'water'] - value[Type == 'soil'])[1])