3

我想在 R 中创建从 2004 年 1 月 1 日到 2010 年 12 月 31 日的每日死亡率数据的时间序列。我现在拥有的原始数据（.csv 文件）具有列日-月-年和每一行是一个死亡案例。因此，如果某一天的死亡率例如等于四，则该日期有四行。如果在特定日期没有报告死亡病例，则该日期在数据集中被省略。

我需要的是一个包含 2557 行的时间序列（从 01/01/2004 到 31/12/2010），其中列出了每天的死亡病例总数。如果某天没有死亡病例，我仍然需要将那一天放在列表中，并为其分配一个“0”。

有谁知道如何做到这一点？

谢谢，戈西亚

原始数据示例：

day month   year
1   1   2004
3   1   2004
3   1   2004
3   1   2004
6   1   2004
7   1   2004

我需要的：

day month   year    deaths
1   1   2004    1
2   1   2004    0
3   1   2004    3
4   1   2004    0
5   1   2004    0
6   1   2004    1
4

1 回答 1

3 
df <- read.table(text="day month   year
1   1   2004
3   1   2004
3   1   2004
3   1   2004
6   1   2004
7   1   2004",header=TRUE)

#transform to dates
dates <- as.Date(with(df,paste(year,month,day,sep="-")))

#contingency table
tab <- as.data.frame(table(dates))
names(tab)[2] <- "deaths"
tab$dates <- as.Date(tab$dates)

#sequence of dates
res <- data.frame(dates=seq(from=min(dates),to=max(dates),by="1 day"))
#merge
res <- merge(res,tab,by="dates",all.x=TRUE)
res[is.na(res$deaths),"deaths"] <- 0
res
#       dates deaths
#1 2004-01-01      1
#2 2004-01-02      0
#3 2004-01-03      3
#4 2004-01-04      0
#5 2004-01-05      0
#6 2004-01-06      1
#7 2004-01-07      1
于 2013-03-07T16:25:02.747 回答