我正在尝试重新编码每 2 小时收集的一些数据,以便找到每个 ID 的起点(即当 obs 不等于零时,即该时间点有数据(,将其称为时间 0,然后对于每个后续时间点称为 2、4、6 等。
例如
ID <- c("f1", "f1", "f1", "f1", "f2", "f2", "f2", "f2", "f3", "f3", "f3", "f3")
time <- rep(c(66, 68, 70, 72), 3)
obs <- c(1, 3, 5, 6, 0, 0, 3, 4, 0, 1, 3, 3)
new.time <- c(0, 2, 4, 6, NA, NA, 0, 2, NA, 0, 2, 4)
data <- as.data.frame(cbind(ID, time, obs, new.time))
希望数据框有效
我有ID,时间和obs,但我想创建"新时间" - 任何帮助,特别是dplyr解决方案
1(我们在最后的注释中将data
定义为data.frame而不是矩阵,然后使用ave
设置new.time
:
不使用包装。
make_no <- function(obs) c(rep(NA, sum(obs == 0)), seq(0, length = sum(obs != 0), by = 2))
transform(data, new.time = ave(obs, ID, FUN = make_no))
给:
ID time obs new.time
1 f1 66 1 0
2 f1 68 3 2
3 f1 70 5 4
4 f1 72 6 6
5 f2 66 0 NA
6 f2 68 0 NA
7 f2 70 3 0
8 f2 72 4 2
9 f3 66 0 NA
10 f3 68 1 0
11 f3 70 3 2
12 f3 72 3 4
2( 或使用 DPLYR:
data %>%
group_by(ID) %>%
mutate(new.time = make_no(obs)) %>%
ungroup
注意
ID <- c("f1", "f1", "f1", "f1", "f2", "f2", "f2", "f2", "f3", "f3", "f3", "f3")
time <- rep(c(66, 68, 70, 72), 3)
obs <- c(1, 3, 5, 6, 0, 0, 3, 4, 0, 1, 3, 3)
data <- data.frame(ID, time, obs)
我们可以创建一个自定义函数并按组应用它,即
f1 <- function(x) {
x1 <- length(x[x != 0])
i1 <- seq(0, length.out = x1, by = 2)
i2 <- c(rep(NA, (length(x) - x1)),i1)
return(i2)
}
#Using `dplyr` to apply it,
library(dplyr)
df %>%
group_by(ID) %>%
mutate(new = f1(obs))
这给了,
# A tibble: 12 x 4 # Groups: ID [3] ID time obs new <fct> <fct> <fct> <dbl> 1 f1 66 1 0 2 f1 68 3 2 3 f1 70 5 4 4 f1 72 6 6 5 f2 66 0 NA 6 f2 68 0 NA 7 f2 70 3 0 8 f2 72 4 2 9 f3 66 0 NA 10 f3 68 1 0 11 f3 70 3 2 12 f3 72 3 4