R用于循环计算行均值比

Question

我有一个像这样的数据框（df）：

Rif   A1  A2  A3  A4  A5  B1  B2  B3  B4  B5
a1   1   1   0   3   7   10  10  7   7   7 
b1   3   3   3   0   0   5   5   5   5   0
c1   1   3   0   0   2   3   3   3   3   3

我将需要计算Vs的列表，如下所示：

V1: mean(A5:A5)/mean(B5:B5)
V2: mean (A4:A5) / mean (B4:B5)
V3: mean (A3:A5) / mean (B3:B5)
....

我认为我应该使用for循环。 我已经尝试过下面编写的循环，但是它不起作用：

ref <- c(1, 2, 3, 4, 5)
 dat1 <- NULL

 for (i in 2:length(ref)) {
   for (j in 1:nrow(df)) {
     dat1[j] <- rowMeans(df[j,i-1:6])/rowMeans(df[j,i+6-1:11])
   }
   df[,paste0("V", i)] <- dat1
 }

我期望的输出是这样的：

Rif   A1  A2  A3  A4  A5  B1  B2  B3  B4  B5  V1  V2  V3  V4  V5
a1    1   1   0   3   7   10  10  7   7   7 
b1    3   3   3   0   0   5   5   5   5   0
c1    1   3   0   0   2   3   3   3   3   3

Answer 1

继续您的步骤，您可以尝试：

for (i in 1:5) {
  df[,paste0("V", i)] <- rowMeans(df[,(i+1):6,drop=FALSE])/rowMeans(df[,(i+6):11,drop=FALSE])
}
df
#  Rif A1 A2 A3 A4 A5 B1 B2 B3 B4 B5        V1        V2        V3        V4        V5
#1  a1  1  1  0  3  7 10 10  7  7  7 0.2926829 0.3548387 0.4761905 0.7142857 1.0000000
#2  b1  3  3  3  0  0  5  5  5  5  0 0.4500000 0.4000000 0.3000000 0.0000000       NaN
#3  c1  1  3  0  0  2  3  3  3  3  3 0.4000000 0.4166667 0.2222222 0.3333333 0.6666667

您可以跳过j循环（对每行进行子设置），然后使用rowMeans 。 另一方面，当您仅对一列进行子集化时，需要drop=FALSE 。 索引当前是硬编码的，但是在数据结构保持不变时将起作用。

或者您使用sapply而不是for-loop例如：

ref <- 1:5
names(ref)  <- paste0("V", ref)
df <- cbind(df, sapply(ref, function(i) {rowMeans(df[,(i+1):6,drop=FALSE])/rowMeans(df[,(i+6):11,drop=FALSE])}))
df
#  Rif A1 A2 A3 A4 A5 B1 B2 B3 B4 B5        V1        V2        V3        V4        V5
#1  a1  1  1  0  3  7 10 10  7  7  7 0.2926829 0.3548387 0.4761905 0.7142857 1.0000000
#2  b1  3  3  3  0  0  5  5  5  5  0 0.4500000 0.4000000 0.3000000 0.0000000       NaN
#3  c1  1  3  0  0  2  3  3  3  3  3 0.4000000 0.4166667 0.2222222 0.3333333 0.6666667