2015-01-13 4 views
0

У меня есть столбец 84 ежемесячных расходов от 1/2004 - 12/2010, который в Excel выглядит как ...Возникли проблемы с временными рядами АиР объекты

12247815.55 
11812697.14 
13741176.13 
21372260.37 
27412419.28 
42447077.96 
55563235.3 
45130678.8 
54579583.53 
43406197.32 
34318334.64 
25321371.4 
...(74 more entries) 

Я пытаюсь запустить СТЛ () из прогнозного пакета на этой серии, и поэтому я загрузить данные:

d <- ts(read.csv("deseason_vVectForTS.csv", 
header = TRUE), 
start=c(2004,1), 
end=c(2010,12), 
frequency = 12) 

(Если я заголовок = FALSE он будет поглощать первую запись - 122 ...- как заголовок для второго столбца, и назовите заголовок первой колонны 'X')

Но вместо моей среды, заполненной Time Series Object from 2004 to 2011 (как было сказано ранее), она просто говорит ts[1:84, 1].

вероятно, связано это тот факт, что,

fit <- stl(d) 

бросает

Error in stl(d) : only univariate series are allowed. 

несмотря на то, что

head(d) 
[1] 12247816 11812697 13741176 21372260 27412419 42447078 

и

d 
      Jan  Feb  Mar  Apr  May  Jun  Jul  Aug  Sep  Oct 
2004 12247816 11812697 13741176 21372260 27412419 42447078 55563235 45130679 54579584 43406197 
("years 2005-2010 look exactly the same, and all rows have columns for Jan-Dec; it just doesn't fit on here neatly - just trying to show the object has taken the ts labeling structure.") 

Что я делаю неправильно? Насколько я знаю, это точно так же, как в прошлом я строил объекты временного ряда ...

+0

Что? http://www.inside-r.org/packages/cran/forecast/docs/stlf Этот фрагмент - это то, на что он похож в Excel. – d8aninja

+0

Просьба на консоль и посмотреть функцию stl в этом пакете, прежде чем отвечать на запросы - определенно используя функцию stl. – d8aninja

+0

stl() не находится в пакете прогноза. Это часть пакета статистики. stlf() находится в пакете прогноза. –

ответ

1

read.csv читает в матрице. Если он имеет только один столбец, он все еще является матрицей. Для его использования в векторе

d <- ts(read.csv("deseason_vVectForTS.csv", 
header = TRUE)[,1], 
start=c(2004,1), 
end=c(2010,12), 
frequency = 12) 

Также проверьте свои факты. stl находится в пакете stats, а не в пакете forecast. Это легко проверить, используя help(stl).

Смежные вопросы