Я новичок в ts и xts object.
При работе с данными временных рядов я столкнулся с проблемой
require(quantmod)
require(forecast)
ticker <- "^GSPC"
getSymbols(ticker, src="yahoo", to = "2013-12-31")
prices <- GSPC[,6] # First get data using package quantmod
# then forecasting using package forecast
prices.ts <- as.ts(prices)
prices.ets <- ets(prices.ts)
prices.fore <- forecast(prices.ets, h=10)
# then plot
plot(prices.fore, xaxt = "n")
Мои проблемы:
1. Когда я пытался сохранить GSPC с датой в файле csv. Я искал и пробовал это
write.zoo((GSPC, file = "GSPC.csv", sep = ",", qmethod = "double"))
Сообщение об ошибке: Error: unexpected ',' in "write.zoo((GSPC,"
, Я проверил синтаксис, он кажется правильным, и попробовал другие комбинации. Все не удалось с аналогичным сообщением об ошибке.
также я попытался index(GSPC)
узнать дату.
а затем cbind(index(GSPC), GSPC[, 6])
. Это тоже не удалось ..
Сообщение об ошибке: Error in merge.xts(..., all = all, fill = fill, suffixes = suffixes) :
dims [product 1762] do not match the length of object [3524]
но когда я проверил длину
> length(GSPC[,6])
[1] 1762
> length(index(GSPC))
[1] 1762
2. сюжет такой
нет x-lab и y-lab. Я пробовал методы принятого ответа, размещенного здесь,. но не удалось.
В частности, я не понимаю цель следующего кода. Это вроде меняет внешний вид сюжета, но никак не меняет внешний вид. Не знаю, теряю ли я какие-то очки.
a = seq(as.Date("2011-11-01"), by="weeks", length=11)
axis(1, at = decimal_date(a), labels = format(a, "%Y %b %d"), cex.axis=0.6)
abline(v = decimal_date(a), col='grey', lwd=0.5)
Кроме того, я хочу построить график из as.Date («2013-01-01»).
Не могли бы вы дать несколько предложений?
Большое спасибо!