使用tsoutliers识别异常值后如何标记异常记录并替换数值

实现步骤与代码示例 #

先将tsoutliers()的运行结果保存为独立对象，避免重复调用函数产生不必要的计算，再基于返回的$index和$replacements字段完成标记和替换操作即可：

1. 存储异常值检测结果

# outlier_res$index 为异常值对应的行位置，outlier_res$replacements 为对应推荐替换值
outlier_res <- tsoutliers(df$Y)

2. 新建异常值标记列

# 初始化标记列为0（代表非异常）
df$outlier <- 0
# 将异常值对应行的标记值设为1
df$outlier[outlier_res$index] <- 1

3. 替换异常值为推荐值

# 按位置一一对应替换异常值
df$Y[outlier_res$index] <- outlier_res$replacements

原代码错误原因 #

which()函数的作用是返回满足条件的元素位置索引，不能直接用于向量级的条件赋值；且$index字段返回的是异常值的位置数值向量，无异常时为空向量，不会返回NA，因此is.na()的判断逻辑也不成立。

内容的提问来源于stack exchange，提问作者Francesca