• QR Code (size=2)
• QR Code (size=3)
• QR Code (size=4)
• QR Code (size=5)
• QR Code (size=6)

0.36	105.43	2.2	11	-0.3	1.16
-0.37	85.27	-9.23	12	-0.3	-11.67
23.72	314.56	20.34	11	1.8	13.27
5.05	87.46	-6.56	20	0.9	-4.33
1.27	195.26	11.9	12	0.0	1.73
4.91	212.21	13.26	16	0.3	3.52
12.79	217.77	14.91	15	0.7	6.77
0.74	130.73	-3.96	14	0.4	-9.73
-11.82	71.85	-12.95	9	-0.9	3.37
7.16	269.85	14.99	12	1.0	8.69
-10.66	107.35	-8.35	18	-0.7	-8.32
13.08	209.47	14.97	17	1.3	2.71
6.03	167.70	10.58	16	1.3	5.72
-1.73	152.65	6.04	9	0.6	1.89
8.09	117.22	-6.66	14	-0.5	-4.52
7.51	155.51	4.38	14	0.8	-0.99
-5.4	90.35	-9.34	11	0.1	-12.3
9.54	192.45	2.55	17	-0.2	5.62
-18.85	143.11	1.98	9	0.5	-2.73
0.59	109.68	-5.84	12	0.4	-6.8
1.19	166.67	2.42	12	0.3	-3.1
18.88	253.33	11.63	11	1.4	5.06
-3.51	163.66	-3.63	16	-1.2	4.26
-5.09	105.27	5.47	13	1.4	-3.74
-3.45	130.23	-15.72	14	-0.7	-4.32
2.94	163.24	6.64	12	0.1	2.75
-1.27	137.29	2.35	11	1.4	-2.26
-18.16	95.87	-19.82	13	-2.2	-15.63
9.65	268.54	8.12	10	0.6	8.58
-11.49	128.07	-9.68	12	-0.8	-10.74
-8.79	115.78	-6.95	16	-1.1	-12.35
13.59	280.66	14.48	11	0.9	11.58
-12.12	126.88	-5.52	12	-1.1	-2.2
8.98	271.65	18.84	12	2.0	9.87
5.27	76.02	1.89	12	1.3	0.92
-7.61	71.94	-20.41	17	-1.2	-14.4
1.43	170.41	15.92	13	0.4	8.79
-14.63	115.71	-9.08	13	-0.9	-3.27
12.32	153.88	0.85	13	0.8	3.26
-13.24	86.94	-19.67	13	-1.7	-8.48
1.18	168.01	-2.53	14	-0.6	6.28
2.26	179.02	-5.32	12	-0.1	1.93
10.8	176.14	1.37	13	0.2	6.54
-2.14	162.12	-8.64	15	-1.6	-0.21
-15.22	158.13	-6.86	8	-1.7	1.23
1.83	162.17	12.73	12	0.7	-0.57
-11.29	100.31	-13.04	16	-1.0	-12.83
-4.09	199.50	1.55	14	-0.1	3.71
15.04	277.45	10.75	15	0.7	9.75
1.83	125.14	5.46	12	1.0	-3.92
-12.46	121.18	-1.57	11	0.0	-9.58
-16.72	90.79	-9.82	14	-0.6	-11.78
6.69	192.63	10.52	9	0.7	-1.01
1.84	155.01	12.55	11	2.9	-9.21
-2.08	191.35	1.58	10	-0.2	2.45
-10.76	87.26	-12.61	19	0.0	-9.03
6.55	129.68	7.22	11	0.8	-1.51
5.69	246.73	8.04	11	1.9	8.94
-12.02	113.94	-13.99	12	-0.5	-4.05
16.34	248.33	1.82	17	0.5	13.72
15.36	210.28	16.07	13	1.4	9.25
-2.71	115.63	-3.52	13	0.0	-6.88
2.46	163.59	6.74	13	0.0	8.32
14.42	237.51	16.8	16	1.3	4.82
13.59	218.64	15.04	9	1.1	3.19
0.2	202.67	12.02	13	0.6	4.73
11.14	95.77	2.19	11	0.6	4.09
6.56	210.17	10.26	12	0.8	10.98
-6.26	125.72	-0.38	12	-0.5	-3.65
-10.42	108.12	-16.18	11	-1.4	-9.34
-6.61	114.78	-7.8	14	-0.3	-6.51
-15.79	97.08	-18.48	17	-1.1	-13.1
12.73	206.74	0.14	17	-0.9	4.82
5.06	187.36	4.17	14	0.8	1.5
-20.99	109.12	-11.91	12	-0.6	-8.17
1.55	113.85	-1.06	14	0.3	-8.49
7.3	176.89	0.8	13	-0.8	2.83
14.33	263.99	8.56	19	0.5	14.67
-0.46	109.70	-9.79	15	-0.3	-9.35
20.73	289.70	20.44	15	2.0	14.87
17.84	247.85	10.21	13	1.3	12.21
8.5	204.71	1.1	16	-0.4	2.3
-15.4	95.63	-6.1	13	0.3	-6.79
11.05	238.13	22.22	12	1.8	5.38
-7.66	177.93	-1.42	11	-0.2	1.97
6.46	214.23	2.46	15	0.3	8.99
6.32	189.69	2.74	15	0.6	7.37
-5.92	166.41	-6.33	15	-0.7	5.62
-11.63	119.08	-0.28	10	-0.1	0.13
-3.57	185.25	3.19	10	-0.4	-3.06
6.26	148.86	-2.1	14	1.0	-0.06
-4.98	119.51	-9.49	15	-0.3	-5.43
-11.31	103.15	-7.63	5	-0.4	-8.49
-2.18	236.28	6.7	12	-0.1	2.05
5.88	176.61	-10.9	12	-0.8	0.51
11.53	274.79	14.39	14	0.8	14.31
-12.22	141.73	-18.17	15	-2.5	2.05
2.13	135.65	-10.91	13	-1.3	0.57
18.81	229.36	12.42	13	1.7	6.91
-1.51	164.01	-4.84	11	-0.1	4.34
6.45	143.27	-2.23	19	0.0	-2.91
13.94	253.27	8.13	18	-0.4	10.31
1.09	260.72	3.98	13	0.2	11.87
9.4	265.43	9.96	15	0.4	7.84
11.39	196.17	18.96	14	2.0	2.58
-14.24	113.10	-6.43	12	0.4	-11.62
3.04	197.78	-3.56	17	0.3	4.03
-15.45	80.84	-3.18	12	-1.0	-8.89
-14.27	92.73	-8.62	5	-0.7	-7.46
10.16	132.36	-1.79	13	0.8	-4.69
-7.26	102.78	-14.67	13	-1.1	-4.95
-13.59	126.54	-11.38	15	-0.8	0.29
-4.74	124.34	-15.23	17	-1.7	-1.31
8.73	256.97	14.4	18	0.9	9.53
11.39	181.67	9.84	14	0.0	1.98
-1.23	115.24	3.63	15	1.1	-4.98
-2.27	191.52	-6.97	18	-1.0	3.88
-0.95	201.67	4.83	10	0.0	0.87
5.04	219.78	6.71	16	-0.2	3.2
-3.7	159.37	-7.64	16	-1.8	-2.82
9.73	217.20	6.52	10	0.0	8.03
6.55	175.97	-1.54	14	-0.1	2.45
7.59	199.34	8.07	14	0.5	3.94
9.62	115.28	1.99	16	0.3	-4.28
4.77	128.92	-1.79	16	1.0	-6.89
10.85	185.98	12.7	13	1.9	8.06
-17.03	98.16	-6.69	10	-0.9	-4.68
2.2	83.02	-14.88	15	-1.0	1.16
-9.23	68.33	-9.86	13	-0.6	-11.67
20.34	311.91	20.06	12	2.3	13.27
-0.3	86.84	-15.82	14	-1.4	-7.58
11.9	190.38	12.6	12	0.5	1.73
13.26	190.27	17.5	10	0.3	3.52
14.91	218.75	-1.32	15	-0.6	6.77
-3.96	103.08	-1.86	12	-0.2	-9.73
-12.95	47.56	-12.36	17	-1.0	3.37
14.99	270.74	18.8	14	0.9	8.69
-8.35	89.26	2.29	12	0.9	-8.32
14.97	203.11	16.87	9	2.7	2.71
10.58	164.82	3.88	14	0.7	5.72
6.04	156.13	1.28	10	0.1	1.89
-6.66	98.96	4.89	13	1.4	-4.52
4.38	150.20	10.23	13	0.6	-0.99
-9.34	72.56	-1.48	10	-0.3	-12.3
2.55	204.97	-8.58	14	-1.9	5.62
1.98	143.22	8.42	15	1.2	-2.73
-5.84	89.13	-11.97	19	-0.5	-6.8
2.42	167.71	0.41	13	1.4	-3.1
11.63	251.31	14.98	12	1.4	5.06
-3.63	160.06	-1.99	15	-1.1	4.26
5.47	109.84	-0.23	13	0.5	-3.74
-15.72	127.69	-11.97	14	-0.9	-4.32
6.64	156.61	7.86	14	0.2	2.75
2.35	140.09	4.3	9	1.0	-2.26
-19.82	77.36	-22.84	15	-2.4	-15.63
8.12	276.42	4.67	14	-0.5	8.58
-9.68	103.32	-3.53	10	0.5	-10.74
-6.95	99.52	-23.53	18	-3.0	-12.35
14.48	283.34	23.36	13	3.6	11.58
-5.52	96.32	5.04	13	0.9	-2.2
18.84	270.99	12.82	15	1.0	9.87
-20.41	48.64	-20.2	10	-0.9	-14.4
15.92	167.69	7.21	11	0.6	8.79
-9.08	102.86	-11.13	13	-0.7	-3.27
0.85	155.24	4.93	17	0.5	3.26
-19.67	67.11	-19.3	15	-2.1	-8.48
-2.53	168.34	-3.86	16	-0.8	6.28
-5.32	174.54	4.41	18	-0.1	1.93
1.37	186.17	7.31	13	0.7	6.54
-8.64	159.59	-3.57	16	-1.3	-0.21
-6.86	164.99	-8.58	16	-1.5	1.23
12.73	159.79	10.01	10	1.2	-0.57
-13.04	73.01	-9.6	14	-0.5	-12.83
1.55	197.10	-3.51	16	-1.0	3.71
10.75	275.37	14.59	12	1.8	9.75
5.46	93.69	-15.66	13	-0.2	-3.92
-1.57	107.49	-2.01	17	0.3	-9.58
-9.82	66.96	-6.81	13	-0.8	-11.78
10.52	194.42	9.99	13	2.5	-1.01
12.55	155.70	3.07	13	1.5	-9.21
1.58	186.49	0.64	17	0.7	2.45
-12.61	57.22	-20.4	12	-2.3	-9.03
7.22	123.34	-7.41	17	-0.6	-1.51
8.04	242.36	6.59	13	1.0	8.94
-13.99	82.18	-23.48	15	-2.7	-4.05
1.82	252.35	5.53	15	0.1	13.72
16.07	204.42	14.74	11	1.4	9.25
-3.52	93.72	-3.71	11	0.4	-6.88
6.74	161.75	5.27	14	-0.5	8.32
16.8	232.98	10.42	15	0.7	4.82
15.04	207.30	7.72	16	0.4	3.19
12.02	187.78	19.97	9	1.2	4.73
2.19	72.99	4.26	15	0.1	4.09
10.26	218.53	5.22	11	0.4	10.98
-0.38	97.73	-3.6	18	-1.8	-3.65
-16.18	85.04	-3.9	15	-0.3	-9.34
-7.8	97.62	-14.1	14	-0.9	-6.51
-18.48	74.57	-19.36	11	-2.3	-13.1
0.14	211.04	6.22	15	0.6	4.82
4.17	169.67	-5.66	15	-0.8	1.5
-11.91	83.29	3.05	9	0.5	-8.17
-1.06	89.14	3.56	15	1.4	-8.49
0.8	170.72	1.2	17	-0.1	2.83
8.56	264.70	9.5	11	1.0	14.67
-9.79	81.07	-5.77	11	-0.1	-9.35
20.44	289.77	15.65	12	1.8	14.87
10.21	249.34	13.19	14	0.4	12.21
1.1	210.07	16.83	9	1.1	2.3
22.22	241.14	21.27	15	2.9	5.38
-1.42	179.88	8.11	14	0.1	1.97
2.46	202.03	3.4	15	0.4	8.99
2.74	181.17	3.72	14	0.0	7.37
-6.33	157.63	-12.43	13	-1.6	5.62
-0.28	94.00	0.15	12	0.1	0.13
3.19	193.96	-5.68	16	-0.4	-3.06
-2.1	137.42	-4.68	11	0.4	-0.06
-9.49	79.51	-3.28	14	0.3	-5.43
-7.63	76.29	-2.44	14	0.2	-8.49
6.7	238.05	16.06	14	0.9	2.05
-10.9	171.42	-13.56	14	-1.1	0.51
14.39	269.46	13.72	14	1.1	14.31
-18.17	141.93	-7.44	14	-0.1	2.05
-10.91	117.37	-21.13	15	-1.3	0.57
12.42	229.70	18.96	11	1.4	6.91
-4.84	161.16	2.81	15	-0.2	4.34
-2.23	131.09	-10.3	13	-0.7	-2.91
8.13	245.30	4.29	10	-0.8	10.31
3.98	263.66	7.45	13	-0.2	11.87
9.96	246.79	13.38	15	0.9	7.84
18.96	201.16	1.94	15	0.2	2.58
-6.43	84.70	-10.53	11	-0.3	-11.62
-3.56	205.61	5.32	12	0.4	4.03
-3.18	63.51	-18.17	16	-2.3	-8.89
-8.62	62.59	-3.01	19	0.6	-7.46
-1.79	131.37	-0.03	16	0.6	-4.69
-14.67	80.17	-5.33	10	-0.7	-4.95
-11.38	122.07	-3.22	13	-0.3	0.29
-15.23	124.11	-11.77	16	-0.8	-1.31
14.4	232.51	17.1	16	1.0	9.53
9.84	181.38	6.96	11	0.1	1.98
3.63	79.77	6.7	7	0.8	-4.98
-6.97	186.64	8.01	10	0.3	3.88
4.83	212.04	-4.94	15	-1.3	0.87
6.71	214.32	6.41	12	0.5	3.2
-7.64	151.13	0	10	-0.6	-2.82
6.52	218.21	15.12	12	1.3	8.03
-1.54	164.89	5.09	14	-0.4	2.45
8.07	191.19	-3.74	14	-0.4	3.94
1.99	102.85	-0.52	11	1.0	-4.28
-1.79	115.25	-19.16	12	-1.1	-6.89
12.7	184.54	14.2	9	2.2	8.06
-6.69	75.90	-9.36	16	0.2	-4.68
-14.88	64.56	-7.35	12	-0.1	1.16
-9.86	73.21	-15.31	14	-0.6	-11.67
20.06	302.35	24.51	13	2.7	13.27
-15.82	95.16	-11.66	12	-0.6	-7.58
12.6	182.14	9.38	17	0.3	1.73
17.5	183.88	9	14	1.2	3.52
-1.32	212.88	2.4	11	0.4	6.77
-1.86	106.70	4.94	12	1.0	-9.73
-12.36	51.07	-14.23	17	-1.0	3.37
18.8	270.00	-0.92	15	-0.7	8.69
2.29	93.67	-1.6	14	-0.6	-8.32
16.87	197.19	11.22	13	0.9	2.71
3.88	151.05	7.85	13	1.2	5.72
1.28	157.62	-6.74	14	-0.5	1.89
4.89	106.16	-0.27	19	0.6	-4.52
10.23	161.53	8.77	14	0.8	-0.99
-1.48	72.98	-8.78	16	0.0	-12.3
-8.58	217.44	-0.65	9	0.0	5.62
8.42	145.14	6.65	11	1.0	-2.73
-11.97	95.55	-7.91	14	0.1	-6.8
0.41	170.14	8.71	14	1.3	-3.1
14.98	243.53	14.75	13	0.7	5.06
-1.99	165.56	-12.66	17	-2.7	4.26
-0.23	115.01	-12.54	16	-0.1	-3.74
-11.97	138.08	-5.35	13	0.1	-4.32
7.86	154.00	-0.28	14	-1.3	2.75
4.3	138.09	-4.23	15	-0.3	-2.26
-22.84	81.63	-16.61	13	-1.4	-15.63
4.67	290.40	17.35	11	0.9	8.58
-3.53	89.91	-9.83	12	-1.1	-10.74
-23.53	106.75	-9.1	8	-0.2	-12.35
23.36	277.99	13.56	10	3.1	11.58
5.04	102.97	6.33	13	1.5	-2.2
12.82	265.42	17.93	13	1.9	9.87
7.21	172.40	5.18	16	0.5	8.79
-11.13	108.80	-13.75	18	-1.1	-3.27
4.93	160.76	-5.54	14	0.2	3.26
-19.3	81.68	-17.38	11	-1.9	-8.48
-3.86	171.91	-10.35	12	-1.3	6.28
4.41	167.56	-5.53	19	-0.5	1.93
7.31	192.48	-2.62	14	-0.8	6.54
-3.57	157.46	5.4	9	-0.6	-0.21
-8.58	162.69	-6.7	11	-1.8	1.23
10.01	152.84	16.92	8	0.8	-0.57
-9.6	85.36	2.75	11	0.2	-12.83
-3.51	185.47	-5.29	13	-0.7	3.71
14.59	266.47	14.76	13	1.0	9.75
-15.66	95.75	6.33	6	0.1	-3.92
-2.01	103.02	3.45	13	0.9	-9.58
-6.81	66.77	-0.03	17	-0.2	-11.78
9.99	194.64	5.88	15	0.6	-1.01
3.07	151.05	-6.04	15	0.2	-9.21
0.64	187.41	-6.59	12	-0.6	2.45
-20.4	34.67	-20.97	9	-2.5	-9.03
-7.41	131.23	-11.09	13	-0.5	-1.51
6.59	238.91	4.07	19	0.4	8.94
-23.48	91.94	-11.7	13	-0.8	-4.05
5.53	251.99	10.01	12	1.2	13.72
14.74	205.17	6.21	15	0.5	9.25
-3.71	85.67	-2.45	17	-0.3	-6.88
5.27	161.06	-3.25	16	-0.4	8.32
10.42	227.24	8.43	18	0.3	4.82
7.72	208.26	6.15	13	0.4	3.19
19.97	193.73	5.54	12	-0.6	4.73
4.26	58.43	-1.76	11	-0.4	4.09
5.22	221.86	9.79	13	1.0	10.98
-3.6	98.29	-3.6	11	0.3	-3.65
-3.9	93.79	-14.57	18	-0.8	-9.34
-14.1	105.73	-9.91	11	-1.1	-6.51
-19.36	82.84	-19	14	-1.6	-13.1
6.22	212.09	4.52	15	1.3	4.82
-5.66	168.92	1.46	11	-0.2	1.5
3.05	77.61	-8.41	16	-0.6	-8.17
3.56	90.76	4.64	12	1.7	-8.49
1.2	165.57	8.97	15	0.1	2.83
9.5	257.70	16.98	14	1.2	14.67
-5.77	78.80	-2.75	11	0.4	-9.35
15.65	274.16	13.81	13	1.0	14.87
13.19	257.67	14.32	11	0.7	12.21
16.83	207.19	13.5	14	1.7	2.3
21.27	243.87	23.43	15	1.7	5.38
8.11	180.64	12.14	15	0.7	1.97
3.4	208.73	9.7	13	0.5	8.99
3.72	186.34	2.05	14	0.8	7.37
-12.43	161.66	-1.58	13	-0.2	5.62
0.15	95.74	-3.52	17	-0.9	0.13
-5.68	194.09	3.91	10	0.6	-3.06
-4.68	138.05	2.69	16	0.6	-0.06
-3.28	80.57	7.9	13	1.2	-5.43
-2.44	51.83	-15.26	17	-1.1	-8.49
16.06	231.51	15.97	12	1.1	2.05
-13.56	166.49	-6.13	10	-0.3	0.51
13.72	276.71	10.53	15	1.3	14.31
-7.44	147.73	1.64	9	-0.2	2.05
-21.13	127.29	-18.04	13	-1.2	0.57
18.96	220.73	15.58	15	0.8	6.91
2.81	156.99	6.27	13	0.3	4.34
-10.3	124.55	-7.52	15	-1.2	-2.91
4.29	242.29	4.07	13	0.3	10.31
7.45	281.36	11.18	19	0.4	11.87
13.38	234.04	20.05	10	1.9	7.84
1.94	196.76	7.49	15	0.6	2.58
-10.53	57.30	-8.03	13	-0.6	-11.62
5.32	207.36	6.33	13	1.1	4.03
-18.17	74.34	-9.29	11	-0.7	-8.89
-3.01	42.15	-7.49	19	0.2	-7.46
-0.03	133.65	-0.18	13	-0.5	-4.69
-5.33	92.60	-3.45	7	-0.1	-4.95
-3.22	122.71	-14.23	17	-1.7	0.29
-11.77	128.95	6.24	10	0.6	-1.31
17.1	238.17	8.75	13	0.4	9.53
6.96	187.85	2.51	13	-0.4	1.98
6.7	82.73	10.71	15	2.5	-4.98
8.01	183.71	9.88	13	0.9	3.88
-4.94	201.56	-3.84	14	0.2	0.87
6.41	211.02	3.13	13	0.4	3.2
0	151.10	-6.2	16	-1.4	-2.82
15.12	226.46	4.9	17	-0.3	8.03
5.09	164.06	-6.89	15	-0.8	2.45
-3.74	193.76	4.03	9	0.2	3.94
-0.52	99.43	-1.5	13	0.7	-4.28
-19.16	107.94	-9.54	10	0.2	-6.89
14.2	179.45	9.77	15	1.1	8.06
-9.36	86.16	-10.78	14	-0.7	-4.68
4.09	102.52	0.36	15	0.2	1.16
-6.52	91.64	-0.37	8	-0.1	-11.67
25.62	312.39	23.72	11	2.5	13.27
8.83	104.45	5.05	15	1.1	-4.33
-4.63	192.31	1.27	16	0.0	1.73
-2.6	215.84	4.91	10	-1.9	3.52
4.35	218.12	12.79	14	-0.7	6.77
-0.56	122.76	0.74	13	0.4	-9.73
-2.41	92.38	-11.82	16	0.8	3.37
11.97	274.09	7.16	11	1.1	8.69
-9.11	108.49	-10.66	14	-0.7	-8.32
3.01	202.46	13.08	10	0.7	2.71
8.01	165.67	6.03	13	1.6	5.72
0.58	153.92	-1.73	15	-0.7	1.89
-10.58	115.48	8.09	11	-1.2	-4.52
6.5	164.07	7.51	15	-0.1	-0.99
-18.24	96.94	-5.4	11	-0.8	-12.3
2.42	192.73	9.54	10	-0.8	5.62
-5.6	106.17	0.59	14	0.1	-6.8
-12.51	97.53	-19.53	17	-1.1	-16.02
-9.17	158.77	1.19	12	-0.1	-3.1
20.08	259.54	18.88	12	1.6	5.06
12.18	153.83	-3.51	17	0.7	4.26
2.52	127.71	-5.09	10	0.6	-3.74
-13.3	125.39	-3.45	16	-0.9	-4.32
-1.36	173.55	2.94	11	-1.1	2.75
-9.69	141.68	-1.27	11	-0.6	-2.26
-4.3	98.84	-18.16	15	0.1	-15.63
10.98	262.03	9.65	11	0.5	8.58
-14.6	128.62	-11.49	14	-1.2	-10.74
-12.82	113.30	-8.79	14	-0.9	-12.35
15.01	282.59	13.59	17	1.2	11.58
-14.26	140.46	-12.12	12	-1.0	-2.2
3.64	273.29	8.98	14	0.2	9.87
-4.4	108.21	5.27	14	-0.7	0.92
-9.37	91.26	-7.61	17	0.4	-14.4
7.13	167.08	1.43	11	0.8	8.79
-6.44	116.00	-14.63	17	-0.4	-3.27
5.44	161.66	12.32	11	0.5	3.26
-1.19	88.70	-13.24	14	-0.6	-8.48
-7.75	166.47	1.18	12	-0.2	6.28
1.38	178.58	2.26	13	0.5	1.93
6.81	171.93	10.8	12	0.1	6.54
-4.43	167.13	-2.14	14	-0.1	-0.21
-10.65	160.57	-15.22	13	-1.1	1.23
8.09	165.68	1.83	12	0.0	-0.57
-9.12	95.36	-11.29	18	-0.4	-12.83
2.3	204.46	-4.09	14	0.2	3.71
3.39	130.25	1.83	9	1.7	-3.92
-6.18	126.49	-12.46	15	-0.3	-9.58
-11.97	95.46	-16.72	12	-1.5	-11.78
11.87	193.62	6.69	18	2.5	-1.01
15.93	288.32	15.04	18	2.0	9.75
-13.59	163.16	1.84	12	-1.3	-9.21
0.07	188.32	-2.08	12	-0.1	2.45
-11.94	118.33	-10.76	9	-0.8	-9.03
2.96	141.43	6.55	12	0.5	-1.51
13.36	239.37	5.69	16	1.7	8.94
-6.79	116.21	-12.02	15	0.4	-4.05
17.56	250.15	16.34	13	1.6	13.72
-1.79	218.34	15.36	10	0.1	9.25
-5.53	127.36	-2.71	12	1.5	-6.88
8.65	167.08	2.46	13	0.4	8.32
1.68	216.47	13.59	13	-0.2	3.19
-2.5	197.31	0.2	13	0.3	4.73
3.59	110.24	11.14	8	0.5	4.09
7.55	212.18	6.56	14	-0.2	10.98
-9.7	123.51	-6.26	14	-1.1	-3.65
-10.05	117.18	-10.42	13	0.1	-9.34
-1.21	119.59	-6.61	14	0.6	-6.51
-12.88	90.61	-15.79	14	-1.1	-13.1
8.25	208.46	12.73	14	1.3	4.82
7.49	187.23	5.06	14	0.6	1.5
-10.33	112.67	-20.99	15	-1.0	-8.17
-4.39	113.99	1.55	12	0.1	-8.49
5.46	176.27	7.3	11	-0.1	2.83
4.27	265.27	14.33	10	0.7	14.67
-1.99	108.58	-0.46	11	0.3	-9.35
18.82	291.56	20.73	6	1.7	14.87
14.98	244.56	17.84	13	1.7	12.21
11.16	198.22	8.5	16	0.8	2.3
1.11	115.28	-15.4	17	1.0	-6.79
4.74	260.73	14.42	10	0.0	4.82
0.7	229.19	11.05	11	0.2	5.38
1.53	175.28	-7.66	14	-0.2	1.97
15.72	226.36	6.46	14	1.5	8.99
8.12	192.94	6.32	15	0.7	7.37
-8.96	156.51	-5.92	16	-1.1	5.62
-13.67	112.89	-11.63	16	-2.3	0.13
-6.97	169.16	-3.57	15	-2.2	-3.06
5.97	149.19	6.26	13	1.6	-0.06
-10.03	127.38	-4.98	15	-1.2	-5.43
-4.84	142.51	-12.22	16	-0.5	2.05
-5.8	120.30	-11.31	12	0.1	-8.49
-0.91	229.77	-2.18	13	-0.5	2.05
7.83	172.36	5.88	14	1.4	0.51
-3.73	130.59	2.13	13	0.3	0.57
12.54	223.00	18.81	12	0.4	6.91
-4.51	163.48	-1.51	15	-1.5	4.34
8.97	150.73	6.45	12	1.4	-2.91
10.59	256.98	13.94	17	0.6	10.31
1.75	256.84	1.09	14	0.3	11.87
8.77	264.16	9.4	13	0.8	7.84
2.63	201.22	11.39	10	0.6	2.58
-18.86	118.12	-14.24	9	-1.9	-11.62
1.34	188.99	3.04	12	-0.5	4.03
4.94	130.55	10.16	10	1.4	-4.69
4.1	106.91	-7.26	15	0.7	-4.95
-6.54	130.06	-13.59	16	-0.5	0.29
6.91	126.91	-4.74	15	0.7	-1.31
-0.01	156.67	-18.85	16	-0.1	-2.73
2.63	261.67	8.73	11	0.6	9.53
17.27	275.37	11.53	15	1.0	14.31
8.65	181.35	11.39	14	0.2	1.98
-5.99	123.87	-1.23	11	0.5	-4.98
-13.33	86.31	-15.45	14	-0.1	-8.89
-3.92	106.65	-14.27	10	-0.3	-7.46
3.21	187.69	-2.27	15	-0.7	3.88
-8.25	187.87	-0.95	12	-1.5	0.87
11	210.93	5.04	16	0.7	3.2
0.73	162.82	-3.7	16	-0.8	-2.82
19.98	214.31	9.73	15	2.2	8.03
9.03	175.95	6.55	14	0.1	2.45
9.57	197.13	7.59	13	0.9	3.94
8.68	120.52	9.62	11	2.5	-4.28
12.03	140.42	4.77	14	0.8	-6.89
15.3	197.83	10.85	13	0.5	8.06
1.37	102.99	-17.03	18	-0.5	-4.68

library(lattice)
library(lmtest)
n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
mywarning <- ''
par1 <- as.numeric(par1)
if(is.na(par1)) {
par1 <- 1
mywarning = 'Warning: you did not specify the column number of the endogenous series! The first column was selected by default.'
}
if (par4=='') par4 <- 0
par4 <- as.numeric(par4)
if (par5=='') par5 <- 0
par5 <- as.numeric(par5)
x <- na.omit(t(y))
k <- length(x[1,])
n <- length(x[,1])
x1 <- cbind(x[,par1], x[,1:k!=par1])
mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
colnames(x1) <- mycolnames #colnames(x)[par1]
x <- x1
if (par3 == 'First Differences'){
(n <- n -1)
x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-B)',colnames(x),sep='')))
for (i in 1:n) {
for (j in 1:k) {
x2[i,j] <- x[i+1,j] - x[i,j]
}
}
x <- x2
}
if (par3 == 'Seasonal Differences (s=12)'){
(n <- n - 12)
x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-B12)',colnames(x),sep='')))
for (i in 1:n) {
for (j in 1:k) {
x2[i,j] <- x[i+12,j] - x[i,j]
}
}
x <- x2
}
if (par3 == 'First and Seasonal Differences (s=12)'){
(n <- n -1)
x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-B)',colnames(x),sep='')))
for (i in 1:n) {
for (j in 1:k) {
x2[i,j] <- x[i+1,j] - x[i,j]
}
}
x <- x2
(n <- n - 12)
x2 <- array(0, dim=c(n,k), dimnames=list(1:n, paste('(1-B12)',colnames(x),sep='')))
for (i in 1:n) {
for (j in 1:k) {
x2[i,j] <- x[i+12,j] - x[i,j]
}
}
x <- x2
}
if(par4 > 0) {
x2 <- array(0, dim=c(n-par4,par4), dimnames=list(1:(n-par4), paste(colnames(x)[par1],'(t-',1:par4,')',sep='')))
for (i in 1:(n-par4)) {
for (j in 1:par4) {
x2[i,j] <- x[i+par4-j,par1]
}
}
x <- cbind(x[(par4+1):n,], x2)
n <- n - par4
}
if(par5 > 0) {
x2 <- array(0, dim=c(n-par5*12,par5), dimnames=list(1:(n-par5*12), paste(colnames(x)[par1],'(t-',1:par5,'s)',sep='')))
for (i in 1:(n-par5*12)) {
for (j in 1:par5) {
x2[i,j] <- x[i+par5*12-j*12,par1]
}
}
x <- cbind(x[(par5*12+1):n,], x2)
n <- n - par5*12
}
if (par2 == 'Include Monthly Dummies'){
x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
for (i in 1:11){
x2[seq(i,n,12),i] <- 1
}
x <- cbind(x, x2)
}
if (par2 == 'Include Quarterly Dummies'){
x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
for (i in 1:3){
x2[seq(i,n,4),i] <- 1
}
x <- cbind(x, x2)
}
(k <- length(x[n,]))
if (par3 == 'Linear Trend'){
x <- cbind(x, c(1:n))
colnames(x)[k+1] <- 't'
}
x
(k <- length(x[n,]))
head(x)
df <- as.data.frame(x)
(mylm <- lm(df))
(mysum <- summary(mylm))
if (n > n25) {
kp3 <- k + 3
nmkm3 <- n - k - 3
gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
numgqtests <- 0
numsignificant1 <- 0
numsignificant5 <- 0
numsignificant10 <- 0
for (mypoint in kp3:nmkm3) {
j <- 0
numgqtests <- numgqtests + 1
for (myalt in c('greater', 'two.sided', 'less')) {
j <- j + 1
gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
}
if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
}
gqarr
}
bitmap(file='test0.png')
plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
points(x[,1]-mysum$resid)
grid()
dev.off()
bitmap(file='test1.png')
plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
grid()
dev.off()
bitmap(file='test2.png')
hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
grid()
dev.off()
bitmap(file='test3.png')
densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
dev.off()
bitmap(file='test4.png')
qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
qqline(mysum$resid)
grid()
dev.off()
(myerror <- as.ts(mysum$resid))
bitmap(file='test5.png')
dum <- cbind(lag(myerror,k=1),myerror)
dum
dum1 <- dum[2:length(myerror),]
dum1
z <- as.data.frame(dum1)
z
plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
lines(lowess(z))
abline(lm(z))
grid()
dev.off()
bitmap(file='test6.png')
acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
grid()
dev.off()
bitmap(file='test7.png')
pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
grid()
dev.off()
bitmap(file='test8.png')
opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
plot(mylm, las = 1, sub='Residual Diagnostics')
par(opar)
dev.off()
if (n > n25) {
bitmap(file='test9.png')
plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
grid()
dev.off()
}
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
a<-table.row.end(a)
myeq <- colnames(x)[1]
myeq <- paste(myeq, '[t] = ', sep='')
for (i in 1:k){
if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
myeq <- paste(myeq, signif(mysum$coefficients[i,1],6), sep=' ')
if (rownames(mysum$coefficients)[i] != '(Intercept)') {
myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
}
}
myeq <- paste(myeq, ' + e[t]')
a<-table.row.start(a)
a<-table.element(a, myeq)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, mywarning)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,hyperlink('ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',header=TRUE)
a<-table.element(a,'Parameter',header=TRUE)
a<-table.element(a,'S.D.',header=TRUE)
a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
a<-table.element(a,'2-tail p-value',header=TRUE)
a<-table.element(a,'1-tail p-value',header=TRUE)
a<-table.row.end(a)
for (i in 1:k){
a<-table.row.start(a)
a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
a<-table.element(a,formatC(signif(mysum$coefficients[i,1],5),format='g',flag='+'))
a<-table.element(a,formatC(signif(mysum$coefficients[i,2],5),format='g',flag=' '))
a<-table.element(a,formatC(signif(mysum$coefficients[i,3],4),format='e',flag='+'))
a<-table.element(a,formatC(signif(mysum$coefficients[i,4],4),format='g',flag=' '))
a<-table.element(a,formatC(signif(mysum$coefficients[i,4]/2,4),format='g',flag=' '))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Multiple R',1,TRUE)
a<-table.element(a,formatC(signif(sqrt(mysum$r.squared),6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'R-squared',1,TRUE)
a<-table.element(a,formatC(signif(mysum$r.squared,6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Adjusted R-squared',1,TRUE)
a<-table.element(a,formatC(signif(mysum$adj.r.squared,6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'F-TEST (value)',1,TRUE)
a<-table.element(a,formatC(signif(mysum$fstatistic[1],6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
a<-table.element(a, signif(mysum$fstatistic[2],6))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
a<-table.element(a, signif(mysum$fstatistic[3],6))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'p-value',1,TRUE)
a<-table.element(a,formatC(signif(1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]),6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
a<-table.element(a,formatC(signif(mysum$sigma,6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
a<-table.element(a,formatC(signif(sum(myerror*myerror),6),format='g',flag=' '))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable3.tab')
if(n < 200) {
a<-table.start()
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Time or Index', 1, TRUE)
a<-table.element(a, 'Actuals', 1, TRUE)
a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
a<-table.row.end(a)
for (i in 1:n) {
a<-table.row.start(a)
a<-table.element(a,i, 1, TRUE)
a<-table.element(a,formatC(signif(x[i],6),format='g',flag=' '))
a<-table.element(a,formatC(signif(x[i]-mysum$resid[i],6),format='g',flag=' '))
a<-table.element(a,formatC(signif(mysum$resid[i],6),format='g',flag=' '))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable4.tab')
if (n > n25) {
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-values',header=TRUE)
a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'breakpoint index',header=TRUE)
a<-table.element(a,'greater',header=TRUE)
a<-table.element(a,'2-sided',header=TRUE)
a<-table.element(a,'less',header=TRUE)
a<-table.row.end(a)
for (mypoint in kp3:nmkm3) {
a<-table.row.start(a)
a<-table.element(a,mypoint,header=TRUE)
a<-table.element(a,formatC(signif(gqarr[mypoint-kp3+1,1],6),format='g',flag=' '))
a<-table.element(a,formatC(signif(gqarr[mypoint-kp3+1,2],6),format='g',flag=' '))
a<-table.element(a,formatC(signif(gqarr[mypoint-kp3+1,3],6),format='g',flag=' '))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable5.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Description',header=TRUE)
a<-table.element(a,'# significant tests',header=TRUE)
a<-table.element(a,'% significant tests',header=TRUE)
a<-table.element(a,'OK/NOK',header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'1% type I error level',header=TRUE)
a<-table.element(a,signif(numsignificant1,6))
a<-table.element(a,formatC(signif(numsignificant1/numgqtests,6),format='g',flag=' '))
if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'5% type I error level',header=TRUE)
a<-table.element(a,signif(numsignificant5,6))
a<-table.element(a,signif(numsignificant5/numgqtests,6))
if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'10% type I error level',header=TRUE)
a<-table.element(a,signif(numsignificant10,6))
a<-table.element(a,signif(numsignificant10/numgqtests,6))
if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable6.tab')
}
}