我有一个时间索引矩阵(xts对象)，我只希望每年11月的第四个星期三。
require(quantmod)
getSymbols("^GSPC", from="1900-01-01") #returns GSPC
GSPC$WED <- weekdays(time(GSPC)) == "Wednesday"
GSPC$NOV <- months(time(GSPC)) == "November"
G <- GSPC[GSPC$WED==1 & GSPC$NOV==1]
那是我进入R所能达到的。为了解决我的问题，我猛烈地抨击了bash。
write.zoo(G, "wen_in_nov")
我做了以下修改:
cat wen_in_nov |
grep -v IND |
cut -c 1-10 |
sed 's/-/ /g' |
awk '{if($3 >= 22 && $3 < 29) print $1, $2, $3, "winner"}' |
sed 's/ /-/g' > fourth_wen
fourth_wen文件需要将-与字符串“winner”分开，所以我只是在vi中做到了。导入回R:fourth_wen <- read.zoo("fourth_wen", format="%Y-%m-%d")
这实际上是自1950年以来11月的第四个星期三。有没有办法用更少的代码在R中完成所有这些工作？
最佳答案
使用.indexmon等直接访问POSIXlt值
GSPC[.indexmon(GSPC)==10 & .indexmday(GSPC) > 22 & .indexmday(GSPC) < 29
&.indexwday(GSPC) == 3]
GSPC.Open GSPC.High GSPC.Low GSPC.Close GSPC.Volume GSPC.Adjusted
2007-11-28 1432.95 1471.62 1432.95 1469.02 4508020000 1469.02
2008-11-26 852.90 887.68 841.37 887.68 5793260000 887.68
2009-11-25 1106.49 1111.18 1104.75 1110.63 3036350000 1110.63
2010-11-24 1183.70 1198.62 1183.70 1198.35 3384250000 1198.35
2011-11-23 1187.48 1187.48 1161.79 1161.79 3798940000 1161.79