c n h u oy w q cav y - r · pdf fileq di o v b o w x q ovdu yx od qow t xv tr u qo m v u w q...

62

Upload: lykien

Post on 31-Jan-2018

216 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

CNwt� CQmi �Q u�H xm u� s�v x�

�QwBU��a O�yW� jQ� w �� CavY x�oWv�O

QO xS� w EL��t

������ Q�DWwv CU�Q�w�

�vWwOv�wUwt O�aUO�U

s�mousavi�pwut�ac�ir

���� R���B

Page 2: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R
Page 3: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�r��t CUQyi

O Q�DioV�B

� �v�tR ��y�QU pw�pYi� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � �v�tR �QU ���

� � � � � � � � � � � � � � � � � �rYi C�Q��eD w OvwQ ��v�tR �QU s�UQD �����

� � � � � � � � � � � � � � � � � � � � �y�QU s�UQD �������

� � � � � � � � � � � � � � � � � � �v�tR �QU OvJ ������

� � � � � � � � � � � � � � � � � �rYi C�Q��eD w OvwQ ������

� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � �v�tR ��y�QU x RHD ��

� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � �ypOt ����

� � � � � � � � � � � � � � � � � � � � � � � � � � R QO �yxir�wt x RHD ���

� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � �oDU�ty ��

� � � � � � � � � � � � � � � � � � � � � � � Tv� Q�w w u�ov��t h QaD ����

� � � � � � � � � � � � � � � � � � � � � �oDU�ty w Tv� Q�wwm h QaD ���

� � � � � � � � � � � � � � � � � �oDU�tyOwN w Tv� Q�wmwD� ���wD h QaD ���

�� � � � � � � � � � � � � � � � � �oDU�tyOwN x�U�Lt �����

�� � � � � � � � � � � � � � � � � � � � Q�ov�oDU�ty ����

hr�

Page 4: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� x��B �iO�YD ��ypOt swOpYi�� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � O�iU xiwv ��

�� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � xtOkt ����

�� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � h QaD ���

�� � � � � � � � � � � � � � � � � � � � � � � � � � � R QO �R�Ux��W ���

�� � � � � � � � � � � � � � � � � � � Q�ov�oDU�ty w swO x�DQt X�wN ����

�� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � �iO�YD uORs�o �

�� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � xtOkt ���

�� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � h QaD ��

� � � � � � � � � � � � � � � � � � � � � � � � � � wQUB p�kDv� Qorta ��

� � � � � � � � � � � � � � � � � � � � �iO�YD uORs�o swO x�DQt X�wN ���

� � � � � � � � � � � � � � � � � � � � � � � � � � � � p��iD Qorta ���

� � � � � � � � � � � � � � � � � � � � � � � � � � � � � �R�Ux��W ���

� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � w�UQoQwD� ��ypOt �

� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � h QaD ���

� � � � � � � � � � � � � � � � � � � � w�UQoQwD� QO � �DU ��v w � �DU � ��

� � � � � � � � � � � � � � � � � � � � � AR��� pOt swO x�DQt X�wN ��

� � � � � � � � � � � � � � � � � � � � � � AR��� Ov �Qi Q�ov�oDU�ty ���

� � � � � � � � � � � � � � � � � � � � � � � � RH �oDU�tyOwN ���D ���

� � � � � � � � � � � � � � � � � � � � � � � � � � � � � �R�Ux��W ���

� � � � � � � � � � � � � � � � � � � � � � � � � xDi� VR�Q� ��ypOt ���

� � � � � � � � � � � � � xOW �R�Ux��W ��yxO�O Q� xDi� VR�Q� ��ypOt ���

� � � � � � � � � � � � � � � � xDi� VR�Q� AR pOt ��v�yH ��tO �QU ���

�� �DU�� ��ypOt swUpYi� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � lQLDt u�ov��t ��ypOt ��

� � � � � � � � � � � � � � � � � � � � X�wN w h QaD �MA�q� Ov �Qi ����

� � � � � � � � � � � � � � � �R�Ux��W w Q�ov�oDU�ty �R ��yp�Ft ���

� � � � � � � � � � � � � � � � � � � � � � � � � � � � � xDi� VR�Q� MA ��ypOt �

� � � � � � � � � � � � � � � � � xOW �R�Ux��W �QU x� xDi� VR�Q� pOt ���

� � � � � � � � � � � � � � � � � � � � � � � � � � ARMA Ov �Qi ����mQD ��ypOt �

� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � h QaD ���

� � � � � � � � � � � � � � � � � � � � � � � � � pOt swO x�DQt X�wN ��

� � � � � � � � � � � � � � � � � � � � � � � � � � �� QHDp�rLD �ARMA ��ypOt ��

� � � � � � � � � � � � � � � � � � � � � � � � � � VR�Q� w �R�Ux��W ����

Page 5: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� � � � � � � � � � � � � � � � � � � � � � � � � � � � p O�D MQv �QU ���

� �DU���v ��ypOt sQ�yJpYi�� � � � � � � � � � � � � � � � � � � � � � � � � � � � � �rYi Q�e ARIMA pOt ����� � � � � � � � � � � � � � � � � � � � � � � � jQ� O�rwD �QU w p��iD ������� � � � � � � � � � � � � � � � � � � � � � � � � � � � � �k�irD pOt ����� � � � � � � � � � � � � � � � � � � � � � � � � � � �yp�Ft w h QaD ����� � � � � � � � � � � � � � � � � � � � � � � � � � VR�Q� w �R�Ux��W ������� � � � � � � � � � � � � � � � � � � � � � � � � � � � � �rYi ARIMA ��ypOt ���� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � h QaD ������ � � � � � � � � � � � � � � � � � � � � � � � � � � � � VR�Q� x wQ ����� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � ARCH ��ypOt ���� � � � � � � � � � � � � � � � � � � � � � � � � � SP��� ��y�QU ������ � � � � � � � � � � � � � � � ARCH pOt h QaD ��Q�O �B�v uOQm pOt ����� � � � � � � � � � � � � � � � � GARCH ��ypOt x� �w� Qt ��y�U� ����� � � � � � � � � � � � � � � � xDi� VR�Q� GARCH pOt w �R�Ux��W ����� � � � � � � � � � � � � � � � � � � � � � � � SP��� �QU Q� VR�Q� ����� � � � � � � � � � � � � � � � � � � � � � � � s�rk� �QU QO �Q�O �B�v ������ � � � � � � � � � � � � � � � � � �R�Ux��W w �v��V�B QO GARCH ����

� �H�Qt

G

Page 6: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

Q�DioV�B

OW x��Q� OvDUy x�Q ��OD�� QO w �OD�t xv�tR u�� QO xm �v�Um ���vW Qw�vt x� R Q�Ri�sQv ��Q� �twta Q�DWwvCU� �R�UxO��B p��k u�o��Q CQwYx� xm CiQo Q�Qk Q� R TQO x� R C��U QO w

http���cran�r�project�org�doc�contrib�Mousavi�R�lang�in�Farsi�pdf

j� Q� R� �Q Q�kL sy ��xOa w CU� xDiQo Q�Qk u�Ovtxkqa xO�iDU� OQwt xOW xDio xDWwv xm OUQ�t Q�v x�OwW�t �Q�RoU�BU �yv R� �H u�ty R� xm OvO�O Q�Qk V� wN C�Lt w h�r OQwt l�vwQDmr� CUB

EL��t� u�wva Q� R Q �L Q�DWwv u��Q��v� sO�Di� R OOt x� X�NEL��t R� ��xQ�B x�yD Qmi x� x�Q xt�O� QO xOv�OOQo�t KQ�t xOv� lOv� Ca��� OL QO � � � � w Q�tu�tR ��OL Q�O�kt ��v�tR ��y�QU Ovv�t� �R QO xS� wx� EL��t u�� w OwW�t s�OkD pYi OvJ QO u KQW xm OQ�O �v�tR ��y�QU x� X�YDN� u EL�t u�rw�Q� ZQi xm OwW�t u�Wv Q��N �Hv�� QO OwW pY�L �QDKkvt xDWwv �D OOQo�t V��Q� w w O����t xt�O� QwQt

CU� xOwtv Um �Q R Q�Ri�sQv �� x�rw� ���vW �xS� w EL��t �kaD ��Q� �t�Qo xOvv�wN xm CU� u��u ��kDQ� w KqY� ��Q� �O�yvW�B w Q�v Q�y�� xvwo Qy R� xOvU� wv �CU�v prN R� �r�N Q �L xR�Hw ��trUt

Owtv Oy�wN xO�iDU� w p��kDU�

�vWwOv�wUwt O�aUO�U

���� R���B � u�QyD

O

Page 7: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

pw� pYi

�v�tR ��y�QU

�v�tR �QU ��

x� Ov�wD�t C�Oy�Wt C�F pta Ovvm�t Q��eD u�tR QO xm CU� ���yOv�Qi C�F R� OvDQ��a �v�tR ��y�QU�� OOQo�t O�m �D �w�Ut �v�tR pY�wi �� xDUUo C�F �wQ �Hv�� QO OQ�PB s�Hv� xDUw�B �� w xDUUo CQwYx� C�Oy�Wt xawtHt w ����� � � �n CQwY x� �v�tR xawtHt �U�vt �v�tR T��kt w �v�tR ��O�t ��NDv�

OwW�t xO�O u�Wv x�� x�� � � � � xn �� w fxt � t � �� �� � � � � ng pmW

�rYi C�Q��eD w OvwQ ��v�tR QU s�UQD ���

�y�QU s�UQD �������

u� Qi�Ut O�OaD xm OwW�t xO�iDU� AirPassengers u�wva Q� R R OwN QO ���yxO�O p��i R� � � � w s�UQD ��Q�CU� xOW C�F ���� xQwO QO Q�t u�� OyO�t u�Wv �Q �Qiv Q�Ry UL Q�� xv�y�t �rrtr�u��

�p�Ft

� data�AirPassengers�

Page 8: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt

� AP �� AirPassengers

� plot�AP ylab �Passengers ����� s���

Ov�wD�t xm OQ�O OwHw window�� s�v x� R u�� R QO �Qo�O ��D CU� q�� ��yO�m ��QH� pY�L ��� pmW

Time

Pas

seng

ers

(100

0’s)

1950 1952 1954 1956 1958 1960

100

200

300

400

500

600

u Qi�Ut xv�y�t �v�tR �QU V �tv ���� pmW

O�vm xHwD Q� R p�Ft x� OyO CUOx� �Q �v�tR �QU R� ��xawtHt Q� R l�

� x �� AirPassengers

� y �� window�x ���� �����

CU� Q� R CQwY x� u xH�Dv xm

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

�� ��� ��� �� ��� ��� � �� �� ��� ��� �� �

��� ��

O�vm xHwD Q� R p�Ft x� x�t �qFt �O�W�� xDW�O pQDvm p�U R� �rYi �wQ O�y�wN� Qo �

� x �� AirPassengers

� y �� window�x ���� c��������

CU� Q� R CQwY x� u xH�Dv xm

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

�� ��� ��� �� ��� ��� � �� �� ��� ��� �� �

O�vm xHwD Q� R p�Ft x� ��

� x �� AirPassengers

� y �� window�x c������� c��������

Page 9: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�v�tR ��y�QU � pYi

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

�� ��� �� ��� ��� � �� �� ��� ��� �� �

��� �� �� ���

OQ�Ov M� Q�D xm O�vm xHwD Q� R �xv�y�t �v�tR �QU x� uwvm �

V� V� V� V V� V� V� V� V V� V�� V��

� ��� ��� ��� �� � � �� ���� ��� ��� �� ���� ���� ���

� �� ���� ���� ��� ����� ���� ����� ���� ���� � ���� ����

� ���� ��� ����� ����� ��� ����� ����� ���� ��� ���� �� � ���

���� ��� ���� � ��� ��� ���� ���� ���� ��� ��� ���� ��

� ���� ���� ���� � � ����� � � ���� ���� ���� ��� ���� ������

������

������

������

������

������

x� Owtv sUQ �Q �yu TBU w xO�O Q�Qk sy QU CWB �Q �Q�t ��yp�U O��� �OD�� xv�y�t ��yxO�O s�UQD ��Q�CU� Q� R KQW x� pw� CQwY ��yO�m O�O s�Hv� �Q Q�m u�� u�wD�t CQwY wO

� x �� matrix�scan��F��R��les�data�ghar�txt�� ncol�� byrowT�

� y �� t�x�

� y �� matrix�yncol� byrowT�

� ts�plot�y�

CU� Q� R KQW x� swO CQwY ��yO�m

� x �� read�table��F��R��les�data�ghar�txt��

� y �� as�matrix�y�

� y �� t�x�

� y �� matrix�y ncol� byrowT�

� ts�plot�y�

�v�tR �QU OvJ �������

O�rwD p��i � p�U �D ��� p�U R� ��r�QDU� QO �qFt O�O u�Wv s��wD CQwYx� �Q �v�tR �QU OvJ u�wD�tu�� OQ�O OwHw �u�D UL Q�� CqmW w �QD�r uw�r�t UL Q�� Q�aWr���t ��Ca�U C�wwr�m UL Q�� jQ�

OwW�t xOv�wN R O�m �UwD online CQwYx� �QU xU

� �http���www�massey�ac�nz��pscowper�ts�cbe�dat�

� CBE �� read�table�www header T�

CBE���� �

CU� Q� R CQwYx� q�� ��yO�m �HwQN

Page 10: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

choc beer elec

� ���� ���� ���

�� ���� ����

� � ��� ����

� �� ���� ����

��Q�O w OvQ�Ov M� Q�D Qo�O CQ��a x� OvDUy x�t w p�U Ok�i q�� ��yxO�O OwW�t x�Lqt xm Qw�u�ty��Q� Owtv p�O�D �v�tR �QU x� ts�� ��D R� xO�iDU� �� �Q ���W� u�� O��� u��Q��v� OvDU�v �v�tR �QU Q�DN�U

O�vm xHwD Q� R ��yO�m x� Q�m u��

� �http���www�massey�ac�nz��pscowper�ts�cbe�dat�

� CBE �� read�table�www header T�

� Elec�ts �� ts�CBE� �� start ���� freq ���

� Beer�ts �� ts�CBE� �� start ���� freq ���

� Choc�ts �� ts�CBE� �� start ���� freq ���

� plot�cbind�Elec�ts Beer�ts Choc�ts��

CU� q�� ��yO�m ��QH� pY�L ��� Q�Owtv

2000

6000

1000

014

000

Ele

c.ts

100

150

200

Bee

r.ts

2000

6000

1960 1965 1970 1975 1980 1985 1990

Cho

c.ts

Time

cbind(Elec.ts, Beer.ts, Choc.ts)

�v�tR �QU OvJ V �tv ��� pmW

�� trend �� seasonal variation

Page 11: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� �v�tR ��y�QU � pYi

�rYi C�Q��eD w OvwQ �������

��v��L� w �CQB Q�O�kt ���rYi C�Q��eD ��OvwQ xmr� OyO�t CUOx� �Q �QU �wor� �yvD xv �v�tR �QU sUQ��rm Cr�L QO w CU� �rYi C�Q��eD p�U Qy QO wor� Q�QmD q�� p�Ft QO �OR�U�t Qy�� R�v �Q ��N Q�O�ktu��D ��Q� pOt u� QDxO�U Ov� wo OvwQ �OvwW�tv Qy�� �� w�vD CQwYx� xm �v�tR �QU QO l�D�tDU�U C�Q��eDQ�t �tHD �� �rYi C�Q��eD Q�F��D u�wD�t �OvwQ K �w uOQw CUOx� ��Q� CU� ��NVy�m w V��Ri� �OvwQQy Q�O�kt xYqN OwW�t s�Hv� aggregate�� ��D R� xO�iDU� �� R QO Q�m u�� OQ� u�� R� xvq�U Q�t x� xv�y�tOR�U�t u�at �yxO�O R�l� Qy ��Q� �Q pYi cycle�� ��D w OwW�tXNWt boxplot�� ��D �UwD pYi

� data�AirPassengers�

� AP �� AirPassengers

� layout�����

� plot�aggregate�AP��

� boxplot�AP � cycle�AP��

CU� q�� ��yO�m ��QH� pY�L ��� Q�OwtvOOQo�t x�Lqt xm Qw�u�ty OwW�t xOv�wN R O�m �UwD online CQwYx� xv�y�t �Q�m�� ��yxO�O �p�Ft

Time

aggr

egat

e(A

P)

1950 1952 1954 1956 1958 1960

2000

5000

1 2 3 4 5 6 7 8 9 10 11 12

100

400

�yv� boxplot w xOW ��tHD xvq�U �QU V �tv ��� pmW

ts�� ��D R� �v�tR �QU x� �yv p�O�D ��Q� w OvDU�v �v�tR �QU �yxO�O u�� CU� unemploy Q�eDt s�v�tHD ��D �UwD xvq�U �UwDt OOQo�t x�Lqt u QO freq �w�vD xQwO w �wQW QDt�Q�B xm �OwW�t xO�iDU�

OwW�t u�at s�UkD Qorta w

� www �� �http���www�massey�ac�nz��pscowper�ts�Maine�dat�

� Maine�month �� read�table�www header TRUE��� outliers

Page 12: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

� Maine�month�ts �� ts�unemploy start c����� �� freq ���

� Maine�annual�ts �� aggregate�Maine�month�ts����

� layout�����

� plot�Maine�month�ts ylab �unemployed �����

� plot�Maine�annual�ts ylab �unemployed �����

CU� q�� ��yO�m ��QH� pY�L ��� Q�Owtv

Time

unem

ploy

ed (

%)

1996 1998 2000 2002 2004 2006

34

56

Time

unem

ploy

ed (

%)

1996 1998 2000 2002 2004

3.5

4.5

�yv� �UwDt w �Q�m�� �QU V �tv ���� pmW

�v�tR ��y�QU x RHD ��

�ypOt ���

�ypOt u��Q��v� OwW�t xOy�Wt �rYi C�Q��eD �� w OvwQ �y�QU R� �r�N QO �OW x�Lqt �q�k xm Qw�u�tyOW���t Q� R CQwYx� xm �CU� ��atH pOt ��v�tR �QU xO�U x� RHD l� OvDUy ���yxir�wt ��Q�O

xt � mt � st � zt

xir�wt CU� �rYi C�Q��eD xOv��tv st w OvwQ xv�Wv mt xir�wt CU� xOW xOy�Wt �QU xt w u�tR t u QO xmOyO�t u�Wv �Q ��N R�v zt

u�Wv Q� R CQwYx� xm �CU� ��� Q pOtl�U�vt pOt �OwW Qy�� OvwQl� CQwYx� �rYi C�Q��eD Qo �OwW�t xO�O

xt � m�t�s

�t�z

�t

�� additive �� multiplicative

Page 13: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� �v�tR ��y�QU � pYi

OOQo�t p�O�D �atH pOt x� �� Q pOt �Q�N� x���Q R� sD� Q�or uDiQo �� OwW�t x�Lqt xm Qw�u�ty

yt � ln�xt� � ln�m�t� � ln�s�t� � ln�z�t�

� mt � st � zt

R QO �yxir wt x� RHD ���

Qo � Ovm�t OQwQ� lQLDt u�ov��t VwQ R� xO�iDU� �� �Q �rYi C�Q��eD w OvwQ decompose�� ��D �R QOx� �w� Qt xO�O p�Ft ��Q� OOQo�t O�H�� Q�Owtv xU �� pmW l� x�ov �OQ�o Q�Qk plot�� ��D pN�O jwi ��D��yDv� QO w OwW�t s�UQD ���� pmW� �� Q w �atH pOt CQwYx� Q� R ��yO�m QO ��q�� p�Ft� jQ� O�rwDpmW� OOQo�t x�Lqt �Hm� �rYi C�Q��eD w OvwQ Q�Owtvl� �wQ xQNq�� w OwW�t p��vO �� Q pOt �yO�m

����

2000

8000

1400

0

obse

rved

2000

6000

1200

0

tren

d

0.90

1.00

1.10

seas

onal

0.94

1.00

1.06

1960 1965 1970 1975 1980 1985 1990

rand

om

Time

Decomposition of multiplicative time series

�QU ��yxir�wt V �tv ���� pmW

�oDU�ty ��

�oDU�ty u�vJ O��� �yv p�rLD w x� RHD ��Q� Ow� Ovy�wN xDU�ty �r�wDt ��yQ�eDt �Cq�L R� �r�N QO��yxO�O j� Q� R� w OwW�t h� QaD �oDU�ty ��D �UwD �v�tR ��y�QU QO �oDU�ty Q�DN�U OOQo u�at

OOQo�t OQwQ� xOW xOy�Wt

Page 14: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

Time

1960 1965 1970 1975 1980 1985 1990

2000

4000

6000

8000

1000

014

000

�rYi C�Q��eD w OvwQ V �tv ���� pmW

Tv�� Q�w w u�ov��t h� QaD ���

Q�eDt u�ov��t E�x� u��Q��v� CU� xat�H u�ov��t k�w QO �OwW�t xO�O u�Wv E hQL �� xm � �� Q O�t�x� xm CU� � pwL C�i�QLv� � Qt u�ov��t E��x � ���� CQ��a OwW�t xO�O V��tv R�v � �� xm CU� x

OwW�t xO�O V��tv R�v �� �� xm OwW�t xDN�vW x Q�eDt Tv�� Q�w

�oDU�ty w Tv�� Q�wwm h� QaD ���

CU� h� QaD p��k Q� R CQwYx� Tv�� Q�wwm s�v x� �twyit �O�W�� xDW�O �x� y� Q�eDt wO Qo �

��x� y� � E��x� �x��y � �y��

R� n xR�Ov� x� ��xvwtv l� Qo � OQ�o�t xR�Ov� �Q �x� y� Q�eDt wO u�� ��N ���DQ� u�R�t k�w QO Tv�� Q�wwmCU� Q� R CQwYx� u x���Q xm Owtv OQwQ� �Q xvwtv Tv�� Q�wwm u�wD�t O�W�� xDW�O �xi� yi�

Cov�x� y� �X

�xi � �x��yi � �y��n

w �� u�� � Q u�� Ovm�t u�at Oa�� uwO� CQwYx� �Q �x� y� Q�eDt wO u�� ���DQ� u�R�t �oDU�ty � Q OQ�Ov OwHw �x� y� ��yQ�eDt u�� ��N ���DQ� xm CU� �vat u�O� OW�� QiY � Q u�� Qo � OQ�O Q�Qk ��

OwW�t h� QaD Q� R CQwYx� �x� y� ��yQ�eDt u�� � �va� xat�H �oDU�ty

��x� y� �E��x� �x��y � �y��

�x�y�

��x� y�

�x�y

Page 15: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� �v�tR ��y�QU � pYi

OOQo�t x�U�Lt Q� R CQwYx� xvwtv �oDU�ty � Q

Cor�x� y� �Cov�x� y�

sd�x� sd�y�

�oDU�tyOwN w Tv�� Q�wmwD� ���wD h� QaD ���

CQwYx� �v�tR �QUl� u�ov��t ��D OwW�t xDN�OQB u�ov��t x� �OD���v�tR ��y�QU �oDU�tyh� QaD ��Q�OW���t t R� �a��D w CU� Q� R

��t� � E�xt�

OQwQ� w CU� ��DU�� u�ov��t QO �v�tR �QU x�ov �OW��v u�tR R� �a��D �va� �OW�� C��F u�ov��t ��D Qo �CU� Q� R CQwYx� u ��xvwtv

�x �nXi��

xt�n

OW���t Q� R CQwYx� �CU� �DU�� u�ov��t QO xm �v�tR �QU l� Tv�� Q�w ��D

���t� � E��xt � ����

CQwYx� u ��xvwtv OQwQ� OOQo�t �� C��F x� p�O�D ���t� x�ov �OW�� �DU�� R�vTv�� Q�w QO �v�tR �QU Qo �CU� Q� R

Var�x� �

P�xt � �x��

n� �

xDU�ty CU� umtt �yQ�eDt �t� OvOw� �DU�� Tv�� Q�w w u�ov��t QO xm O�Owtv x�Lqt �Q �v�tR �QU uwvm �D�oDU� �yv ��ys�o Cw�iD x� �ki �yQ�eDt u�� �oDU�ty xm CU� swO x�DQt ��DU�� �v�tR �QU OvW�� R�v

Ov� wo �Q�N��D �Q �yQ�eDt u�� ��ys�o O�OaD OW�� xDW�O�QU Qo � Ovt�v ��B��B �oDU�ty �� w �oDU�tyOwN �Q hrDNt ��yu�tR QO VOwN �� Q�eDtl� �oDU�tyCQwYx� �QCU� k R� �a��D xm �k �va� �acvf�Tv�� Q�wwmwD� ��D u�wD�t x�ov �OW�� swO x�DQt ��DU�� �v�tR

Owtv u��� Q� R�k � E��xt � ���xt�k � ���

��D CU� E�xt� � E�xt�k� � � �va� �CU�v u�tR x� xDU��w u�ov��t �Q� R OQ�Ov u�tR x� �oDU� �k ��DOwW�t h� QaD Q� R CQwYx� k Q�N��D �� �acf� �oDU�tyOwN

�k ��k��

OW���t �� � �� �Q� R �CU� �� � � xm CiQo xH�Dv u�wD�t q�� h� QaD R� xO�iDU� ���� stationary �� lag � Autocorrelation � Serial correlation

Page 16: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt ��

�oDU�tyOwN x�U�Lt �������

xOy�Wt CiH n� � u�wD�t OW�� xDUUo �v�tR �QU l� x�� x�� � � � � xn xOy�Wt n O�vm ZQi uwvm ��va�

�x�� x��� �x�� x��� � � � � �xn��� xn�

ZQi swO Q�eDt u�wva x� �Q xOy�Wt u�twO w pw� Q�eDt u�wva x� GwR Qy QO �Q xOy�Wt u�rw� Qo � CN�U �QCU� Q� R CQwY x� �oDU�ty � Q �O�vm

r� �

n��Pt��

�xt � �x���

� �xt�� � �x���

�s

n��Pt��

�xt � �x���

�� n��Pt��

�xt�� � �x���

��

Ovt�v C�Oy�Wt u�� �oDU�tyOwN � Q �Q u w CU� �x� � �n��

nPt��

xt w �x� � �n��

n��Pt��

xt u QO xm

CU� hrDNt ��yu�tR QO x Q�eDt u�� �oDU�ty k�w QO xmOwW�t xO�iDU� Q� R pwtQi R� r� OW�� nQR� �i�m xR�Ov� x� n Qo �

r� �

n��Pt��

�xt � �x� �xt�� � �x�

nPt��

�xt � �x��

Owtv x�U�Lt k Q�N��D �� �Q �oDU�tyOwN � Q u�wD�t j� Q� u�ty x� CU� �x � �n

nPt��

xt u QO xm

rk �

n�kPt��

�xt � �x� �xt�k � �x�

nPt��

�xt � �x��k � � �� �� � � �

R� �Q �v�tR �QU C�aq�� R� �r�N CQwY u�� QO �Q� R� OQm x�U�Lt �Q �yrk O���v k � nQ�O�kt ��Q� �qwtat�O�yO�t CUO

OwW�t xDiQo Q�v QO � �pnQ�Okt � uOw� Q�O�vat K�U �� rk ��xrY�i OQwQ� ��Q�

�Q�ov�oDU�ty �������

OwW�t xO�iDU� ��ypw� QwLt� k ULQ� ��yZQa QwLt� rk Q�Owtv R� ��v�tR ��y�QU Q�UiD w Q��aD ��Q�QwLt CU� Oa�� uwO� OQ�O Q�Qk �wQ �oDU�tyOwN Q�O�kt xm �yZQa QwLt Ovt�v Q�ov�oDU�ty �Q Q�Owtv u����� correlogram

Page 17: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� �v�tR ��y�QU � pYi

OQ�O �Q�oxvwtv �v�tR xrY�i x� �oDU� OyO�t u�Wv �Q �yQ�N��D xm �ypw�p�Ft u�wvax� OwW�t xO�iDU� R QO lag�plot�� ��D R� sy x� C�Uv hrDNt ��yQ�N��D s�UQD ��Q�

� lag�plot�LakeHuron lag� labelsF do�linesF diag�col �red��

��Q� � �� pmW� OOQo�t s�UQD Q�N��D Q�yJ �D LakeHuron ��yxO�O �CU� lag� uwJ �q�� p�Ft QO

lag 1

Lake

Hur

on57

657

757

857

958

058

158

2

576 577 578 579 580 581 582

lag 2

Lake

Hur

on

lag 3

Lake

Hur

on

lag 4

Lake

Hur

on

576

577

578

579

580

581

582

576 577 578 579 580 581 582

Q�N��D ��yQ�Owtv ���� pmW

CU� Q� R CQwYx� u �rm pmW OwW�t xO�iDU� acf�� ��D R� R QO Q�ov�oDU�ty s�UQD w rk x�U�Lt

acf�x lag�maxNULL typec��correlation��covariance��partial�� plotTRUE�

u QO xmCU� �yxO�O Q�OQ� !x

�x���Q R� xm CU� ZQiV�B ��Q�O R �OOQo hPL u�twoQ u�� Qo � CU� Q�v OQwt Q�N��D O�OaD !lag�maxOW���t � ��v�t QO u sD� Q�or w �yxvwtv xR�Ov� N xm �OOQo�t x�U�Lt � �log�N�

�correlation� �OOQo hPL u�twoQ u�� Qo � Ovm�t u�at �Q acf ��D �wv w CU� �yQDm �Q�m R� ��xDWQ !typeCU� ��D u�� ZQiV�B

OOQo�t sUQ Q�ov�oDU�ty Q�Owtv CQwY u�� QO wCU�TRUEZQiV�B ��Q�O �k�vt u�twoQ u�� !plotOwW�t A�J �yQ�N��D ��R� x� acf Q�O�kt w OwW�tv sUQ xOW xDio Q�Owtv �OW�� FALSE u�twoQ u�� Qo �R� �yv x� �UQDUO ��Q� w OOQo�t xQ�NP acf�x��acf CQwYx� �Q�OQ� QO Q�O�kt u�� xm O�W�� xDW�O xHwDk� Qo � CU� acf�acf�k��� CQwYx� q�� Q�OQ� u��Q��v� OwW�t xO�iDU� Q�OQ�T�Ov� u�wva x� Q�N��D u�tR

CU� l� Q��Q� acf�acf��� Q�Okt �OW��O���tv x�Lqt �Q Q� R ��yO�m acf ��D p�Ft u�wvax� uwvm �

Page 18: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

� www �� �http���www�massey�ac�nz��pscowper�ts�wave�dat�

� wave�dat �� read�table �www headerT�

� par�mfrowc���� marc�������

� plot�ts�waveht� ylab wave height �mm� �

� acf�waveht main���

CU� q�� ��yO�m ��QH� pY�L "�� Q�Owtv

wav

e he

ight

(m

m)

0 100 200 300 400

−50

00

500

0 5 10 15 20 25

−0.

50.

00.

51.

0

AC

F

u� Q�ov�oDU�ty w �v�tR �QU V �tv ���� pmW

OwW�t x�U�Lt Q� R CQwYx� k�� QO acf Q�Okt p�k p�Ft xt�O� QO

� acf�waveht��acf���

��� ���������

OwW�t x�U�Lt Q� R CQwYx� �OW�� k�� QO Tv�� Q�wwmwD� x�U�Lt Qw�vt Qo � ��OOHt

� acf�waveht type c��covariance����acf���

��� ��������

xrY�i xv�Wv xm OvQ�O OwHw �yrk u���B w q�� QO u�J�N CQwYx� �N wO O�vm�t x�Lqt xm Qw�u�tyOW���t �yxvwtv xR�Ov� N u QO xm �CU� � �p

NQ��Q� � uOw� Q�O�vat K�U �� w OvDUy u�v�t��

u�� R� rk Q�O�kt Qo � OQ�o�t Q�Qk xO�iDU� OQwt �QU �oDU�ty u��aD ��Q� u R� xm CU� �Q�DN�U rk Q�O�ktx� O��� �Hv�� QO xD�r� OQm OQ �Q � K�U QO �k � �va� QiY ZQi u�wD�t x�ov �Ot uwQ�� Q�N� �w�N�N wO xrY�i uwQO Q�O�kt xm u swO CU� �vFDUt xOa�k u�� R� r� � � xm u�� pw� CW�O xHwD xDmv xUQ�N��D ��Q� �qFt � K�U QO xm u swU OvQ�Ov �Q�O�vat Cw�iD QiY �� �t� OvW��v QiY ��k�kOCU� umttO�vmv �Qw�O u T�U� Q� w O�v�Ov ��N pwtWt �Q � r� �RH�� uO�Di� uwQ�� �D wO ���yxO�O O�OaD �� U�vDt�Q�ov�oDU�ty QO OyO�t u�Wv xDUy Vy�m l� CQwYx� Q�ov�oDU�ty QO �Q OwN �yxO�O OvwQ �qwtat

Page 19: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� �v�tR ��y�QU � pYi

w CU� �QU QO x��Wt Q�O�kt R� �W�v xm � C�Ft w nQR� Q�O�kt ��� pmW acf�AirPassengers� ��yxO�OO�vm�t x�Lqt �Q �OvDUy syx� l�ORv ��yu�tR QO

0.0 0.5 1.0 1.5

−0.

20.

00.

20.

40.

60.

81.

0

Lag

AC

F

Series AirPassengers

AirPassengers ��yxO�O Q�ov�oDU�ty ���� pmW

QO ��yxO�O �QU xvq�U �w�vD OwW�t Qy�� Q�ov�oDU�ty �wQ �OW�� xDW�O OwHw R�v �rYi C�Q��eD Qo �x���Q xmCU� p�Um� Q�N��D QO �oDU�ty Q�Okt QFm �OL OwW�t x�Lqt �w�vD u�ty �� R�v u Q�ov�oDU�tyQ�DiQ uwJ Tmar�� OyO�t u�Wv �Q Ov�xOW �OH sy R� x�t �� �UwD xm �xt� xt���� GwR u�� ��N C�Ft�� �ivt x���Q xm Ov�xOW �OH sy R� �p�U ��� x�t � CQwYx� Q�O�kt �CU� �Uwv�U ���� QkD �rYi C�Q��eD

CU� u�DUtR lJwm Q�O�kt p��vO x� OyO�t MQ u�DU��D QO xm q�� Q�O�kt p�Ft ��Q� OvQ�O syxm �rYi C�Q��eD w OvwQ hPL R� TB� �v�tR �QU QO �oDU�tyOwN X�NWD Q�ov�oDU�ty �rY� hOyhPL OvwQ xir�wt w K�LYD �rYi C�Q��eD decompose�� ��D R� xO�iDU� �� Q� R O�m QO CU� �Ov�xOW OQwQ�

OOQo�t s�UQD ���� pmW QO Q�ov�oDU�ty w �#�� pmW QO �iO�YD xir�wt �QU TBU OOQo�t

� data�AirPassengers�

� AP �� AirPassengers

� AP�decom �� decompose�AP �multiplicative��

� plot�ts�AP�decom�random���������

� acf�AP�decom�random��������

�rYi C�Q��eD K�LYD w CU� pmW �Uwv�Um �yrk ���Q�t xm OyO�t u�Wv ���� pmW QO Q�ov�oDU�typkDUt CQwYx� �Q x�t xOR�wO ��yxir�wt x� x� RHD �Q� R OUQ�tv CUQO Q�v x� Q��aD u�� CU� xOw�v QF�wt �qt�mJuly R� �rY� ��yxO�O �QU Q��at h�QLv� xm O�vm�t x�Lqt �OwW s�Hv� �QDW�� �UQQ� Qo � Ovm�t OQwQ�C�Q��eD K�LYD R� TB Q��at h�QLv� w CU� � OvwQ hPL R� TB �Q��at h�QLv� CU� � Q��Q� June �D

Page 20: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt ��

Time

ts(A

P.d

ecom

$ran

dom

[7:1

38])

0 20 40 60 80 100 120

0.90

0.95

1.00

1.05

1.10

AirPassengers �iO�YD xir�wt ����� pmW

O�vm xHwD Q� R ��yO�m x� �yQ��at h�QLv� x�U�Lt ��Q� CU� � � Q��Q� �rYi

� AP �� AirPassengers

� AP�decom �� decompose�AP �multiplicative��

� sd�AP��������

��� ��������

� sd�AP������� � AP�decom�trend��������

��� ��������

� sd�AP�decom�random��������

��� ���������

Page 21: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� �v�tR ��y�QU � pYi

0 5 10 15 20

−0.

20.

00.

20.

40.

60.

81.

0

Lag

AC

F

Series AP.decom$random[7:138]

AirPassengers �iO�YD xir�wt Q�ov�oDU�ty ����� pmW

Page 22: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

swO pYi

x��B �iO�YD ��ypOt

���yQorta ��vt OwW�tEL� �OQ�O �Oa� pwYi QO �O�� R OQ� Q�m xm�iO�YD ��ypOt R� ��xQ�B pYi u�� QOpOt QDU�vtX�NWD ��Q� �O�OH Q�R�� w Ovvm�t xO�iDU� �r�N �yv R� �v�tR ��y�QU xm OOQo�t �iQatOOQo�t QDK �w w O����t �QDW�� �U� R�v xDWPo ��yEL� R� ��xQ�B u�vJty OQ�o�t Q�Qk �UQQ� OQwt

�O�iU xiwv ��

xtOkt ���

Ovm�t u��� t u�tR QO �Q pOt w xOy�Wt ��yxO�O �QU u�� Cw�iD R� �W�v ���N �pOt l� QO xOv�t�k��xt �va� �yxOv�t�k�� �QU CQwY u�� QO CU� ZwQit pOt �UwD �yt �va� u xOW OQwQ� Q�O�kt w yt �QU

CU� Q� R CQwYx�

xt � yt � �yt

�� White noise

��

Page 23: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� x �B �iO�YD ��ypOt pYi

Tmavt Q�ov�oDU�ty QO �rYi C�Q��eD �� w OvwQ Ovv�t �v�tR �QU ��yxir�wt xm OW x�Lqt p�k pYi QOOvW�� xDU�ty O���v �yxOv�t�k�� �QU �t� �OvDUy xDU�ty �v�tR ��y�QU ��awv xm OW xO�O u�vJty OwW�tOOQo�t Q� R h� QaD Hwt xO�� u�� OyO u�Wv �Q �YNWt �wor� �wv I�y O���v �yv Q�ov�oDU�ty �Pr w

h� QaD ���

w�� w�� � � � � wn��yQ�eDt Qo � �Ov� wo ��DWN� xDUUo O�iU xiwv �Q fwt � t � �� �� � � � � ng�v�tR�QUxH�Dv u�wD�t pqkDU� C�Y�N R� OvW�� �� Tv�� Q�w w QiY u�ov��t �� u�Um� � RwD ��Q�O w pkDUt sy R�u�wD�t x�ov �Ovm �wQ�B p�tQv p�tDL� �r�oJ ��D R� �yQ�eDt Qo � CU� Cor�wi� wj� � i �� j xm CiQo

CU� wt � N� � ��� xm CW�O �Uw�o O�iU xiwv �QU

R QO R�Ux��W ���

C�kw� ��a� �Q pOt �UwD xOW �R�Ux��W �QU OwW�t xO�iDU� �yxO�O ��R�Ux��W ��Q� xDi�� VR�Q� pOtOOQo R��tDt xOW xOy�Wt ��yxO�O R� �D Ov� wo R�v ��oDN�U ��yxO�O

��Q� u R� u�wD�t w O��tv O�H�� xOv� ��Q� �Q ���yw� Q�vU Ov�wD�t �p�Ft ��Q� CU� �O�it Ov�Qi �R�Ux��Wrnorm�� � p�Ft ��Q� O����t u�t�U �� wN� Q�m u�� R QO Owtv xO�iDU� pOt ��yQDt�Q�B u�v�t�� xR��

Ovm O�H�� �Q OvDUy OQ�Ov�DU� p�tQv u�Um� � RwD ��Q�O xm pkDUt Q�eDt � Ov�wD�tO�vm xHwD Q� R ��yO�m x� uwvm �

� set�seed���

� w �� rnorm�����

� plot�ts�w type �l��

Q�� Qy QO xm �vat u�O� �OW�� u�Um� �R�Ux��W ��QH� Q�� Qy �wQW x�kv xm CU� u ��Q� set�seed�� ��D��y���D� hrDNt ��y�QH� �� OwW hPL ��D u�� Qo � �t� OOQo�t O�H�� u�Um� OOa � xt�vQ� ��QH�

CU� ��� pmW q�� ��yO�m ��QH� xH�Dv OOQo�t pY�L Cw�iDt

Q�ov�oDU�ty w swO x�DQt X�wN ���

Tv�� Q�wwmwD� ��D ��Q�O �OW�� �� Tv�� Q�w w QiY u�ov��t ��Q�O w xDU�ty�v w fwtgCQwY x��iO�YD Q�eDt Qo �CU� Q� R

��k� �

������� k �

k ��

�Q� RCov�wt� wt�k� � E��wt � ���wt�k � ���

�� Discrete white noise �� Simulation �� Synthetic

Page 24: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt ��

Time

w

0 20 40 60 80 100

−2

−1

01

2

O�iU xiwv �v�tR �QU V �tv ��� pmW

OwW�t xH�Dv x�ov �CU� � � E�wt� � E�wt�k� � uwJ

E�wt� wt�k� � E�wt�E�wt�k� �

�����E�w�

t � � ��a k �

� k ��

!R� CUDQ��a �oDU�tyOwN ��D w

�k �ckc�

������ k �

k ��

x� pta QO �t� OvW�� QiY Q��Q� O��� k �� ��R� x� xOW �R�Ux��W O�iU xiwv ��yxO�O �oDU�tyOwN ��D�vat u�O� u�� w OvwW�t k�w u�v�t�� xR�� pN�O �t� �CU�v QiY ��k�kO �yrk Q�O�kt �Q�oxvwtv C�Q��eD CrauOt uwQ�� �Q�N��D � w � uOw� Q�O�vat K�U x� xHwD �� �DL OvQ�Ov QiY �� �Q�O�vat Cw�iD �yv xm CU�

O�vm xHwD Q� R p�Ft x� uwvm � CU� p�mW� q� R�v �yrk R� OQwt l�

� set�seed���

� acf�rnorm������

OQ�O OwHw �a��D R QO xD�r� OQ�O �oORuwQ�� l� � Q�N��D QO xm �CU� ��� pmW q�� ��yO�m ��QH� xH�DvCU� ARMAacf�� ��D u�� Ovm�t x�U�Lt �Q�v CQwYx� �Q O�iU xiwv l� rk �oDU�tyOwN ��D xm

O�vm xHwD Q� R p�Ft x� Owtv s�UQD plot�� ��D �� u�wD�t �Q ��D u�� pY�L

� wn �� ARMAacf�ar� ma� lag�max���

Page 25: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� x �B �iO�YD ��ypOt pYi

0 5 10 15 20

−0.

20.

00.

20.

40.

60.

81.

0

Lag

AC

F

Series rnorm(100)

O�iU xiwv �oDU�tyOwN ���D V �tv �� pmW

� plot�wn type h �

QiY Q��Q� Q�O�kt �yQ�N��D Q��U ��Q� w CU� l� Q��Q� k � QO xm �CU� ��� pmW q�� ��yO�m ��QH� xH�DvCU�

��iO�YD uORs�o ��

xtOkt ���

�Q�o �Qi pOt R� u VR�Q� �DL �qwtat CU� �iO�YD ��yOvwQ ��Q� ��U�vt VR�Q� ���r�e �iO�YD uORs�oCU� QDU�vt ARIMA uwJ

h� QaD ���

!Qo � CU� �iO�YD uORs�o fxtg x�ov CU� �v�tR �QU l� fxtg xm O�vm ZQi

xt � xt�� � wt

�� u�vJty w q�� xrO�at QO xt�� � xt�� � wt�� �v� Ro��H �� CU� O�iU xiwv �QU fwtg u QO xm!xm OwW�t xH�Dv ��wQUB �v� Ro��H� VwQ u�� xt�O� w �r�k QO xt�� �v� Ro��H

xt � wt � wt�� � wt�� � � � ��� Random walk �� Back substitution

Page 26: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

5 10 15 20

0.0

0.2

0.4

0.6

0.8

1.0

Index

wn

O�iU xiwv �Q�v �oDU�tyOwN ���D V �tv �� pmW

!�Pr OwW�t �wQW t � � u�tR R� w CU�v C��yv�� jwi �QU pta QO

xt � w� � w� � � � �� wt

CQwYx� Qortal� u�wD�t u��Q��v� OOQo�t p�ta� R�v xO�J�B �v�tR ��y�QU ��ypOt ��Q� wQUB �v� Ro��HOwtv h� QaD Q� R

�wQUB p�kDv� Qorta ���

CU� Q� R CQwYx� xm B Qorta

B�xt� � xt��

���DQt xOW xDio Qorta Qo � Ov� wo R�v Q�N��D Qorta u x� OQ�wt R� ��xQ�B QO xm �OwW�t xO�t�v wQUB p�kDv� Qorta!x�ov �OOQo Q�QmD

Bn�xt� � xt�n

OwW�t xDWwv Qorta u�� �� �iO�YD uORs�o uwvm �

xt � B�xt� � wt � ���B�xt � wt � xt � ���B���wt

� �� �B �B� � � � � �wt � xt � wt � wt�� � wt�� � � � ��� Backward shift operator

Page 27: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� x �B �iO�YD ��ypOt pYi

�iO�YD uORs�o swO x�DQt X�wN ���

CU� Q� R CQwYx� �x � x� xHwD �� �iO�YD uORs�o Ov�Qi swO x�DQt X�wN

�k�t� � Cov�xt� xt�k� � Cov

�tX

i��

wi�t�kXj��

wj

�Xi�j

Cov�wi� wj�

CU� Q� R CQwYx� �oDU�tyOwN � Q uwvm � CU� �DU���v Ov�Qi u��Q��v� �CU� u�tR R� �a��D Tv�� Q�wwm

�k�t� �Cov�xt� xt�k�pVar�xt� Var�xt�k�

�t��p

t���t� k����

�p� � k�t

l�ORv �k Q�Okt u��Q��v� CU� QDmJwm ��x�Lqt p��k CQwYx� t x� C�Uv k Q�Okt p�LQy� nQR� ��yt ��Q�R� Ov�m �r�N CQwYx� w Ow� Oy�wN C�Ft Q�O�kt ��Q�O �iO�YD uORs�o ��Q� Q�ov�oDU�ty w CU� l� x�

O����t Vy�m l� Q�Okt

p��iD Qorta ���

xm OW�� �iO�YD uORs�o fxtg �QU Qo � p�Ft ��Q� Ovm p�O�D �DU�� �QU x� �Q �DU���v �QU Ov�wD�t Qorta u��xt� xt�� � wt CQwYx� u x���Q xm Ovm�t O�H�� �Q O�iU xiwv �QU fxtg pw� x�DQt p �iD �CU� �DU���v

CU�OwW�t h� QaD Q� R CQwYx� r �va� p �iD Qorta

rxt � xt � xt��

OOQo u��� wQUB p�kDv� Qorta UL Q� Ov�wD�t r u��Q��v� CU� rxt � ���B�xt xm O�W�� xDW�O xHwDCU� Q� R CQwYx� p �iD Qorta QDq�� ��yx�DQt

rn � ���B�n

R�Ux��W ���

OwW�t Qy�� �yQ�Owtv QO pOt �rY� C��YwYN wCU� O�it �v�tR �QUl� xar��t ��Q� �R�Ux��W ���r�exwkr�� u�mt� l� u�wvax� Ov�wD�t pOt �O�O u�Wv OwN R� �Q x��Wt C��YwYN �N� Q�D xO�O xm �Dkw u��Q��v�

OOQo �iQat u ��Q�Ovm �R�Ux��W �Q �iO�YD uORs�o l� Ov�wD�t Q� R ��yO�m

� set�seed���

� x �� w �� rnorm������

� Di�erence operator

Page 28: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt

� for �t in ������� x�t� �� x�t � �� � w�t�

� plot�x type �l��

uORs�o for xkrL OOQo�t �yOQ�Okt x x� ��vt w OyO�t C�Uv w x� �Q O�iU xiwv �QU QwDUO u�rw�xkrL R� u QO xm Owtv O�H�� ���yO�m �� �Q �iO�YD uORs�o �R�Ux��W u�wD�t xD�r� Ovm�t O�H�� �Q �iO�YD

OW�� xOWv xO�iDU�

� set�seed��� � so you can reproduce the results

� v �� rnorm������ � v contains ���� iid N���� variates

� x �� cumsum�v� � x is a random walk

� plot�x type �l��

�QU Q�ov�oDU�ty uOQw CUOx� ��Q� CU� xOW xO�O u�Wv ��� pmW QO xOW O�H�� �QU �wQ Qy�

0 200 400 600 800 1000

020

4060

Index

x

�iO�YD uORs�o V �tv ��� pmW

OwW�t xOy�Wt ��� pmW QO xH�DvTBU OOQo�t �QH� Q� R QwDUO �p�k ��yxt�vQ� xt�O� QO �xOW �R�Ux��W

� acf�x�

OOQo�t O�iU xiwv ��yxO�O �QU x� p�O�D p �iD l� QO �iO�YD uORs�o �OW x�Lqt �q�k xm Qw�u�ty�p �iD Qorta p�ta� ��Q� CU� jO�Y R�v xOW �R�Ux��W �QU OQwt QO xOW xDio r�t xm O�O O��� uwvm ���� pmW QO xH�Dv TBU OwW�t �QH� Q� R QwDUO �Q�N� ��yxt�vQ� xt�O� QO CU� OwHwt R QO di��� ��D

OwW�t xOy�Wt

� acf�di��x��

u�v�t�� xR�� R� �oORuwQ�� OQwt wO OwW�tv x�Lqt Q�ov�oDU�ty ��Q� �YNWt �wor� �Q�N� pmW QO�R�Ux��W �QUCiQ�t Q��Dv� xm Qw�uty u��Q��v� CU�v �vDa� p��k � uOw� Q�O�vat K�U QO xm OQ�O OwHw

Page 29: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

x �B �iO�YD ��ypOt pYi

0 5 10 15 20 25 30

0.0

0.2

0.4

0.6

0.8

1.0

Lag

AC

F

Series x

�iO�YD uORs�o Q�ov�oDU�tyOwN V �tv ��� pmW

Ovm�t �wQ�B �iO�YD uORs�o Ov�Qi l� R� x xOW

w�UQoQwD� ��ypOt ��

h� QaD ���

Qo � �Ov� wo AR�p� Q�YDN� x� p x�DQt w�UQoQwD� Ov�Qi �Q fxtg �v�tR �QU

xt � �xt�� � �xt�� � � � �� pxt�p � wt

p�kDv� Qorta UL Q� Q�N� xrO�at Qo � OvDUy p �� �� pOt ��yQDt�Q�B i w O�iU xiwv fwtg u QO xmOw� Oy�wN Q� R CQwYx� �OOQo u��� wQUB

p�B�xt � ��� �B � �B� � � � � � pB

p�xt � wt

O�vm xHwD Q� R C�mv x� uwvm �

OW���t � � � u QO xm CU� AR��� X�N Cr�L �iO�YD uORs�o w�UQoQwD� u x� �U u�ty x� w CU� xDWPo ��yu�tR QO �yx x� C�Uv xt uw�UQoQ k�w QO pOt u��

Ov� wo

!�� CU� Q��Q� t u�tR QO OQwQ�

�xt � �xt�� � �xt�� � � � �� t�pxt

Page 30: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

0 5 10 15 20 25 30

0.0

0.2

0.4

0.6

0.8

1.0

Lag

AC

F

Series diff(x)

p��iD p�ta� �� �iO�YD uORs�o Q�ov�oDU�tyOwN V �tv ��� pmW

OOQo�t OQwQ� ��N C�a� Qt pk�OL VwQ x� pOt ��yQDt�Q�B

w�UQoQwD� QO ���DU���v w ���DU�� ���

QOk Qo � OW�� �rDNt �� w �k�kL ��yxW� Q ��Q�O Ov�wD�t xm Ovt�v xYNWt xrO�at �Q p�B� � xrO�atuORs�o Ov�Qi QO CU� �DU�� Ov�Qi x�ov �OOQo OL�w R� QDW�� p�B� ��xrtHOvJ ��yxW� Q s�tD jr�tu��aD ��Q� Q� R p�Ft Q�yJ x� uwvm � CU� �DU���v Ov�Qi w OwW�t B � � xm CU� � � � B �iO�YD

O�vm xHwD AR Ov�Qi ���DU���v w ���DU��

B � � Q��Q� �� ��B � xrO�at xW� Q �Q� R �CU� �DU�� xt � �

�xt�� �wt CQwYx� AR��� pOt ��

CU� QDoQR� OL�w R� xm CU�

wQUB p�kDv� Qorta R� xO�iDU� �� �Q� R CU� �DU�� xt � xt�� � �xt�� � wt CQwYx� AR��� pOt ��

��xrtHOvJ u�� xW� Q OwW�t xH�Dv ��B� ���xt � wt �va� ��B

�� B��xt � wt xrO�atOW���t OL�w R� QDoQR� w CU� B � � xm O��t CUOx� �B� � �

�B � ��� xrO�at pL ��

CU� OL�w �yxW� Q R� �m� �Q� R �CU� �DU���v xt � ��xt�� �

��xt�� � wt CQwYx� AR��� pOt ��

����B � ���B � ��xt � wt �va� ��

��B� �B � ��xt � wt wQUB p�kDv� Qorta R� xO�iDU� ��

CU� B � ���� ��yxW� Q ��Q�O �B� � ����B � ���B � �� ��xrtHOvJ OwW�t xH�Dv

CU� OL�w R� QDoQR� B � �� jr�t QOk �Q� R OOQo�t ���DU���v �U B � � Q�Okt �ki

�va� � � �B� � ��yxW� Q �Q� R �CU� �DU�� xt � ��

xt�� � wt CQwYx� AR��� pOt ��

Page 31: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� x �B �iO�YD ��ypOt pYi

� Q��Q� �yxW� Q R� l� Qy jr�t QOk OW���t i �p�� w CU� �rDNt O�Oa� xm OvDUy B � ��i

CU� OL�w R� QDoQR� xm OW���t

xrO�at ��yxW� Q u�wD�t xm OQ�O OwHw polyroot�� s�v x� �a��D R QO ��xrtHOvJ ��yxW� Q x�U�Lt ��Q�xOW xDio ��D �� sQ�yJ w swU OQ�wt �p�Ft u�wvax� OQm �UQQ� �Q �y�QU ���DU�� w Owtv x�U�Lt �Q xYNWt

OwW�t x�U�Lt!swU p�Ft

� polyroot�c�������

��� ���i ����i

!R� OvDQ��a �yxW� Q jr�t QOk

� Mod�polyroot�c��������

��� � �

!sQ�yJ p�Ft

� polyroot�c��������

��� ���i ���i

!R� OvDQ��a �yxW� Q jr�t QOk

� Mod�polyroot�c���������

��� � �

AR��� pOt swO x�DQt X�wN ���

CU� Q� R CQwYx� AR��� pOt �OW xOy�Wt �q�k xm Qw�u�ty

xt � xt�� � wt

xm O�O u�Wv u�wD�t OW���t �� Tv�� Q�w w QiY u�ov��t �� O�iU xiwv fwtg u QO xm

�x �

�k � k������ ��

xm CWwv u�wD�t jj � � Ov�Qi ���DU�� �QW w B Qorta R� xO�iDU� ��

��� B�xt � wt

Owtv x�U�Lt �Q xt u�wD�t uwvm �

xt � ��� B���wt

� wt � wt�� � �wt�� � � � � ��Xi��

iwt�i

Page 32: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

O�vm xHwD u�ov��t x�U�Lt x� uwvm �

E�xt� � E

� �Xi��

iwt�i

�Xi��

iE�wt�i� �

CU� Q� R KQW x� Tv�� Q�wwmwD�

�k � Cov�xt� xt�k� � Cov

� �Xi��

iwt�i��Xj��

jwt�k�j

�X

j�k�iij Cov�wt�i� wt�k�j�

� k��X�

i���i � k������ ��

AR��� Ov��Qi Q�ov�oDU�ty ���

O�vm s�UkD Tv�� Q�w Q� �Q Tv�� Q�wwmwD� xm CU� �i�m �oDU�tyOwN ��D x�U�Lt ��Q�

�k � k �k �

p�t QiY CtU x� �QDW�� CaQU �� QDmJwm ��y ��R� x� Q�ov�oDU�ty u��Q��v� CU� jj � � u QO xmO��t �� pmW QO �ivt w C�Ft Q�O�kt ��Q� Q�ov�oDU�ty wO p�Ft xt�O� QO Ovm�t

� plot����� rho����� ����� type �b� xlab�lag k��

� rho �� function�k alpha� alpha�k

� par�mfrowc���� marc�������

� plot����� rho����� ���� type �b� xlab�lag k�

� ylabexpression�r�k�� mainexpression�alpha�����

� plot����� rho����� ����� type �b� xlab�lag k�

� ylabexpression�r�k�� mainexpression�alpha������

� abline�h� lty��

� �� RH �oDU�tyOwN ���D ���

x� �ki xt Qo � �DL �OvDUy QiY hr�Nt �yQ�N��D s�tD ��Q� �y�oDU�tyOwN Q�O�kt ��k Q�N� xrO�at j���t�y�oDU�ty QF� xm CU� �oDU�ty R� �awv �k Q�N��D QO �� RH �oDU�tyOwN ��D OW�� xDW�O �oDU� xt��Ow� Oy�wN QiY k � � s�tD ��Q� AR��� �� RH �oDU�tyOwN ��D p�Ft ��Q� OOQo hPL QDmJwm Q�N��D ���CU� AR�k� pOt xDi�� VR�Q� � Q u�t�k �� Q��Q� k Q�N��D QO �� RH �oDU�tyOwN ��D ��rm Cr�L QO�oDU�tyOwN ��D Q�ov�oDU�ty u��Q��v� CU� QiY Q��Q� k � p ��R� x� k ��Q AR�p� Ov�Qi QO �va�w x�U�Lt ��Q� pacf�� s�v x� �a��D R QO OQ�o Q�Qk xO�iDU� OQwt Ov�wD�t AR Ov�Qi x�DQt u��aD ��Q� �� RH

OQ�O OwHw �� RH �oDU�tyOwN ��D s�UQD� Partial autocorrelation function

Page 33: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� x �B �iO�YD ��ypOt pYi

0 2 4 6 8 10

0.0

0.4

0.8

α = 0.7

lag kr k

0 2 4 6 8 10

−0.

50.

00.

51.

0

α = − 0.7

lag k

r k

� � �������� ��Q� AR��� Ov �Qi Q�ov�oDU�ty V �tv ��� pmW

R�Ux��W ���

QO �� RH �oDU�tyOwN w �oDU�tyOwN ��wD w Ov�Qi w xOW �R�Ux��W Q� R CQwYx� R QO AR��� pOt�oDU�tyOwN Q�ov�oDU�ty QO k � � Q�O�kt ��Q� O�vm�t x�Lqt xm Qw�u�ty OOQo�t s�UQD "�� pmW

OQ�Ov OwHw �Q�O�vat �oDU�ty ��� RH

xDi�� VR�Q� �ypOt ���

xOW R�Ux��W �yxO�O Q� xDi�� VR�Q� �ypOt ���

�QU ��Q� x�ar uw�UQoQwD� pOt �Q� R ��yO�m QO O����t VR�Q� �yxO�O Q� ar�� ��D �UwD R QO AR�p� pOtR� xO�iDU� �� u ���v�HtTv�� Q�w xm �O����t VR�Q� xOW xO�O QDt�Q�B w � u�v�t�� xR�� �� xOW �R�Ux��W

OwW�t G�QNDU� x�ar�asy�var

� set�seed���

� x �� w �� rnorm�����

� for �t in ������ x�t� �� ��� � x�t � �� � w�t�

� x�ar �� ar�x method �mle��

� x�ar�order

pOt x�DQ

��� �

� x�ar�ar

pOt QDt�Q�B

��� asymtotic

Page 34: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

0 20 40 60 80 100

−1

01

23

Indexx

0 5 10 15 20

−0.

20.

20.

61.

0

Lag

AC

F

5 10 15 20

−0.

20.

20.

40.

6

Lag

Par

tial A

CF

V �y�oDU�tyOwN w AR��� Ov �Qi V �tv ��� pmW

��� ���������

� x�ar�ar � c��� �� � sqrt�x�ar�asy�var�

pOt QDt�Q�B u�v�t�� xR��

��� ��������� ���������

��NDv� CU� Q�wDU� ���tvDUQO QFm �OL ��D T�U� Q� �OW xO�iDU� �yO�m QO VR�Q� Ov�Qi QO xm ��mle VwQx��� x�U�Lt CU� pOt QDt�Q�B pk�OL Qov��� xm OwW�t s�Hv� ��AIC x��� R� �Ov�Qi x� �w� Qt p Q�Okt

CU� Q� R CQwYx� xOW xDio

AIC � ��� log�likelihood � �� number of parameters

CUQO x�DQt �q�� ��yO�m xm O�vm xHwD Ovm�t Q��DN� �Q VR�Q� u� QDy� AIC u� QDmJwm �� pOt �ar�� ��D QOQ�Okt R� QDmJwm xm �Ot CUOx� � � �� Q��Q� AR��� pOt QDt�Q�B OQwQ� Owtv ���� R�� �Q p � � �va� pOtOQ�Ov OwHw E�L u�� R� �O�Q�� w OOQo�t pOt QDt�Q�B Q�Okt pt�W � u�v�t�� xR�� �t� CU� � ��

xDi�� VR�Q� AR pOt ��v�yH �tO QU ���

�aO� u�wD�t O�OQD �� CU� ��xv�Nro QF� R� �W�v xm OyO�t u�Wv �Q �tO OvwQ V��Ri� �� p�U R� �QU u���QU �� �aO� u�� �Q�oR�U ��Q� CU� xOW hPL �iO�YD xO�OB l� p�kDv� R� �W�v ��xOv��Ri OvwQ xm Owtv��� Maximum likelihood estimation ��� Akaike Information Criterion

Page 35: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� x �B �iO�YD ��ypOt pYi

Owtv xO�iDU� u�at ��wD R� xm u uwO� xOwtv �R�UpOt �Q OvwQ xir�wt CU� umtt ��yxO�OO�vm xHwD xvq�U ��tO �QU u�ov��t x� xDi�� VR�Q� AR pOt x�

� www �� �http���www�massey�ac�nz��pscowper�ts�global�dat�

� Global �� scan�www�

� Global�ts �� ts�Global st c����� �� end c����� ��� fr ���

� Global�ar �� ar�aggregate�Global�ts FUN mean� method �mle��

� mean�aggregate�Global�ts FUN mean��

��� ����������

� Global�ar�order

��� �

� Global�ar�ar

��� ���������� ���������� ���������� ����������

� acf�Global�ar�res�����Global�ar�order�� lag �� main���

0 10 20 30 40 50

−0.

20.

00.

20.

40.

60.

81.

0

Lag

AC

F

�yxOv�t�k�� Q�ov�oDU�ty V �tv ��� pmW

Q�N��D QO rk Q�Okt RH x� �Q� R OyO�t u�Wv �Q O�iU xiwv �QU l� ��yxOv�t�k�� ���� pmW� Q�ov�oDU�tyOUQ�t Q�v x� U�vt ��yxO�O Q� AR�� pOt VR�Q� �Pr OvQ�O Q�Qk u�v�t�� xR�� QO �yv x�k� ��

xO�O CUOx� �Q pOt ��yQDt�Q�B w xvq�U �UwDt ��tO u�ov��t Q�Okt xm Q�N� ��yO�m G��Dv x� xHwD �� uwvm �CWwv u�wD�t �CU�

�xt � � ��� ���xt��� ���� � ���xt��� ���� ����xt��� ����� ����xt�� ����

Page 36: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

swU pYi

�DU�� ��ypOt

��O�m � �Q fxtg �v�tR �QU OwW�t �Q�o�B EL� u xt�O� uwvm � OW KQ�t ��r��t ���DU�� x� H�Q �q�k�va� �Ovmv Q��eD u�tR QO p�kDv� QF� QO �QU s�wD � RwD ��D Qo � Ov� wo ��DU��

f �xt� � xt� � � � � � xtk� � f �xt��h� xt��h� � � � � xtk�h�

Cov�xt� xs�Tv�� Q�wwm w OwW�t pt�W R�v �Q u�tR QO C��FTv�� Q�w w u�ov��t �Om�wt ���DU�� xm O�vm xHwDOW�� xDW�O �oDU� k � jt� sj Q�N��D x� �ki

lQLDt u�ov��t ��ypOt ��

X�wN w h� QaD �MA�q� Ov��Qi ���

xrtH q u� QN u�vJty w O�iU xiwv �Q�H xrtH R� ��N �mQD q x�DQt R� �MA� lQLDt u�ov��t Ov�QiCU� Q� R CQwYx� u x���Q OW���t O�iU xiwv �r�k

xt � wt � ��wt�� � � � �� �qwt�q�� Strictly stationary

Page 37: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� �DU � ��ypOt pYi

p�kDv� Qorta UL Q� u�wD�t �Q q�� xrO�at CU� ��w Tv�� Q�w w QiY u�ov��t �� O�iU xiwv fwtg u QO xm

Owtv �U� wvR�� �B wQUB

xt � �� � ��B � ��B� � � � �� �qB

q�wt � q�B�wt

��DU�� CqtH R� �y�vDt tH pt�WMA ��yOv�Qi xm u�� x� Q�v CU� q xHQO R� ��xrtHOvJ q u QO xmOQ�O OwHw u�tR x� C�Uv Tv�� Q�wwmwD� w u�ov��t QO C��F w CU� Q�QkQ� ���DU�� �Pr �CU� O�iU xiwv

CU� x�U�Lt p��k �oO�U x� fxtgTv�� Q�w w u�ov��t

E�xt� � E�wt � ��wt�� � � � �� �qwt�q�

� E�wt� � ��E�wt��� � � � �� �qE�wt�q�

OvDUy QiY Q��Q� �yu�ov��t R� l� Qy �Q� R

Var�xt� � Var�wt � ��wt�� � � � �� �qwt�q�

� Var�wt� � ��� Var�wt��� � � � �� ��

q Var�wt�q�

� ��w�� � ��

� � � � �� ��q �

OvDUy pkDUt Qo�Om� R� w ��w u�Um� Tv�� Q�w ��Q�O CqtH s�tD �Q� R

CU� Q� R QwYx� k ��Q� �oDU�tyOwN ��D

��k� �

�����������������

� k � q�kP

i��

�i�i�k

qP

i��

��i

k � �� � � � � q

k � q

O�iU xiwv pkDUt CqtH pt�W xt�� w xt �Q� R �CU� QiY Q��Q� ��D k � q ��Q� CU� �� � � u QO xmCU� QiY �yv Tv�� Q�wwm u�vJty w OvDUy

xrtH uwO� �y�vDt�v x�DQt �� �DU�� w�UQoQwD� Ov�Qi UL Q� �Q u u�wD� Qo � Ov� wo �Q�PBuwQ�w �Q MA Ov�QiOwW u��� Q� R CQwYx� Ov�wD�t xt � ��� �B�wt x���Q �� MA��� Ov�Qi �p�Ft ��Q� CWwv ��N

wt � ��� �B���xt � xt � �xt�� � ��xt�� � � � �

Ovm jOY ���Qoty �QW �D CU� j�j � � u QO xmOL�w R� QDoQR� �oty �B� ��yxW� Q jr�t QOk Qo � CU� Q�PBuwQ�w Ov�Qi l� MA�q� Ov�Qi ��rm Qw�x�

OvW���� invertible

Page 38: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt

R�Ux��W w Q�ov�oDU�ty �R �yp�Ft ���

w xkrL R� xO�iDU� �� �a��D Q�m u�� ��Q� Owtv xO��B R QO u�wD�t �Q MA�q� Ov�Qi ��Q� �oDU�tyOwN ��DOwW�t h� QaD �QW

� rho �� function�k beta� f

� q �� length�beta� � �

� if �k � q� ACF �� � else f

� s� �� �� s� �� �

� for �i in ���q�k���� s� �� s� � beta�i� � beta�i�k�

� for �i in ���q���� s� �� s� � beta�i���

� ACF �� s� � s�g

� ACFg

Ov�Qi ��Q� Q� R ��yO�m CU� x�U�Lt p��k MA�q� Ov�Qi ��Q� �oDU�tyOwN ��D �q�� ��yO�m R� xO�iDU� ��CU� ��� pmW CQwY x� �� � �� w �� �� ��� � �� ��yQDt�Q�B �� MA���

� beta �� c�� ��� ��� ����

� rho�k �� rep�� ���

� for �k in ����� rho�k�k� �� rho�k beta�

� plot����� c�� rho�k� pch � ylab expression�rho�k�� xlab�lag k��

� abline�� ��

0 2 4 6 8 10

0.0

0.2

0.4

0.6

0.8

1.0

lag k

ρ k

MA�� pOt Q�ov�oDU�ty V �tv ��� pmW

Page 39: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�DU � ��ypOt pYi

w �v�tR �QU OQ�o�t Q�Qk xO�iDU� OQwt Q�ov�oDU�ty s�UQD w MA��� Ov�Qi �R�Ux��W ��Q� Q� R ��yO�m�oDU�tyOwN ��D Q�Okt xU u�rw� �CiQ�t Q��Dv� xm Qw�u�ty CU� xOW s�UQD ��� pmW QO Q�ov�oDU�tyOvQ�O Q�Qk u�v�t�� xR�� QO �yrk Q�O�kt �� R� QDW�� ��yQ�N��D ��Q� w OvDUy QiY R� �Q�O�vat Cw�iD ��Q�O

CU� �Q�oxvwtv C�Q��eD R� �W�v xm OW�� xDW�O OwHw �oORuwQ�� CU� umtt sy � xD�r�

� set�seed���

� b �� c���� ��� ����

� x �� w �� rnorm������

� for �t in ������� f

� for �j in ���� x�t� �� x�t� � b�j� � w�t � j�

� g

� par�mfrowc���� marc�������

� plot�x type �l� xlab�time t��

� acf�x main���

0 200 400 600 800 1000

−4

−2

02

4

time t

x

0 5 10 15 20 25 30

0.0

0.4

0.8

Lag

AC

F

u� Q�ov�oDU�ty w MA�� �R�Ux��W Ov �Qi V �tv �� pmW

Page 40: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

xDi� VR�Q�MA ��ypOt ��

xOW R�Ux��W QU x� xDi�� VR�Q� pOt ���

u�twoQ �� pOt x�DQt xm OQ�O OwHw arima�� s�v x� �a��D R QO O��� VR�Q� �yxO�O Q� Ov�wD�t MA�q� pOtQUm �Q u�ov��t ZQiV�B CQwYx� xOW xDio ��D �ar�� ��D hqN Q� OwW�t s��vD order�c� � �q�

OOQo�t OQwQ� ���O�t R� ZQa xrtH QO w Ovm�tvC�a� Qt �wtHt ��yQDt�Q�B OQwQ� ��Q� arima�� ��D OvwW�t OQwQ� �OOa sD� Qwor�l��UwD pOt ��yQDt�Q�B

OW�� method�c��CSS�� u�twoQ xm �DQwY QO xD�r� OR�U�t pk�OL �Q ��QW��yQDt�Q�B ��R� x� xm CU� CQwY u�� x� MA�q� Ov�Qi VR�Q� ��Q� ��QW C�a� Qt �wtHt sD� Qwor� K� wDQ� R CQwYx� �wt �va� VOQwQ� w wt �yxOv�t�k�� OOQo�t x�U�Lt QQmt Qw�x� �yxOv�t�k�� C�a� Qt �wtHt

CU�

S� ���� � � � � ��q� �nXt��

�w�t �

nXt��

nxt � � ��� �wt�� � � � �� ��q �wt�q�

o�

R� �Q�O�kt ��OOa �wHDUH OvW�� QiY �� Q��Q� ��yv Q�QmD �wQW QO �w�� � � � � �wt�q Q�O�kt xm u Q� �wQWtOOQo pk�OL q�� C�a� Qt �wtHt xm Ovm�t u�at �Q �yQDt�Q�B

��yQDt�Q�B O����t VR�Q� Q�N� xOW �R�Ux��W ��yxO�O x� x�ma �lQLDt u�ov��t pOt �Q� R ��yO�m QOpt�W �coeff�� � s�e� of coeff� ���� QkD� OQ�o�t Q�Qk � u�v�t�� xR�� QO xt�vQ� �HwQN QO xOW OQwQ�Q��Dv� xm u�vJty ��vt CU� xOW xO�iDU� �yv R� �R�Ux��W QO xm OwW�t � ��� ��� �� ��yQDt�Q�B Q�O�kt

OQ�Ov QiY �� �Q�O�vat Cw�iD ��O�t R� ZQa CiQ�t��D QO include�mean�FALSE u�twoQ R� Q�m u�� ��Q� Owtv s��vD QiY Q��Q� �Q u�ov��t Q�Okt u�wD�t

CU� include�mean�TRUE ZQiV�B u��Q��v� OwW�t xO�iDU� arima��

� set�seed���

� b �� c���� ��� ����

� x �� w �� rnorm������

� for �t in ������� f

� for �j in ���� x�t� �� x�t� � b�j� � w�t � j�

� g

� x�ma �� arima�x order c�� � ���

� x�ma

Call�

arima�x � x� order � c��� �� ���

�� Intercept

Page 41: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� �DU � ��ypOt pYi

Coefficients�

ma� ma ma� intercept

����� ������ ������ �����

s�e� ����� ������ ����� ������

sigma� estimated as ������ log likelihood � �������� aic � ������

ARMA Ov�Qi ���mQD ��ypOt ��

h� QaD ���

CU� Q� R CQwYx� p x�DQt R� w�UQoQwD� Ov�Qi fxtg �QU �OW xOy�Wt �q�k xm Qw�u�ty

xt � �xt�� � �xt�� � � � �� pxt�p � wt

MA w AR Ov�Qi xm �Dkw �ypOt R� �O�it xDUO OvDUy pOt ��yQDt�Q�B �yi w O�iU xiwv fwtg u QO xm�ARMA� lQLDt u�ov��t w�UQoQwD� Ov�Qi ��Q�O fxtg �v�tR �QU OOQo�t pY�L �OvwW�t xOwRi� sy x�

CU� Q� R CQwYx� u x���Q OwW�t xO�O u�Wv ARMA�p�q� �� w CU� �p�q� x�DQt R�

xt � �xt�� � �xt�� � � � �� pxt�p � wt � ��wt�� � ��wt�� � � � �� �qwt�q

CWwv wQUB p�kDv� Qorta CQwYx� u�wD�t �Q jwi x���Q CU� O�iU xiwv fwtg u QO xm

p�B�xt � q�B�wt

CU� xHwD QwNQO ARMA�p�q� OQwt QO Q� R OQ�wt

OvW�� OL�w R� QDoQR� ��yxW� Q s�tD jr�t QOk Qo � CU� �DU�� Ov�Qi ��

OvW�� OL�w R� QDoQR� ��yxW� Q s�tD jr�t QOk Qo � CU� Q�PBuwQ�w Ov�Qi ��

CU� ARMA�p� � pOt X�N Cr�L AR�p� pOt ��

CU� ARMA� �q� pOt X�N Cr�L MA�q� pOt ��

�yQDt�Q�B XwYN QO �QDW�� ��Q�m ���r�e ARMA pOt �OwW�t s�Hv� VR�Q� �Dkw !QDt�Q�B pk�OL ��OQ�O ���yvD x� AR �� w MA ��ypOt x� C�Uv

Page 42: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

pOt swO x�DQt X�wN ���

CqtHUL Q� �Q u �OD�� QO xm CU� QDy� �fxtg �v�tR �QU ARMA�p�q� swO x�DQt X�wN u��aD ��Q�O� Q�o� Q�v QO �Q ARMA����� pOt p�Ft ��Q� OvDUy pkDUt sy R� �yv �Q� R O�U� wv� O�iU xiwv

xt � xt�� � wt � �wt��

Q� �Q Q�N� xrO�at CU� Var�wt� � ��w w E�wt� u Tv�� Q�w w u�ov��t CU� O�iU xiwv wt u QO xm

O� QO O�iU xiwv ��yxir�wt UL Q� xm xOwtv DQt �Qw� xt T�U�

xt � ��� B����� � �B�wt

CU� Q� R CQwYx� q�� x���Q CU�Q hQ� �U�

xt � �� � B � �B� � � � � ��� � �B�wt

� �Xi��

iBi

�� � �B�wt

�� �

�Xi��

i��Bi�� ��Xi��

i�Bi��

wt

� wt � � � ���Xi��

i��wt�i

Tv�� Q�w w CU� E�wt�i� � �yi s�tD ��R� x� �Q� R �CU� QiY Q��Q� E�xt� u�ov��t �q�� x���Q x� xHwD ��

Var�xt� � Var

wt � � � ��

�Xi��

i��wt�i

� ��w � ��

w�� ���� � ����

OwW�t xH�Dv k � ��Q� �k Tv�� Q�wwmwD� ��Q�

Cov�xt� xt�k� � � � ��k����w � � � �����

wk�Xi��

�i��

� � � ��k����w � � � �����

wk��� ����

OwW�t xH�Dv �k �oDU�tyOwN ��D ��Q�

�k � �k��� � Cov�xt� xt�k��Var�xt�

�k��� � ���� � ��

� � � � ��

Page 43: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� �DU � ��ypOt pYi

�k � �k�� CiQo xH�Dv u�wD�t Q�N� x���Q R�

�� QHDp�rLD ARMA ��ypOt ��

VR�Q� w R�Ux��W ���

��D �UwD u�wD�t �Q OwW�t EL� �Oa� pYi QO xm ARIMA ��yOv�Qi u R� QD�twta w ARMA Ov�Qic�p� �q� x�DQt �� arima�� ��D �UwD Ov�wD�t ARMA�p�q� pOt Owtv �R�Ux��W R QO arima�sim��

O�O VR�Q�pOt TBU �OwW�t �R�Ux��W � � �� w � � �� �� ARMA����� Ov�Qi ��Q� �yxO�O Q� R ��yO�m QOl�ORv xvwtv R� � w ��yQDt�Q�B OQwQ� CiQ�t Q��Dv� xm Qw�u�ty O����t VR�Q� �yv Q� ARMA�����

CU� q�� QO �yv Q�O�kt

� set�seed���

� x �� arima�sim�n ����� list�ar ���� ma �����

� coef�arima�x order c�� � ����

ar� ma� intercept

����������� ��������� �����������

p�O�D MQv QU ���

x��� w O����t VR�Q� p�O�D MQv �QU Q� ARMA����� w AR��� �MA��� ��ypOt �Q� R ��yO�m QO��� pmW OW���t �yxO�O ��Q� �QDU�vt pOt ARMA����� Ov�Qi OvOQo�t xU��kt sy �� �yv AICw �CU� lJwm ��y�oDU�tyOwN ��Q�O xm OyO�t u�Wv �Q ARMA����� pOt xOv�t�k�� Q�ov�oDU�ty

Ovm�t O���D �Q pOt u�� R� xO�iDU� w OvDUy O�iU xiwv �� Q�oR�U

� www �� �http���www�massey�ac�nz��pscowper�ts�pounds�nz�dat�

� x �� read�table�www header T�

� x�ts �� ts�x st ���� fr ��

� x�ma �� arima�x�ts order c�� � ���

� x�ar �� arima�x�ts order c�� � ���

� x�arma �� arima�x�ts order c�� � ���

� AIC�x�ma�

��� ���������

� AIC�x�ar�

��� ���������

� AIC�x�arma�

��� ���������

� x�arma

Page 44: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

Call�

arima�x � x�ts� order � c��� �� ���

Coefficients�

ar� ma� intercept

���� ���� ����

s�e� ���� ��� ����

sigma� estimated as ������� log likelihood � ���� aic � ����

� acf�resid�x�arma��

0 1 2 3

−0.

20.

00.

20.

40.

60.

81.

0

Lag

AC

F

p O�D MQv �QU ��Q� ARMA���� �QU xOv�t�k�� Q�ov�oDU�ty �� pmW

Page 45: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

sQ�yJ pYi

�DU���v ��ypOt

OvwQ w �rYi C�Q�eD R� QF��Dt �Q� R �OvDUy �DU���v �y�QU R� �Q��U� Ot xDWPo ��ypYi QO xmQw�u�tyu�� QO OW p�O�D �DU�� �QU x� p �iD Q��l� �� xm Ow� �yv R� �m� �iO�YD uORs�o p�Ft u�wvax� OvW���t

CU� lQLDt u�ov��t w w�UQoQwD� CqtH pt�W xm O����t �U� �iO�YD uORs�o pOt pYiOwW�t xO�t�v �ARIMA Q�YDN� x� xm �CU� ��oJQ�Bm� �� w ��tHD G�DLt �r �iD �QU

SARIMA CQwYx� �Q u u�wD�t w OR�U�t �DU���v �Q u xm CU� �rYi CqtH pt�W ARIMA Ov�QiOW���t �v�tR ��y�QU R� �Q��U� p�rLD QO �Ovtv�wv Q�R�� �rYi ARIMA ��ypOt O�O u�Wv

��W �r�t �v�tR ��y�QU QO r�t u�� Ovm�t Q�eD Kw w Qw�x� Tv�� Q�w �DU���v ��y�QU R� ��a� QOw�UQoQwD� pOt ��y�QU �wv u�� �R�UpOt ��Q� �yOQm� wQ R� �m� OQ�O RwQ� u�mt� R�v �t�rk� ��yOQwmQ QO CU��GARCH Q�YDN� x� u xDi�� s�taD �wv OwW�t xO�t�v �ARCH pOt Q�YDN� x� w CU� Tv�� Q�w ��Q�

OwW�t xO�t�v

�� aggregated �� integrated �� AutoRegressive Integreted Moving Avrage �� AutoRegressive

Conditional Heteroskedastic �� Generalised AutoRegressive Conditional Heteroskedastic

Page 46: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt ��

�rYi Q�e ARIMA pOt ��

jQ� O�rwD QU w p��iD ���

x�wN w �iO�YD uORs�o pFt OW�� �iO�YD x�wN �Ovm hPL �Q OvwQ Ov�wD�t fxtg �QU Q� p �iD p�ta�p �iD u�rw� �xt � Xt�� � wt �va� �iO�YD uORs�o QO CU� ��N OvwQ xm �OQwt pFt OW�� u�at

OW���t �DU�� xm CU� rxt � xt � xt�� � wt

rxt � xt�xt�� x�ov �CU� O�iU xiwv ��y��N ����N OvwQl� QO xm xt � a�bt�wt Ov�Qi p��kt QOVR�Q� R QO arima�� ��D O�O u�Wv MA��� �� �Q u u�wD�t xm CU� �DU�� lQLDt u�ov��t Ov�Qi l��O�yO VR�Q� R �UwD u�at OvwQ Q� �r �iD pOt O�y�wN� �tW Qo � OQ�tW�tv R�Ht C��F Q�Okt �� �Q �ypOtxO�iDU� �� �r �iD �QU Q� �Q ARMA pOt u�wD�t xD�r� OwW�t s��vD QiY u�ov��t �� �r �iD �QU x�ov

O�O VR�Q� d� w include�mean�FALSE w arima�� ��D R�rxt � yt � b�wt�wt�� u��Q��v� �OW�� xt � a� bt�wt xm O�vm ZQi ARIMA pOt �� ���DQ� QOQ� MA��� Qo � O��t CUOx� xt � x� �

Pti�� yi x���Q yt � rxt Twmat h� QaD R� xO�iDU� �� CU�

xOW �R�Ux��W fxtg Tv�� Q�w u��Q��v� CU� �� Q��Q� wt�� � Q �O��� VR�Q� fytg �r �iD �v�tR �QUCi�� Oy�wN V��Ri� s�kDUt �N pwL

O�vm xOy�Wt Q� R QO O�v�wD�t �Q Q�N� h�Qo �Q�B C��� RH

yt � rxt� xt � xt��

� a� bt � wt � fa� b�t� �� � wt��g� b� wt � wt��

� x� �tX

i��

yi � x� �tX

i��

�b� wi � wt���

� x� � bt � wt � xt

Page 47: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� �DU ��v ��ypOt � pYi

OOQo�t x� � a w CU� w� � xm O�W�� xDW�O xHwDOwW�t xH�Dv xt�� ��Q� wQUB �v� Ro��H ��

xt � xt�� � b � wt � �wt��

� xt�� � b � wt�� � �wt�� � b� wt � �wt��

� xt�� � �b� wt � �� � ��wt�� � �wt��

� xt�� � �b� wt � �� � ��wt�� � �� � ��wt�� � �wt��

���

� x� � bt � wt � �� � ��tX

i��

wt

!R� CUDQ��a Tv�� Q�wVar�xt� � ��

wf� � �� � ����t� ��gQ��Q� Tv�� Q�w CQwY u�� QO xm OW�� � � �� xm u�� Qot �O����t V��Ri� R�v Vv�� Q�w t Q�Okt V��Ri� �� xm

CU� ��w

sD� Q�or ��yxO�O p�Ft ��Q� OQ�o�t CQwY di� �UwD R QO p �iD u�rw� �OW x�Lqt �q�k xm Qw�u�ty�QU ��v�tR �QU � ��� pmW� Q�Owtv xU Q� R ��yO�m O� Q�o� Q�v QO �Q xDU� QDmr� O�rwD x� �w� Qt �a���p�Ft �r �iD �QU wO QO OwW�t xOy�Wt xm Qw�u�ty OyO�t u�Wv �Q sD� Q�or �r �iD �QU w �r �iD

OwW�tv Qy�� Qo�O OvwQ xir�wt �xOW xDio

� www �� �http���www�massey�ac�nz��pscowper�ts�cbe�dat�

� CBE �� read�table�www he T�

� Elec�ts �� ts�CBE� �� start ���� freq ���

� par�mfrowc���� marc�������

� plot�Elec�ts�

� plot�di��Elec�ts��

� plot�di��log�Elec�ts���

��k�irD pOt ���

�� w Ov� wo d x�DQt R� �k�irD pOt �Q u x�ov �OOQo fwtg O�iU xiwv p �iD u�t�d R� TB fxtg �QU Qo ��QU OW���t wQUB p�kDv� Qorta B u QO xm CU� rd � ��� B�d xm s�O�O OwW�t xO�O u�Wv I�d�

Qo � �CU� d x�DQt R� �k�irD fxtg���B�d � wt

p�Ft u�wvax� CW�O R�v q�� x�DQt �� �Q di��� ��D u�wD�t CU� I��� X�N Cr�L �iO�YD uORs�o OW��u�wva Q� R ��D u�� QO �Qo�O u�twoQ CU� QmP p��k di��x� d��� �� w di��di��x�� p �iD Q�� wO ��Q��� Integrated model

Page 48: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

Time

Ele

c.ts

1960 1965 1970 1975 1980 1985 1990

2000

8000

1400

0

Time

diff(

Ele

c.ts

)

1960 1965 1970 1975 1980 1985 1990

−15

000

1000

Time

diff(

log(

Ele

c.ts

))

1960 1965 1970 1975 1980 1985 1990

−0.

150.

000.

100.

20

�tD Q�or �r��iD �QU w �r��iD �QU ��v�tR �QU V �tv ���� pmW

u�� CU� OL�w Q��Q� ZQiV�B CQwYx� u Q�Okt OvmXNWt �Q p �iD Q�N��D Ov�wD�t xm CU� OwHwt lagxir�wt di��x� lag���� p�Ft ��Q� OOQo �atH pOt QO �rYi C�Q��eD xir�wt hPL �U Ov�wD�t u�twoQCw�iDCU� umtt �Hv�� QO O��tv�t hPL xv�y�t ��y�QU QO �Q �atH pOt �rYi C�Q��eD QF� w ��N OvwQu O�m TBU �O�vm xHwD wt � vt � vt�� p�Ft x� s�y�� iQ ��Q� OW�� xDW�O s�y�� lag w d u�twoQ wO u��

CU� Q� R CQwYx� di��� ��D QO

� w �� di��v lag� d��

�yp�Ft w h� QaD ���

Ov�Qi x� p�O�D fxtg �QU p �iD u�t�d Qo � �CU� ARIMA�p�d�q� Ov�Qi ��Q�O fxtg �v�tR �QU��Hx� Qo � uwvm � CU� p�B�yt � q�B�wt x�ov �OW�� yt � ��� B�dxt Qo � OOQo ARMA�p�q�

xm OwW�t xH�Dv �OQ�o Q�Qk xt Q�eDt yt

p�B���� B�dxt � q�B�wt

ARIMA ��ypOt R� p�Ft OvJ x� uwvm � OvDUy q w p x�DQt R��DQD x� ���y��xrtH OvJ q w p u QO xmO�vm xHwD

Page 49: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� �DU ��v ��ypOt � pYi

wQUB p�kDv� Qorta CQwYx� uwvm � CU� pOt QDt�Q�B � u QO xm xt � xt�� � wt � �wt�� Ov�Qi OwW�t xDWwv

xt � xt�� � wt � �wt�� � ��� B�xt � �� � �B�wt

lQLDt u�ov��t �k�irD pOt u x� w �CU� ARIMA� ����� CQwYx� fxtg �QU �xU��kt �� u��Q��v�ARIMA� �d�q��IMA�d�q� �rm Cr�L QO OwW�t xO�O u�Wv IMA����� CQwYx� xm Ov� wo

CU�

Qorta CQwYx� uwvm � CU� pOt QDt�Q�B u QO xm xt � xt�� � xt�� � xt�� � wt Ov�Qi OwW�t xDWwv wQUB p�kDv�

xt � xt�� � xt�� � xt�� � wt � ��� B���� B�xt � wt

�rm Cr�L QO OwW�t xO�O u�Wv ARI����� CQwYx� xm Ov� wo �k�irD w�UQoQwD� pOt u x� wCU� ARI�p�d��ARIMA�p�d� �

VR�Q� w R�Ux��W ���

u��aD c�p�d�q� �� pOt x�DQt xm O����t VR�Q� R QO arima�� ��D �� �yxO�O Q� ARIMA�p�d�q� Ov�Qi�R�Ux��W xt � ��xt��� xt�� �wt � ��wt�� x� �w� Qt ��yxO�O Q� R ��yO�m QO p�Ft ��Q� OOQo�t

OwW�t xO�O VR�Q� �yv Q� ARIMA������� pOt TBU w xOW

� set�seed���

� x �� w �� rnorm������

� for �i in ������� x�i� �� ��� � x�i � �� � x�i � �� � ��� �

� x�i � �� � w�i� � ��� � w�i � ��

� arima�x order c�� � ���

Call�

arima�x x order c�� � ���

Coe�cients�

ar� ma�

������ ������

s�e� ������ ������

sigma�� estimated as ������ log likelihood �������� aic �������

OQ�O s�v arima�sim�� ��D u�� OyO�t s�Hv� �Q q�� O�m u�ty Q�m xm OQ�O OwHw ��xv�N��Dm ��D R QO xD�r�CU� Q� R CQwYx� xOW �R�Ux��W sDU�U uwvm �

� x �� arima�sim�model list�order c�� � �� ar ���

Page 50: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt ��

� ma ���� n �����

� arima�x order c�� � ���

Call�

arima�x x order c�� � ���

Coe�cients�

ar� ma�

������ ������

s�e� ������ ������

sigma�� estimated as ������ log likelihood �������� aic �������

�rYi ARIMA ��ypOt ��

h� QaD ���

�atH pOt �rYi C�Q��eD QF� �D Ovm�t xO�iDU� �s� pYi xR�Ov� x� �Q�N��D �� p �iD R� ARIMA �rYi pOtu�ov��t xrtH Qov��� p �iD s Q�N��D w Ovm�t hPL �Q OvwQ xir�wt p �iD Qorta OL�w Q�N��D �� OOQo hPLpOt OW���t s Q�N��D �� lQLDt u�ov��t w w�UQoQwD� CqtH pt�W ARIMA �rYi pOt CU� lQLDt

CWwv Q� R CQwYx� wQUB p�kDv� Qorta UL Q� u�wD�t �Q ARIMA�p� d� q��P�D�Q�s �rYi

�P �Bs�p�B���� Bs�D��� B�dxt � �Q�B

s� q�B�wt

�twta Cr�L QO OvW���t q w Q �p �P x�DQt R� ��yxrtH OvJ �DQD x� q w �Q �p ��P u QO xm��xrtH OvJ CqtH� xYNWt xrO�at ��yxW� Q jr�t QOk w D � d � Qo � �t� CU� �DU���v pOt u���rYi pOt p�Ft OvJ x� uwvm � Ow� Oy�wN �DU�� x�ov �OvW�� OL�w R� QDoQR� �oty �Q�N� xrO�at AJ hQ�

O�vm xHwD ARIMA

xt � x�� � wt CQwYx� xm ARIMA� � � ���� � ��� ��� �rYi �w�vD �� AR xO�U pOt p�U x�t �wQ xDWPo p�U x�t xm u Q� �wQWt CU� U�vt xv�y�t ��yxO�O ��Q� �rOt u�vJ CU�

OW�� j�����j � � Qo � CU� �DU�� pOt u�� OW�� xDW�Po QF� �Q�H

w ���H��H �� O�O �U� xt � xt�� � xt��� � xt��� � wt CQwYx� u�wD�t �Q q�� Ov� pOt �O��t CUOx� ���B

������B�xt � wt �� w ��� B������B�xt � wt x���Q �Q�oQwDm �iOwW�t pY�L ARIMA� � �� ���� � ��� Ov�Qi ��rYi ARIMA �rm x���Q �� xU��kt �� xmr�t u�� x� xHwD �� CWwv u�wD�t R�v rxt � rxt��� � wt CQwYx� �Q pOt u�� xm O�vm xHwD�CU� �DU���v pOt u�� OQ�O xDWPo p�U QO �x�t �va�� u�tR u�ty QO Q��eD x� �oDU� t u�tR QO Q��eD xm

OW���t B � � xW� Q ��Q�O xm CU� ��� B� xrtH pt�W AJ hQ� ��xrtH OvJ �Q� R

Page 51: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� �DU ��v ��ypOt � pYi

u�� CU� xt � �� � �B�wt � wt � �wt� CQwYx� pYi Q�yJ lQLDt u�ov��t pOt l� �OvW�� �iO�YD OvwQ ��Q�O �yxO�O Qo � CU� U�vt OvwQ uwO� ��yxO�O ��Q� �ki w CU� �DU�� pOtCQwYx� u�wD�t w �CU� xt � xt���wt��wt� xm O����t�U� pw� x�DQt p �iD CQwYx� pOt�w�vD QO p �iD �OW�� �iO�YD OvwQ ��Q�O �rYi CqtH Qo � CWwv R�vARIMA� � �� �� � � ��

CQwYx� �Q u w �OyO�t CUOx� �Q xt � xt� � wt � �wt� x���Q w OOQo�t p�ta� R�v �rYiCWwv u�wD�t R�v ARIMA� � � �� � �� ��

CW�Ol� Q�N��D �� p �iD u�wD�t u��Q��v� �Ovm�t hPL �Q ��N OvwQ �s Q�N��D �� p �iD xm CW�O xHwD O���lQLDt u�ov��t CqtH x�ov �OQ�O OwHw ��N OvwQ xm �r�L QO OOQo p�ta� p �iD QOl� Q�N��D Qo � $Q�N ��

OOQo�t �iQat O�iU xiwv CQwYx�pOtl� CQwYx� O�iU xiwv w �rYi C�Q��eD w ��N OvwQ u QO xm �w�vD �� �v�tR �QU x� �p�Ft u�wvax�

O�vm xHwD OQ�O OwHw �atHxt � a � bt� st� � wt

O�vm xHwD Q�N��D �� p �iD u�rw� x� uwvm � s� � st�� u��Q��v� OyO�t u�Wv �Q Q� t xOv�t�k�� �t� xm

���B�xt � xt � xt�

� a� bt� �at� b�t� �� � st� � st�� � wt � wt�

� b� wt � wt�

u�rw� xm O�vm ZQi uwvm � Owtv u��� b C��F xrtH �� ARIMA� � � �� ����� CQwYx� �Q u u�wD�t xmx�ov �OOQo p�ta� Q�N��D �� p �iD R� p�k l� Q�N��D �� p �iD x�DQt

��� B���� B�xt � ��� B��b � st� � st��� � wt � wt���

� wt � wt�� � wt� � wt��

CU� C��F Q�Okt Ok�i xm CWwv ARIMA� � �� ��� � �� �� CQwYx� �Q u u�wD �t xm

VR�Q� x�wQ ���

u��Q��v� OW�� CqtH R� ���mQD w OOaDt ��yQDt�Q�B pt�W Ov�wD�t xwkr�� CQwYx� ��rYi ARIMA ��ypOtR� QDU�vt pOt ��NDv� ��Q� w OwW u�LDt� �yxO�O Q� VR�Q� ��Q� pOt R� �a�Uw xOwOLt xm CU� xDU��Wxiwv CQwYx� xOv�t�k�� Q�ov�oDU�ty ��Q�O ��k�vt �xOW xDi�� VR�Q�U�vt pOt OOQo xO�iDU� AIC x���

CU� O�iUxm OQ�O OwHw seasonal u�wva Q� R �v�twoQ ��D u�� QO OwW�t xO�iDU� VR�Q� ��Q� arima�� ��D R� �R QOxOW xDiQo Q�v QO xDU� QDmr� O�rwD ��yxO�O sD� Q�or �QU Q� R p�Ft QO OOQo�t GQO u QO �rYi ��yxir�wt

Page 52: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt ��

�Q �QDy� VR�Q� ARI pOt Ovm�t hPL �Q ��N OvwQ xm CU� d � � OQwt wO Qy QO p �iD QDt�Q�B CU�OQ�O �QDmJwm AIC �Q� R �OyO�t u�Wv

� www �� �http���www�massey�ac�nz��pscowper�ts�cbe�dat�

� CBE �� read�table�www he T�

� Elec�ts �� ts�CBE� �� start ���� freq ���

� AIC �arima�log�Elec�ts� order c�����

� seas list�order c����� �����

��� ���������

� AIC �arima�log�Elec�ts� order c�����

� seas list�order c����� �����

��� ���������

��D Qt� u�� QO CrwyU ��Q� OwW�t xO�iDU� ��N w �aU R� QDU�vt pOt uDi�� w �ypOt xvt�O pQDvm ��Q�VwQ �Dkw OQm� wQ u�� OyO CUOx� �Q QDU�vt pOt AIC x��� T�U� Q� �D CWwv u�wD�t �Q �mJwm

O�vm xHwD Q� R ��yO�m x� uwvm � CU� �QDQ�wDU� sD� Qwor� w OyO�t ��wH QDy� OW�� CSS

� www �� �http���www�massey�ac�nz��pscowper�ts�cbe�dat�

� CBE �� read�table�www he T�

� Elec�ts �� ts�CBE� �� start ���� freq ���

� get�best�arima �� function�x�ts maxord c���������

� f

� best�aic �� �e�

� n �� length�x�ts�

� for �p in ��maxord���� for�d in ��maxord���� for�q in ��maxord����

� for �P in ��maxord���� for�D in ��maxord���� for�Q in ��maxord����

� f

� �t �� arima�x�ts order c�pdq�

� seas list�order c�PDQ�

� frequency�x�ts�� method �CSS��

� �t�aic �� �� � �t�loglik � �log�n� � �� � length��t�coef�

� if ��t�aic � best�aic�

� f

� best�aic �� �t�aic

� best��t �� �t

� best�model �� c�pdqPDQ�

� g

� g

� list�best�aic best��t best�model�

� g

� best�arima�elec �� get�best�arima� log�Elec�ts�

� maxord c���������

� best��t�elec �� best�arima�elec�����

� acf� resid�best��t�elec� �

� best�arima�elec �����

��� � � � � � �

Page 53: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� �DU ��v ��ypOt � pYi

� ts�plot� cbind� window�Elec�tsstart �����

� exp�predict�best��t�elec����pred� � lty ����

QDq�� D�Qt uwtR CU� ARIMA� � �� ����� � ���� swO x�DQt R� QDC�U�vt pOt �q�� ��yO�m ��QH� ��pmW OvDUy ���� pmW� O�iU xiwv ���� QkD xOW xDio pOt ��yxOv�t�k�� �Q� R �OUQ�tv Q�v x� �QwQ pOt

OyO�t u�Wv �Q xOW �v��V�B Q�O�kt ���

Time

1982 1984 1986 1988 1990 1992

8000

1000

012

000

1400

0

jQ� O�rwD xv�y�t �QU �v��V�B V �tv ��� pmW

ARCH ��ypOt ��

SP��� �yQU ���

OQ�O OwHw R QO MASS xv�N��Dm QO McGraw�Hill �v�Btm x� �w� Qt SP� ��yxO�O

� library�MASS�

� data�SP����

� par�mfrowc���� marc������� � plot�SP��� type l �

� acf�SP��� main���

�� OyO�t u�Wv �Q �DU�� Ov�Qi pw� x�ov QO OyO�t u�Wv �Q xOW xDio C�aq�� �v�tR �QU ��� pmWQ�Okt u� QDoQR� �QN ErF QO u Q�Okt w u� QDmJwm �QU �Uw ErF QO Tv�� Q�w xm OwW�t swrat QDW�� CkO

Page 54: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt ��

0.0 0.5 1.0 1.5 2.0

0.0

0.2

0.4

0.6

0.8

1.0

Lag

AC

F

Series resid(best.fit.elec)

jQ� O�rwD xv�y�t �QU ��yxOv�t�k�� V �tv ��� pmW

w OW��v C��F u�tR x� C�Uv Tv�� Q�w Qo � OwW�tv s�Hv� �t�vt Q�Ut QO V��Ri� Qo�O CQ��a x� OQ�O �Qx� x�ov �O����t V��Ri� OvwQ �� Tv�� Q�w xm � jQ� O�rwD w ���t�B�wy ��yxO�O Ovv�t �OQ�o s�Hv� s�vt C�Q��eDxO�O pFt� OOQo u�tR x� xDU��w w OW�� �W��Ri� Ow� QB ��Q�O �QU Tv�� Q�w Qo � Ov� wo heteroskedastic uQ�eDt �QU Q�ov�oDU�ty xm O�W�� xDW�O xHwD Ov� wo ��QW heteroskedastic �QU x� x�ov �SP� Cw�iDt u�tR QO xm u�� p�rO x�CU� �DU���v �QU �t� �OQ�Ov O�iU xiwv �� �Q�O�vat Cw�iD �SP� � �yCt�k� Qt Q�ov�oDU�ty x�ov �OyO�t u�Wv �Q O�iU xiwv Q�ov�oDU�ty �� pmW j���t Qo � CU� Q�eDtTv�� Q�wu�ov��t xm u �QW x�� OvDUy Tv�� Q�w pO�at C�a� Qt u�� OyO�t u�Wv �Q u�tR x� C�Uv C�Q��eD Q�O�kt � � Q��Q� � Q�t�UO �� �D � x� wv�S � R� SP� xO�O u�ov��t �OOQo s��vD QiY CQwYx� �yvOwW�t xDiQo Q�v QO �� �D �� R� T�Ov� u��Q��v� CU� �mJwm Q�Okt Tv�� Q�w �� xU��kt QO xm �CU�

CU� Q� R CQwYx� SP� xOW K�LYD u�ov��t �� � Qt Q�O�kt Q�ov�oDU�ty

library�MASS�

data�SP����

acf��SP��� � mean�SP��������

OyO�t u�Wv �Q �Q�O��B�v w ��QW heteroskedastic Kw w x� ��� pmW

Page 55: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� �DU ��v ��ypOt � pYi

0 500 1000 1500 2000 2500

−6

−2

24

Index

SP

500

0 5 10 15 20 25 30 35

0.0

0.4

0.8

Lag

AC

F

SP��� ��yxO�O �QU V �tv ���� pmW

0 5 10 15 20 25 30 35

0.0

0.2

0.4

0.6

0.8

1.0

Lag

AC

F

Series (SP500 − mean(SP500))^2

xOW K�LYD u�ov��t �� Q O�kt �� Qt V �tv ���� pmW

Page 56: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt ��

ARCH pOt h� QaD �Q�O��B�v uOQm pOt ���

x� uO�UQ ��Q� Ovm u��D �Q Tv�� Q�w QO ��QW C�Q��eD xm O�O CUOx� �rOt O��� �Q�DiQ u�� uOQm pOt ��Q�xO�O u�Wv ARCH��� �� xm CU� heteroskedastic ��QW w�UQoQwD� x�DQt u�rw� f�g �QU �r�t u��

Qo � �OwW�t�t � wt

q� � ���t��

OvDUy pOt ��yQDt�Q�B � �� w CU� OL�wTv�� Q�w w QiY u�ov��t �� O�iU xiwv fwtg q�� x���Q QO OW��O� Q�o� Tv�� Q�w Q�N� x���Q R� volality u��� ��Q�

Var��t� � E���t �

� E���t �E�� � ���t���

� E�� � ���t���

� � � �Var��t���

AR��� Ov�Qi �� �Q Q�N� xrO�at Qo � �CU� QiY u�ov��t ��Q�O f�g w OL�wTv�� Q�w ��Q�O fwtg xm u�� x� Q�vO�vm xU��kt

xt � � � �xt�� � wt

Ovm�t Q� DiQ AR��� pFt CUQO ARCH��� Ov �Qi T v� � Q�w xm O�vm�t x�Lqt OwW�t pY�L$CU� U�vt ARCH pOt VR�Q� xm Ovm�t XNWt �y��N � Qt ��y�oDU�tyOwN

GARCH �ypOt x� �w� Qt �y�U� ���

CQwYx� ARCH�p� Ov�Qi O����U� q�� ��yQ�N��D �� p x�DQt Ov�Qil� x� Ov�wD�t ARCH pOt x�DQt u�rw�CU� Q� R

�t � wt

vuut� �

pXi��

p��t�i

CU� OL�w Tv�� Q�w w QiY u�ov��t �� O�iU xiwv fwtg u QO xmARCH�p� xm �CU� GARCH�q�p� CQwYx� ARCH pOt R� �r�t ��yxO�O OQ� Q�m �� QD�Uw ��y�U�

Qo � �OW���t GARCH�q�p� pOt ��Q�O f�tg �QU CU� GARCH� �p� X�N Cr�L

�t � wt

pht

u QO xm

ht � � �

pXi��

i��t�i �

qXj��

�jht�j

OvDUy pOt ��yQDt�Q�B �j �j � � �� � � � � q� w i �i � � �� � � � � p� w CU�

Page 57: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� �DU ��v ��ypOt � pYi

xDi�� VR�Q� GARCH pOt w R�Ux��W ���

ht � � � t�� � �t�� xm at � wt

pht CQwYx� GARCH����� pOt �� ���yxO�O Q� R ��yO�m QO

CU� ��� pmW CQwYx� Q�ov�oDU�ty w OW�� � � �� � � O��� �Q�O��B ��Q� OwW�t �R�Ux��W CU�

� set�seed���

� alpha� �� ���

� alpha� �� ���

� beta� �� ���

� w �� rnorm�������

� a �� rep�� ������

� h �� rep�� ������

� for �i in �������� f

� h�i� �� alpha� � alpha� � �a�i � ����� � beta� � h�i ���

� a�i� �� w�i� � sqrt�h�i��

� g

� par�mfrowc���� marc�������

� acf�a�

� acf�a�� main���

0 10 20 30 40

0.0

0.4

0.8

Lag

AC

F

Series a

0 10 20 30 40

0.0

0.4

0.8

Lag

AC

F

�yv� �� Qt w Q O�kt Q�ov�oDU�ty V �tv ���� pmW

Qy�� �oDU�ty Q�O�kt � Qt QO xm OQ�O xDU�ty�v Q�O�kt xm OyO�t u�Wv �Q GARCH pOt X�wN a �QUOOQo�t

Page 58: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt �

garch�� ��D R� xO�iDU� �� xm O����t VR�Q� xOW �R�Ux��W ��yxO�O �wQ Q� GARCH pOt �Q� R p�Ft QO���� R�� xQ�� wO pOt ��yQDt�Q�B xm OOQo�t x�Lqt OwW�t Ci�� tseries �Q�Ri�sQv xDU� QO xm OwW�t s�Hv���awv xm CU� GARCH����� Q��Q� Q�N� ��D ZQiV�B OQ�o�t Q�Qk � u�v�t�� xR�� xOwOLt QO w OwW�t

OwW�t xO�iDU� order�c�p�q� u�twoQ R� QDq�� ��yx�DQ ��Q� �t� �CU� woNU�B

� set�seed���

� alpha� �� ���

� alpha� �� ���

� beta� �� ���

� w �� rnorm�������

� a �� rep�� ������

� h �� rep�� ������

� for �i in �������� f

� h�i� �� alpha� � alpha� � �a�i � ����� � beta� � h�i ���

� a�i� �� w�i� � sqrt�h�i��

� g

� library�tseries�

� library�quadprog�

� library�zoo�

� a�garch �� garch�a grad �numerical� trace FALSE�

� con�nt�a�garch�

OOQo�t pY�L Q� R G��Dv �q�� ��yO�m ��QH� ��

��� � ���� �

a� ��������� ���������

a� ��������� ���������

b� ��������� ���������

CQwYx� �Q �OOa C��U�Lt pL�Qt w OyO�t CUOx� �Q C��U�Lt xYqN trace�F u�twoQ p�Ft QOOyO�t CUOx� �Q �QDQ�wDU� �OOa VwQ grad��numerical� u�twoQ OyO�tv u�Wv �r��iD

SP��� QU Q� VR�Q� ���

Q� R x���Q R� fwtg �va� GARCH pOt xOv�t�k�� �QU O��� �t VR�Q� SP� �QU Q� GARCH pOtOwW�t x�U�Lt

�wt ��tp�ht

QiY u�ov��t �� O�iU xiwv CQwYx� O��� xH�Dv x�ov �OW�� U�vt �yxOv�t�k�� �QU ��Q� GARCH pOt Qo �GARCH����� Cr�L QO OwW�t Qy�� OL�w Tv�� Q�w w

�ht � �� � ����t�� � ���ht��

Page 59: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

� �DU ��v ��ypOt � pYi

CU� t � �� � � � � n ��R� x��h� � u QO xm�� Q�Okt u u��Q��v� CU�v OwHwt xOv�t�k�� �QU Q�Okt u�rw� OQ�o�t s�Hv� garch ��D �UwD C��U�LtO��t CUOx� �yxOv�t�k�� C�a� Qt �QU w xOv�t�k�� �QU Q�ov�oDU�ty w OwW�t hPL ���� R� xO�iDU�

s�rk� QU QO Q�O��B�v ���

C�k�kLD OL�w x� �w� Qt ���� �� �� ��yxO�O CU� xOW s�Hv� �D�ar��t s�rk� QO �Q�O��B�v x� H�Q ��Q�N�xOW s�UQD �� pmW QO w xOW xOv�wN Q�v OQwt �QU CU� xDiQo Q�Qk pta lqt East Anglia x�oWv�OCUOx� �CU� xOt �����CtUk ��yO�m QO xm get�best�arima ��D �UwDU�vt ARIMA pOt CU�CU� xOW p�O�D �vwDU x� �Q�U Cr�L R� t�� ��D R� xO�iDU� �� pOt ��yQDt�Q�B ��Q� u�v�t�� xR�� O��t

� stemp �� scan��http���www�massey�ac�nz��pscowper�ts�stemp�dat��

� stemp�ts �� ts�stemp start ���� freq ���

� plot�stemp�ts�

� stemp�best �� get�best�arima�stemp�ts maxord rep�����

� stemp�best�����

��� � � � � � �

Time

stem

p.ts

1850 1900 1950 2000

−1.

0−

0.5

0.0

0.5

s�rk� �Q�O �B�v �QU V �tv ���� pmW

OOQo�t OQwQ� pOt ��yQDt�Q�B jwi ��yO�m xt�O� QO

Page 60: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt ��

� stemp�arima �� arima�stemp�ts order c�����

� seas list�order c����� ����

� t� con�nt�stemp�arima� �

ar� ma� ma� sar� sar� sma�

��� � ��������� ��������� ��������� ��������� ����������� ����������

���� � ��������� ��������� ��������� ��������� ���������� ����������

VR�Q� ��OOHt Q� R ��yO�m �� pOt �Pr �OQ�Ov QiY �� �Q�O�vat Cw�iD sar� �va� AR �rYi xir�wt u�twOO����t

� stemp�arima �� arima�stemp�ts order c�����

� seas list�order c����� ����

� t� con�nt�stemp�arima� �

ar� ma� ma� sar� sma�

��� � ��������� ��������� ��������� ��������� ����������

���� � ��������� ��������� ��������� ��������� ����������

��vt CU� xOW s�UQD "�� pmW QO ARIMA pOt ��yxOv�t�k�� Q�ov�oDU�ty �VR�Q� �� wmv pQDvm ��Q�OwW�t x�Lqt Q�N� pmW QO R�v �yxOv�t�k�� � Qt Q�ov�oDU�ty �UQQ� ��Q�

� stemp�res �� resid�stemp�arima�

� par�mfrowc���� marc�������

� acf�stemp�res�

� acf�stemp�res�� main���

O����t VR�Q� �yxOv�t�k�� �QU Q� GARCH pOt u��Q��v� �OQ�O OwHw �Q�O��B�v xm CU� uWwQ �qt�m

R�Ux��W w �v��V�B QO GARCH ���

OyO�t u�Wv �yxOv�t�k�� QO �v�tR �QU p�rLD QO �Q OwN QF� GARCH pOt OW x�Lqt �q�k xm Qw�u�tyw xt�� ��r�t �v�tR ��y�QU C�ar��t QO pOt u�� �rY� OQ� Q�m u��Q��v� OQ�Ov ��kv �v��V�B QO �Q�F��D �Pr

CU� �R�Ux��W Ov�Qi ��Q� s�rk�

Page 61: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

�� �DU ��v ��ypOt � pYi

0.0 0.5 1.0 1.5 2.0 2.5

0.0

0.4

0.8

Lag

AC

F

Series stemp.res

0.0 0.5 1.0 1.5 2.0 2.5

0.0

0.4

0.8

Lag

AC

F

u� C�a� Qt w �yxOv�t�k�� Q�ov�oDU�ty V �tv ���� pmW

Page 62: C N H U Oy W Q Cav Y - R · PDF fileQ Di o V B O W x Q OvDU yx OD QOw t xv tR u QO m v U W Q wv R is wt aQ D w C U R UxO B p ku o QCQ w Y x m i Q oQ k R TQO R

��� ��vWwOv�wUwt ��

�H�Qt

��� Cowpertwait� Paul SP� Metcalf Andrew V ����� Introductory Time

Series with R� Springer� ��p�

��� Ihaka Ross� ����� Time Series Analysis� University of Auckland� � �p�

��� Shumway Robert H� Sto�er David S� ����� Time Series Analysis and

Its Applications �With R Examples�� Third edition� Spriger� ��p�