在计算 R 中两个数据帧之间的相关性之前如何忽略双零

Question

我有两个具有相同列数（100：样本）和行数（9600：基因）的数据框。这两个数据帧是从两个不同的程序输出的，我想计算它们之间的相关性。

我的示例数据集：

df1 <-data.frame(Sample1 =c(0.52,2.5,8.3,10.5,5.3),Sample2=c(0,0,2,1,0), Sample3=c(0,12,13,14,0))
rownames(df1)<-c("KO1","KO2","KO3","KO4","KO5")
df2<- data.frame(Sample1=c(1,2,3,4,5),Sample2=c(0,0,8,9,0),Sample3=c(0,12,13,14,0))
rownames(df2)<-c("KO1","KO2","KO3","KO4","KO5")
df<-data.frame(df1,df2)

>df1
      Sample1 Sample2 Sample3
KO1    0.52       0       0
KO2    2.50       0      12
KO3    8.30       2      13
KO4   10.50       1      14
KO5    5.30       0       0

>df2
      Sample1 Sample2 Sample3
KO1       1       0       0
KO2       2       0      12
KO3       3       8      13
KO4       4       9      14
KO5       5       0       0

在计算相关性时，我想删除两个数据框中都为零的条目。例如，对于样本 1，应包含每一行，但对于样本 2，不应包含 KO1、KO2 和 KO5，同样对于样本 3 、 KO1 和 KO5 不应包括在内。在这里，我计算两个数据帧之间的列相关性。
我尝试了以下代码：

output_without_zero<- with(subset(df, !(df1 == 0 & df2 == 0)), cor(df1,df2,method = "spearman"))
output_with_zero<- cor(df1,df2,method = "spearman")

我预计从相关性中删除零应该不同于包含它们。但我得到了它们两个相同的相关矩阵。 如何获得欲望输出？

提前谢谢您

Answer 1

将零替换为

NA

df1[df1 == 0] <- NA
df2[df2 == 0] <- NA

作者：

complete.obs

cor(df1, df2, method = "spearman",  use = "complete.obs")
        Sample1 Sample2 Sample3
Sample1       1       1       1
Sample2      -1      -1      -1
Sample3       1       1       1

作者：

pairwise.complete.obs

cor(df1, df2, method = "spearman",  use = "pairwise.complete.obs")
        Sample1 Sample2 Sample3
Sample1     0.7       1       1
Sample2    -1.0      -1      -1
Sample3     1.0       1       1

Answer 2

我认为你的意思是不应该包括KO1、KO2和KO5。

另外，我不确定当没有

subset(df, !(df1 == 0 & df2 == 0)

时

df

是否会起作用。

尝试-

inds <- rowSums(df1 == 0 & df2 == 0) == 0
#If you want to remove rows which has any occurrence of 0 in either dataframe
#inds <- rowSums(df1 == 0 | df2 == 0) == 0

df1_1 <- df1[inds, ]
df2_1 <- df2[inds, ]
output_with_zero<- cor(df1,df2,method = "spearman")
output_without_zero<- cor(df1_1,df2_1,method = "spearman")

Answer 3

我认为以下回答了您的问题，当两个值都为零时删除对。

# construct vectors
set.seed(1234)
A <- sample(0:10, 100, replace=T)
B <- sample(0:10, 100, replace=T)

# remove elements of A and B where both are 0
cor(A[A != 0 | B != 0], B[A !=0 | B != 0])

取自https://stackoverflow.com/a/37396482