时间戳

首页 > 解决方案 > 基于R中的标识符将数据提取到一系列新数据帧中

r - 基于R中的标识符将数据提取到一系列新数据帧中

问题描述

我正在尝试做一些我确信很容易的事情，但我什至不知道要正确搜索什么才能找到灵魂。

我有很长的时间序列数据。一列表示状态，其值为 0、1、2 或 -1。我想提取并处理 state = 2 的数据。提取所有这些数据没有问题；我想要做的是分别处理 state=2 的每个时期。也就是说，我想通过数据，并且 state=2，做一些分析（例如，与此相关的其他变量的初始时期线性拟合时间段），然后转到下一个 state=2 的实例并重复此分析。或者，将每个时期的数据提取到自己的数据帧中并进行分析（但这会创建数百个小数据帧）。

我什至无法弄清楚如何识别起点（即 i=2 和 i-1=1 的位置）。

任何指向我应该查看的命令或包的指针都将不胜感激。

标签： r

解决方案

您可以将apply函数与MARGIN = 1（按行）一起使用，然后您将能够按索引访问每个结果，因为应用函数的值是一个列表。请看下图：

# Simulation of the data ----
library(lubridate)

set.seed(123)
n <- 20
df <- data.frame(
  id = factor(sample(-1:2, 20, replace = TRUE)),
  t_stamp = sort(sample(seq(dmy_hm("01-01-2018 00:00"), dmy_hm("01-01-2019 00:00"), by = 10), 20)),
  matrix(rnorm(100 * n), nrow = n, dimnames = list(NULL, c(paste0("X", 1:50), paste0("Y", 1:50)))))

# Calculation & output -----
# filter out all id == 2
df_m <- df[df$id == 2, ]

# apply to all the data linear regression ----
df_res <- apply(df_m, 1, function(xs) {
  rw <- as.numeric(xs[-(1:2)])
  x <- rw[1:50]
  y <- rw[51:100]
  smry <- summary(lm(y ~ x))
  list(xs[2], smry)

})

# access to the second incident
df_res[[2]]
# [[1]]
# t_stamp 
# "2018-03-26 13:01:40" 
# 
# [[2]]
# 
# Call:
#  lm(formula = y ~ x)
#
# Residuals:
#   Min       1Q   Median       3Q      Max 
# -2.44448 -0.79877  0.09427  0.88524  3.11190 
# 
# Coefficients:
#   Estimate Std. Error t value Pr(>|t|)
# (Intercept)   0.1731     0.1592   1.088    0.282
# x             0.1595     0.1544   1.033    0.307
# 
# Residual standard error: 1.125 on 48 degrees of freedom
# Multiple R-squared:  0.02173, Adjusted R-squared:  0.001349 
# F-statistic: 1.066 on 1 and 48 DF,  p-value: 0.307

# check the second incident
df_m[2, ]
# id             t_stamp         X1       X2        X3        X4   ...        X50        Y1         Y2       Y3        Y4
# 4  2 2018-03-26 13:01:40 -0.7288912 2.168956 -1.018575 0.6443765 ... -0.1321751 0.8983962 -0.2608322 1.036548 -1.691862

推荐阅读