我只是试一试,因为我知道这里有一些聪明人可能有一个r码。 我不能自己编码。
所以我得到了一个数据集,其中包含2000-01和2008-12之间的名称和年月。看起来像这样:
Name Date
A 2000-01
A 2000-02
A ...
A 2008-12
A 2000-01
B 2000-01
B ...
B 2008-12
C and so on..
对于我的键列中的每个名称,每年都会有一个值。这是我能要求的最好的。不幸的是,有些年份我的关键专栏没有价值。 在我的数据集中进一步查看名称A:
因此,如果我在2000 - 2008年之间每年都没有1个观察结果,那么我希望根据下一年的月份,从年份和月份中得到一个没有我的关键列值的行。观察。 在这个例子中:
2003-02有我的keycolumn的值,而2002-02没有,我想从2002-02和Name A的日期取回行。 简而言之:"根据明年的关键列保留前一年的行数"
有一些简单的方法来编码吗? 谢谢:))
答案 0 :(得分:1)
没有直接简单的方法来编码您所描述的内容,但是将问题分解为更容易的部分当然是可能的。问题的核心部分如下。给定具有非NA值的行的数据帧,例如
year month
1 2002 12
2 2005 11
3 2006 01
4 2008 07
对于每一行,检查数据帧以查看前一年是否存在;如果是,则返回行,如果不是,则返回前一年和同一月的附加行。这是一个可能看起来像
的函数check_ym <- function(y, m, dat) {
if ((y - 1) %in% dat$year) {
return(data.frame(Date = paste(y, m, sep = "-"), stringsAsFactors = FALSE))
} else {
return(data.frame(Date = paste(c(y - 1, y), c(m, m), sep = "-"), stringsAsFactors = FALSE))
}
}
现在,让我们制作一些假数据。
library(dplyr)
library(tidyr)
library(purrr)
# Simulate data
set.seed(123)
x <- data.frame(Date = paste(sample(2000:2008, 4),
sprintf("%02d", sample(1:12, 4, replace = TRUE)),
sep = "-"),
KeyColumn = floor(runif(4, 1, 10)))
d <- data.frame(Date = paste(rep(2000:2008, each = 12),
sprintf("%02d", rep(1:12, times = 9)),
sep = "-")) %>%
left_join(x)
识别非NA行:
dd <- d %>%
na.omit() %>%
separate(Date, into = c("year", "month")) %>%
mutate(year = as.numeric(year))
dd
# year month KeyColumn
# 1 2002 12 5
# 2 2005 11 5
# 3 2006 01 5
# 4 2008 07 9
然后,我们运行上面的函数,遍历year
和month
列。这给了我们
out <- map2_df(dd$year, dd$month, .f = check_ym, dat = dd)
out
# Date
# 1 2001-12
# 2 2002-12
# 3 2004-11
# 4 2005-11
# 5 2006-01
# 6 2007-07
# 7 2008-07
最后,我们将此与原始数据相结合:
inner_join(out, d)
# Joining, by = "Date"
# Date KeyColumn
# 1 2001-12 NA
# 2 2002-12 5
# 3 2004-11 NA
# 4 2005-11 5
# 5 2006-01 5
# 6 2007-07 NA
# 7 2008-07 9
这只适用于一个Name
。我们也可以为许多Name
执行此操作。首先创建一些假数据:
# Simulate data
set.seed(123)
d <- map_df(setNames(1:3, LETTERS[1:3]), function(...) {
x <- data.frame(Date = paste(sample(2000:2008, 4),
sprintf("%02d", sample(1:12, 4, replace = TRUE)),
sep = "-"),
KeyColumn = floor(runif(4, 1, 10)))
data.frame(Date = paste(rep(2000:2008, each = 12),
sprintf("%02d", rep(1:12, times = 9)),
sep = "-")) %>%
left_join(x)
}, .id = "Name")
dd <- d %>%
na.omit() %>%
separate(Date, into = c("year", "month")) %>%
mutate(year = as.numeric(year))
dd
# Name year month KeyColumn
# 1 A 2002 12 5
# 2 A 2005 11 5
# 3 A 2006 01 5
# 4 A 2008 07 9
# 5 B 2000 04 6
# 6 B 2004 01 7
# 7 B 2005 12 9
# 8 B 2006 03 9
# 9 B 2000 04 6
# 10 C 2003 12 1
# 11 C 2005 04 7
# 12 C 2006 11 5
# 13 C 2008 02 8
现在,使用split
将数据框拆分为三个数据框Name
;对于每个子数据帧,我们应用check_ym()
,然后我们将结果组合在一起并将其与原始数据相结合:
lapply(split(dd, dd$Name), function(dat) {
map2_df(dat$year, dat$month, .f = check_ym, dat = dat)
}) %>%
bind_rows(.id = "Name") %>%
inner_join(d)
# Joining, by = c("Name", "Date")
# Name Date KeyColumn
# 1 A 2001-12 NA
# 2 A 2002-12 5
# 3 A 2004-11 NA
# 4 A 2005-11 5
# 5 A 2006-01 5
# 6 A 2007-07 NA
# 7 A 2008-07 9
# 8 B 2000-04 6
# 9 B 2003-01 NA
# 10 B 2004-01 7
# 11 B 2005-12 9
# 12 B 2006-03 9
# 13 C 2002-12 NA
# 14 C 2003-12 1
# 15 C 2004-04 NA
# 16 C 2005-04 7
# 17 C 2006-11 5
# 18 C 2007-02 NA
# 19 C 2008-02 8