Tworzenie obiektów szeregów czasowych w R
Szereg czasowy można traktować jako wektor lub macierz liczb uzupełnioną o informacje o tym, kiedy zostały one zarejestrowane. W R informacje te są przechowywane w obiekcie ts. W większości ćwiczeń będziesz korzystać z szeregów czasowych dostępnych w gotowych pakietach. Jeśli jednak chcesz pracować z własnymi danymi, musisz wiedzieć, jak utworzyć obiekt ts w R.
Spójrz na poniższy przykład usnim_2002, zawierający marże odsetkowe netto banków w USA za rok 2002 (źródło: FFIEC).
> usnim_2002
usnim
1 2002-01-01 4.08
2 2002-04-01 4.10
3 2002-07-01 4.06
4 2002-10-01 4.04
> # ts(data, start, frequency, ...)
> usnim_ts = ts(usnim_2002[, 2], start = c(2002, 1), frequency = 4)
Funkcja ts() przyjmuje trzy argumenty:
data– wszystkie dane zusnim_2002z wyjątkiem kolumny z datą; nie jest ona potrzebna, ponieważ obiekttsprzechowuje informacje o czasie osobno.start– data pierwszej obserwacji w postacic(rok, okres). Styczeń odpowiada okresowi 1; analogicznie, data startowa w kwietniu to 2, w lipcu – 3, a w październiku – 4. Innymi słowy, okres odpowiada kwartałowi roku.frequency– wartość 4, ponieważ dane mają częstotliwość kwartalną.
W tym ćwiczeniu wczytasz dane szeregu czasowego z pliku xlsx za pomocą read_excel() – funkcji z pakietu readxl – i zapiszesz je jako obiekt ts. Zarówno plik xlsx, jak i pakiet zostały już załadowane do twojego środowiska pracy.
To ćwiczenie jest częścią kursu
Prognozowanie w R
Instrukcje do ćwiczenia
- Użyj funkcji
read_excel(), aby wczytać dane z pliku"exercise1.xlsx"do zmiennejmydata. - Zastosuj funkcję
head()namydataw konsoli R, aby sprawdzić pierwsze kilka wierszy danych. Zwróć uwagę na daty – w 1981 roku znajdują się cztery obserwacje, co wskazuje na dane kwartalne z częstotliwością czterech wierszy na rok. Pierwsza obserwacja, czyli data startowa, toMar-81– pierwszy z czterech wierszy dla roku 1981, co oznacza, że marzec odpowiada pierwszemu okresowi. - Utwórz obiekt
tso nazwiemytsza pomocą funkcjits(). Ustaw argumentydata,startifrequencyna podstawie poczynionych obserwacji.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# Read the data from Excel into R
___ <- ___("exercise1.xlsx")
# Look at the first few lines of mydata
___
# Create a ts object called myts
myts <- ts(___[___], start = c(___, ___), frequency = ___)