在R中提取超过月平均值的每小时值(及其索引)

# 加载数据集
data <- read.csv("data.csv")

# 转换日期时间列格式
data$datetime <- as.POSIXct(data$datetime, format = "%Y-%m-%d %H:%M:%S")

# 转换为时间序列对象
ts_data <- xts(data$value, order.by = data$datetime)

# 计算每小时平均值
hourly_avg <- aggregate(ts_data, FUN = mean, by = as.POSIXlt(index(ts_data))$hour)

# 筛选超过月平均值的每小时值及其索引
hourly_values <- hourly_avg[hourly_avg > mean(hourly_avg)]
hourly_indices <- index(hourly_values)

# 输出结果
hourly_values
hourly_indices