我仍然是傻瓜新手,请原谅我,如果这是一个简单的问题,但我无法在stackoverflow上找到答案。我想为每个主要的美国假期创建一组指标变量,只需将一个函数应用到我的日期字段,该函数可以检测哪些天是假日,然后我可以将Model.matrix等转换为一组指标变量。
例如,我有从2012年1月1日到2013年9月15日的每日数据,我想为复活节创建一个变量指标。
我目前正在使用timeDate
软件包将一年时间传递给他们的函数Easter()
以查找日期。然后,我将日期键入以下代码以创建指标变量。
Easter(2012)
EasterInd2012<-as.numeric(DATASET$Date=="2012-04-08")
答案 0 :(得分:10)
获得一般假期指标变量的最简单方法是创建您感兴趣的所有假期的向量,然后在数据框中匹配这些日期。这样的事情应该有效:
library(timeDate)
# Sample data
Date <- seq(as.Date("2012-01-01"), as.Date("2013-09-15"), by="1 day")
DATASET <- data.frame(rnorm(624), Date)
# Vector of holidays
holidays <- c(as.Date("2012-01-01"),
as.Date(Easter(2013)),
as.Date("2012-12-25"),
as.Date("2012-12-31"))
# 1 if holiday, 0 if not. Could also be a factor, like c("Yes", "No")
DATASET$holiday <- ifelse(DATASET$Date %in% holidays, 1, 0)
您可以手动输入日期,也可以使用timeDate
的一些内置假日功能(listHolidays()
功能显示所有这些功能)。所以你也可以这样构建holidays
:
holidays <- c(as.Date("2012-01-01"),
as.Date(Easter(2013)),
as.Date(USLaborDay(2012)),
as.Date(USThanksgivingDay(2012)),
as.Date(USMemorialDay(2012)),
as.Date("2012-12-25"),
as.Date("2012-12-31"))
要获得每个假期的具体指标,您需要一次执行一个:
EasterInd2012 <- ifelse(DATASET$Date==as.Date(Easter(2012)), 1, 0)
LaborDay2012 <- ifelse(DATASET$Date==as.Date(UsLaborDay(2012)), 1, 0)
# etc.