Mil*_*der 3 awk grep cut r sed
我有一个时间索引矩阵(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文件需要-将字符串"赢家" 分开,所以我只是在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