Question

今天第一篇文章。

目前，我正在从事研究，我必须分析心率数据。现在，我已经完成了大部分代码，我发现我正在使用的算法（RHRV）利用心跳的出现以及相匹配的秒数。

我的简单txt文件以秒为单位，因此看起来像这样（1列，433行）： 0.905 1.943 2.941 4.005 4.998 5.921

这些是心跳，第二次发生。

我现在的问题是：如何删除不可能发生的情况，例如： 2次事件/心跳的发生速度超过0.2秒（200毫秒）

是否有一种简单的方法或功能可以让我更改数据？

期待您的回答！

干杯

桑德

Answer 1

# create some random values
set.seed(3)
heartBeats <- sort(runif(10))

# which of them are reasonable?
keepThese <- heartBeats - lag(heartBeats) > 0.02
keepThese[is.na(keepThese)] <- TRUE # remove NAs

# and now subset the original vector
heartBeats[keepThese]

Answer 2

我想不出一种矢量化的方法来做到这一点，但这可行。请注意，这对于长向量来说将是昂贵的，因为表达式y <- c(y, x[i])将继续重建向量y。

x <- c(171.1156, 171.1345, 171.256, 171.309, 171.465, 171.4986, 171.986, 172.563)
m <- 1
y <- numeric(0)
y[m] <- x[1]
for (i in 2:NROW(x)) {
    if (y[m] + 0.2 < x[i]) {
        m <- m + 1
        y <- c(y, x[i])
    }
}
y
# output
[1] 171.1156 171.4650 171.9860 172.5630

我遍历x的元素，每次遇到与y中最新接受的值相距0.2以上的值时，我用{{ 1}}。

当变量靠近时如何消除变量（即小于200毫秒->删除行）

2 个答案: