Question

我有一个data.frame，其中包含多个列（即V1 ... Vn+1），其值为1或0，每列都是一个时间步长。

我想知道值为1的平均time（列数）。序列1 1 1 1 1 1的值为1。

目前我可以想到计算这个的方法是计算1s之间0的平均计数（+1），但它是有缺陷的。

例如，具有这些值1 0 0 1 0 1的行的结果为2.5（2 + 1 = 3; 3/2 = 1.5 ; 1.5 + 1 = 2.5）。

但是，如果序列以0开始或结束，则应在没有它们的情况下计算此结果的结果。例如，0 1 0 0 1 1将计算为1 0 0 1 1，结果为3。

有缺陷的，例如1 0 1 1 0 0计算为1 0 1 1，结果为2，但这不是理想的结果（1.5）

考虑到以零开头或结尾的问题，有没有办法按行计算值1之间的列数？

# example data.frame with desired result
df <- structure(list(Trial = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), Location = c(1L, 
1L, 1L, 1L, 2L, 2L, 2L, 2L), Position = c(1L, 2L, 3L, 4L, 1L, 
2L, 3L, 4L), V1 = c(1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L), V2 = c(1L, 
1L, 1L, 0L, 1L, 0L, 0L, 0L), V3 = c(1L, 1L, 1L, 0L, 1L, 0L, 0L, 
1L), V4 = c(1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L), V5 = c(1L, 0L, 0L, 
0L, 1L, 0L, 0L, 0L), V6 = c(1L, 1L, 1L, 0L, 1L, 1L, 0L, 0L), 
    Result = c(1, 3, 2, NA, 1, 2.5, 3, 1.5)), .Names = c("Trial", 
"Location", "Position", "V1", "V2", "V3", "V4", "V5", "V6", "Result"
), class = "data.frame", row.names = c(NA, -8L))

df1 <- df[,4:9]

#This code `apply(df1,1,function(x) which(rev(x)==1)[1])) calculates the number of columns back until a value of 1, or forward without `rev`. But this doesn't quite help with the flaw.

Answer 1

如果第一个和最后一个值之间的范围是k，并且该范围内的总数为n，那么平均差距为(k-1)/(n-1)。您可以使用以下方法计算：

apply(df1, 1, function(x) {
  w <- which(x == 1)
  if (length(w) <= 1) NA
  else diff(range(w)) / (length(w)-1)
})
# [1] 1.0 2.0 2.0  NA 1.0 2.5 3.0 1.5

按行，获取x

1 个答案: