我有两个数据帧:
df1<- data.frame(POS=c(1, 2, 3, 4, 5),
A=c(0.2,0,0,0.8,0),
G=c(0,0.3,0,0.2,0.5),
T=c(0.4,0.7,1,0,0.5),
C=c(0.4,0,0,0,0))
df2<- data.frame(POS=c(1, 3, 4, 5, 6),
A=c(0.3,0,0.6,0,0),
G=c(0,0,0.4,0.7,1),
T=c(0.2,1,0,0.3,0),
C=c(0.5,0,0,0,0))
我想得到每个POS的总平方差
我需要的第一件事是两个相同维度的矩阵。如何为用0填充的缺失位置添加行?
一旦我有了两个矩阵,我将执行(df1-df2)^2,然后是rowsums方法。
总结一下:如何为每个data.frame?中丢失的行添加一个充满零的新行
这样地:
POS A G T C
1 0.2 0.0 0.4 0.4
2 0.0 0.3 0.7 0.0
3 0.0 0.0 1.0 0.0
4 0.8 0.2 0.0 0.0
5 0.0 0.5 0.5 0.0
6 0.0 0.0 0.0 0.0
POS A G T C
1 0.3 0.0 0.2 0.5
2 0.0 0.0 0.0 0.0
3 0.0 0.0 1.0 0.0
4 0.6 0.4 0.0 0.0
5 0.0 0.7 0.3 0.0
6 0.0 1.0 0.0 0.0
正如李哲远在评论中所说,您可以使用每个df的merge()
,将另一个df的POS列放在第一列(POS)上。这会添加用NA
填充的缺失行,因此您必须用零填充它们:
df1 <- merge(df1, df2$POS, all = TRUE, by=1)
df1[is.na(df1)] <- 0
df1
POS A G T C
1 1 0.2 0.0 0.4 0.4
2 2 0.0 0.3 0.7 0.0
3 3 0.0 0.0 1.0 0.0
4 4 0.8 0.2 0.0 0.0
5 5 0.0 0.5 0.5 0.0
6 6 0.0 0.0 0.0 0.0
df2 <- merge(df2, df1$POS, all = TRUE, by=1)
df2[is.na(df2)] <- 0
df2
POS A G T C
1 1 0.3 0.0 0.2 0.5
2 2 0.0 0.0 0.0 0.0
3 3 0.0 0.0 1.0 0.0
4 4 0.6 0.4 0.0 0.0
5 5 0.0 0.7 0.3 0.0
6 6 0.0 1.0 0.0 0.0