R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(0
+ ,17.080
+ ,9
+ ,0
+ ,17.240
+ ,8
+ ,0
+ ,17.200
+ ,7
+ ,0
+ ,15.580
+ ,9
+ ,0
+ ,18.950
+ ,9
+ ,0
+ ,23.360
+ ,8
+ ,0
+ ,19.190
+ ,6
+ ,0
+ ,14.150
+ ,6
+ ,0
+ ,19.870
+ ,8
+ ,0
+ ,19.420
+ ,9
+ ,0
+ ,16.020
+ ,6
+ ,0
+ ,17.350
+ ,8
+ ,0
+ ,21.930
+ ,8
+ ,0
+ ,21.180
+ ,8
+ ,0
+ ,22.580
+ ,8
+ ,0
+ ,19.320
+ ,7
+ ,0
+ ,14.720
+ ,5
+ ,0
+ ,18.780
+ ,6
+ ,0
+ ,23.520
+ ,9
+ ,0
+ ,26.040
+ ,9
+ ,0
+ ,14.000
+ ,5
+ ,0
+ ,12.560
+ ,5
+ ,0
+ ,0.839
+ ,4
+ ,0
+ ,25.780
+ ,7
+ ,0
+ ,29.880
+ ,9
+ ,0
+ ,14.040
+ ,3
+ ,0
+ ,23.480
+ ,9
+ ,0
+ ,17.550
+ ,5
+ ,0
+ ,29.800
+ ,8
+ ,0
+ ,21.000
+ ,9
+ ,0
+ ,12.820
+ ,5
+ ,0
+ ,30.000
+ ,9
+ ,0
+ ,26.730
+ ,8
+ ,0
+ ,20.930
+ ,7
+ ,0
+ ,16.120
+ ,5
+ ,0
+ ,21.750
+ ,8
+ ,0
+ ,27.250
+ ,9
+ ,0
+ ,20.710
+ ,8
+ ,0
+ ,15.470
+ ,8
+ ,0
+ ,20.040
+ ,8
+ ,0
+ ,31.350
+ ,9
+ ,0
+ ,24.200
+ ,8
+ ,0
+ ,17.760
+ ,5
+ ,0
+ ,19.310
+ ,8
+ ,0
+ ,13.430
+ ,5
+ ,0
+ ,20.760
+ ,9
+ ,0
+ ,16.240
+ ,7
+ ,0
+ ,13.440
+ ,8
+ ,0
+ ,16.500
+ ,6
+ ,0
+ ,27.320
+ ,8
+ ,0
+ ,20.170
+ ,5
+ ,0
+ ,27.970
+ ,9
+ ,0
+ ,35.560
+ ,9
+ ,0
+ ,17.030
+ ,8
+ ,0
+ ,16.340
+ ,6
+ ,0
+ ,25.700
+ ,9
+ ,0
+ ,30.160
+ ,9
+ ,0
+ ,24.190
+ ,7
+ ,0
+ ,15.690
+ ,4
+ ,0
+ ,16.980
+ ,8
+ ,0
+ ,21.230
+ ,8
+ ,0
+ ,24.810
+ ,8
+ ,0
+ ,14.810
+ ,6
+ ,0
+ ,15.770
+ ,4
+ ,0
+ ,19.400
+ ,8
+ ,0
+ ,17.470
+ ,6
+ ,0
+ ,20.690
+ ,9
+ ,0
+ ,16.310
+ ,7
+ ,0
+ ,15.360
+ ,5
+ ,0
+ ,25.600
+ ,9
+ ,0
+ ,19.620
+ ,8
+ ,0
+ ,25.310
+ ,8
+ ,0
+ ,27.150
+ ,9
+ ,0
+ ,24.570
+ ,9
+ ,0
+ ,20.900
+ ,9
+ ,0
+ ,17.890
+ ,7
+ ,0
+ ,18.580
+ ,5
+ ,0
+ ,14.520
+ ,5
+ ,0
+ ,38.420
+ ,9
+ ,0
+ ,17.190
+ ,6
+ ,0
+ ,21.110
+ ,7
+ ,0
+ ,16.950
+ ,6
+ ,0
+ ,22.110
+ ,8
+ ,0
+ ,17.940
+ ,8
+ ,0
+ ,19.170
+ ,7
+ ,0
+ ,21.440
+ ,8
+ ,0
+ ,12.530
+ ,7
+ ,0
+ ,26.590
+ ,9
+ ,0
+ ,15.800
+ ,5
+ ,0
+ ,21.260
+ ,9
+ ,0
+ ,30.290
+ ,9
+ ,0
+ ,29.640
+ ,9
+ ,0
+ ,16.110
+ ,7
+ ,0
+ ,22.150
+ ,8
+ ,0
+ ,23.880
+ ,8
+ ,0
+ ,21.960
+ ,9
+ ,0
+ ,17.510
+ ,9
+ ,0
+ ,21.650
+ ,9
+ ,0
+ ,16.820
+ ,7
+ ,0
+ ,15.230
+ ,8
+ ,0
+ ,12.920
+ ,8
+ ,0
+ ,16.490
+ ,7
+ ,0
+ ,25.880
+ ,9
+ ,0
+ ,0.796
+ ,4
+ ,0
+ ,25.740
+ ,9
+ ,0
+ ,19.790
+ ,6
+ ,0
+ ,23.540
+ ,8
+ ,0
+ ,17.180
+ ,6
+ ,0
+ ,17.420
+ ,7
+ ,0
+ ,16.030
+ ,7
+ ,0
+ ,26.390
+ ,8
+ ,0
+ ,18.290
+ ,7
+ ,0
+ ,20.840
+ ,7
+ ,0
+ ,22.200
+ ,7
+ ,0
+ ,14.730
+ ,7
+ ,0
+ ,23.410
+ ,8
+ ,0
+ ,16.980
+ ,7
+ ,0
+ ,11.960
+ ,5
+ ,0
+ ,18.720
+ ,8
+ ,0
+ ,22.190
+ ,7
+ ,0
+ ,24.200
+ ,9
+ ,0
+ ,18.270
+ ,7
+ ,0
+ ,14.610
+ ,7
+ ,0
+ ,13.380
+ ,6
+ ,0
+ ,20.900
+ ,8
+ ,0
+ ,16.970
+ ,8
+ ,0
+ ,15.620
+ ,8
+ ,0
+ ,20.400
+ ,9
+ ,0
+ ,16.090
+ ,7
+ ,0
+ ,24.580
+ ,8
+ ,0
+ ,26.500
+ ,9
+ ,0
+ ,14.290
+ ,8
+ ,0
+ ,16.750
+ ,8
+ ,0
+ ,19.470
+ ,9
+ ,0
+ ,20.690
+ ,8
+ ,0
+ ,15.720
+ ,6
+ ,0
+ ,13.480
+ ,6
+ ,0
+ ,22.880
+ ,8
+ ,0
+ ,17.730
+ ,9
+ ,0
+ ,0.791
+ ,5
+ ,0
+ ,19.050
+ ,7
+ ,0
+ ,24.630
+ ,9
+ ,0
+ ,14.310
+ ,6
+ ,0
+ ,26.310
+ ,9
+ ,0
+ ,31.140
+ ,9
+ ,0
+ ,21.350
+ ,9
+ ,0
+ ,15.270
+ ,6
+ ,0
+ ,22.930
+ ,8
+ ,0
+ ,30.420
+ ,9
+ ,0
+ ,29.270
+ ,8
+ ,0
+ ,26.650
+ ,8
+ ,0
+ ,23.010
+ ,9
+ ,0
+ ,24.600
+ ,9
+ ,0
+ ,25.920
+ ,9
+ ,0
+ ,17.500
+ ,7
+ ,0
+ ,17.590
+ ,8
+ ,0
+ ,15.360
+ ,6
+ ,0
+ ,22.590
+ ,9
+ ,0
+ ,20.480
+ ,9
+ ,0
+ ,25.710
+ ,9
+ ,0
+ ,20.460
+ ,7
+ ,0
+ ,17.800
+ ,8
+ ,0
+ ,15.520
+ ,5
+ ,0
+ ,19.530
+ ,8
+ ,0
+ ,28.930
+ ,9
+ ,0
+ ,17.130
+ ,6
+ ,0
+ ,28.510
+ ,9
+ ,0
+ ,16.240
+ ,6
+ ,0
+ ,26.310
+ ,8
+ ,0
+ ,18.190
+ ,5
+ ,0
+ ,16.580
+ ,7
+ ,0
+ ,21.580
+ ,7
+ ,0
+ ,17.890
+ ,4
+ ,0
+ ,30.040
+ ,9
+ ,0
+ ,25.030
+ ,8
+ ,0
+ ,19.330
+ ,9
+ ,0
+ ,20.910
+ ,9
+ ,0
+ ,23.160
+ ,9
+ ,0
+ ,17.040
+ ,5
+ ,0
+ ,16.060
+ ,9
+ ,0
+ ,11.650
+ ,7
+ ,0
+ ,21.020
+ ,6
+ ,0
+ ,23.200
+ ,9
+ ,0
+ ,22.300
+ ,9
+ ,0
+ ,17.160
+ ,9
+ ,0
+ ,17.900
+ ,7
+ ,0
+ ,11.460
+ ,5
+ ,0
+ ,21.870
+ ,8
+ ,0
+ ,27.170
+ ,9
+ ,0
+ ,17.960
+ ,7
+ ,1
+ ,19.530
+ ,9
+ ,0
+ ,13.350
+ ,8
+ ,0
+ ,21.190
+ ,9
+ ,0
+ ,16.660
+ ,6
+ ,0
+ ,18.260
+ ,6
+ ,0
+ ,27.090
+ ,8
+ ,0
+ ,28.710
+ ,9
+ ,0
+ ,10.920
+ ,5
+ ,0
+ ,22.620
+ ,6
+ ,0
+ ,21.040
+ ,6
+ ,0
+ ,21.660
+ ,9
+ ,0
+ ,16.900
+ ,7
+ ,0
+ ,29.730
+ ,9
+ ,0
+ ,21.450
+ ,8
+ ,0
+ ,19.710
+ ,5
+ ,0
+ ,20.950
+ ,7
+ ,0
+ ,16.970
+ ,6
+ ,0
+ ,24.550
+ ,9
+ ,0
+ ,19.200
+ ,7
+ ,0
+ ,21.640
+ ,9
+ ,0
+ ,21.300
+ ,9
+ ,0
+ ,29.930
+ ,8
+ ,0
+ ,25.290
+ ,9
+ ,0
+ ,17.260
+ ,7
+ ,0
+ ,24.420
+ ,9
+ ,0
+ ,11.020
+ ,4
+ ,0
+ ,20.560
+ ,9
+ ,0
+ ,18.080
+ ,5
+ ,0
+ ,23.050
+ ,8
+ ,0
+ ,19.690
+ ,9
+ ,0
+ ,15.560
+ ,8
+ ,0
+ ,10.720
+ ,3
+ ,0
+ ,20.420
+ ,9
+ ,0
+ ,15.120
+ ,8
+ ,0
+ ,14.230
+ ,6
+ ,0
+ ,36.810
+ ,9
+ ,0
+ ,19.910
+ ,8
+ ,0
+ ,18.970
+ ,8
+ ,0
+ ,13.700
+ ,7
+ ,0
+ ,13.380
+ ,6
+ ,0
+ ,20.160
+ ,8
+ ,0
+ ,26.390
+ ,9
+ ,0
+ ,13.890
+ ,4
+ ,0
+ ,16.120
+ ,7
+ ,0
+ ,21.350
+ ,8
+ ,0
+ ,26.810
+ ,8
+ ,0
+ ,32.230
+ ,9
+ ,0
+ ,17.960
+ ,6
+ ,0
+ ,20.100
+ ,8
+ ,0
+ ,15.230
+ ,6
+ ,0
+ ,17.440
+ ,8
+ ,0
+ ,24.850
+ ,9
+ ,0
+ ,23.350
+ ,8
+ ,0
+ ,14.150
+ ,7
+ ,0
+ ,20.760
+ ,9
+ ,0
+ ,24.350
+ ,8
+ ,0
+ ,17.280
+ ,7
+ ,0
+ ,28.500
+ ,9
+ ,0
+ ,18.440
+ ,8
+ ,0
+ ,17.540
+ ,9
+ ,0
+ ,13.430
+ ,6
+ ,0
+ ,23.030
+ ,8
+ ,0
+ ,22.460
+ ,9
+ ,0
+ ,24.760
+ ,8
+ ,0
+ ,32.390
+ ,9
+ ,0
+ ,24.570
+ ,9
+ ,0
+ ,23.820
+ ,8
+ ,0
+ ,16.400
+ ,7
+ ,0
+ ,15.890
+ ,5
+ ,0
+ ,20.560
+ ,7
+ ,0
+ ,22.260
+ ,8
+ ,0
+ ,18.860
+ ,9
+ ,0
+ ,28.330
+ ,9
+ ,0
+ ,17.150
+ ,6
+ ,0
+ ,26.310
+ ,8
+ ,0
+ ,25.500
+ ,7
+ ,0
+ ,19.120
+ ,9
+ ,0
+ ,18.770
+ ,7
+ ,0
+ ,19.350
+ ,7
+ ,0
+ ,15.390
+ ,5
+ ,0
+ ,28.030
+ ,9
+ ,0
+ ,29.230
+ ,9
+ ,0
+ ,23.580
+ ,8
+ ,0
+ ,20.940
+ ,8
+ ,0
+ ,18.550
+ ,9
+ ,0
+ ,15.350
+ ,6
+ ,0
+ ,21.350
+ ,7
+ ,0
+ ,19.300
+ ,5
+ ,0
+ ,21.820
+ ,9
+ ,0
+ ,13.590
+ ,5
+ ,0
+ ,20.020
+ ,7
+ ,0
+ ,16.990
+ ,6
+ ,0
+ ,25.000
+ ,8
+ ,0
+ ,23.660
+ ,7
+ ,0
+ ,20.690
+ ,8
+ ,0
+ ,14.180
+ ,4
+ ,0
+ ,23.330
+ ,8
+ ,0
+ ,15.140
+ ,5
+ ,0
+ ,17.580
+ ,8
+ ,0
+ ,25.350
+ ,7
+ ,0
+ ,25.640
+ ,7
+ ,0
+ ,24.870
+ ,9
+ ,0
+ ,15.910
+ ,9
+ ,0
+ ,16.240
+ ,8
+ ,0
+ ,27.980
+ ,9
+ ,0
+ ,16.910
+ ,6
+ ,0
+ ,19.990
+ ,8
+ ,0
+ ,18.690
+ ,9
+ ,0
+ ,10.040
+ ,4
+ ,0
+ ,14.270
+ ,6
+ ,0
+ ,18.260
+ ,7
+ ,0
+ ,26.880
+ ,9
+ ,0
+ ,16.570
+ ,8
+ ,0
+ ,16.720
+ ,6
+ ,0
+ ,20.150
+ ,8
+ ,0
+ ,23.710
+ ,7
+ ,0
+ ,21.150
+ ,5
+ ,0
+ ,23.280
+ ,8
+ ,0
+ ,14.950
+ ,7
+ ,0
+ ,28.840
+ ,11
+ ,0
+ ,23.280
+ ,10
+ ,0
+ ,33.810
+ ,14
+ ,0
+ ,21.700
+ ,11
+ ,0
+ ,34.700
+ ,11
+ ,0
+ ,30.580
+ ,12
+ ,0
+ ,18.110
+ ,10
+ ,0
+ ,25.240
+ ,11
+ ,0
+ ,26.420
+ ,10
+ ,0
+ ,37.410
+ ,14
+ ,0
+ ,43.360
+ ,13
+ ,0
+ ,48.420
+ ,14
+ ,0
+ ,45.500
+ ,12
+ ,0
+ ,28.410
+ ,12
+ ,0
+ ,31.660
+ ,10
+ ,0
+ ,38.160
+ ,13
+ ,0
+ ,25.610
+ ,10
+ ,0
+ ,36.540
+ ,11
+ ,0
+ ,24.810
+ ,10
+ ,0
+ ,26.650
+ ,11
+ ,0
+ ,32.030
+ ,10
+ ,0
+ ,35.490
+ ,13
+ ,1
+ ,22.360
+ ,14
+ ,0
+ ,32.220
+ ,11
+ ,0
+ ,31.110
+ ,10
+ ,0
+ ,34.900
+ ,11
+ ,0
+ ,31.470
+ ,13
+ ,0
+ ,25.200
+ ,10
+ ,0
+ ,22.920
+ ,10
+ ,0
+ ,28.890
+ ,12
+ ,0
+ ,22.460
+ ,10
+ ,0
+ ,19.370
+ ,10
+ ,0
+ ,26.460
+ ,10
+ ,0
+ ,29.570
+ ,11
+ ,0
+ ,40.070
+ ,11
+ ,0
+ ,23.860
+ ,11
+ ,0
+ ,32.510
+ ,10
+ ,0
+ ,27.620
+ ,11
+ ,0
+ ,30.110
+ ,11
+ ,0
+ ,43.050
+ ,13
+ ,0
+ ,39.060
+ ,13
+ ,0
+ ,35.830
+ ,11
+ ,0
+ ,32.360
+ ,11
+ ,0
+ ,34.360
+ ,14
+ ,0
+ ,30.580
+ ,11
+ ,0
+ ,30.070
+ ,10
+ ,0
+ ,34.890
+ ,10
+ ,0
+ ,28.640
+ ,10
+ ,1
+ ,34.280
+ ,14
+ ,0
+ ,28.190
+ ,13
+ ,0
+ ,22.500
+ ,10
+ ,0
+ ,46.830
+ ,14
+ ,0
+ ,23.520
+ ,10
+ ,0
+ ,31.080
+ ,11
+ ,0
+ ,39.940
+ ,13
+ ,0
+ ,43.930
+ ,12
+ ,1
+ ,32.080
+ ,13
+ ,0
+ ,25.920
+ ,10
+ ,0
+ ,31.930
+ ,13
+ ,1
+ ,16.940
+ ,11
+ ,1
+ ,39.570
+ ,14
+ ,0
+ ,23.460
+ ,11
+ ,1
+ ,47.890
+ ,13
+ ,0
+ ,35.150
+ ,11
+ ,0
+ ,27.540
+ ,11
+ ,0
+ ,27.200
+ ,10
+ ,0
+ ,24.630
+ ,11
+ ,0
+ ,26.330
+ ,11
+ ,0
+ ,30.480
+ ,10
+ ,0
+ ,31.110
+ ,11
+ ,0
+ ,37.450
+ ,13
+ ,1
+ ,23.840
+ ,12
+ ,0
+ ,20.940
+ ,10
+ ,0
+ ,31.830
+ ,10
+ ,1
+ ,30.740
+ ,14
+ ,0
+ ,39.770
+ ,11
+ ,0
+ ,33.540
+ ,10
+ ,0
+ ,34.110
+ ,11
+ ,1
+ ,23.870
+ ,10
+ ,0
+ ,31.710
+ ,11
+ ,0
+ ,38.870
+ ,13
+ ,0
+ ,26.460
+ ,13
+ ,0
+ ,25.040
+ ,10
+ ,0
+ ,35.870
+ ,11
+ ,0
+ ,38.450
+ ,11
+ ,0
+ ,29.710
+ ,12
+ ,0
+ ,28.910
+ ,10
+ ,0
+ ,18.230
+ ,10
+ ,0
+ ,24.170
+ ,11
+ ,0
+ ,21.750
+ ,10
+ ,0
+ ,27.350
+ ,11
+ ,0
+ ,42.730
+ ,14
+ ,0
+ ,29.760
+ ,13
+ ,1
+ ,38.350
+ ,12
+ ,0
+ ,40.650
+ ,11
+ ,0
+ ,23.180
+ ,11
+ ,0
+ ,35.960
+ ,11
+ ,0
+ ,33.950
+ ,14
+ ,0
+ ,27.510
+ ,12
+ ,0
+ ,26.730
+ ,10
+ ,0
+ ,25.560
+ ,12
+ ,0
+ ,25.420
+ ,11
+ ,0
+ ,26.080
+ ,10
+ ,0
+ ,23.540
+ ,11
+ ,1
+ ,25.990
+ ,13
+ ,0
+ ,14.580
+ ,10
+ ,0
+ ,37.950
+ ,10
+ ,0
+ ,24.910
+ ,11
+ ,0
+ ,30.600
+ ,13
+ ,0
+ ,25.450
+ ,10
+ ,0
+ ,29.930
+ ,11
+ ,0
+ ,33.050
+ ,10
+ ,1
+ ,47.560
+ ,13
+ ,0
+ ,37.740
+ ,11
+ ,0
+ ,28.550
+ ,10
+ ,0
+ ,29.880
+ ,11
+ ,0
+ ,24.980
+ ,11
+ ,0
+ ,31.690
+ ,14
+ ,0
+ ,28.870
+ ,11
+ ,0
+ ,27.040
+ ,13
+ ,0
+ ,35.150
+ ,11
+ ,0
+ ,34.250
+ ,11
+ ,0
+ ,22.870
+ ,10
+ ,0
+ ,24.340
+ ,13
+ ,0
+ ,23.650
+ ,10
+ ,1
+ ,30.860
+ ,13
+ ,0
+ ,26.960
+ ,10
+ ,0
+ ,28.680
+ ,12
+ ,0
+ ,28.130
+ ,10
+ ,1
+ ,43.090
+ ,14
+ ,0
+ ,32.550
+ ,12
+ ,1
+ ,34.130
+ ,10
+ ,0
+ ,45.930
+ ,11
+ ,0
+ ,41.110
+ ,14
+ ,0
+ ,19.160
+ ,12
+ ,0
+ ,18.580
+ ,10
+ ,1
+ ,29.750
+ ,10
+ ,0
+ ,33.500
+ ,10
+ ,0
+ ,29.010
+ ,10
+ ,0
+ ,22.410
+ ,12
+ ,0
+ ,42.250
+ ,13
+ ,0
+ ,32.230
+ ,11
+ ,0
+ ,52.240
+ ,12
+ ,0
+ ,40.730
+ ,11
+ ,0
+ ,40.800
+ ,12
+ ,0
+ ,26.060
+ ,11
+ ,1
+ ,31.690
+ ,11
+ ,0
+ ,44.110
+ ,12
+ ,0
+ ,37.910
+ ,12
+ ,0
+ ,30.890
+ ,13
+ ,0
+ ,24.650
+ ,11
+ ,1
+ ,33.430
+ ,12
+ ,0
+ ,32.000
+ ,10
+ ,0
+ ,29.130
+ ,12
+ ,0
+ ,48.770
+ ,13
+ ,0
+ ,23.580
+ ,10
+ ,0
+ ,32.790
+ ,12
+ ,0
+ ,25.810
+ ,10
+ ,0
+ ,23.470
+ ,12
+ ,0
+ ,26.910
+ ,10
+ ,0
+ ,28.270
+ ,11
+ ,0
+ ,18.730
+ ,10
+ ,1
+ ,37.510
+ ,12
+ ,0
+ ,25.380
+ ,14
+ ,0
+ ,27.580
+ ,10
+ ,0
+ ,30.500
+ ,10
+ ,0
+ ,30.790
+ ,12
+ ,0
+ ,22.010
+ ,10
+ ,0
+ ,18.580
+ ,10
+ ,1
+ ,22.160
+ ,13
+ ,0
+ ,34.030
+ ,12
+ ,0
+ ,35.010
+ ,12
+ ,0
+ ,25.780
+ ,11
+ ,1
+ ,30.780
+ ,13
+ ,1
+ ,31.860
+ ,12
+ ,0
+ ,16.650
+ ,10
+ ,0
+ ,20.810
+ ,11
+ ,0
+ ,29.740
+ ,11
+ ,1
+ ,32.970
+ ,13
+ ,0
+ ,40.730
+ ,12
+ ,0
+ ,44.480
+ ,13
+ ,0
+ ,39.840
+ ,13
+ ,0
+ ,22.500
+ ,10
+ ,0
+ ,27.520
+ ,12
+ ,1
+ ,23.040
+ ,12
+ ,0
+ ,36.800
+ ,14
+ ,1
+ ,31.020
+ ,11
+ ,0
+ ,28.620
+ ,10
+ ,1
+ ,26.770
+ ,13
+ ,0
+ ,30.230
+ ,11
+ ,0
+ ,36.810
+ ,11
+ ,0
+ ,32.550
+ ,13
+ ,0
+ ,36.920
+ ,12
+ ,0
+ ,23.560
+ ,10
+ ,0
+ ,45.910
+ ,10
+ ,0
+ ,30.820
+ ,12
+ ,1
+ ,32.970
+ ,13
+ ,0
+ ,32.580
+ ,11
+ ,0
+ ,22.160
+ ,10
+ ,0
+ ,32.470
+ ,11
+ ,0
+ ,43.240
+ ,11
+ ,0
+ ,23.620
+ ,11
+ ,0
+ ,25.630
+ ,11
+ ,0
+ ,32.060
+ ,11
+ ,0
+ ,35.850
+ ,14
+ ,0
+ ,47.200
+ ,12
+ ,0
+ ,33.310
+ ,13
+ ,0
+ ,50.830
+ ,13
+ ,1
+ ,34.980
+ ,10
+ ,0
+ ,24.170
+ ,12
+ ,0
+ ,23.640
+ ,10
+ ,0
+ ,23.410
+ ,10
+ ,1
+ ,27.590
+ ,12
+ ,1
+ ,29.530
+ ,11
+ ,0
+ ,32.310
+ ,12
+ ,0
+ ,30.780
+ ,11
+ ,0
+ ,33.690
+ ,11
+ ,0
+ ,35.290
+ ,12
+ ,0
+ ,28.660
+ ,12
+ ,0
+ ,28.910
+ ,14
+ ,0
+ ,30.220
+ ,11
+ ,0
+ ,31.270
+ ,10
+ ,0
+ ,28.660
+ ,11
+ ,0
+ ,26.050
+ ,12
+ ,0
+ ,30.560
+ ,13
+ ,0
+ ,25.690
+ ,12
+ ,0
+ ,25.010
+ ,11
+ ,0
+ ,33.200
+ ,11
+ ,0
+ ,21.230
+ ,11
+ ,0
+ ,37.800
+ ,14
+ ,0
+ ,38.470
+ ,11
+ ,1
+ ,37.850
+ ,13
+ ,0
+ ,39.240
+ ,12
+ ,0
+ ,21.320
+ ,10
+ ,0
+ ,27.520
+ ,12
+ ,0
+ ,24.490
+ ,13
+ ,0
+ ,34.560
+ ,10
+ ,0
+ ,30.730
+ ,10
+ ,0
+ ,26.880
+ ,10
+ ,0
+ ,33.290
+ ,10
+ ,0
+ ,42.710
+ ,14
+ ,0
+ ,35.300
+ ,12
+ ,0
+ ,29.280
+ ,11
+ ,0
+ ,26.890
+ ,11
+ ,0
+ ,23.320
+ ,12
+ ,0
+ ,29.340
+ ,14
+ ,1
+ ,22.760
+ ,14
+ ,0
+ ,31.100
+ ,10
+ ,0
+ ,28.940
+ ,11
+ ,1
+ ,46.370
+ ,11
+ ,0
+ ,24.350
+ ,10
+ ,0
+ ,28.380
+ ,10
+ ,0
+ ,30.350
+ ,12
+ ,0
+ ,48.310
+ ,12
+ ,0
+ ,28.120
+ ,11
+ ,0
+ ,27.140
+ ,12
+ ,0
+ ,30.860
+ ,10
+ ,0
+ ,35.190
+ ,12
+ ,0
+ ,42.320
+ ,13
+ ,0
+ ,27.700
+ ,10
+ ,0
+ ,33.410
+ ,12
+ ,0
+ ,30.900
+ ,10
+ ,0
+ ,25.310
+ ,13
+ ,0
+ ,28.220
+ ,12
+ ,1
+ ,30.380
+ ,10
+ ,0
+ ,29.350
+ ,12
+ ,0
+ ,25.680
+ ,10
+ ,0
+ ,23.870
+ ,11
+ ,0
+ ,24.990
+ ,12
+ ,0
+ ,41.300
+ ,11
+ ,0
+ ,30.010
+ ,12
+ ,0
+ ,31.320
+ ,10
+ ,0
+ ,35.770
+ ,13
+ ,0
+ ,32.220
+ ,12
+ ,0
+ ,32.800
+ ,11
+ ,0
+ ,26.590
+ ,11
+ ,0
+ ,28.220
+ ,11
+ ,0
+ ,21.400
+ ,11
+ ,0
+ ,42.030
+ ,12
+ ,0
+ ,29.970
+ ,14
+ ,1
+ ,31.200
+ ,11
+ ,0
+ ,25.620
+ ,11
+ ,0
+ ,30.820
+ ,12
+ ,0
+ ,38.060
+ ,14
+ ,1
+ ,33.390
+ ,11
+ ,1
+ ,31.520
+ ,13
+ ,0
+ ,24.580
+ ,11
+ ,0
+ ,23.910
+ ,10
+ ,0
+ ,31.410
+ ,13
+ ,0
+ ,25.790
+ ,12
+ ,1
+ ,31.040
+ ,11
+ ,1
+ ,40.450
+ ,13
+ ,1
+ ,47.630
+ ,14
+ ,0
+ ,21.000
+ ,10
+ ,1
+ ,30.690
+ ,11
+ ,0
+ ,27.850
+ ,11
+ ,0
+ ,42.840
+ ,15
+ ,1
+ ,45.060
+ ,15
+ ,0
+ ,29.060
+ ,18
+ ,0
+ ,51.020
+ ,19
+ ,1
+ ,35.190
+ ,19
+ ,1
+ ,36.880
+ ,16
+ ,0
+ ,44.290
+ ,17
+ ,0
+ ,42.790
+ ,15
+ ,0
+ ,45.000
+ ,15
+ ,0
+ ,26.350
+ ,15
+ ,1
+ ,26.790
+ ,15
+ ,1
+ ,21.980
+ ,15
+ ,1
+ ,33.450
+ ,19
+ ,0
+ ,30.820
+ ,18
+ ,0
+ ,33.870
+ ,16
+ ,1
+ ,30.820
+ ,17
+ ,1
+ ,29.030
+ ,16
+ ,1
+ ,30.040
+ ,15
+ ,0
+ ,57.930
+ ,15
+ ,0
+ ,39.850
+ ,15
+ ,0
+ ,42.200
+ ,18
+ ,0
+ ,47.240
+ ,17
+ ,0
+ ,37.310
+ ,15
+ ,1
+ ,34.060
+ ,17
+ ,0
+ ,35.000
+ ,17
+ ,0
+ ,36.740
+ ,16
+ ,0
+ ,56.330
+ ,17
+ ,1
+ ,31.220
+ ,15
+ ,1
+ ,33.300
+ ,15
+ ,1
+ ,26.080
+ ,16
+ ,0
+ ,36.450
+ ,16
+ ,1
+ ,37.990
+ ,15
+ ,1
+ ,40.860
+ ,18
+ ,0
+ ,28.870
+ ,15
+ ,1
+ ,40.700
+ ,16
+ ,0
+ ,39.600
+ ,17
+ ,0
+ ,42.990
+ ,16
+ ,0
+ ,29.810
+ ,16
+ ,1
+ ,22.640
+ ,15
+ ,1
+ ,44.040
+ ,18
+ ,1
+ ,22.780
+ ,15
+ ,0
+ ,45.040
+ ,16
+ ,0
+ ,56.380
+ ,17
+ ,1
+ ,48.720
+ ,16
+ ,1
+ ,42.700
+ ,16
+ ,1
+ ,37.270
+ ,15
+ ,0
+ ,28.530
+ ,18
+ ,1
+ ,27.950
+ ,16
+ ,0
+ ,32.110
+ ,15)
+ ,dim=c(3
+ ,654)
+ ,dimnames=list(c('Smoke'
+ ,'FEV'
+ ,'Age')
+ ,1:654))
> y <- array(NA,dim=c(3,654),dimnames=list(c('Smoke','FEV','Age'),1:654))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '2'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
Attaching package: 'zoo'
The following object(s) are masked from package:base :
as.Date.numeric
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- 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'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> 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[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
FEV Smoke Age
1 17.080 0 9
2 17.240 0 8
3 17.200 0 7
4 15.580 0 9
5 18.950 0 9
6 23.360 0 8
7 19.190 0 6
8 14.150 0 6
9 19.870 0 8
10 19.420 0 9
11 16.020 0 6
12 17.350 0 8
13 21.930 0 8
14 21.180 0 8
15 22.580 0 8
16 19.320 0 7
17 14.720 0 5
18 18.780 0 6
19 23.520 0 9
20 26.040 0 9
21 14.000 0 5
22 12.560 0 5
23 0.839 0 4
24 25.780 0 7
25 29.880 0 9
26 14.040 0 3
27 23.480 0 9
28 17.550 0 5
29 29.800 0 8
30 21.000 0 9
31 12.820 0 5
32 30.000 0 9
33 26.730 0 8
34 20.930 0 7
35 16.120 0 5
36 21.750 0 8
37 27.250 0 9
38 20.710 0 8
39 15.470 0 8
40 20.040 0 8
41 31.350 0 9
42 24.200 0 8
43 17.760 0 5
44 19.310 0 8
45 13.430 0 5
46 20.760 0 9
47 16.240 0 7
48 13.440 0 8
49 16.500 0 6
50 27.320 0 8
51 20.170 0 5
52 27.970 0 9
53 35.560 0 9
54 17.030 0 8
55 16.340 0 6
56 25.700 0 9
57 30.160 0 9
58 24.190 0 7
59 15.690 0 4
60 16.980 0 8
61 21.230 0 8
62 24.810 0 8
63 14.810 0 6
64 15.770 0 4
65 19.400 0 8
66 17.470 0 6
67 20.690 0 9
68 16.310 0 7
69 15.360 0 5
70 25.600 0 9
71 19.620 0 8
72 25.310 0 8
73 27.150 0 9
74 24.570 0 9
75 20.900 0 9
76 17.890 0 7
77 18.580 0 5
78 14.520 0 5
79 38.420 0 9
80 17.190 0 6
81 21.110 0 7
82 16.950 0 6
83 22.110 0 8
84 17.940 0 8
85 19.170 0 7
86 21.440 0 8
87 12.530 0 7
88 26.590 0 9
89 15.800 0 5
90 21.260 0 9
91 30.290 0 9
92 29.640 0 9
93 16.110 0 7
94 22.150 0 8
95 23.880 0 8
96 21.960 0 9
97 17.510 0 9
98 21.650 0 9
99 16.820 0 7
100 15.230 0 8
101 12.920 0 8
102 16.490 0 7
103 25.880 0 9
104 0.796 0 4
105 25.740 0 9
106 19.790 0 6
107 23.540 0 8
108 17.180 0 6
109 17.420 0 7
110 16.030 0 7
111 26.390 0 8
112 18.290 0 7
113 20.840 0 7
114 22.200 0 7
115 14.730 0 7
116 23.410 0 8
117 16.980 0 7
118 11.960 0 5
119 18.720 0 8
120 22.190 0 7
121 24.200 0 9
122 18.270 0 7
123 14.610 0 7
124 13.380 0 6
125 20.900 0 8
126 16.970 0 8
127 15.620 0 8
128 20.400 0 9
129 16.090 0 7
130 24.580 0 8
131 26.500 0 9
132 14.290 0 8
133 16.750 0 8
134 19.470 0 9
135 20.690 0 8
136 15.720 0 6
137 13.480 0 6
138 22.880 0 8
139 17.730 0 9
140 0.791 0 5
141 19.050 0 7
142 24.630 0 9
143 14.310 0 6
144 26.310 0 9
145 31.140 0 9
146 21.350 0 9
147 15.270 0 6
148 22.930 0 8
149 30.420 0 9
150 29.270 0 8
151 26.650 0 8
152 23.010 0 9
153 24.600 0 9
154 25.920 0 9
155 17.500 0 7
156 17.590 0 8
157 15.360 0 6
158 22.590 0 9
159 20.480 0 9
160 25.710 0 9
161 20.460 0 7
162 17.800 0 8
163 15.520 0 5
164 19.530 0 8
165 28.930 0 9
166 17.130 0 6
167 28.510 0 9
168 16.240 0 6
169 26.310 0 8
170 18.190 0 5
171 16.580 0 7
172 21.580 0 7
173 17.890 0 4
174 30.040 0 9
175 25.030 0 8
176 19.330 0 9
177 20.910 0 9
178 23.160 0 9
179 17.040 0 5
180 16.060 0 9
181 11.650 0 7
182 21.020 0 6
183 23.200 0 9
184 22.300 0 9
185 17.160 0 9
186 17.900 0 7
187 11.460 0 5
188 21.870 0 8
189 27.170 0 9
190 17.960 0 7
191 19.530 1 9
192 13.350 0 8
193 21.190 0 9
194 16.660 0 6
195 18.260 0 6
196 27.090 0 8
197 28.710 0 9
198 10.920 0 5
199 22.620 0 6
200 21.040 0 6
201 21.660 0 9
202 16.900 0 7
203 29.730 0 9
204 21.450 0 8
205 19.710 0 5
206 20.950 0 7
207 16.970 0 6
208 24.550 0 9
209 19.200 0 7
210 21.640 0 9
211 21.300 0 9
212 29.930 0 8
213 25.290 0 9
214 17.260 0 7
215 24.420 0 9
216 11.020 0 4
217 20.560 0 9
218 18.080 0 5
219 23.050 0 8
220 19.690 0 9
221 15.560 0 8
222 10.720 0 3
223 20.420 0 9
224 15.120 0 8
225 14.230 0 6
226 36.810 0 9
227 19.910 0 8
228 18.970 0 8
229 13.700 0 7
230 13.380 0 6
231 20.160 0 8
232 26.390 0 9
233 13.890 0 4
234 16.120 0 7
235 21.350 0 8
236 26.810 0 8
237 32.230 0 9
238 17.960 0 6
239 20.100 0 8
240 15.230 0 6
241 17.440 0 8
242 24.850 0 9
243 23.350 0 8
244 14.150 0 7
245 20.760 0 9
246 24.350 0 8
247 17.280 0 7
248 28.500 0 9
249 18.440 0 8
250 17.540 0 9
251 13.430 0 6
252 23.030 0 8
253 22.460 0 9
254 24.760 0 8
255 32.390 0 9
256 24.570 0 9
257 23.820 0 8
258 16.400 0 7
259 15.890 0 5
260 20.560 0 7
261 22.260 0 8
262 18.860 0 9
263 28.330 0 9
264 17.150 0 6
265 26.310 0 8
266 25.500 0 7
267 19.120 0 9
268 18.770 0 7
269 19.350 0 7
270 15.390 0 5
271 28.030 0 9
272 29.230 0 9
273 23.580 0 8
274 20.940 0 8
275 18.550 0 9
276 15.350 0 6
277 21.350 0 7
278 19.300 0 5
279 21.820 0 9
280 13.590 0 5
281 20.020 0 7
282 16.990 0 6
283 25.000 0 8
284 23.660 0 7
285 20.690 0 8
286 14.180 0 4
287 23.330 0 8
288 15.140 0 5
289 17.580 0 8
290 25.350 0 7
291 25.640 0 7
292 24.870 0 9
293 15.910 0 9
294 16.240 0 8
295 27.980 0 9
296 16.910 0 6
297 19.990 0 8
298 18.690 0 9
299 10.040 0 4
300 14.270 0 6
301 18.260 0 7
302 26.880 0 9
303 16.570 0 8
304 16.720 0 6
305 20.150 0 8
306 23.710 0 7
307 21.150 0 5
308 23.280 0 8
309 14.950 0 7
310 28.840 0 11
311 23.280 0 10
312 33.810 0 14
313 21.700 0 11
314 34.700 0 11
315 30.580 0 12
316 18.110 0 10
317 25.240 0 11
318 26.420 0 10
319 37.410 0 14
320 43.360 0 13
321 48.420 0 14
322 45.500 0 12
323 28.410 0 12
324 31.660 0 10
325 38.160 0 13
326 25.610 0 10
327 36.540 0 11
328 24.810 0 10
329 26.650 0 11
330 32.030 0 10
331 35.490 0 13
332 22.360 1 14
333 32.220 0 11
334 31.110 0 10
335 34.900 0 11
336 31.470 0 13
337 25.200 0 10
338 22.920 0 10
339 28.890 0 12
340 22.460 0 10
341 19.370 0 10
342 26.460 0 10
343 29.570 0 11
344 40.070 0 11
345 23.860 0 11
346 32.510 0 10
347 27.620 0 11
348 30.110 0 11
349 43.050 0 13
350 39.060 0 13
351 35.830 0 11
352 32.360 0 11
353 34.360 0 14
354 30.580 0 11
355 30.070 0 10
356 34.890 0 10
357 28.640 0 10
358 34.280 1 14
359 28.190 0 13
360 22.500 0 10
361 46.830 0 14
362 23.520 0 10
363 31.080 0 11
364 39.940 0 13
365 43.930 0 12
366 32.080 1 13
367 25.920 0 10
368 31.930 0 13
369 16.940 1 11
370 39.570 1 14
371 23.460 0 11
372 47.890 1 13
373 35.150 0 11
374 27.540 0 11
375 27.200 0 10
376 24.630 0 11
377 26.330 0 11
378 30.480 0 10
379 31.110 0 11
380 37.450 0 13
381 23.840 1 12
382 20.940 0 10
383 31.830 0 10
384 30.740 1 14
385 39.770 0 11
386 33.540 0 10
387 34.110 0 11
388 23.870 1 10
389 31.710 0 11
390 38.870 0 13
391 26.460 0 13
392 25.040 0 10
393 35.870 0 11
394 38.450 0 11
395 29.710 0 12
396 28.910 0 10
397 18.230 0 10
398 24.170 0 11
399 21.750 0 10
400 27.350 0 11
401 42.730 0 14
402 29.760 0 13
403 38.350 1 12
404 40.650 0 11
405 23.180 0 11
406 35.960 0 11
407 33.950 0 14
408 27.510 0 12
409 26.730 0 10
410 25.560 0 12
411 25.420 0 11
412 26.080 0 10
413 23.540 0 11
414 25.990 1 13
415 14.580 0 10
416 37.950 0 10
417 24.910 0 11
418 30.600 0 13
419 25.450 0 10
420 29.930 0 11
421 33.050 0 10
422 47.560 1 13
423 37.740 0 11
424 28.550 0 10
425 29.880 0 11
426 24.980 0 11
427 31.690 0 14
428 28.870 0 11
429 27.040 0 13
430 35.150 0 11
431 34.250 0 11
432 22.870 0 10
433 24.340 0 13
434 23.650 0 10
435 30.860 1 13
436 26.960 0 10
437 28.680 0 12
438 28.130 0 10
439 43.090 1 14
440 32.550 0 12
441 34.130 1 10
442 45.930 0 11
443 41.110 0 14
444 19.160 0 12
445 18.580 0 10
446 29.750 1 10
447 33.500 0 10
448 29.010 0 10
449 22.410 0 12
450 42.250 0 13
451 32.230 0 11
452 52.240 0 12
453 40.730 0 11
454 40.800 0 12
455 26.060 0 11
456 31.690 1 11
457 44.110 0 12
458 37.910 0 12
459 30.890 0 13
460 24.650 0 11
461 33.430 1 12
462 32.000 0 10
463 29.130 0 12
464 48.770 0 13
465 23.580 0 10
466 32.790 0 12
467 25.810 0 10
468 23.470 0 12
469 26.910 0 10
470 28.270 0 11
471 18.730 0 10
472 37.510 1 12
473 25.380 0 14
474 27.580 0 10
475 30.500 0 10
476 30.790 0 12
477 22.010 0 10
478 18.580 0 10
479 22.160 1 13
480 34.030 0 12
481 35.010 0 12
482 25.780 0 11
483 30.780 1 13
484 31.860 1 12
485 16.650 0 10
486 20.810 0 11
487 29.740 0 11
488 32.970 1 13
489 40.730 0 12
490 44.480 0 13
491 39.840 0 13
492 22.500 0 10
493 27.520 0 12
494 23.040 1 12
495 36.800 0 14
496 31.020 1 11
497 28.620 0 10
498 26.770 1 13
499 30.230 0 11
500 36.810 0 11
501 32.550 0 13
502 36.920 0 12
503 23.560 0 10
504 45.910 0 10
505 30.820 0 12
506 32.970 1 13
507 32.580 0 11
508 22.160 0 10
509 32.470 0 11
510 43.240 0 11
511 23.620 0 11
512 25.630 0 11
513 32.060 0 11
514 35.850 0 14
515 47.200 0 12
516 33.310 0 13
517 50.830 0 13
518 34.980 1 10
519 24.170 0 12
520 23.640 0 10
521 23.410 0 10
522 27.590 1 12
523 29.530 1 11
524 32.310 0 12
525 30.780 0 11
526 33.690 0 11
527 35.290 0 12
528 28.660 0 12
529 28.910 0 14
530 30.220 0 11
531 31.270 0 10
532 28.660 0 11
533 26.050 0 12
534 30.560 0 13
535 25.690 0 12
536 25.010 0 11
537 33.200 0 11
538 21.230 0 11
539 37.800 0 14
540 38.470 0 11
541 37.850 1 13
542 39.240 0 12
543 21.320 0 10
544 27.520 0 12
545 24.490 0 13
546 34.560 0 10
547 30.730 0 10
548 26.880 0 10
549 33.290 0 10
550 42.710 0 14
551 35.300 0 12
552 29.280 0 11
553 26.890 0 11
554 23.320 0 12
555 29.340 0 14
556 22.760 1 14
557 31.100 0 10
558 28.940 0 11
559 46.370 1 11
560 24.350 0 10
561 28.380 0 10
562 30.350 0 12
563 48.310 0 12
564 28.120 0 11
565 27.140 0 12
566 30.860 0 10
567 35.190 0 12
568 42.320 0 13
569 27.700 0 10
570 33.410 0 12
571 30.900 0 10
572 25.310 0 13
573 28.220 0 12
574 30.380 1 10
575 29.350 0 12
576 25.680 0 10
577 23.870 0 11
578 24.990 0 12
579 41.300 0 11
580 30.010 0 12
581 31.320 0 10
582 35.770 0 13
583 32.220 0 12
584 32.800 0 11
585 26.590 0 11
586 28.220 0 11
587 21.400 0 11
588 42.030 0 12
589 29.970 0 14
590 31.200 1 11
591 25.620 0 11
592 30.820 0 12
593 38.060 0 14
594 33.390 1 11
595 31.520 1 13
596 24.580 0 11
597 23.910 0 10
598 31.410 0 13
599 25.790 0 12
600 31.040 1 11
601 40.450 1 13
602 47.630 1 14
603 21.000 0 10
604 30.690 1 11
605 27.850 0 11
606 42.840 0 15
607 45.060 1 15
608 29.060 0 18
609 51.020 0 19
610 35.190 1 19
611 36.880 1 16
612 44.290 0 17
613 42.790 0 15
614 45.000 0 15
615 26.350 0 15
616 26.790 1 15
617 21.980 1 15
618 33.450 1 19
619 30.820 0 18
620 33.870 0 16
621 30.820 1 17
622 29.030 1 16
623 30.040 1 15
624 57.930 0 15
625 39.850 0 15
626 42.200 0 18
627 47.240 0 17
628 37.310 0 15
629 34.060 1 17
630 35.000 0 17
631 36.740 0 16
632 56.330 0 17
633 31.220 1 15
634 33.300 1 15
635 26.080 1 16
636 36.450 0 16
637 37.990 1 15
638 40.860 1 18
639 28.870 0 15
640 40.700 1 16
641 39.600 0 17
642 42.990 0 16
643 29.810 0 16
644 22.640 1 15
645 44.040 1 18
646 22.780 1 15
647 45.040 0 16
648 56.380 0 17
649 48.720 1 16
650 42.700 1 16
651 37.270 1 15
652 28.530 0 18
653 27.950 1 16
654 32.110 0 15
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Smoke Age
3.409 -2.148 2.330
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-16.8175 -3.4785 -0.3986 3.5367 20.8720
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 3.40917 0.82212 4.147 3.82e-05 ***
Smoke -2.14803 0.81515 -2.635 0.00861 **
Age 2.32990 0.08262 28.199 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 5.704 on 651 degrees of freedom
Multiple R-squared: 0.5767, Adjusted R-squared: 0.5754
F-statistic: 443.5 on 2 and 651 DF, p-value: < 2.2e-16
> 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
+ }
[,1] [,2] [,3]
[1,] 2.240053e-01 4.480106e-01 0.7759947
[2,] 1.067402e-01 2.134805e-01 0.8932598
[3,] 1.003764e-01 2.007529e-01 0.8996236
[4,] 5.610237e-02 1.122047e-01 0.9438976
[5,] 2.710777e-02 5.421554e-02 0.9728922
[6,] 1.298102e-02 2.596203e-02 0.9870190
[7,] 5.741566e-03 1.148313e-02 0.9942584
[8,] 5.568981e-03 1.113796e-02 0.9944310
[9,] 3.679815e-03 7.359630e-03 0.9963202
[10,] 3.652987e-03 7.305974e-03 0.9963470
[11,] 1.763253e-03 3.526506e-03 0.9982367
[12,] 9.527218e-04 1.905444e-03 0.9990473
[13,] 4.772629e-04 9.545258e-04 0.9995227
[14,] 4.363683e-04 8.727367e-04 0.9995636
[15,] 9.461192e-04 1.892238e-03 0.9990539
[16,] 5.170542e-04 1.034108e-03 0.9994829
[17,] 3.441240e-04 6.882480e-04 0.9996559
[18,] 1.079865e-02 2.159730e-02 0.9892014
[19,] 2.779234e-02 5.558468e-02 0.9722077
[20,] 4.530268e-02 9.060536e-02 0.9546973
[21,] 4.657194e-02 9.314387e-02 0.9534281
[22,] 3.312640e-02 6.625280e-02 0.9668736
[23,] 2.763733e-02 5.527467e-02 0.9723627
[24,] 6.109014e-02 1.221803e-01 0.9389099
[25,] 4.685239e-02 9.370479e-02 0.9531476
[26,] 3.485327e-02 6.970655e-02 0.9651467
[27,] 4.600920e-02 9.201839e-02 0.9539908
[28,] 4.896015e-02 9.792031e-02 0.9510398
[29,] 3.788722e-02 7.577444e-02 0.9621128
[30,] 2.872116e-02 5.744232e-02 0.9712788
[31,] 2.076280e-02 4.152560e-02 0.9792372
[32,] 1.755010e-02 3.510021e-02 0.9824499
[33,] 1.254061e-02 2.508122e-02 0.9874594
[34,] 1.447556e-02 2.895111e-02 0.9855244
[35,] 1.048599e-02 2.097197e-02 0.9895140
[36,] 1.644110e-02 3.288220e-02 0.9835589
[37,] 1.301540e-02 2.603079e-02 0.9869846
[38,] 1.085886e-02 2.171772e-02 0.9891411
[39,] 8.314486e-03 1.662897e-02 0.9916855
[40,] 5.968345e-03 1.193669e-02 0.9940317
[41,] 4.812733e-03 9.625465e-03 0.9951873
[42,] 3.792656e-03 7.585311e-03 0.9962073
[43,] 6.500997e-03 1.300199e-02 0.9934990
[44,] 4.624180e-03 9.248360e-03 0.9953758
[45,] 5.457054e-03 1.091411e-02 0.9945429
[46,] 5.984355e-03 1.196871e-02 0.9940156
[47,] 5.475276e-03 1.095055e-02 0.9945247
[48,] 1.917817e-02 3.835635e-02 0.9808218
[49,] 1.839950e-02 3.679900e-02 0.9816005
[50,] 1.401304e-02 2.802609e-02 0.9859870
[51,] 1.082465e-02 2.164929e-02 0.9891754
[52,] 1.187321e-02 2.374642e-02 0.9881268
[53,] 1.140385e-02 2.280770e-02 0.9885962
[54,] 9.787986e-03 1.957597e-02 0.9902120
[55,] 9.606214e-03 1.921243e-02 0.9903938
[56,] 7.244244e-03 1.448849e-02 0.9927558
[57,] 5.948148e-03 1.189630e-02 0.9940519
[58,] 4.692311e-03 9.384623e-03 0.9953077
[59,] 3.961586e-03 7.923172e-03 0.9960384
[60,] 3.136204e-03 6.272407e-03 0.9968638
[61,] 2.291847e-03 4.583694e-03 0.9977082
[62,] 1.935734e-03 3.871467e-03 0.9980643
[63,] 1.592899e-03 3.185797e-03 0.9984071
[64,] 1.150221e-03 2.300442e-03 0.9988498
[65,] 8.481607e-04 1.696321e-03 0.9991518
[66,] 6.421303e-04 1.284261e-03 0.9993579
[67,] 5.412869e-04 1.082574e-03 0.9994587
[68,] 4.332715e-04 8.665431e-04 0.9995667
[69,] 3.043456e-04 6.086913e-04 0.9996957
[70,] 2.488936e-04 4.977871e-04 0.9997511
[71,] 1.795929e-04 3.591857e-04 0.9998204
[72,] 1.527798e-04 3.055597e-04 0.9998472
[73,] 1.055340e-04 2.110680e-04 0.9998945
[74,] 1.188926e-03 2.377851e-03 0.9988111
[75,] 8.697350e-04 1.739470e-03 0.9991303
[76,] 6.488329e-04 1.297666e-03 0.9993512
[77,] 4.692708e-04 9.385416e-04 0.9995307
[78,] 3.366838e-04 6.733676e-04 0.9996633
[79,] 3.015989e-04 6.031979e-04 0.9996984
[80,] 2.154303e-04 4.308606e-04 0.9997846
[81,] 1.531417e-04 3.062835e-04 0.9998469
[82,] 2.160937e-04 4.321873e-04 0.9997839
[83,] 1.645868e-04 3.291737e-04 0.9998354
[84,] 1.181857e-04 2.363713e-04 0.9998818
[85,] 9.517793e-05 1.903559e-04 0.9999048
[86,] 1.100988e-04 2.201976e-04 0.9998899
[87,] 1.127240e-04 2.254480e-04 0.9998873
[88,] 9.470139e-05 1.894028e-04 0.9999053
[89,] 6.677775e-05 1.335555e-04 0.9999332
[90,] 4.908316e-05 9.816632e-05 0.9999509
[91,] 3.768945e-05 7.537890e-05 0.9999623
[92,] 5.258641e-05 1.051728e-04 0.9999474
[93,] 4.088610e-05 8.177220e-05 0.9999591
[94,] 3.179547e-05 6.359094e-05 0.9999682
[95,] 4.192467e-05 8.384935e-05 0.9999581
[96,] 8.721229e-05 1.744246e-04 0.9999128
[97,] 6.967241e-05 1.393448e-04 0.9999303
[98,] 5.179989e-05 1.035998e-04 0.9999482
[99,] 1.878889e-04 3.757778e-04 0.9998121
[100,] 1.408561e-04 2.817121e-04 0.9998591
[101,] 1.148581e-04 2.297162e-04 0.9998851
[102,] 8.648950e-05 1.729790e-04 0.9999135
[103,] 6.248800e-05 1.249760e-04 0.9999375
[104,] 4.692151e-05 9.384302e-05 0.9999531
[105,] 3.855043e-05 7.710085e-05 0.9999614
[106,] 3.652845e-05 7.305691e-05 0.9999635
[107,] 2.637748e-05 5.275496e-05 0.9999736
[108,] 1.933321e-05 3.866642e-05 0.9999807
[109,] 1.533797e-05 3.067594e-05 0.9999847
[110,] 1.447295e-05 2.894590e-05 0.9999855
[111,] 1.056359e-05 2.112717e-05 0.9999894
[112,] 7.991947e-06 1.598389e-05 0.9999920
[113,] 6.019359e-06 1.203872e-05 0.9999940
[114,] 4.804952e-06 9.609904e-06 0.9999952
[115,] 3.791250e-06 7.582501e-06 0.9999962
[116,] 2.623161e-06 5.246322e-06 0.9999974
[117,] 1.834441e-06 3.668882e-06 0.9999982
[118,] 1.749247e-06 3.498494e-06 0.9999983
[119,] 1.420512e-06 2.841023e-06 0.9999986
[120,] 9.831981e-07 1.966396e-06 0.9999990
[121,] 9.539544e-07 1.907909e-06 0.9999990
[122,] 1.149693e-06 2.299386e-06 0.9999989
[123,] 9.773907e-07 1.954781e-06 0.9999990
[124,] 7.713555e-07 1.542711e-06 0.9999992
[125,] 6.046051e-07 1.209210e-06 0.9999994
[126,] 4.491061e-07 8.982122e-07 0.9999996
[127,] 7.066295e-07 1.413259e-06 0.9999993
[128,] 6.925867e-07 1.385173e-06 0.9999993
[129,] 6.522770e-07 1.304554e-06 0.9999993
[130,] 4.523252e-07 9.046505e-07 0.9999995
[131,] 3.132089e-07 6.264178e-07 0.9999997
[132,] 2.479710e-07 4.959420e-07 0.9999998
[133,] 1.734467e-07 3.468934e-07 0.9999998
[134,] 2.179012e-07 4.358024e-07 0.9999998
[135,] 1.963036e-06 3.926072e-06 0.9999980
[136,] 1.384310e-06 2.768620e-06 0.9999986
[137,] 9.742127e-07 1.948425e-06 0.9999990
[138,] 7.255460e-07 1.451092e-06 0.9999993
[139,] 5.436875e-07 1.087375e-06 0.9999995
[140,] 7.746724e-07 1.549345e-06 0.9999992
[141,] 6.034055e-07 1.206811e-06 0.9999994
[142,] 4.290713e-07 8.581427e-07 0.9999996
[143,] 3.061649e-07 6.123297e-07 0.9999997
[144,] 3.761883e-07 7.523767e-07 0.9999996
[145,] 6.076400e-07 1.215280e-06 0.9999994
[146,] 6.043698e-07 1.208740e-06 0.9999994
[147,] 4.322349e-07 8.644698e-07 0.9999996
[148,] 3.010577e-07 6.021154e-07 0.9999997
[149,] 2.165783e-07 4.331565e-07 0.9999998
[150,] 1.553353e-07 3.106706e-07 0.9999998
[151,] 1.370728e-07 2.741456e-07 0.9999999
[152,] 9.635540e-08 1.927108e-07 0.9999999
[153,] 6.892953e-08 1.378591e-07 0.9999999
[154,] 5.852539e-08 1.170508e-07 0.9999999
[155,] 4.124679e-08 8.249358e-08 1.0000000
[156,] 2.906511e-08 5.813022e-08 1.0000000
[157,] 2.481317e-08 4.962633e-08 1.0000000
[158,] 1.778964e-08 3.557929e-08 1.0000000
[159,] 1.287936e-08 2.575872e-08 1.0000000
[160,] 1.232294e-08 2.464587e-08 1.0000000
[161,] 8.415109e-09 1.683022e-08 1.0000000
[162,] 7.525075e-09 1.505015e-08 1.0000000
[163,] 5.079516e-09 1.015903e-08 1.0000000
[164,] 4.809004e-09 9.618009e-09 1.0000000
[165,] 4.413236e-09 8.826472e-09 1.0000000
[166,] 3.282183e-09 6.564366e-09 1.0000000
[167,] 2.425086e-09 4.850173e-09 1.0000000
[168,] 3.252044e-09 6.504088e-09 1.0000000
[169,] 3.700097e-09 7.400194e-09 1.0000000
[170,] 2.939577e-09 5.879154e-09 1.0000000
[171,] 2.940306e-09 5.880612e-09 1.0000000
[172,] 2.363015e-09 4.726030e-09 1.0000000
[173,] 1.607532e-09 3.215065e-09 1.0000000
[174,] 1.237423e-09 2.474847e-09 1.0000000
[175,] 2.422577e-09 4.845153e-09 1.0000000
[176,] 3.999106e-09 7.998211e-09 1.0000000
[177,] 3.634005e-09 7.268009e-09 1.0000000
[178,] 2.483433e-09 4.966866e-09 1.0000000
[179,] 1.759947e-09 3.519894e-09 1.0000000
[180,] 2.533195e-09 5.066391e-09 1.0000000
[181,] 1.754863e-09 3.509726e-09 1.0000000
[182,] 1.345198e-09 2.690396e-09 1.0000000
[183,] 9.026601e-10 1.805320e-09 1.0000000
[184,] 7.008059e-10 1.401612e-09 1.0000000
[185,] 4.796587e-10 9.593174e-10 1.0000000
[186,] 3.271411e-10 6.542822e-10 1.0000000
[187,] 6.667960e-10 1.333592e-09 1.0000000
[188,] 5.089296e-10 1.017859e-09 1.0000000
[189,] 3.414016e-10 6.828031e-10 1.0000000
[190,] 2.373705e-10 4.747410e-10 1.0000000
[191,] 2.573343e-10 5.146687e-10 1.0000000
[192,] 2.411486e-10 4.822971e-10 1.0000000
[193,] 1.940259e-10 3.880518e-10 1.0000000
[194,] 2.281388e-10 4.562776e-10 1.0000000
[195,] 2.047769e-10 4.095538e-10 1.0000000
[196,] 1.500946e-10 3.001893e-10 1.0000000
[197,] 1.091410e-10 2.182819e-10 1.0000000
[198,] 1.203590e-10 2.407180e-10 1.0000000
[199,] 7.976895e-11 1.595379e-10 1.0000000
[200,] 8.487321e-11 1.697464e-10 1.0000000
[201,] 5.878695e-11 1.175739e-10 1.0000000
[202,] 3.890484e-11 7.780968e-11 1.0000000
[203,] 2.555698e-11 5.111395e-11 1.0000000
[204,] 1.675842e-11 3.351685e-11 1.0000000
[205,] 1.216418e-11 2.432836e-11 1.0000000
[206,] 9.082824e-12 1.816565e-11 1.0000000
[207,] 1.862691e-11 3.725382e-11 1.0000000
[208,] 1.241987e-11 2.483974e-11 1.0000000
[209,] 8.733826e-12 1.746765e-11 1.0000000
[210,] 5.675321e-12 1.135064e-11 1.0000000
[211,] 3.795122e-12 7.590244e-12 1.0000000
[212,] 3.081162e-12 6.162324e-12 1.0000000
[213,] 2.508443e-12 5.016886e-12 1.0000000
[214,] 1.666106e-12 3.332213e-12 1.0000000
[215,] 1.521589e-12 3.043178e-12 1.0000000
[216,] 1.918526e-12 3.837051e-12 1.0000000
[217,] 1.275767e-12 2.551534e-12 1.0000000
[218,] 1.040851e-12 2.081703e-12 1.0000000
[219,] 1.441461e-12 2.882921e-12 1.0000000
[220,] 1.070114e-12 2.140227e-12 1.0000000
[221,] 1.075027e-11 2.150054e-11 1.0000000
[222,] 7.492032e-12 1.498406e-11 1.0000000
[223,] 5.602160e-12 1.120432e-11 1.0000000
[224,] 6.292004e-12 1.258401e-11 1.0000000
[225,] 5.162965e-12 1.032593e-11 1.0000000
[226,] 3.537620e-12 7.075240e-12 1.0000000
[227,] 2.511296e-12 5.022591e-12 1.0000000
[228,] 1.753208e-12 3.506415e-12 1.0000000
[229,] 1.375260e-12 2.750519e-12 1.0000000
[230,] 8.976779e-13 1.795356e-12 1.0000000
[231,] 8.994439e-13 1.798888e-12 1.0000000
[232,] 1.734311e-12 3.468623e-12 1.0000000
[233,] 1.158580e-12 2.317160e-12 1.0000000
[234,] 7.935020e-13 1.587004e-12 1.0000000
[235,] 5.492219e-13 1.098444e-12 1.0000000
[236,] 4.959271e-13 9.918541e-13 1.0000000
[237,] 3.222633e-13 6.445267e-13 1.0000000
[238,] 2.165537e-13 4.331074e-13 1.0000000
[239,] 2.282531e-13 4.565063e-13 1.0000000
[240,] 1.811969e-13 3.623939e-13 1.0000000
[241,] 1.304385e-13 2.608771e-13 1.0000000
[242,] 9.156733e-14 1.831347e-13 1.0000000
[243,] 8.045689e-14 1.609138e-13 1.0000000
[244,] 6.349158e-14 1.269832e-13 1.0000000
[245,] 8.971571e-14 1.794314e-13 1.0000000
[246,] 7.377194e-14 1.475439e-13 1.0000000
[247,] 4.862556e-14 9.725113e-14 1.0000000
[248,] 3.286259e-14 6.572518e-14 1.0000000
[249,] 2.444541e-14 4.889082e-14 1.0000000
[250,] 4.989394e-14 9.978789e-14 1.0000000
[251,] 3.192042e-14 6.384083e-14 1.0000000
[252,] 2.180279e-14 4.360559e-14 1.0000000
[253,] 1.663080e-14 3.326161e-14 1.0000000
[254,] 1.114214e-14 2.228428e-14 1.0000000
[255,] 7.275717e-15 1.455143e-14 1.0000000
[256,] 4.625471e-15 9.250942e-15 1.0000000
[257,] 5.020982e-15 1.004196e-14 1.0000000
[258,] 4.240752e-15 8.481505e-15 1.0000000
[259,] 2.725819e-15 5.451639e-15 1.0000000
[260,] 2.439674e-15 4.879349e-15 1.0000000
[261,] 2.949467e-15 5.898934e-15 1.0000000
[262,] 3.052817e-15 6.105633e-15 1.0000000
[263,] 1.968140e-15 3.936280e-15 1.0000000
[264,] 1.252195e-15 2.504390e-15 1.0000000
[265,] 8.131297e-16 1.626259e-15 1.0000000
[266,] 6.513661e-16 1.302732e-15 1.0000000
[267,] 6.268912e-16 1.253782e-15 1.0000000
[268,] 4.112578e-16 8.225155e-16 1.0000000
[269,] 2.627753e-16 5.255505e-16 1.0000000
[270,] 3.098079e-16 6.196159e-16 1.0000000
[271,] 2.109186e-16 4.218371e-16 1.0000000
[272,] 1.402076e-16 2.804153e-16 1.0000000
[273,] 1.290651e-16 2.581302e-16 1.0000000
[274,] 9.052780e-17 1.810556e-16 1.0000000
[275,] 6.000216e-17 1.200043e-16 1.0000000
[276,] 3.761910e-17 7.523819e-17 1.0000000
[277,] 2.378131e-17 4.756263e-17 1.0000000
[278,] 1.735404e-17 3.470808e-17 1.0000000
[279,] 1.457673e-17 2.915346e-17 1.0000000
[280,] 9.318117e-18 1.863623e-17 1.0000000
[281,] 6.244802e-18 1.248960e-17 1.0000000
[282,] 3.945971e-18 7.891941e-18 1.0000000
[283,] 2.491737e-18 4.983473e-18 1.0000000
[284,] 2.288355e-18 4.576711e-18 1.0000000
[285,] 2.602343e-18 5.204685e-18 1.0000000
[286,] 3.136442e-18 6.272884e-18 1.0000000
[287,] 1.911871e-18 3.823743e-18 1.0000000
[288,] 4.824789e-18 9.649579e-18 1.0000000
[289,] 5.823705e-18 1.164741e-17 1.0000000
[290,] 4.539162e-18 9.078325e-18 1.0000000
[291,] 2.882887e-18 5.765773e-18 1.0000000
[292,] 1.947285e-18 3.894569e-18 1.0000000
[293,] 2.277023e-18 4.554047e-18 1.0000000
[294,] 1.781459e-18 3.562919e-18 1.0000000
[295,] 1.418913e-18 2.837825e-18 1.0000000
[296,] 9.427096e-19 1.885419e-18 1.0000000
[297,] 6.417888e-19 1.283578e-18 1.0000000
[298,] 7.527626e-19 1.505525e-18 1.0000000
[299,] 4.929595e-19 9.859190e-19 1.0000000
[300,] 3.346547e-19 6.693094e-19 1.0000000
[301,] 2.753113e-19 5.506226e-19 1.0000000
[302,] 3.360050e-19 6.720099e-19 1.0000000
[303,] 2.114676e-19 4.229352e-19 1.0000000
[304,] 2.252829e-19 4.505658e-19 1.0000000
[305,] 1.350046e-19 2.700091e-19 1.0000000
[306,] 1.046653e-19 2.093305e-19 1.0000000
[307,] 6.657550e-20 1.331510e-19 1.0000000
[308,] 1.141087e-19 2.282175e-19 1.0000000
[309,] 1.321068e-19 2.642137e-19 1.0000000
[310,] 7.925220e-20 1.585044e-19 1.0000000
[311,] 2.108444e-19 4.216888e-19 1.0000000
[312,] 1.654201e-19 3.308402e-19 1.0000000
[313,] 1.007548e-19 2.015096e-19 1.0000000
[314,] 6.459766e-20 1.291953e-19 1.0000000
[315,] 2.667491e-19 5.334982e-19 1.0000000
[316,] 3.013311e-18 6.026621e-18 1.0000000
[317,] 5.586698e-17 1.117340e-16 1.0000000
[318,] 4.264728e-17 8.529456e-17 1.0000000
[319,] 3.772526e-17 7.545052e-17 1.0000000
[320,] 3.079071e-17 6.158143e-17 1.0000000
[321,] 2.015829e-17 4.031659e-17 1.0000000
[322,] 2.895620e-17 5.791240e-17 1.0000000
[323,] 2.016457e-17 4.032914e-17 1.0000000
[324,] 1.459479e-17 2.918958e-17 1.0000000
[325,] 1.344102e-17 2.688205e-17 1.0000000
[326,] 8.503756e-18 1.700751e-17 1.0000000
[327,] 1.921525e-17 3.843050e-17 1.0000000
[328,] 1.335480e-17 2.670960e-17 1.0000000
[329,] 1.067512e-17 2.135023e-17 1.0000000
[330,] 1.062332e-17 2.124664e-17 1.0000000
[331,] 7.629286e-18 1.525857e-17 1.0000000
[332,] 5.116640e-18 1.023328e-17 1.0000000
[333,] 4.475244e-18 8.950488e-18 1.0000000
[334,] 3.237084e-18 6.474169e-18 1.0000000
[335,] 3.042167e-18 6.084334e-18 1.0000000
[336,] 5.767240e-18 1.153448e-17 1.0000000
[337,] 3.618650e-18 7.237300e-18 1.0000000
[338,] 2.216833e-18 4.433666e-18 1.0000000
[339,] 9.446372e-18 1.889274e-17 1.0000000
[340,] 1.043691e-17 2.087381e-17 1.0000000
[341,] 1.034735e-17 2.069470e-17 1.0000000
[342,] 6.854741e-18 1.370948e-17 1.0000000
[343,] 4.233547e-18 8.467094e-18 1.0000000
[344,] 9.974402e-18 1.994880e-17 1.0000000
[345,] 9.040189e-18 1.808038e-17 1.0000000
[346,] 1.047916e-17 2.095831e-17 1.0000000
[347,] 7.224567e-18 1.444913e-17 1.0000000
[348,] 5.043269e-18 1.008654e-17 1.0000000
[349,] 3.128196e-18 6.256391e-18 1.0000000
[350,] 2.169585e-18 4.339170e-18 1.0000000
[351,] 3.484332e-18 6.968665e-18 1.0000000
[352,] 2.194485e-18 4.388970e-18 1.0000000
[353,] 1.921281e-18 3.842562e-18 1.0000000
[354,] 2.375049e-18 4.750098e-18 1.0000000
[355,] 2.246567e-18 4.493135e-18 1.0000000
[356,] 7.929885e-18 1.585977e-17 1.0000000
[357,] 6.480048e-18 1.296010e-17 1.0000000
[358,] 4.084536e-18 8.169072e-18 1.0000000
[359,] 4.166509e-18 8.333017e-18 1.0000000
[360,] 2.539135e-17 5.078269e-17 1.0000000
[361,] 1.922697e-17 3.845394e-17 1.0000000
[362,] 1.256639e-17 2.513278e-17 1.0000000
[363,] 8.942245e-18 1.788449e-17 1.0000000
[364,] 1.753282e-17 3.506563e-17 1.0000000
[365,] 3.051472e-17 6.102945e-17 1.0000000
[366,] 3.823691e-17 7.647382e-17 1.0000000
[367,] 2.316111e-15 4.632222e-15 1.0000000
[368,] 2.221447e-15 4.442893e-15 1.0000000
[369,] 1.557565e-15 3.115131e-15 1.0000000
[370,] 1.002609e-15 2.005217e-15 1.0000000
[371,] 9.912987e-16 1.982597e-15 1.0000000
[372,] 7.733021e-16 1.546604e-15 1.0000000
[373,] 5.583079e-16 1.116616e-15 1.0000000
[374,] 3.608494e-16 7.216988e-16 1.0000000
[375,] 2.627916e-16 5.255832e-16 1.0000000
[376,] 2.649934e-16 5.299868e-16 1.0000000
[377,] 3.428508e-16 6.857016e-16 1.0000000
[378,] 2.842435e-16 5.684870e-16 1.0000000
[379,] 2.032187e-16 4.064373e-16 1.0000000
[380,] 5.609866e-16 1.121973e-15 1.0000000
[381,] 6.000663e-16 1.200133e-15 1.0000000
[382,] 4.901115e-16 9.802229e-16 1.0000000
[383,] 3.405888e-16 6.811776e-16 1.0000000
[384,] 2.243140e-16 4.486280e-16 1.0000000
[385,] 1.894119e-16 3.788237e-16 1.0000000
[386,] 3.427800e-16 6.855600e-16 1.0000000
[387,] 2.428881e-16 4.857762e-16 1.0000000
[388,] 2.592016e-16 5.184033e-16 1.0000000
[389,] 4.827301e-16 9.654603e-16 1.0000000
[390,] 3.412420e-16 6.824840e-16 1.0000000
[391,] 2.207030e-16 4.414059e-16 1.0000000
[392,] 5.558724e-16 1.111745e-15 1.0000000
[393,] 5.913022e-16 1.182604e-15 1.0000000
[394,] 6.491265e-16 1.298253e-15 1.0000000
[395,] 4.579279e-16 9.158558e-16 1.0000000
[396,] 5.127331e-16 1.025466e-15 1.0000000
[397,] 4.659205e-16 9.318410e-16 1.0000000
[398,] 1.004719e-15 2.009437e-15 1.0000000
[399,] 3.374884e-15 6.749767e-15 1.0000000
[400,] 4.260957e-15 8.521914e-15 1.0000000
[401,] 4.552499e-15 9.104997e-15 1.0000000
[402,] 3.402120e-15 6.804241e-15 1.0000000
[403,] 3.016904e-15 6.033808e-15 1.0000000
[404,] 1.969706e-15 3.939413e-15 1.0000000
[405,] 2.394960e-15 4.789921e-15 1.0000000
[406,] 2.064183e-15 4.128367e-15 1.0000000
[407,] 1.380510e-15 2.761020e-15 1.0000000
[408,] 1.605349e-15 3.210697e-15 1.0000000
[409,] 1.660324e-15 3.320648e-15 1.0000000
[410,] 1.297971e-14 2.595942e-14 1.0000000
[411,] 3.546742e-14 7.093484e-14 1.0000000
[412,] 3.286566e-14 6.573132e-14 1.0000000
[413,] 2.602906e-14 5.205812e-14 1.0000000
[414,] 1.839908e-14 3.679815e-14 1.0000000
[415,] 1.190362e-14 2.380724e-14 1.0000000
[416,] 1.137849e-14 2.275698e-14 1.0000000
[417,] 2.018682e-13 4.037364e-13 1.0000000
[418,] 2.935215e-13 5.870430e-13 1.0000000
[419,] 1.962649e-13 3.925299e-13 1.0000000
[420,] 1.289634e-13 2.579268e-13 1.0000000
[421,] 1.180865e-13 2.361730e-13 1.0000000
[422,] 1.070410e-13 2.140819e-13 1.0000000
[423,] 7.113451e-14 1.422690e-13 1.0000000
[424,] 9.354320e-14 1.870864e-13 1.0000000
[425,] 8.762785e-14 1.752557e-13 1.0000000
[426,] 7.316996e-14 1.463399e-13 1.0000000
[427,] 6.670905e-14 1.334181e-13 1.0000000
[428,] 1.595290e-13 3.190581e-13 1.0000000
[429,] 1.328925e-13 2.657849e-13 1.0000000
[430,] 8.885544e-14 1.777109e-13 1.0000000
[431,] 5.915276e-14 1.183055e-13 1.0000000
[432,] 4.480245e-14 8.960489e-14 1.0000000
[433,] 2.940201e-14 5.880402e-14 1.0000000
[434,] 5.324699e-14 1.064940e-13 1.0000000
[435,] 3.458305e-14 6.916610e-14 1.0000000
[436,] 5.306194e-14 1.061239e-13 1.0000000
[437,] 9.499042e-13 1.899808e-12 1.0000000
[438,] 8.575067e-13 1.715013e-12 1.0000000
[439,] 4.509175e-12 9.018350e-12 1.0000000
[440,] 8.752934e-12 1.750587e-11 1.0000000
[441,] 6.793054e-12 1.358611e-11 1.0000000
[442,] 6.738078e-12 1.347616e-11 1.0000000
[443,] 4.633098e-12 9.266197e-12 1.0000000
[444,] 9.566243e-12 1.913249e-11 1.0000000
[445,] 1.419576e-11 2.839151e-11 1.0000000
[446,] 1.017123e-11 2.034246e-11 1.0000000
[447,] 7.200108e-10 1.440022e-09 1.0000000
[448,] 1.737336e-09 3.474671e-09 1.0000000
[449,] 2.740246e-09 5.480492e-09 1.0000000
[450,] 2.266682e-09 4.533365e-09 1.0000000
[451,] 1.781930e-09 3.563860e-09 1.0000000
[452,] 5.747965e-09 1.149593e-08 1.0000000
[453,] 5.822174e-09 1.164435e-08 1.0000000
[454,] 4.635110e-09 9.270219e-09 1.0000000
[455,] 4.378252e-09 8.756504e-09 1.0000000
[456,] 3.391545e-09 6.783091e-09 1.0000000
[457,] 2.866033e-09 5.732065e-09 1.0000000
[458,] 2.200800e-09 4.401599e-09 1.0000000
[459,] 1.443774e-08 2.887549e-08 1.0000000
[460,] 1.242333e-08 2.484667e-08 1.0000000
[461,] 8.908097e-09 1.781619e-08 1.0000000
[462,] 6.618694e-09 1.323739e-08 1.0000000
[463,] 9.907631e-09 1.981526e-08 1.0000000
[464,] 7.124038e-09 1.424808e-08 1.0000000
[465,] 5.183360e-09 1.036672e-08 1.0000000
[466,] 8.847582e-09 1.769516e-08 1.0000000
[467,] 1.032958e-08 2.065916e-08 1.0000000
[468,] 2.370867e-08 4.741734e-08 1.0000000
[469,] 1.703253e-08 3.406507e-08 1.0000000
[470,] 1.296402e-08 2.592803e-08 1.0000000
[471,] 9.330875e-09 1.866175e-08 1.0000000
[472,] 9.576963e-09 1.915393e-08 1.0000000
[473,] 1.701060e-08 3.402119e-08 1.0000000
[474,] 3.541853e-08 7.083705e-08 1.0000000
[475,] 2.623987e-08 5.247973e-08 1.0000000
[476,] 2.032753e-08 4.065507e-08 1.0000000
[477,] 1.742556e-08 3.485111e-08 1.0000000
[478,] 1.279566e-08 2.559132e-08 1.0000000
[479,] 9.233798e-09 1.846760e-08 1.0000000
[480,] 2.525780e-08 5.051561e-08 1.0000000
[481,] 4.323203e-08 8.646407e-08 1.0000000
[482,] 3.097237e-08 6.194475e-08 1.0000000
[483,] 2.227355e-08 4.454710e-08 1.0000000
[484,] 3.272455e-08 6.544910e-08 1.0000000
[485,] 6.771350e-08 1.354270e-07 0.9999999
[486,] 6.746752e-08 1.349350e-07 0.9999999
[487,] 6.754793e-08 1.350959e-07 0.9999999
[488,] 5.948113e-08 1.189623e-07 0.9999999
[489,] 7.261469e-08 1.452294e-07 0.9999999
[490,] 5.307143e-08 1.061429e-07 0.9999999
[491,] 4.004364e-08 8.008728e-08 1.0000000
[492,] 2.867299e-08 5.734598e-08 1.0000000
[493,] 2.779011e-08 5.558023e-08 1.0000000
[494,] 1.961644e-08 3.923288e-08 1.0000000
[495,] 2.127007e-08 4.254013e-08 1.0000000
[496,] 1.534137e-08 3.068273e-08 1.0000000
[497,] 1.356642e-08 2.713283e-08 1.0000000
[498,] 1.212256e-08 2.424512e-08 1.0000000
[499,] 1.811129e-07 3.622258e-07 0.9999998
[500,] 1.316127e-07 2.632255e-07 0.9999999
[501,] 9.463028e-08 1.892606e-07 0.9999999
[502,] 7.118843e-08 1.423769e-07 0.9999999
[503,] 7.326089e-08 1.465218e-07 0.9999999
[504,] 5.453345e-08 1.090669e-07 0.9999999
[505,] 2.024274e-07 4.048547e-07 0.9999998
[506,] 2.163318e-07 4.326636e-07 0.9999998
[507,] 1.881231e-07 3.762463e-07 0.9999998
[508,] 1.386536e-07 2.773071e-07 0.9999999
[509,] 1.002790e-07 2.005579e-07 0.9999999
[510,] 6.343515e-07 1.268703e-06 0.9999994
[511,] 4.620787e-07 9.241575e-07 0.9999995
[512,] 4.578784e-06 9.157569e-06 0.9999954
[513,] 6.238220e-06 1.247644e-05 0.9999938
[514,] 7.623590e-06 1.524718e-05 0.9999924
[515,] 6.740315e-06 1.348063e-05 0.9999933
[516,] 6.099415e-06 1.219883e-05 0.9999939
[517,] 4.804586e-06 9.609171e-06 0.9999952
[518,] 3.565127e-06 7.130254e-06 0.9999964
[519,] 2.614083e-06 5.228166e-06 0.9999974
[520,] 1.911298e-06 3.822597e-06 0.9999981
[521,] 1.558160e-06 3.116320e-06 0.9999984
[522,] 1.251844e-06 2.503688e-06 0.9999987
[523,] 9.929755e-07 1.985951e-06 0.9999990
[524,] 1.095647e-06 2.191293e-06 0.9999989
[525,] 7.859857e-07 1.571971e-06 0.9999992
[526,] 6.148098e-07 1.229620e-06 0.9999994
[527,] 4.450772e-07 8.901545e-07 0.9999996
[528,] 4.372931e-07 8.745862e-07 0.9999996
[529,] 3.434385e-07 6.868770e-07 0.9999997
[530,] 3.531257e-07 7.062514e-07 0.9999996
[531,] 3.220551e-07 6.441103e-07 0.9999997
[532,] 2.465349e-07 4.930697e-07 0.9999998
[533,] 3.790859e-07 7.581718e-07 0.9999996
[534,] 2.820160e-07 5.640320e-07 0.9999997
[535,] 3.791097e-07 7.582194e-07 0.9999996
[536,] 3.687689e-07 7.375378e-07 0.9999996
[537,] 4.302548e-07 8.605096e-07 0.9999996
[538,] 4.881928e-07 9.763855e-07 0.9999995
[539,] 4.162518e-07 8.325036e-07 0.9999996
[540,] 6.792687e-07 1.358537e-06 0.9999993
[541,] 6.799691e-07 1.359938e-06 0.9999993
[542,] 5.041105e-07 1.008221e-06 0.9999995
[543,] 3.625883e-07 7.251767e-07 0.9999996
[544,] 3.186158e-07 6.372316e-07 0.9999997
[545,] 3.686766e-07 7.373531e-07 0.9999996
[546,] 2.856536e-07 5.713071e-07 0.9999997
[547,] 1.981753e-07 3.963506e-07 0.9999998
[548,] 1.509833e-07 3.019667e-07 0.9999998
[549,] 2.196512e-07 4.393024e-07 0.9999998
[550,] 2.255875e-07 4.511751e-07 0.9999998
[551,] 5.221898e-07 1.044380e-06 0.9999995
[552,] 3.861646e-07 7.723292e-07 0.9999996
[553,] 2.703253e-07 5.406506e-07 0.9999997
[554,] 4.163077e-06 8.326154e-06 0.9999958
[555,] 3.432100e-06 6.864199e-06 0.9999966
[556,] 2.427430e-06 4.854861e-06 0.9999976
[557,] 1.743676e-06 3.487353e-06 0.9999983
[558,] 1.412541e-05 2.825082e-05 0.9999859
[559,] 1.043210e-05 2.086419e-05 0.9999896
[560,] 9.088961e-06 1.817792e-05 0.9999909
[561,] 6.855168e-06 1.371034e-05 0.9999931
[562,] 5.348457e-06 1.069691e-05 0.9999947
[563,] 7.618201e-06 1.523640e-05 0.9999924
[564,] 5.373104e-06 1.074621e-05 0.9999946
[565,] 3.861757e-06 7.723514e-06 0.9999961
[566,] 2.911943e-06 5.823886e-06 0.9999971
[567,] 3.815681e-06 7.631362e-06 0.9999962
[568,] 2.972394e-06 5.944788e-06 0.9999970
[569,] 2.463327e-06 4.926654e-06 0.9999975
[570,] 1.782893e-06 3.565785e-06 0.9999982
[571,] 1.291229e-06 2.582458e-06 0.9999987
[572,] 1.287965e-06 2.575931e-06 0.9999987
[573,] 1.435833e-06 2.871667e-06 0.9999986
[574,] 3.201148e-06 6.402295e-06 0.9999968
[575,] 2.254873e-06 4.509746e-06 0.9999977
[576,] 1.677506e-06 3.355012e-06 0.9999983
[577,] 1.184183e-06 2.368366e-06 0.9999988
[578,] 7.910574e-07 1.582115e-06 0.9999992
[579,] 5.720867e-07 1.144173e-06 0.9999994
[580,] 4.191123e-07 8.382246e-07 0.9999996
[581,] 2.821445e-07 5.642890e-07 0.9999997
[582,] 4.189034e-07 8.378067e-07 0.9999996
[583,] 7.719783e-07 1.543957e-06 0.9999992
[584,] 6.963970e-07 1.392794e-06 0.9999993
[585,] 5.154507e-07 1.030901e-06 0.9999995
[586,] 4.152325e-07 8.304649e-07 0.9999996
[587,] 2.713336e-07 5.426672e-07 0.9999997
[588,] 1.852116e-07 3.704232e-07 0.9999998
[589,] 1.700575e-07 3.401151e-07 0.9999998
[590,] 1.095945e-07 2.191889e-07 0.9999999
[591,] 9.695117e-08 1.939023e-07 0.9999999
[592,] 8.161652e-08 1.632330e-07 0.9999999
[593,] 5.587082e-08 1.117416e-07 0.9999999
[594,] 6.487308e-08 1.297462e-07 0.9999999
[595,] 4.226475e-08 8.452950e-08 1.0000000
[596,] 6.968886e-08 1.393777e-07 0.9999999
[597,] 6.746310e-07 1.349262e-06 0.9999993
[598,] 1.000682e-06 2.001365e-06 0.9999990
[599,] 7.060150e-07 1.412030e-06 0.9999993
[600,] 5.659589e-07 1.131918e-06 0.9999994
[601,] 4.287853e-07 8.575707e-07 0.9999996
[602,] 1.346728e-06 2.693456e-06 0.9999987
[603,] 8.130887e-06 1.626177e-05 0.9999919
[604,] 7.743994e-06 1.548799e-05 0.9999923
[605,] 7.062113e-06 1.412423e-05 0.9999929
[606,] 5.198588e-06 1.039718e-05 0.9999948
[607,] 3.575714e-06 7.151429e-06 0.9999964
[608,] 2.865227e-06 5.730454e-06 0.9999971
[609,] 3.293815e-06 6.587629e-06 0.9999967
[610,] 8.314640e-06 1.662928e-05 0.9999917
[611,] 7.947691e-06 1.589538e-05 0.9999921
[612,] 1.959485e-05 3.918969e-05 0.9999804
[613,] 1.940568e-05 3.881135e-05 0.9999806
[614,] 6.683553e-05 1.336711e-04 0.9999332
[615,] 6.347597e-05 1.269519e-04 0.9999365
[616,] 6.276405e-05 1.255281e-04 0.9999372
[617,] 6.026350e-05 1.205270e-04 0.9999397
[618,] 4.209552e-05 8.419103e-05 0.9999579
[619,] 2.051784e-03 4.103568e-03 0.9979482
[620,] 1.628384e-03 3.256769e-03 0.9983716
[621,] 1.165487e-03 2.330973e-03 0.9988345
[622,] 1.005421e-03 2.010842e-03 0.9989946
[623,] 6.855702e-04 1.371140e-03 0.9993144
[624,] 5.673310e-04 1.134662e-03 0.9994327
[625,] 5.766558e-04 1.153312e-03 0.9994233
[626,] 3.650191e-04 7.300382e-04 0.9996350
[627,] 2.182819e-03 4.365638e-03 0.9978172
[628,] 1.404696e-03 2.809393e-03 0.9985953
[629,] 8.657444e-04 1.731489e-03 0.9991343
[630,] 1.574490e-03 3.148979e-03 0.9984255
[631,] 9.467942e-04 1.893588e-03 0.9990532
[632,] 7.209183e-04 1.441837e-03 0.9992791
[633,] 4.417472e-04 8.834944e-04 0.9995583
[634,] 3.205196e-04 6.410391e-04 0.9996795
[635,] 2.148641e-04 4.297282e-04 0.9997851
[636,] 1.116257e-04 2.232514e-04 0.9998884
[637,] 7.707622e-05 1.541524e-04 0.9999229
[638,] 8.384015e-05 1.676803e-04 0.9999162
[639,] 1.788827e-04 3.577653e-04 0.9998211
[640,] 8.106760e-05 1.621352e-04 0.9999189
[641,] 4.620747e-04 9.241495e-04 0.9995379
[642,] 2.731672e-04 5.463344e-04 0.9997268
[643,] 6.059864e-02 1.211973e-01 0.9394014
> postscript(file="/var/www/html/freestat/rcomp/tmp/1tfvw1292683540.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> 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()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/2tfvw1292683540.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/346cz1292683540.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/446cz1292683540.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/546cz1292683540.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 654
Frequency = 1
1 2 3 4 5 6
-7.29831264 -4.80840833 -2.51850401 -8.79831264 -5.42831264 1.31159167
7 8 9 10 11 12
1.80140030 -3.23859970 -2.17840833 -4.95831264 -1.36859970 -4.69840833
13 14 15 16 17 18
-0.11840833 -0.86840833 0.53159167 -0.39850401 -0.33869539 1.39140030
19 20 21 22 23 24
-0.85831264 1.66168736 -1.05869539 -2.49869539 -11.88979108 6.06149599
25 26 27 28 29 30
5.50168736 3.64111323 -0.89831264 2.49130461 7.75159167 -3.37831264
31 32 33 34 35 36
-2.23869539 5.62168736 4.68159167 1.21149599 1.06130461 -0.29840833
37 38 39 40 41 42
2.87168736 -1.33840833 -6.57840833 -2.00840833 6.97168736 2.15159167
43 44 45 46 47 48
2.70130461 -2.73840833 -1.62869539 -3.61831264 -3.47850401 -8.60840833
49 50 51 52 53 54
-0.88859970 5.27159167 5.11130461 3.59168736 11.18168736 -5.01840833
55 56 57 58 59 60
-1.04859970 1.32168736 5.78168736 4.47149599 2.96120892 -5.06840833
61 62 63 64 65 66
-0.81840833 2.76159167 -2.57859970 3.04120892 -2.64840833 0.08140030
67 68 69 70 71 72
-3.68831264 -3.40850401 0.30130461 1.22168736 -2.42840833 3.26159167
73 74 75 76 77 78
2.77168736 0.19168736 -3.47831264 -1.82850401 3.52130461 -0.53869539
79 80 81 82 83 84
14.04168736 -0.19859970 1.39149599 -0.43859970 0.06159167 -4.10840833
85 86 87 88 89 90
-0.54850401 -0.60840833 -7.18850401 2.21168736 0.74130461 -3.11831264
91 92 93 94 95 96
5.91168736 5.26168736 -3.60850401 0.10159167 1.83159167 -2.41831264
97 98 99 100 101 102
-6.86831264 -2.72831264 -2.89850401 -6.81840833 -9.12840833 -3.22850401
103 104 105 106 107 108
1.50168736 -11.93279108 1.36168736 2.40140030 1.49159167 -0.20859970
109 110 111 112 113 114
-2.29850401 -3.68850401 4.34159167 -1.42850401 1.12149599 2.48149599
115 116 117 118 119 120
-4.98850401 1.36159167 -2.73850401 -3.09869539 -3.32840833 2.47149599
121 122 123 124 125 126
-0.17831264 -1.44850401 -5.10850401 -4.00859970 -1.14840833 -5.07840833
127 128 129 130 131 132
-6.42840833 -3.97831264 -3.62850401 2.53159167 2.12168736 -7.75840833
133 134 135 136 137 138
-5.29840833 -4.90831264 -1.35840833 -1.66859970 -3.90859970 0.83159167
139 140 141 142 143 144
-6.64831264 -14.26769539 -0.66850401 0.25168736 -3.07859970 1.93168736
145 146 147 148 149 150
6.76168736 -3.02831264 -2.11859970 0.88159167 6.04168736 7.22159167
151 152 153 154 155 156
4.60159167 -1.36831264 0.22168736 1.54168736 -2.21850401 -4.45840833
157 158 159 160 161 162
-2.02859970 -1.78831264 -3.89831264 1.33168736 0.74149599 -4.24840833
163 164 165 166 167 168
0.46130461 -2.51840833 4.55168736 -0.25859970 4.13168736 -1.14859970
169 170 171 172 173 174
4.26159167 3.13130461 -3.13850401 1.86149599 5.16120892 5.66168736
175 176 177 178 179 180
2.98159167 -5.04831264 -3.46831264 -1.21831264 1.98130461 -8.31831264
181 182 183 184 185 186
-8.06850401 3.63140030 -1.17831264 -2.07831264 -7.21831264 -1.81850401
187 188 189 190 191 192
-3.59869539 -0.17840833 2.79168736 -1.75850401 -2.70027897 -8.69840833
193 194 195 196 197 198
-3.18831264 -0.72859970 0.87140030 5.04159167 4.33168736 -4.13869539
199 200 201 202 203 204
5.23140030 3.65140030 -2.71831264 -2.81850401 5.35168736 -0.59840833
205 206 207 208 209 210
4.65130461 1.23149599 -0.41859970 0.17168736 -0.51850401 -2.73831264
211 212 213 214 215 216
-3.07831264 7.88159167 0.91168736 -2.45850401 0.04168736 -1.70879108
217 218 219 220 221 222
-3.81831264 3.02130461 1.00159167 -4.68831264 -6.48840833 0.32111323
223 224 225 226 227 228
-3.95831264 -6.92840833 -3.15859970 12.43168736 -2.13840833 -3.07840833
229 230 231 232 233 234
-6.01850401 -4.00859970 -1.88840833 2.01168736 1.16120892 -3.59850401
235 236 237 238 239 240
-0.69840833 4.76159167 7.85168736 0.57140030 -1.94840833 -2.15859970
241 242 243 244 245 246
-4.60840833 0.47168736 1.30159167 -5.56850401 -3.61831264 2.30159167
247 248 249 250 251 252
-2.43850401 4.12168736 -3.60840833 -6.83831264 -3.95859970 0.98159167
253 254 255 256 257 258
-1.91831264 2.71159167 8.01168736 0.19168736 1.77159167 -3.31850401
259 260 261 262 263 264
0.83130461 0.84149599 0.21159167 -5.51831264 3.95168736 -0.23859970
265 266 267 268 269 270
4.26159167 5.78149599 -5.25831264 -0.94850401 -0.36850401 0.33130461
271 272 273 274 275 276
3.65168736 4.85168736 1.53159167 -1.10840833 -5.82831264 -2.03859970
277 278 279 280 281 282
1.63149599 4.24130461 -2.55831264 -1.46869539 0.30149599 -0.39859970
283 284 285 286 287 288
2.95159167 3.94149599 -1.35840833 1.45120892 1.28159167 0.08130461
289 290 291 292 293 294
-4.46840833 5.63149599 5.92149599 0.49168736 -8.46831264 -5.80840833
295 296 297 298 299 300
3.60168736 -0.47859970 -2.05840833 -5.68831264 -2.68879108 -3.11859970
301 302 303 304 305 306
-1.45850401 2.50168736 -5.47840833 -0.66859970 -1.89840833 3.99149599
307 308 309 310 311 312
6.09130461 1.23159167 -4.76850401 -0.19812126 -3.42821695 -2.21783419
313 314 315 316 317 318
-7.33812126 5.66187874 -0.78802557 -8.59821695 -3.79812126 -0.28821695
319 320 321 322 323 324
1.38216581 9.66207012 12.39216581 14.13197443 -2.95802557 4.95178305
325 326 327 328 329 330
4.46207012 -1.09821695 7.50187874 -1.89821695 -2.38812126 5.32178305
331 332 333 334 335 336
1.79207012 -11.51980052 3.18187874 4.40178305 5.86187874 -2.22792988
337 338 339 340 341 342
-1.50821695 -3.78821695 -2.47802557 -4.24821695 -7.33821695 -0.24821695
343 344 345 346 347 348
0.53187874 11.03187874 -5.17812126 5.80178305 -1.41812126 1.07187874
349 350 351 352 353 354
9.35207012 5.36207012 6.79187874 3.32187874 -1.66783419 1.54187874
355 356 357 358 359 360
3.36178305 8.18178305 1.93178305 0.40019948 -5.50792988 -4.20821695
361 362 363 364 365 366
10.80216581 -3.18821695 2.04187874 6.24207012 12.56197443 0.53010379
367 368 369 370 371 372
-0.78821695 -1.76792988 -9.95008759 5.69019948 -5.57812126 16.34010379
373 374 375 376 377 378
6.11187874 -1.49812126 0.49178305 -4.40812126 -2.70812126 3.77178305
379 380 381 382 383 384
2.07187874 3.75207012 -5.37999190 -5.76821695 5.12178305 -3.13980052
385 386 387 388 389 390
10.73187874 6.83178305 5.07187874 -0.69018328 2.67187874 5.17207012
391 392 393 394 395 396
-7.23792988 -1.66821695 6.83187874 9.41187874 -1.65802557 2.20178305
397 398 399 400 401 402
-8.47821695 -4.86812126 -4.95821695 -1.68812126 6.70216581 -3.93792988
403 404 405 406 407 408
9.13000810 11.61187874 -5.85812126 6.92187874 -2.07783419 -3.85802557
409 410 411 412 413 414
0.02178305 -5.80802557 -3.61812126 -0.62821695 -5.49812126 -5.55989621
415 416 417 418 419 420
-12.12821695 11.24178305 -4.12812126 -3.09792988 -1.25821695 0.89187874
421 422 423 424 425 426
6.34178305 16.01010379 8.70187874 1.84178305 0.84187874 -4.05812126
427 428 429 430 431 432
-4.33783419 -0.16812126 -6.65792988 6.11187874 5.21187874 -3.83821695
433 434 435 436 437 438
-9.35792988 -3.05821695 -0.68989621 0.25178305 -2.68802557 1.42178305
439 440 441 442 443 444
9.21019948 1.18197443 9.56981672 16.89187874 5.08216581 -12.20802557
445 446 447 448 449 450
-8.12821695 5.18981672 6.79178305 2.30178305 -8.95802557 8.55207012
451 452 453 454 455 456
3.19187874 20.87197443 11.69187874 9.43197443 -2.97812126 4.79991241
457 458 459 460 461 462
12.74197443 6.54197443 -2.80792988 -4.38812126 4.21000810 5.29178305
463 464 465 466 467 468
-2.23802557 15.07207012 -3.12821695 1.42197443 -0.89821695 -7.89802557
469 470 471 472 473 474
0.20178305 -0.76812126 -7.97821695 8.29000810 -10.64783419 0.87178305
475 476 477 478 479 480
3.79178305 -0.57802557 -4.69821695 -8.12821695 -9.38989621 2.66197443
481 482 483 484 485 486
3.64197443 -3.25812126 -0.76989621 2.64000810 -10.05821695 -8.22812126
487 488 489 490 491 492
0.70187874 1.42010379 9.36197443 10.78207012 6.14207012 -4.20821695
493 494 495 496 497 498
-3.84802557 -6.17999190 0.77216581 4.12991241 1.91178305 -4.77989621
499 500 501 502 503 504
1.19187874 7.77187874 -1.14792988 5.55197443 -3.14821695 19.20178305
505 506 507 508 509 510
-0.54802557 1.42010379 3.54187874 -4.54821695 3.43187874 14.20187874
511 512 513 514 515 516
-5.41812126 -3.40812126 3.02187874 -0.17783419 15.83197443 -0.38792988
517 518 519 520 521 522
17.13207012 10.41981672 -7.19802557 -3.06821695 -3.29821695 -1.62999190
523 524 525 526 527 528
2.63991241 0.94197443 1.74187874 4.65187874 3.92197443 -2.70802557
529 530 531 532 533 534
-7.11783419 1.18187874 4.56178305 -0.37812126 -5.31802557 -3.13792988
535 536 537 538 539 540
-5.67802557 -4.02812126 4.16187874 -7.80812126 1.77216581 9.43187874
541 542 543 544 545 546
6.30010379 7.87197443 -5.38821695 -3.84802557 -9.20792988 7.85178305
547 548 549 550 551 552
4.02178305 0.17178305 6.58178305 6.68216581 3.93197443 0.24187874
553 554 555 556 557 558
-2.14812126 -8.04802557 -6.68783419 -11.11980052 4.39178305 -0.09812126
559 560 561 562 563 564
19.47991241 -2.35821695 1.67178305 -1.01802557 16.94197443 -0.91812126
565 566 567 568 569 570
-4.22802557 4.15178305 3.82197443 8.62207012 0.99178305 2.04197443
571 572 573 574 575 576
4.19178305 -8.38792988 -3.14802557 5.81981672 -2.01802557 -1.02821695
577 578 579 580 581 582
-5.16812126 -6.37802557 12.26187874 -1.35802557 4.61178305 2.07207012
583 584 585 586 587 588
0.85197443 3.76187874 -2.44812126 -0.81812126 -7.63812126 10.66197443
589 590 591 592 593 594
-6.05783419 4.30991241 -3.41812126 -0.54802557 2.03216581 6.49991241
595 596 597 598 599 600
-0.02989621 -4.45812126 -2.79821695 -2.28792988 -5.57802557 4.14991241
601 602 603 604 605 606
8.90010379 13.75019948 -5.70821695 3.79991241 -1.18812126 4.48226150
607 608 609 610 611 612
8.85029517 -16.28745143 3.34264426 -10.33932207 -1.65960914 1.27245288
613 614 615 616 617 618
4.43226150 6.64226150 -12.00773850 -9.41970483 -14.22970483 -12.07932207
619 620 621 622 623 624
-14.52745143 -6.81764281 -10.04951345 -9.50960914 -6.16970483 19.57226150
625 626 627 628 629 630
1.49226150 -3.14745143 4.22245288 -1.04773850 -6.80951345 -8.01754712
631 632 633 634 635 636
-3.94764281 13.31245288 -4.98970483 -2.90970483 -12.45960914 -4.23764281
637 638 639 640 641 642
1.78029517 -2.33941776 -9.48773850 2.16039086 -3.41754712 2.30235719
643 644 645 646 647 648
-10.87764281 -13.56970483 0.84058224 -13.42970483 4.35235719 13.36245288
649 650 651 652 653 654
10.18039086 4.16039086 1.06029517 -16.81745143 -10.58960914 -6.24773850
> postscript(file="/var/www/html/freestat/rcomp/tmp/6110q1292683540.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> dum <- cbind(lag(myerror,k=1),myerror)
> dum
Time Series:
Start = 0
End = 654
Frequency = 1
lag(myerror, k = 1) myerror
0 -7.29831264 NA
1 -4.80840833 -7.29831264
2 -2.51850401 -4.80840833
3 -8.79831264 -2.51850401
4 -5.42831264 -8.79831264
5 1.31159167 -5.42831264
6 1.80140030 1.31159167
7 -3.23859970 1.80140030
8 -2.17840833 -3.23859970
9 -4.95831264 -2.17840833
10 -1.36859970 -4.95831264
11 -4.69840833 -1.36859970
12 -0.11840833 -4.69840833
13 -0.86840833 -0.11840833
14 0.53159167 -0.86840833
15 -0.39850401 0.53159167
16 -0.33869539 -0.39850401
17 1.39140030 -0.33869539
18 -0.85831264 1.39140030
19 1.66168736 -0.85831264
20 -1.05869539 1.66168736
21 -2.49869539 -1.05869539
22 -11.88979108 -2.49869539
23 6.06149599 -11.88979108
24 5.50168736 6.06149599
25 3.64111323 5.50168736
26 -0.89831264 3.64111323
27 2.49130461 -0.89831264
28 7.75159167 2.49130461
29 -3.37831264 7.75159167
30 -2.23869539 -3.37831264
31 5.62168736 -2.23869539
32 4.68159167 5.62168736
33 1.21149599 4.68159167
34 1.06130461 1.21149599
35 -0.29840833 1.06130461
36 2.87168736 -0.29840833
37 -1.33840833 2.87168736
38 -6.57840833 -1.33840833
39 -2.00840833 -6.57840833
40 6.97168736 -2.00840833
41 2.15159167 6.97168736
42 2.70130461 2.15159167
43 -2.73840833 2.70130461
44 -1.62869539 -2.73840833
45 -3.61831264 -1.62869539
46 -3.47850401 -3.61831264
47 -8.60840833 -3.47850401
48 -0.88859970 -8.60840833
49 5.27159167 -0.88859970
50 5.11130461 5.27159167
51 3.59168736 5.11130461
52 11.18168736 3.59168736
53 -5.01840833 11.18168736
54 -1.04859970 -5.01840833
55 1.32168736 -1.04859970
56 5.78168736 1.32168736
57 4.47149599 5.78168736
58 2.96120892 4.47149599
59 -5.06840833 2.96120892
60 -0.81840833 -5.06840833
61 2.76159167 -0.81840833
62 -2.57859970 2.76159167
63 3.04120892 -2.57859970
64 -2.64840833 3.04120892
65 0.08140030 -2.64840833
66 -3.68831264 0.08140030
67 -3.40850401 -3.68831264
68 0.30130461 -3.40850401
69 1.22168736 0.30130461
70 -2.42840833 1.22168736
71 3.26159167 -2.42840833
72 2.77168736 3.26159167
73 0.19168736 2.77168736
74 -3.47831264 0.19168736
75 -1.82850401 -3.47831264
76 3.52130461 -1.82850401
77 -0.53869539 3.52130461
78 14.04168736 -0.53869539
79 -0.19859970 14.04168736
80 1.39149599 -0.19859970
81 -0.43859970 1.39149599
82 0.06159167 -0.43859970
83 -4.10840833 0.06159167
84 -0.54850401 -4.10840833
85 -0.60840833 -0.54850401
86 -7.18850401 -0.60840833
87 2.21168736 -7.18850401
88 0.74130461 2.21168736
89 -3.11831264 0.74130461
90 5.91168736 -3.11831264
91 5.26168736 5.91168736
92 -3.60850401 5.26168736
93 0.10159167 -3.60850401
94 1.83159167 0.10159167
95 -2.41831264 1.83159167
96 -6.86831264 -2.41831264
97 -2.72831264 -6.86831264
98 -2.89850401 -2.72831264
99 -6.81840833 -2.89850401
100 -9.12840833 -6.81840833
101 -3.22850401 -9.12840833
102 1.50168736 -3.22850401
103 -11.93279108 1.50168736
104 1.36168736 -11.93279108
105 2.40140030 1.36168736
106 1.49159167 2.40140030
107 -0.20859970 1.49159167
108 -2.29850401 -0.20859970
109 -3.68850401 -2.29850401
110 4.34159167 -3.68850401
111 -1.42850401 4.34159167
112 1.12149599 -1.42850401
113 2.48149599 1.12149599
114 -4.98850401 2.48149599
115 1.36159167 -4.98850401
116 -2.73850401 1.36159167
117 -3.09869539 -2.73850401
118 -3.32840833 -3.09869539
119 2.47149599 -3.32840833
120 -0.17831264 2.47149599
121 -1.44850401 -0.17831264
122 -5.10850401 -1.44850401
123 -4.00859970 -5.10850401
124 -1.14840833 -4.00859970
125 -5.07840833 -1.14840833
126 -6.42840833 -5.07840833
127 -3.97831264 -6.42840833
128 -3.62850401 -3.97831264
129 2.53159167 -3.62850401
130 2.12168736 2.53159167
131 -7.75840833 2.12168736
132 -5.29840833 -7.75840833
133 -4.90831264 -5.29840833
134 -1.35840833 -4.90831264
135 -1.66859970 -1.35840833
136 -3.90859970 -1.66859970
137 0.83159167 -3.90859970
138 -6.64831264 0.83159167
139 -14.26769539 -6.64831264
140 -0.66850401 -14.26769539
141 0.25168736 -0.66850401
142 -3.07859970 0.25168736
143 1.93168736 -3.07859970
144 6.76168736 1.93168736
145 -3.02831264 6.76168736
146 -2.11859970 -3.02831264
147 0.88159167 -2.11859970
148 6.04168736 0.88159167
149 7.22159167 6.04168736
150 4.60159167 7.22159167
151 -1.36831264 4.60159167
152 0.22168736 -1.36831264
153 1.54168736 0.22168736
154 -2.21850401 1.54168736
155 -4.45840833 -2.21850401
156 -2.02859970 -4.45840833
157 -1.78831264 -2.02859970
158 -3.89831264 -1.78831264
159 1.33168736 -3.89831264
160 0.74149599 1.33168736
161 -4.24840833 0.74149599
162 0.46130461 -4.24840833
163 -2.51840833 0.46130461
164 4.55168736 -2.51840833
165 -0.25859970 4.55168736
166 4.13168736 -0.25859970
167 -1.14859970 4.13168736
168 4.26159167 -1.14859970
169 3.13130461 4.26159167
170 -3.13850401 3.13130461
171 1.86149599 -3.13850401
172 5.16120892 1.86149599
173 5.66168736 5.16120892
174 2.98159167 5.66168736
175 -5.04831264 2.98159167
176 -3.46831264 -5.04831264
177 -1.21831264 -3.46831264
178 1.98130461 -1.21831264
179 -8.31831264 1.98130461
180 -8.06850401 -8.31831264
181 3.63140030 -8.06850401
182 -1.17831264 3.63140030
183 -2.07831264 -1.17831264
184 -7.21831264 -2.07831264
185 -1.81850401 -7.21831264
186 -3.59869539 -1.81850401
187 -0.17840833 -3.59869539
188 2.79168736 -0.17840833
189 -1.75850401 2.79168736
190 -2.70027897 -1.75850401
191 -8.69840833 -2.70027897
192 -3.18831264 -8.69840833
193 -0.72859970 -3.18831264
194 0.87140030 -0.72859970
195 5.04159167 0.87140030
196 4.33168736 5.04159167
197 -4.13869539 4.33168736
198 5.23140030 -4.13869539
199 3.65140030 5.23140030
200 -2.71831264 3.65140030
201 -2.81850401 -2.71831264
202 5.35168736 -2.81850401
203 -0.59840833 5.35168736
204 4.65130461 -0.59840833
205 1.23149599 4.65130461
206 -0.41859970 1.23149599
207 0.17168736 -0.41859970
208 -0.51850401 0.17168736
209 -2.73831264 -0.51850401
210 -3.07831264 -2.73831264
211 7.88159167 -3.07831264
212 0.91168736 7.88159167
213 -2.45850401 0.91168736
214 0.04168736 -2.45850401
215 -1.70879108 0.04168736
216 -3.81831264 -1.70879108
217 3.02130461 -3.81831264
218 1.00159167 3.02130461
219 -4.68831264 1.00159167
220 -6.48840833 -4.68831264
221 0.32111323 -6.48840833
222 -3.95831264 0.32111323
223 -6.92840833 -3.95831264
224 -3.15859970 -6.92840833
225 12.43168736 -3.15859970
226 -2.13840833 12.43168736
227 -3.07840833 -2.13840833
228 -6.01850401 -3.07840833
229 -4.00859970 -6.01850401
230 -1.88840833 -4.00859970
231 2.01168736 -1.88840833
232 1.16120892 2.01168736
233 -3.59850401 1.16120892
234 -0.69840833 -3.59850401
235 4.76159167 -0.69840833
236 7.85168736 4.76159167
237 0.57140030 7.85168736
238 -1.94840833 0.57140030
239 -2.15859970 -1.94840833
240 -4.60840833 -2.15859970
241 0.47168736 -4.60840833
242 1.30159167 0.47168736
243 -5.56850401 1.30159167
244 -3.61831264 -5.56850401
245 2.30159167 -3.61831264
246 -2.43850401 2.30159167
247 4.12168736 -2.43850401
248 -3.60840833 4.12168736
249 -6.83831264 -3.60840833
250 -3.95859970 -6.83831264
251 0.98159167 -3.95859970
252 -1.91831264 0.98159167
253 2.71159167 -1.91831264
254 8.01168736 2.71159167
255 0.19168736 8.01168736
256 1.77159167 0.19168736
257 -3.31850401 1.77159167
258 0.83130461 -3.31850401
259 0.84149599 0.83130461
260 0.21159167 0.84149599
261 -5.51831264 0.21159167
262 3.95168736 -5.51831264
263 -0.23859970 3.95168736
264 4.26159167 -0.23859970
265 5.78149599 4.26159167
266 -5.25831264 5.78149599
267 -0.94850401 -5.25831264
268 -0.36850401 -0.94850401
269 0.33130461 -0.36850401
270 3.65168736 0.33130461
271 4.85168736 3.65168736
272 1.53159167 4.85168736
273 -1.10840833 1.53159167
274 -5.82831264 -1.10840833
275 -2.03859970 -5.82831264
276 1.63149599 -2.03859970
277 4.24130461 1.63149599
278 -2.55831264 4.24130461
279 -1.46869539 -2.55831264
280 0.30149599 -1.46869539
281 -0.39859970 0.30149599
282 2.95159167 -0.39859970
283 3.94149599 2.95159167
284 -1.35840833 3.94149599
285 1.45120892 -1.35840833
286 1.28159167 1.45120892
287 0.08130461 1.28159167
288 -4.46840833 0.08130461
289 5.63149599 -4.46840833
290 5.92149599 5.63149599
291 0.49168736 5.92149599
292 -8.46831264 0.49168736
293 -5.80840833 -8.46831264
294 3.60168736 -5.80840833
295 -0.47859970 3.60168736
296 -2.05840833 -0.47859970
297 -5.68831264 -2.05840833
298 -2.68879108 -5.68831264
299 -3.11859970 -2.68879108
300 -1.45850401 -3.11859970
301 2.50168736 -1.45850401
302 -5.47840833 2.50168736
303 -0.66859970 -5.47840833
304 -1.89840833 -0.66859970
305 3.99149599 -1.89840833
306 6.09130461 3.99149599
307 1.23159167 6.09130461
308 -4.76850401 1.23159167
309 -0.19812126 -4.76850401
310 -3.42821695 -0.19812126
311 -2.21783419 -3.42821695
312 -7.33812126 -2.21783419
313 5.66187874 -7.33812126
314 -0.78802557 5.66187874
315 -8.59821695 -0.78802557
316 -3.79812126 -8.59821695
317 -0.28821695 -3.79812126
318 1.38216581 -0.28821695
319 9.66207012 1.38216581
320 12.39216581 9.66207012
321 14.13197443 12.39216581
322 -2.95802557 14.13197443
323 4.95178305 -2.95802557
324 4.46207012 4.95178305
325 -1.09821695 4.46207012
326 7.50187874 -1.09821695
327 -1.89821695 7.50187874
328 -2.38812126 -1.89821695
329 5.32178305 -2.38812126
330 1.79207012 5.32178305
331 -11.51980052 1.79207012
332 3.18187874 -11.51980052
333 4.40178305 3.18187874
334 5.86187874 4.40178305
335 -2.22792988 5.86187874
336 -1.50821695 -2.22792988
337 -3.78821695 -1.50821695
338 -2.47802557 -3.78821695
339 -4.24821695 -2.47802557
340 -7.33821695 -4.24821695
341 -0.24821695 -7.33821695
342 0.53187874 -0.24821695
343 11.03187874 0.53187874
344 -5.17812126 11.03187874
345 5.80178305 -5.17812126
346 -1.41812126 5.80178305
347 1.07187874 -1.41812126
348 9.35207012 1.07187874
349 5.36207012 9.35207012
350 6.79187874 5.36207012
351 3.32187874 6.79187874
352 -1.66783419 3.32187874
353 1.54187874 -1.66783419
354 3.36178305 1.54187874
355 8.18178305 3.36178305
356 1.93178305 8.18178305
357 0.40019948 1.93178305
358 -5.50792988 0.40019948
359 -4.20821695 -5.50792988
360 10.80216581 -4.20821695
361 -3.18821695 10.80216581
362 2.04187874 -3.18821695
363 6.24207012 2.04187874
364 12.56197443 6.24207012
365 0.53010379 12.56197443
366 -0.78821695 0.53010379
367 -1.76792988 -0.78821695
368 -9.95008759 -1.76792988
369 5.69019948 -9.95008759
370 -5.57812126 5.69019948
371 16.34010379 -5.57812126
372 6.11187874 16.34010379
373 -1.49812126 6.11187874
374 0.49178305 -1.49812126
375 -4.40812126 0.49178305
376 -2.70812126 -4.40812126
377 3.77178305 -2.70812126
378 2.07187874 3.77178305
379 3.75207012 2.07187874
380 -5.37999190 3.75207012
381 -5.76821695 -5.37999190
382 5.12178305 -5.76821695
383 -3.13980052 5.12178305
384 10.73187874 -3.13980052
385 6.83178305 10.73187874
386 5.07187874 6.83178305
387 -0.69018328 5.07187874
388 2.67187874 -0.69018328
389 5.17207012 2.67187874
390 -7.23792988 5.17207012
391 -1.66821695 -7.23792988
392 6.83187874 -1.66821695
393 9.41187874 6.83187874
394 -1.65802557 9.41187874
395 2.20178305 -1.65802557
396 -8.47821695 2.20178305
397 -4.86812126 -8.47821695
398 -4.95821695 -4.86812126
399 -1.68812126 -4.95821695
400 6.70216581 -1.68812126
401 -3.93792988 6.70216581
402 9.13000810 -3.93792988
403 11.61187874 9.13000810
404 -5.85812126 11.61187874
405 6.92187874 -5.85812126
406 -2.07783419 6.92187874
407 -3.85802557 -2.07783419
408 0.02178305 -3.85802557
409 -5.80802557 0.02178305
410 -3.61812126 -5.80802557
411 -0.62821695 -3.61812126
412 -5.49812126 -0.62821695
413 -5.55989621 -5.49812126
414 -12.12821695 -5.55989621
415 11.24178305 -12.12821695
416 -4.12812126 11.24178305
417 -3.09792988 -4.12812126
418 -1.25821695 -3.09792988
419 0.89187874 -1.25821695
420 6.34178305 0.89187874
421 16.01010379 6.34178305
422 8.70187874 16.01010379
423 1.84178305 8.70187874
424 0.84187874 1.84178305
425 -4.05812126 0.84187874
426 -4.33783419 -4.05812126
427 -0.16812126 -4.33783419
428 -6.65792988 -0.16812126
429 6.11187874 -6.65792988
430 5.21187874 6.11187874
431 -3.83821695 5.21187874
432 -9.35792988 -3.83821695
433 -3.05821695 -9.35792988
434 -0.68989621 -3.05821695
435 0.25178305 -0.68989621
436 -2.68802557 0.25178305
437 1.42178305 -2.68802557
438 9.21019948 1.42178305
439 1.18197443 9.21019948
440 9.56981672 1.18197443
441 16.89187874 9.56981672
442 5.08216581 16.89187874
443 -12.20802557 5.08216581
444 -8.12821695 -12.20802557
445 5.18981672 -8.12821695
446 6.79178305 5.18981672
447 2.30178305 6.79178305
448 -8.95802557 2.30178305
449 8.55207012 -8.95802557
450 3.19187874 8.55207012
451 20.87197443 3.19187874
452 11.69187874 20.87197443
453 9.43197443 11.69187874
454 -2.97812126 9.43197443
455 4.79991241 -2.97812126
456 12.74197443 4.79991241
457 6.54197443 12.74197443
458 -2.80792988 6.54197443
459 -4.38812126 -2.80792988
460 4.21000810 -4.38812126
461 5.29178305 4.21000810
462 -2.23802557 5.29178305
463 15.07207012 -2.23802557
464 -3.12821695 15.07207012
465 1.42197443 -3.12821695
466 -0.89821695 1.42197443
467 -7.89802557 -0.89821695
468 0.20178305 -7.89802557
469 -0.76812126 0.20178305
470 -7.97821695 -0.76812126
471 8.29000810 -7.97821695
472 -10.64783419 8.29000810
473 0.87178305 -10.64783419
474 3.79178305 0.87178305
475 -0.57802557 3.79178305
476 -4.69821695 -0.57802557
477 -8.12821695 -4.69821695
478 -9.38989621 -8.12821695
479 2.66197443 -9.38989621
480 3.64197443 2.66197443
481 -3.25812126 3.64197443
482 -0.76989621 -3.25812126
483 2.64000810 -0.76989621
484 -10.05821695 2.64000810
485 -8.22812126 -10.05821695
486 0.70187874 -8.22812126
487 1.42010379 0.70187874
488 9.36197443 1.42010379
489 10.78207012 9.36197443
490 6.14207012 10.78207012
491 -4.20821695 6.14207012
492 -3.84802557 -4.20821695
493 -6.17999190 -3.84802557
494 0.77216581 -6.17999190
495 4.12991241 0.77216581
496 1.91178305 4.12991241
497 -4.77989621 1.91178305
498 1.19187874 -4.77989621
499 7.77187874 1.19187874
500 -1.14792988 7.77187874
501 5.55197443 -1.14792988
502 -3.14821695 5.55197443
503 19.20178305 -3.14821695
504 -0.54802557 19.20178305
505 1.42010379 -0.54802557
506 3.54187874 1.42010379
507 -4.54821695 3.54187874
508 3.43187874 -4.54821695
509 14.20187874 3.43187874
510 -5.41812126 14.20187874
511 -3.40812126 -5.41812126
512 3.02187874 -3.40812126
513 -0.17783419 3.02187874
514 15.83197443 -0.17783419
515 -0.38792988 15.83197443
516 17.13207012 -0.38792988
517 10.41981672 17.13207012
518 -7.19802557 10.41981672
519 -3.06821695 -7.19802557
520 -3.29821695 -3.06821695
521 -1.62999190 -3.29821695
522 2.63991241 -1.62999190
523 0.94197443 2.63991241
524 1.74187874 0.94197443
525 4.65187874 1.74187874
526 3.92197443 4.65187874
527 -2.70802557 3.92197443
528 -7.11783419 -2.70802557
529 1.18187874 -7.11783419
530 4.56178305 1.18187874
531 -0.37812126 4.56178305
532 -5.31802557 -0.37812126
533 -3.13792988 -5.31802557
534 -5.67802557 -3.13792988
535 -4.02812126 -5.67802557
536 4.16187874 -4.02812126
537 -7.80812126 4.16187874
538 1.77216581 -7.80812126
539 9.43187874 1.77216581
540 6.30010379 9.43187874
541 7.87197443 6.30010379
542 -5.38821695 7.87197443
543 -3.84802557 -5.38821695
544 -9.20792988 -3.84802557
545 7.85178305 -9.20792988
546 4.02178305 7.85178305
547 0.17178305 4.02178305
548 6.58178305 0.17178305
549 6.68216581 6.58178305
550 3.93197443 6.68216581
551 0.24187874 3.93197443
552 -2.14812126 0.24187874
553 -8.04802557 -2.14812126
554 -6.68783419 -8.04802557
555 -11.11980052 -6.68783419
556 4.39178305 -11.11980052
557 -0.09812126 4.39178305
558 19.47991241 -0.09812126
559 -2.35821695 19.47991241
560 1.67178305 -2.35821695
561 -1.01802557 1.67178305
562 16.94197443 -1.01802557
563 -0.91812126 16.94197443
564 -4.22802557 -0.91812126
565 4.15178305 -4.22802557
566 3.82197443 4.15178305
567 8.62207012 3.82197443
568 0.99178305 8.62207012
569 2.04197443 0.99178305
570 4.19178305 2.04197443
571 -8.38792988 4.19178305
572 -3.14802557 -8.38792988
573 5.81981672 -3.14802557
574 -2.01802557 5.81981672
575 -1.02821695 -2.01802557
576 -5.16812126 -1.02821695
577 -6.37802557 -5.16812126
578 12.26187874 -6.37802557
579 -1.35802557 12.26187874
580 4.61178305 -1.35802557
581 2.07207012 4.61178305
582 0.85197443 2.07207012
583 3.76187874 0.85197443
584 -2.44812126 3.76187874
585 -0.81812126 -2.44812126
586 -7.63812126 -0.81812126
587 10.66197443 -7.63812126
588 -6.05783419 10.66197443
589 4.30991241 -6.05783419
590 -3.41812126 4.30991241
591 -0.54802557 -3.41812126
592 2.03216581 -0.54802557
593 6.49991241 2.03216581
594 -0.02989621 6.49991241
595 -4.45812126 -0.02989621
596 -2.79821695 -4.45812126
597 -2.28792988 -2.79821695
598 -5.57802557 -2.28792988
599 4.14991241 -5.57802557
600 8.90010379 4.14991241
601 13.75019948 8.90010379
602 -5.70821695 13.75019948
603 3.79991241 -5.70821695
604 -1.18812126 3.79991241
605 4.48226150 -1.18812126
606 8.85029517 4.48226150
607 -16.28745143 8.85029517
608 3.34264426 -16.28745143
609 -10.33932207 3.34264426
610 -1.65960914 -10.33932207
611 1.27245288 -1.65960914
612 4.43226150 1.27245288
613 6.64226150 4.43226150
614 -12.00773850 6.64226150
615 -9.41970483 -12.00773850
616 -14.22970483 -9.41970483
617 -12.07932207 -14.22970483
618 -14.52745143 -12.07932207
619 -6.81764281 -14.52745143
620 -10.04951345 -6.81764281
621 -9.50960914 -10.04951345
622 -6.16970483 -9.50960914
623 19.57226150 -6.16970483
624 1.49226150 19.57226150
625 -3.14745143 1.49226150
626 4.22245288 -3.14745143
627 -1.04773850 4.22245288
628 -6.80951345 -1.04773850
629 -8.01754712 -6.80951345
630 -3.94764281 -8.01754712
631 13.31245288 -3.94764281
632 -4.98970483 13.31245288
633 -2.90970483 -4.98970483
634 -12.45960914 -2.90970483
635 -4.23764281 -12.45960914
636 1.78029517 -4.23764281
637 -2.33941776 1.78029517
638 -9.48773850 -2.33941776
639 2.16039086 -9.48773850
640 -3.41754712 2.16039086
641 2.30235719 -3.41754712
642 -10.87764281 2.30235719
643 -13.56970483 -10.87764281
644 0.84058224 -13.56970483
645 -13.42970483 0.84058224
646 4.35235719 -13.42970483
647 13.36245288 4.35235719
648 10.18039086 13.36245288
649 4.16039086 10.18039086
650 1.06029517 4.16039086
651 -16.81745143 1.06029517
652 -10.58960914 -16.81745143
653 -6.24773850 -10.58960914
654 NA -6.24773850
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -4.80840833 -7.29831264
[2,] -2.51850401 -4.80840833
[3,] -8.79831264 -2.51850401
[4,] -5.42831264 -8.79831264
[5,] 1.31159167 -5.42831264
[6,] 1.80140030 1.31159167
[7,] -3.23859970 1.80140030
[8,] -2.17840833 -3.23859970
[9,] -4.95831264 -2.17840833
[10,] -1.36859970 -4.95831264
[11,] -4.69840833 -1.36859970
[12,] -0.11840833 -4.69840833
[13,] -0.86840833 -0.11840833
[14,] 0.53159167 -0.86840833
[15,] -0.39850401 0.53159167
[16,] -0.33869539 -0.39850401
[17,] 1.39140030 -0.33869539
[18,] -0.85831264 1.39140030
[19,] 1.66168736 -0.85831264
[20,] -1.05869539 1.66168736
[21,] -2.49869539 -1.05869539
[22,] -11.88979108 -2.49869539
[23,] 6.06149599 -11.88979108
[24,] 5.50168736 6.06149599
[25,] 3.64111323 5.50168736
[26,] -0.89831264 3.64111323
[27,] 2.49130461 -0.89831264
[28,] 7.75159167 2.49130461
[29,] -3.37831264 7.75159167
[30,] -2.23869539 -3.37831264
[31,] 5.62168736 -2.23869539
[32,] 4.68159167 5.62168736
[33,] 1.21149599 4.68159167
[34,] 1.06130461 1.21149599
[35,] -0.29840833 1.06130461
[36,] 2.87168736 -0.29840833
[37,] -1.33840833 2.87168736
[38,] -6.57840833 -1.33840833
[39,] -2.00840833 -6.57840833
[40,] 6.97168736 -2.00840833
[41,] 2.15159167 6.97168736
[42,] 2.70130461 2.15159167
[43,] -2.73840833 2.70130461
[44,] -1.62869539 -2.73840833
[45,] -3.61831264 -1.62869539
[46,] -3.47850401 -3.61831264
[47,] -8.60840833 -3.47850401
[48,] -0.88859970 -8.60840833
[49,] 5.27159167 -0.88859970
[50,] 5.11130461 5.27159167
[51,] 3.59168736 5.11130461
[52,] 11.18168736 3.59168736
[53,] -5.01840833 11.18168736
[54,] -1.04859970 -5.01840833
[55,] 1.32168736 -1.04859970
[56,] 5.78168736 1.32168736
[57,] 4.47149599 5.78168736
[58,] 2.96120892 4.47149599
[59,] -5.06840833 2.96120892
[60,] -0.81840833 -5.06840833
[61,] 2.76159167 -0.81840833
[62,] -2.57859970 2.76159167
[63,] 3.04120892 -2.57859970
[64,] -2.64840833 3.04120892
[65,] 0.08140030 -2.64840833
[66,] -3.68831264 0.08140030
[67,] -3.40850401 -3.68831264
[68,] 0.30130461 -3.40850401
[69,] 1.22168736 0.30130461
[70,] -2.42840833 1.22168736
[71,] 3.26159167 -2.42840833
[72,] 2.77168736 3.26159167
[73,] 0.19168736 2.77168736
[74,] -3.47831264 0.19168736
[75,] -1.82850401 -3.47831264
[76,] 3.52130461 -1.82850401
[77,] -0.53869539 3.52130461
[78,] 14.04168736 -0.53869539
[79,] -0.19859970 14.04168736
[80,] 1.39149599 -0.19859970
[81,] -0.43859970 1.39149599
[82,] 0.06159167 -0.43859970
[83,] -4.10840833 0.06159167
[84,] -0.54850401 -4.10840833
[85,] -0.60840833 -0.54850401
[86,] -7.18850401 -0.60840833
[87,] 2.21168736 -7.18850401
[88,] 0.74130461 2.21168736
[89,] -3.11831264 0.74130461
[90,] 5.91168736 -3.11831264
[91,] 5.26168736 5.91168736
[92,] -3.60850401 5.26168736
[93,] 0.10159167 -3.60850401
[94,] 1.83159167 0.10159167
[95,] -2.41831264 1.83159167
[96,] -6.86831264 -2.41831264
[97,] -2.72831264 -6.86831264
[98,] -2.89850401 -2.72831264
[99,] -6.81840833 -2.89850401
[100,] -9.12840833 -6.81840833
[101,] -3.22850401 -9.12840833
[102,] 1.50168736 -3.22850401
[103,] -11.93279108 1.50168736
[104,] 1.36168736 -11.93279108
[105,] 2.40140030 1.36168736
[106,] 1.49159167 2.40140030
[107,] -0.20859970 1.49159167
[108,] -2.29850401 -0.20859970
[109,] -3.68850401 -2.29850401
[110,] 4.34159167 -3.68850401
[111,] -1.42850401 4.34159167
[112,] 1.12149599 -1.42850401
[113,] 2.48149599 1.12149599
[114,] -4.98850401 2.48149599
[115,] 1.36159167 -4.98850401
[116,] -2.73850401 1.36159167
[117,] -3.09869539 -2.73850401
[118,] -3.32840833 -3.09869539
[119,] 2.47149599 -3.32840833
[120,] -0.17831264 2.47149599
[121,] -1.44850401 -0.17831264
[122,] -5.10850401 -1.44850401
[123,] -4.00859970 -5.10850401
[124,] -1.14840833 -4.00859970
[125,] -5.07840833 -1.14840833
[126,] -6.42840833 -5.07840833
[127,] -3.97831264 -6.42840833
[128,] -3.62850401 -3.97831264
[129,] 2.53159167 -3.62850401
[130,] 2.12168736 2.53159167
[131,] -7.75840833 2.12168736
[132,] -5.29840833 -7.75840833
[133,] -4.90831264 -5.29840833
[134,] -1.35840833 -4.90831264
[135,] -1.66859970 -1.35840833
[136,] -3.90859970 -1.66859970
[137,] 0.83159167 -3.90859970
[138,] -6.64831264 0.83159167
[139,] -14.26769539 -6.64831264
[140,] -0.66850401 -14.26769539
[141,] 0.25168736 -0.66850401
[142,] -3.07859970 0.25168736
[143,] 1.93168736 -3.07859970
[144,] 6.76168736 1.93168736
[145,] -3.02831264 6.76168736
[146,] -2.11859970 -3.02831264
[147,] 0.88159167 -2.11859970
[148,] 6.04168736 0.88159167
[149,] 7.22159167 6.04168736
[150,] 4.60159167 7.22159167
[151,] -1.36831264 4.60159167
[152,] 0.22168736 -1.36831264
[153,] 1.54168736 0.22168736
[154,] -2.21850401 1.54168736
[155,] -4.45840833 -2.21850401
[156,] -2.02859970 -4.45840833
[157,] -1.78831264 -2.02859970
[158,] -3.89831264 -1.78831264
[159,] 1.33168736 -3.89831264
[160,] 0.74149599 1.33168736
[161,] -4.24840833 0.74149599
[162,] 0.46130461 -4.24840833
[163,] -2.51840833 0.46130461
[164,] 4.55168736 -2.51840833
[165,] -0.25859970 4.55168736
[166,] 4.13168736 -0.25859970
[167,] -1.14859970 4.13168736
[168,] 4.26159167 -1.14859970
[169,] 3.13130461 4.26159167
[170,] -3.13850401 3.13130461
[171,] 1.86149599 -3.13850401
[172,] 5.16120892 1.86149599
[173,] 5.66168736 5.16120892
[174,] 2.98159167 5.66168736
[175,] -5.04831264 2.98159167
[176,] -3.46831264 -5.04831264
[177,] -1.21831264 -3.46831264
[178,] 1.98130461 -1.21831264
[179,] -8.31831264 1.98130461
[180,] -8.06850401 -8.31831264
[181,] 3.63140030 -8.06850401
[182,] -1.17831264 3.63140030
[183,] -2.07831264 -1.17831264
[184,] -7.21831264 -2.07831264
[185,] -1.81850401 -7.21831264
[186,] -3.59869539 -1.81850401
[187,] -0.17840833 -3.59869539
[188,] 2.79168736 -0.17840833
[189,] -1.75850401 2.79168736
[190,] -2.70027897 -1.75850401
[191,] -8.69840833 -2.70027897
[192,] -3.18831264 -8.69840833
[193,] -0.72859970 -3.18831264
[194,] 0.87140030 -0.72859970
[195,] 5.04159167 0.87140030
[196,] 4.33168736 5.04159167
[197,] -4.13869539 4.33168736
[198,] 5.23140030 -4.13869539
[199,] 3.65140030 5.23140030
[200,] -2.71831264 3.65140030
[201,] -2.81850401 -2.71831264
[202,] 5.35168736 -2.81850401
[203,] -0.59840833 5.35168736
[204,] 4.65130461 -0.59840833
[205,] 1.23149599 4.65130461
[206,] -0.41859970 1.23149599
[207,] 0.17168736 -0.41859970
[208,] -0.51850401 0.17168736
[209,] -2.73831264 -0.51850401
[210,] -3.07831264 -2.73831264
[211,] 7.88159167 -3.07831264
[212,] 0.91168736 7.88159167
[213,] -2.45850401 0.91168736
[214,] 0.04168736 -2.45850401
[215,] -1.70879108 0.04168736
[216,] -3.81831264 -1.70879108
[217,] 3.02130461 -3.81831264
[218,] 1.00159167 3.02130461
[219,] -4.68831264 1.00159167
[220,] -6.48840833 -4.68831264
[221,] 0.32111323 -6.48840833
[222,] -3.95831264 0.32111323
[223,] -6.92840833 -3.95831264
[224,] -3.15859970 -6.92840833
[225,] 12.43168736 -3.15859970
[226,] -2.13840833 12.43168736
[227,] -3.07840833 -2.13840833
[228,] -6.01850401 -3.07840833
[229,] -4.00859970 -6.01850401
[230,] -1.88840833 -4.00859970
[231,] 2.01168736 -1.88840833
[232,] 1.16120892 2.01168736
[233,] -3.59850401 1.16120892
[234,] -0.69840833 -3.59850401
[235,] 4.76159167 -0.69840833
[236,] 7.85168736 4.76159167
[237,] 0.57140030 7.85168736
[238,] -1.94840833 0.57140030
[239,] -2.15859970 -1.94840833
[240,] -4.60840833 -2.15859970
[241,] 0.47168736 -4.60840833
[242,] 1.30159167 0.47168736
[243,] -5.56850401 1.30159167
[244,] -3.61831264 -5.56850401
[245,] 2.30159167 -3.61831264
[246,] -2.43850401 2.30159167
[247,] 4.12168736 -2.43850401
[248,] -3.60840833 4.12168736
[249,] -6.83831264 -3.60840833
[250,] -3.95859970 -6.83831264
[251,] 0.98159167 -3.95859970
[252,] -1.91831264 0.98159167
[253,] 2.71159167 -1.91831264
[254,] 8.01168736 2.71159167
[255,] 0.19168736 8.01168736
[256,] 1.77159167 0.19168736
[257,] -3.31850401 1.77159167
[258,] 0.83130461 -3.31850401
[259,] 0.84149599 0.83130461
[260,] 0.21159167 0.84149599
[261,] -5.51831264 0.21159167
[262,] 3.95168736 -5.51831264
[263,] -0.23859970 3.95168736
[264,] 4.26159167 -0.23859970
[265,] 5.78149599 4.26159167
[266,] -5.25831264 5.78149599
[267,] -0.94850401 -5.25831264
[268,] -0.36850401 -0.94850401
[269,] 0.33130461 -0.36850401
[270,] 3.65168736 0.33130461
[271,] 4.85168736 3.65168736
[272,] 1.53159167 4.85168736
[273,] -1.10840833 1.53159167
[274,] -5.82831264 -1.10840833
[275,] -2.03859970 -5.82831264
[276,] 1.63149599 -2.03859970
[277,] 4.24130461 1.63149599
[278,] -2.55831264 4.24130461
[279,] -1.46869539 -2.55831264
[280,] 0.30149599 -1.46869539
[281,] -0.39859970 0.30149599
[282,] 2.95159167 -0.39859970
[283,] 3.94149599 2.95159167
[284,] -1.35840833 3.94149599
[285,] 1.45120892 -1.35840833
[286,] 1.28159167 1.45120892
[287,] 0.08130461 1.28159167
[288,] -4.46840833 0.08130461
[289,] 5.63149599 -4.46840833
[290,] 5.92149599 5.63149599
[291,] 0.49168736 5.92149599
[292,] -8.46831264 0.49168736
[293,] -5.80840833 -8.46831264
[294,] 3.60168736 -5.80840833
[295,] -0.47859970 3.60168736
[296,] -2.05840833 -0.47859970
[297,] -5.68831264 -2.05840833
[298,] -2.68879108 -5.68831264
[299,] -3.11859970 -2.68879108
[300,] -1.45850401 -3.11859970
[301,] 2.50168736 -1.45850401
[302,] -5.47840833 2.50168736
[303,] -0.66859970 -5.47840833
[304,] -1.89840833 -0.66859970
[305,] 3.99149599 -1.89840833
[306,] 6.09130461 3.99149599
[307,] 1.23159167 6.09130461
[308,] -4.76850401 1.23159167
[309,] -0.19812126 -4.76850401
[310,] -3.42821695 -0.19812126
[311,] -2.21783419 -3.42821695
[312,] -7.33812126 -2.21783419
[313,] 5.66187874 -7.33812126
[314,] -0.78802557 5.66187874
[315,] -8.59821695 -0.78802557
[316,] -3.79812126 -8.59821695
[317,] -0.28821695 -3.79812126
[318,] 1.38216581 -0.28821695
[319,] 9.66207012 1.38216581
[320,] 12.39216581 9.66207012
[321,] 14.13197443 12.39216581
[322,] -2.95802557 14.13197443
[323,] 4.95178305 -2.95802557
[324,] 4.46207012 4.95178305
[325,] -1.09821695 4.46207012
[326,] 7.50187874 -1.09821695
[327,] -1.89821695 7.50187874
[328,] -2.38812126 -1.89821695
[329,] 5.32178305 -2.38812126
[330,] 1.79207012 5.32178305
[331,] -11.51980052 1.79207012
[332,] 3.18187874 -11.51980052
[333,] 4.40178305 3.18187874
[334,] 5.86187874 4.40178305
[335,] -2.22792988 5.86187874
[336,] -1.50821695 -2.22792988
[337,] -3.78821695 -1.50821695
[338,] -2.47802557 -3.78821695
[339,] -4.24821695 -2.47802557
[340,] -7.33821695 -4.24821695
[341,] -0.24821695 -7.33821695
[342,] 0.53187874 -0.24821695
[343,] 11.03187874 0.53187874
[344,] -5.17812126 11.03187874
[345,] 5.80178305 -5.17812126
[346,] -1.41812126 5.80178305
[347,] 1.07187874 -1.41812126
[348,] 9.35207012 1.07187874
[349,] 5.36207012 9.35207012
[350,] 6.79187874 5.36207012
[351,] 3.32187874 6.79187874
[352,] -1.66783419 3.32187874
[353,] 1.54187874 -1.66783419
[354,] 3.36178305 1.54187874
[355,] 8.18178305 3.36178305
[356,] 1.93178305 8.18178305
[357,] 0.40019948 1.93178305
[358,] -5.50792988 0.40019948
[359,] -4.20821695 -5.50792988
[360,] 10.80216581 -4.20821695
[361,] -3.18821695 10.80216581
[362,] 2.04187874 -3.18821695
[363,] 6.24207012 2.04187874
[364,] 12.56197443 6.24207012
[365,] 0.53010379 12.56197443
[366,] -0.78821695 0.53010379
[367,] -1.76792988 -0.78821695
[368,] -9.95008759 -1.76792988
[369,] 5.69019948 -9.95008759
[370,] -5.57812126 5.69019948
[371,] 16.34010379 -5.57812126
[372,] 6.11187874 16.34010379
[373,] -1.49812126 6.11187874
[374,] 0.49178305 -1.49812126
[375,] -4.40812126 0.49178305
[376,] -2.70812126 -4.40812126
[377,] 3.77178305 -2.70812126
[378,] 2.07187874 3.77178305
[379,] 3.75207012 2.07187874
[380,] -5.37999190 3.75207012
[381,] -5.76821695 -5.37999190
[382,] 5.12178305 -5.76821695
[383,] -3.13980052 5.12178305
[384,] 10.73187874 -3.13980052
[385,] 6.83178305 10.73187874
[386,] 5.07187874 6.83178305
[387,] -0.69018328 5.07187874
[388,] 2.67187874 -0.69018328
[389,] 5.17207012 2.67187874
[390,] -7.23792988 5.17207012
[391,] -1.66821695 -7.23792988
[392,] 6.83187874 -1.66821695
[393,] 9.41187874 6.83187874
[394,] -1.65802557 9.41187874
[395,] 2.20178305 -1.65802557
[396,] -8.47821695 2.20178305
[397,] -4.86812126 -8.47821695
[398,] -4.95821695 -4.86812126
[399,] -1.68812126 -4.95821695
[400,] 6.70216581 -1.68812126
[401,] -3.93792988 6.70216581
[402,] 9.13000810 -3.93792988
[403,] 11.61187874 9.13000810
[404,] -5.85812126 11.61187874
[405,] 6.92187874 -5.85812126
[406,] -2.07783419 6.92187874
[407,] -3.85802557 -2.07783419
[408,] 0.02178305 -3.85802557
[409,] -5.80802557 0.02178305
[410,] -3.61812126 -5.80802557
[411,] -0.62821695 -3.61812126
[412,] -5.49812126 -0.62821695
[413,] -5.55989621 -5.49812126
[414,] -12.12821695 -5.55989621
[415,] 11.24178305 -12.12821695
[416,] -4.12812126 11.24178305
[417,] -3.09792988 -4.12812126
[418,] -1.25821695 -3.09792988
[419,] 0.89187874 -1.25821695
[420,] 6.34178305 0.89187874
[421,] 16.01010379 6.34178305
[422,] 8.70187874 16.01010379
[423,] 1.84178305 8.70187874
[424,] 0.84187874 1.84178305
[425,] -4.05812126 0.84187874
[426,] -4.33783419 -4.05812126
[427,] -0.16812126 -4.33783419
[428,] -6.65792988 -0.16812126
[429,] 6.11187874 -6.65792988
[430,] 5.21187874 6.11187874
[431,] -3.83821695 5.21187874
[432,] -9.35792988 -3.83821695
[433,] -3.05821695 -9.35792988
[434,] -0.68989621 -3.05821695
[435,] 0.25178305 -0.68989621
[436,] -2.68802557 0.25178305
[437,] 1.42178305 -2.68802557
[438,] 9.21019948 1.42178305
[439,] 1.18197443 9.21019948
[440,] 9.56981672 1.18197443
[441,] 16.89187874 9.56981672
[442,] 5.08216581 16.89187874
[443,] -12.20802557 5.08216581
[444,] -8.12821695 -12.20802557
[445,] 5.18981672 -8.12821695
[446,] 6.79178305 5.18981672
[447,] 2.30178305 6.79178305
[448,] -8.95802557 2.30178305
[449,] 8.55207012 -8.95802557
[450,] 3.19187874 8.55207012
[451,] 20.87197443 3.19187874
[452,] 11.69187874 20.87197443
[453,] 9.43197443 11.69187874
[454,] -2.97812126 9.43197443
[455,] 4.79991241 -2.97812126
[456,] 12.74197443 4.79991241
[457,] 6.54197443 12.74197443
[458,] -2.80792988 6.54197443
[459,] -4.38812126 -2.80792988
[460,] 4.21000810 -4.38812126
[461,] 5.29178305 4.21000810
[462,] -2.23802557 5.29178305
[463,] 15.07207012 -2.23802557
[464,] -3.12821695 15.07207012
[465,] 1.42197443 -3.12821695
[466,] -0.89821695 1.42197443
[467,] -7.89802557 -0.89821695
[468,] 0.20178305 -7.89802557
[469,] -0.76812126 0.20178305
[470,] -7.97821695 -0.76812126
[471,] 8.29000810 -7.97821695
[472,] -10.64783419 8.29000810
[473,] 0.87178305 -10.64783419
[474,] 3.79178305 0.87178305
[475,] -0.57802557 3.79178305
[476,] -4.69821695 -0.57802557
[477,] -8.12821695 -4.69821695
[478,] -9.38989621 -8.12821695
[479,] 2.66197443 -9.38989621
[480,] 3.64197443 2.66197443
[481,] -3.25812126 3.64197443
[482,] -0.76989621 -3.25812126
[483,] 2.64000810 -0.76989621
[484,] -10.05821695 2.64000810
[485,] -8.22812126 -10.05821695
[486,] 0.70187874 -8.22812126
[487,] 1.42010379 0.70187874
[488,] 9.36197443 1.42010379
[489,] 10.78207012 9.36197443
[490,] 6.14207012 10.78207012
[491,] -4.20821695 6.14207012
[492,] -3.84802557 -4.20821695
[493,] -6.17999190 -3.84802557
[494,] 0.77216581 -6.17999190
[495,] 4.12991241 0.77216581
[496,] 1.91178305 4.12991241
[497,] -4.77989621 1.91178305
[498,] 1.19187874 -4.77989621
[499,] 7.77187874 1.19187874
[500,] -1.14792988 7.77187874
[501,] 5.55197443 -1.14792988
[502,] -3.14821695 5.55197443
[503,] 19.20178305 -3.14821695
[504,] -0.54802557 19.20178305
[505,] 1.42010379 -0.54802557
[506,] 3.54187874 1.42010379
[507,] -4.54821695 3.54187874
[508,] 3.43187874 -4.54821695
[509,] 14.20187874 3.43187874
[510,] -5.41812126 14.20187874
[511,] -3.40812126 -5.41812126
[512,] 3.02187874 -3.40812126
[513,] -0.17783419 3.02187874
[514,] 15.83197443 -0.17783419
[515,] -0.38792988 15.83197443
[516,] 17.13207012 -0.38792988
[517,] 10.41981672 17.13207012
[518,] -7.19802557 10.41981672
[519,] -3.06821695 -7.19802557
[520,] -3.29821695 -3.06821695
[521,] -1.62999190 -3.29821695
[522,] 2.63991241 -1.62999190
[523,] 0.94197443 2.63991241
[524,] 1.74187874 0.94197443
[525,] 4.65187874 1.74187874
[526,] 3.92197443 4.65187874
[527,] -2.70802557 3.92197443
[528,] -7.11783419 -2.70802557
[529,] 1.18187874 -7.11783419
[530,] 4.56178305 1.18187874
[531,] -0.37812126 4.56178305
[532,] -5.31802557 -0.37812126
[533,] -3.13792988 -5.31802557
[534,] -5.67802557 -3.13792988
[535,] -4.02812126 -5.67802557
[536,] 4.16187874 -4.02812126
[537,] -7.80812126 4.16187874
[538,] 1.77216581 -7.80812126
[539,] 9.43187874 1.77216581
[540,] 6.30010379 9.43187874
[541,] 7.87197443 6.30010379
[542,] -5.38821695 7.87197443
[543,] -3.84802557 -5.38821695
[544,] -9.20792988 -3.84802557
[545,] 7.85178305 -9.20792988
[546,] 4.02178305 7.85178305
[547,] 0.17178305 4.02178305
[548,] 6.58178305 0.17178305
[549,] 6.68216581 6.58178305
[550,] 3.93197443 6.68216581
[551,] 0.24187874 3.93197443
[552,] -2.14812126 0.24187874
[553,] -8.04802557 -2.14812126
[554,] -6.68783419 -8.04802557
[555,] -11.11980052 -6.68783419
[556,] 4.39178305 -11.11980052
[557,] -0.09812126 4.39178305
[558,] 19.47991241 -0.09812126
[559,] -2.35821695 19.47991241
[560,] 1.67178305 -2.35821695
[561,] -1.01802557 1.67178305
[562,] 16.94197443 -1.01802557
[563,] -0.91812126 16.94197443
[564,] -4.22802557 -0.91812126
[565,] 4.15178305 -4.22802557
[566,] 3.82197443 4.15178305
[567,] 8.62207012 3.82197443
[568,] 0.99178305 8.62207012
[569,] 2.04197443 0.99178305
[570,] 4.19178305 2.04197443
[571,] -8.38792988 4.19178305
[572,] -3.14802557 -8.38792988
[573,] 5.81981672 -3.14802557
[574,] -2.01802557 5.81981672
[575,] -1.02821695 -2.01802557
[576,] -5.16812126 -1.02821695
[577,] -6.37802557 -5.16812126
[578,] 12.26187874 -6.37802557
[579,] -1.35802557 12.26187874
[580,] 4.61178305 -1.35802557
[581,] 2.07207012 4.61178305
[582,] 0.85197443 2.07207012
[583,] 3.76187874 0.85197443
[584,] -2.44812126 3.76187874
[585,] -0.81812126 -2.44812126
[586,] -7.63812126 -0.81812126
[587,] 10.66197443 -7.63812126
[588,] -6.05783419 10.66197443
[589,] 4.30991241 -6.05783419
[590,] -3.41812126 4.30991241
[591,] -0.54802557 -3.41812126
[592,] 2.03216581 -0.54802557
[593,] 6.49991241 2.03216581
[594,] -0.02989621 6.49991241
[595,] -4.45812126 -0.02989621
[596,] -2.79821695 -4.45812126
[597,] -2.28792988 -2.79821695
[598,] -5.57802557 -2.28792988
[599,] 4.14991241 -5.57802557
[600,] 8.90010379 4.14991241
[601,] 13.75019948 8.90010379
[602,] -5.70821695 13.75019948
[603,] 3.79991241 -5.70821695
[604,] -1.18812126 3.79991241
[605,] 4.48226150 -1.18812126
[606,] 8.85029517 4.48226150
[607,] -16.28745143 8.85029517
[608,] 3.34264426 -16.28745143
[609,] -10.33932207 3.34264426
[610,] -1.65960914 -10.33932207
[611,] 1.27245288 -1.65960914
[612,] 4.43226150 1.27245288
[613,] 6.64226150 4.43226150
[614,] -12.00773850 6.64226150
[615,] -9.41970483 -12.00773850
[616,] -14.22970483 -9.41970483
[617,] -12.07932207 -14.22970483
[618,] -14.52745143 -12.07932207
[619,] -6.81764281 -14.52745143
[620,] -10.04951345 -6.81764281
[621,] -9.50960914 -10.04951345
[622,] -6.16970483 -9.50960914
[623,] 19.57226150 -6.16970483
[624,] 1.49226150 19.57226150
[625,] -3.14745143 1.49226150
[626,] 4.22245288 -3.14745143
[627,] -1.04773850 4.22245288
[628,] -6.80951345 -1.04773850
[629,] -8.01754712 -6.80951345
[630,] -3.94764281 -8.01754712
[631,] 13.31245288 -3.94764281
[632,] -4.98970483 13.31245288
[633,] -2.90970483 -4.98970483
[634,] -12.45960914 -2.90970483
[635,] -4.23764281 -12.45960914
[636,] 1.78029517 -4.23764281
[637,] -2.33941776 1.78029517
[638,] -9.48773850 -2.33941776
[639,] 2.16039086 -9.48773850
[640,] -3.41754712 2.16039086
[641,] 2.30235719 -3.41754712
[642,] -10.87764281 2.30235719
[643,] -13.56970483 -10.87764281
[644,] 0.84058224 -13.56970483
[645,] -13.42970483 0.84058224
[646,] 4.35235719 -13.42970483
[647,] 13.36245288 4.35235719
[648,] 10.18039086 13.36245288
[649,] 4.16039086 10.18039086
[650,] 1.06029517 4.16039086
[651,] -16.81745143 1.06029517
[652,] -10.58960914 -16.81745143
[653,] -6.24773850 -10.58960914
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -4.80840833 -7.29831264
2 -2.51850401 -4.80840833
3 -8.79831264 -2.51850401
4 -5.42831264 -8.79831264
5 1.31159167 -5.42831264
6 1.80140030 1.31159167
7 -3.23859970 1.80140030
8 -2.17840833 -3.23859970
9 -4.95831264 -2.17840833
10 -1.36859970 -4.95831264
11 -4.69840833 -1.36859970
12 -0.11840833 -4.69840833
13 -0.86840833 -0.11840833
14 0.53159167 -0.86840833
15 -0.39850401 0.53159167
16 -0.33869539 -0.39850401
17 1.39140030 -0.33869539
18 -0.85831264 1.39140030
19 1.66168736 -0.85831264
20 -1.05869539 1.66168736
21 -2.49869539 -1.05869539
22 -11.88979108 -2.49869539
23 6.06149599 -11.88979108
24 5.50168736 6.06149599
25 3.64111323 5.50168736
26 -0.89831264 3.64111323
27 2.49130461 -0.89831264
28 7.75159167 2.49130461
29 -3.37831264 7.75159167
30 -2.23869539 -3.37831264
31 5.62168736 -2.23869539
32 4.68159167 5.62168736
33 1.21149599 4.68159167
34 1.06130461 1.21149599
35 -0.29840833 1.06130461
36 2.87168736 -0.29840833
37 -1.33840833 2.87168736
38 -6.57840833 -1.33840833
39 -2.00840833 -6.57840833
40 6.97168736 -2.00840833
41 2.15159167 6.97168736
42 2.70130461 2.15159167
43 -2.73840833 2.70130461
44 -1.62869539 -2.73840833
45 -3.61831264 -1.62869539
46 -3.47850401 -3.61831264
47 -8.60840833 -3.47850401
48 -0.88859970 -8.60840833
49 5.27159167 -0.88859970
50 5.11130461 5.27159167
51 3.59168736 5.11130461
52 11.18168736 3.59168736
53 -5.01840833 11.18168736
54 -1.04859970 -5.01840833
55 1.32168736 -1.04859970
56 5.78168736 1.32168736
57 4.47149599 5.78168736
58 2.96120892 4.47149599
59 -5.06840833 2.96120892
60 -0.81840833 -5.06840833
61 2.76159167 -0.81840833
62 -2.57859970 2.76159167
63 3.04120892 -2.57859970
64 -2.64840833 3.04120892
65 0.08140030 -2.64840833
66 -3.68831264 0.08140030
67 -3.40850401 -3.68831264
68 0.30130461 -3.40850401
69 1.22168736 0.30130461
70 -2.42840833 1.22168736
71 3.26159167 -2.42840833
72 2.77168736 3.26159167
73 0.19168736 2.77168736
74 -3.47831264 0.19168736
75 -1.82850401 -3.47831264
76 3.52130461 -1.82850401
77 -0.53869539 3.52130461
78 14.04168736 -0.53869539
79 -0.19859970 14.04168736
80 1.39149599 -0.19859970
81 -0.43859970 1.39149599
82 0.06159167 -0.43859970
83 -4.10840833 0.06159167
84 -0.54850401 -4.10840833
85 -0.60840833 -0.54850401
86 -7.18850401 -0.60840833
87 2.21168736 -7.18850401
88 0.74130461 2.21168736
89 -3.11831264 0.74130461
90 5.91168736 -3.11831264
91 5.26168736 5.91168736
92 -3.60850401 5.26168736
93 0.10159167 -3.60850401
94 1.83159167 0.10159167
95 -2.41831264 1.83159167
96 -6.86831264 -2.41831264
97 -2.72831264 -6.86831264
98 -2.89850401 -2.72831264
99 -6.81840833 -2.89850401
100 -9.12840833 -6.81840833
101 -3.22850401 -9.12840833
102 1.50168736 -3.22850401
103 -11.93279108 1.50168736
104 1.36168736 -11.93279108
105 2.40140030 1.36168736
106 1.49159167 2.40140030
107 -0.20859970 1.49159167
108 -2.29850401 -0.20859970
109 -3.68850401 -2.29850401
110 4.34159167 -3.68850401
111 -1.42850401 4.34159167
112 1.12149599 -1.42850401
113 2.48149599 1.12149599
114 -4.98850401 2.48149599
115 1.36159167 -4.98850401
116 -2.73850401 1.36159167
117 -3.09869539 -2.73850401
118 -3.32840833 -3.09869539
119 2.47149599 -3.32840833
120 -0.17831264 2.47149599
121 -1.44850401 -0.17831264
122 -5.10850401 -1.44850401
123 -4.00859970 -5.10850401
124 -1.14840833 -4.00859970
125 -5.07840833 -1.14840833
126 -6.42840833 -5.07840833
127 -3.97831264 -6.42840833
128 -3.62850401 -3.97831264
129 2.53159167 -3.62850401
130 2.12168736 2.53159167
131 -7.75840833 2.12168736
132 -5.29840833 -7.75840833
133 -4.90831264 -5.29840833
134 -1.35840833 -4.90831264
135 -1.66859970 -1.35840833
136 -3.90859970 -1.66859970
137 0.83159167 -3.90859970
138 -6.64831264 0.83159167
139 -14.26769539 -6.64831264
140 -0.66850401 -14.26769539
141 0.25168736 -0.66850401
142 -3.07859970 0.25168736
143 1.93168736 -3.07859970
144 6.76168736 1.93168736
145 -3.02831264 6.76168736
146 -2.11859970 -3.02831264
147 0.88159167 -2.11859970
148 6.04168736 0.88159167
149 7.22159167 6.04168736
150 4.60159167 7.22159167
151 -1.36831264 4.60159167
152 0.22168736 -1.36831264
153 1.54168736 0.22168736
154 -2.21850401 1.54168736
155 -4.45840833 -2.21850401
156 -2.02859970 -4.45840833
157 -1.78831264 -2.02859970
158 -3.89831264 -1.78831264
159 1.33168736 -3.89831264
160 0.74149599 1.33168736
161 -4.24840833 0.74149599
162 0.46130461 -4.24840833
163 -2.51840833 0.46130461
164 4.55168736 -2.51840833
165 -0.25859970 4.55168736
166 4.13168736 -0.25859970
167 -1.14859970 4.13168736
168 4.26159167 -1.14859970
169 3.13130461 4.26159167
170 -3.13850401 3.13130461
171 1.86149599 -3.13850401
172 5.16120892 1.86149599
173 5.66168736 5.16120892
174 2.98159167 5.66168736
175 -5.04831264 2.98159167
176 -3.46831264 -5.04831264
177 -1.21831264 -3.46831264
178 1.98130461 -1.21831264
179 -8.31831264 1.98130461
180 -8.06850401 -8.31831264
181 3.63140030 -8.06850401
182 -1.17831264 3.63140030
183 -2.07831264 -1.17831264
184 -7.21831264 -2.07831264
185 -1.81850401 -7.21831264
186 -3.59869539 -1.81850401
187 -0.17840833 -3.59869539
188 2.79168736 -0.17840833
189 -1.75850401 2.79168736
190 -2.70027897 -1.75850401
191 -8.69840833 -2.70027897
192 -3.18831264 -8.69840833
193 -0.72859970 -3.18831264
194 0.87140030 -0.72859970
195 5.04159167 0.87140030
196 4.33168736 5.04159167
197 -4.13869539 4.33168736
198 5.23140030 -4.13869539
199 3.65140030 5.23140030
200 -2.71831264 3.65140030
201 -2.81850401 -2.71831264
202 5.35168736 -2.81850401
203 -0.59840833 5.35168736
204 4.65130461 -0.59840833
205 1.23149599 4.65130461
206 -0.41859970 1.23149599
207 0.17168736 -0.41859970
208 -0.51850401 0.17168736
209 -2.73831264 -0.51850401
210 -3.07831264 -2.73831264
211 7.88159167 -3.07831264
212 0.91168736 7.88159167
213 -2.45850401 0.91168736
214 0.04168736 -2.45850401
215 -1.70879108 0.04168736
216 -3.81831264 -1.70879108
217 3.02130461 -3.81831264
218 1.00159167 3.02130461
219 -4.68831264 1.00159167
220 -6.48840833 -4.68831264
221 0.32111323 -6.48840833
222 -3.95831264 0.32111323
223 -6.92840833 -3.95831264
224 -3.15859970 -6.92840833
225 12.43168736 -3.15859970
226 -2.13840833 12.43168736
227 -3.07840833 -2.13840833
228 -6.01850401 -3.07840833
229 -4.00859970 -6.01850401
230 -1.88840833 -4.00859970
231 2.01168736 -1.88840833
232 1.16120892 2.01168736
233 -3.59850401 1.16120892
234 -0.69840833 -3.59850401
235 4.76159167 -0.69840833
236 7.85168736 4.76159167
237 0.57140030 7.85168736
238 -1.94840833 0.57140030
239 -2.15859970 -1.94840833
240 -4.60840833 -2.15859970
241 0.47168736 -4.60840833
242 1.30159167 0.47168736
243 -5.56850401 1.30159167
244 -3.61831264 -5.56850401
245 2.30159167 -3.61831264
246 -2.43850401 2.30159167
247 4.12168736 -2.43850401
248 -3.60840833 4.12168736
249 -6.83831264 -3.60840833
250 -3.95859970 -6.83831264
251 0.98159167 -3.95859970
252 -1.91831264 0.98159167
253 2.71159167 -1.91831264
254 8.01168736 2.71159167
255 0.19168736 8.01168736
256 1.77159167 0.19168736
257 -3.31850401 1.77159167
258 0.83130461 -3.31850401
259 0.84149599 0.83130461
260 0.21159167 0.84149599
261 -5.51831264 0.21159167
262 3.95168736 -5.51831264
263 -0.23859970 3.95168736
264 4.26159167 -0.23859970
265 5.78149599 4.26159167
266 -5.25831264 5.78149599
267 -0.94850401 -5.25831264
268 -0.36850401 -0.94850401
269 0.33130461 -0.36850401
270 3.65168736 0.33130461
271 4.85168736 3.65168736
272 1.53159167 4.85168736
273 -1.10840833 1.53159167
274 -5.82831264 -1.10840833
275 -2.03859970 -5.82831264
276 1.63149599 -2.03859970
277 4.24130461 1.63149599
278 -2.55831264 4.24130461
279 -1.46869539 -2.55831264
280 0.30149599 -1.46869539
281 -0.39859970 0.30149599
282 2.95159167 -0.39859970
283 3.94149599 2.95159167
284 -1.35840833 3.94149599
285 1.45120892 -1.35840833
286 1.28159167 1.45120892
287 0.08130461 1.28159167
288 -4.46840833 0.08130461
289 5.63149599 -4.46840833
290 5.92149599 5.63149599
291 0.49168736 5.92149599
292 -8.46831264 0.49168736
293 -5.80840833 -8.46831264
294 3.60168736 -5.80840833
295 -0.47859970 3.60168736
296 -2.05840833 -0.47859970
297 -5.68831264 -2.05840833
298 -2.68879108 -5.68831264
299 -3.11859970 -2.68879108
300 -1.45850401 -3.11859970
301 2.50168736 -1.45850401
302 -5.47840833 2.50168736
303 -0.66859970 -5.47840833
304 -1.89840833 -0.66859970
305 3.99149599 -1.89840833
306 6.09130461 3.99149599
307 1.23159167 6.09130461
308 -4.76850401 1.23159167
309 -0.19812126 -4.76850401
310 -3.42821695 -0.19812126
311 -2.21783419 -3.42821695
312 -7.33812126 -2.21783419
313 5.66187874 -7.33812126
314 -0.78802557 5.66187874
315 -8.59821695 -0.78802557
316 -3.79812126 -8.59821695
317 -0.28821695 -3.79812126
318 1.38216581 -0.28821695
319 9.66207012 1.38216581
320 12.39216581 9.66207012
321 14.13197443 12.39216581
322 -2.95802557 14.13197443
323 4.95178305 -2.95802557
324 4.46207012 4.95178305
325 -1.09821695 4.46207012
326 7.50187874 -1.09821695
327 -1.89821695 7.50187874
328 -2.38812126 -1.89821695
329 5.32178305 -2.38812126
330 1.79207012 5.32178305
331 -11.51980052 1.79207012
332 3.18187874 -11.51980052
333 4.40178305 3.18187874
334 5.86187874 4.40178305
335 -2.22792988 5.86187874
336 -1.50821695 -2.22792988
337 -3.78821695 -1.50821695
338 -2.47802557 -3.78821695
339 -4.24821695 -2.47802557
340 -7.33821695 -4.24821695
341 -0.24821695 -7.33821695
342 0.53187874 -0.24821695
343 11.03187874 0.53187874
344 -5.17812126 11.03187874
345 5.80178305 -5.17812126
346 -1.41812126 5.80178305
347 1.07187874 -1.41812126
348 9.35207012 1.07187874
349 5.36207012 9.35207012
350 6.79187874 5.36207012
351 3.32187874 6.79187874
352 -1.66783419 3.32187874
353 1.54187874 -1.66783419
354 3.36178305 1.54187874
355 8.18178305 3.36178305
356 1.93178305 8.18178305
357 0.40019948 1.93178305
358 -5.50792988 0.40019948
359 -4.20821695 -5.50792988
360 10.80216581 -4.20821695
361 -3.18821695 10.80216581
362 2.04187874 -3.18821695
363 6.24207012 2.04187874
364 12.56197443 6.24207012
365 0.53010379 12.56197443
366 -0.78821695 0.53010379
367 -1.76792988 -0.78821695
368 -9.95008759 -1.76792988
369 5.69019948 -9.95008759
370 -5.57812126 5.69019948
371 16.34010379 -5.57812126
372 6.11187874 16.34010379
373 -1.49812126 6.11187874
374 0.49178305 -1.49812126
375 -4.40812126 0.49178305
376 -2.70812126 -4.40812126
377 3.77178305 -2.70812126
378 2.07187874 3.77178305
379 3.75207012 2.07187874
380 -5.37999190 3.75207012
381 -5.76821695 -5.37999190
382 5.12178305 -5.76821695
383 -3.13980052 5.12178305
384 10.73187874 -3.13980052
385 6.83178305 10.73187874
386 5.07187874 6.83178305
387 -0.69018328 5.07187874
388 2.67187874 -0.69018328
389 5.17207012 2.67187874
390 -7.23792988 5.17207012
391 -1.66821695 -7.23792988
392 6.83187874 -1.66821695
393 9.41187874 6.83187874
394 -1.65802557 9.41187874
395 2.20178305 -1.65802557
396 -8.47821695 2.20178305
397 -4.86812126 -8.47821695
398 -4.95821695 -4.86812126
399 -1.68812126 -4.95821695
400 6.70216581 -1.68812126
401 -3.93792988 6.70216581
402 9.13000810 -3.93792988
403 11.61187874 9.13000810
404 -5.85812126 11.61187874
405 6.92187874 -5.85812126
406 -2.07783419 6.92187874
407 -3.85802557 -2.07783419
408 0.02178305 -3.85802557
409 -5.80802557 0.02178305
410 -3.61812126 -5.80802557
411 -0.62821695 -3.61812126
412 -5.49812126 -0.62821695
413 -5.55989621 -5.49812126
414 -12.12821695 -5.55989621
415 11.24178305 -12.12821695
416 -4.12812126 11.24178305
417 -3.09792988 -4.12812126
418 -1.25821695 -3.09792988
419 0.89187874 -1.25821695
420 6.34178305 0.89187874
421 16.01010379 6.34178305
422 8.70187874 16.01010379
423 1.84178305 8.70187874
424 0.84187874 1.84178305
425 -4.05812126 0.84187874
426 -4.33783419 -4.05812126
427 -0.16812126 -4.33783419
428 -6.65792988 -0.16812126
429 6.11187874 -6.65792988
430 5.21187874 6.11187874
431 -3.83821695 5.21187874
432 -9.35792988 -3.83821695
433 -3.05821695 -9.35792988
434 -0.68989621 -3.05821695
435 0.25178305 -0.68989621
436 -2.68802557 0.25178305
437 1.42178305 -2.68802557
438 9.21019948 1.42178305
439 1.18197443 9.21019948
440 9.56981672 1.18197443
441 16.89187874 9.56981672
442 5.08216581 16.89187874
443 -12.20802557 5.08216581
444 -8.12821695 -12.20802557
445 5.18981672 -8.12821695
446 6.79178305 5.18981672
447 2.30178305 6.79178305
448 -8.95802557 2.30178305
449 8.55207012 -8.95802557
450 3.19187874 8.55207012
451 20.87197443 3.19187874
452 11.69187874 20.87197443
453 9.43197443 11.69187874
454 -2.97812126 9.43197443
455 4.79991241 -2.97812126
456 12.74197443 4.79991241
457 6.54197443 12.74197443
458 -2.80792988 6.54197443
459 -4.38812126 -2.80792988
460 4.21000810 -4.38812126
461 5.29178305 4.21000810
462 -2.23802557 5.29178305
463 15.07207012 -2.23802557
464 -3.12821695 15.07207012
465 1.42197443 -3.12821695
466 -0.89821695 1.42197443
467 -7.89802557 -0.89821695
468 0.20178305 -7.89802557
469 -0.76812126 0.20178305
470 -7.97821695 -0.76812126
471 8.29000810 -7.97821695
472 -10.64783419 8.29000810
473 0.87178305 -10.64783419
474 3.79178305 0.87178305
475 -0.57802557 3.79178305
476 -4.69821695 -0.57802557
477 -8.12821695 -4.69821695
478 -9.38989621 -8.12821695
479 2.66197443 -9.38989621
480 3.64197443 2.66197443
481 -3.25812126 3.64197443
482 -0.76989621 -3.25812126
483 2.64000810 -0.76989621
484 -10.05821695 2.64000810
485 -8.22812126 -10.05821695
486 0.70187874 -8.22812126
487 1.42010379 0.70187874
488 9.36197443 1.42010379
489 10.78207012 9.36197443
490 6.14207012 10.78207012
491 -4.20821695 6.14207012
492 -3.84802557 -4.20821695
493 -6.17999190 -3.84802557
494 0.77216581 -6.17999190
495 4.12991241 0.77216581
496 1.91178305 4.12991241
497 -4.77989621 1.91178305
498 1.19187874 -4.77989621
499 7.77187874 1.19187874
500 -1.14792988 7.77187874
501 5.55197443 -1.14792988
502 -3.14821695 5.55197443
503 19.20178305 -3.14821695
504 -0.54802557 19.20178305
505 1.42010379 -0.54802557
506 3.54187874 1.42010379
507 -4.54821695 3.54187874
508 3.43187874 -4.54821695
509 14.20187874 3.43187874
510 -5.41812126 14.20187874
511 -3.40812126 -5.41812126
512 3.02187874 -3.40812126
513 -0.17783419 3.02187874
514 15.83197443 -0.17783419
515 -0.38792988 15.83197443
516 17.13207012 -0.38792988
517 10.41981672 17.13207012
518 -7.19802557 10.41981672
519 -3.06821695 -7.19802557
520 -3.29821695 -3.06821695
521 -1.62999190 -3.29821695
522 2.63991241 -1.62999190
523 0.94197443 2.63991241
524 1.74187874 0.94197443
525 4.65187874 1.74187874
526 3.92197443 4.65187874
527 -2.70802557 3.92197443
528 -7.11783419 -2.70802557
529 1.18187874 -7.11783419
530 4.56178305 1.18187874
531 -0.37812126 4.56178305
532 -5.31802557 -0.37812126
533 -3.13792988 -5.31802557
534 -5.67802557 -3.13792988
535 -4.02812126 -5.67802557
536 4.16187874 -4.02812126
537 -7.80812126 4.16187874
538 1.77216581 -7.80812126
539 9.43187874 1.77216581
540 6.30010379 9.43187874
541 7.87197443 6.30010379
542 -5.38821695 7.87197443
543 -3.84802557 -5.38821695
544 -9.20792988 -3.84802557
545 7.85178305 -9.20792988
546 4.02178305 7.85178305
547 0.17178305 4.02178305
548 6.58178305 0.17178305
549 6.68216581 6.58178305
550 3.93197443 6.68216581
551 0.24187874 3.93197443
552 -2.14812126 0.24187874
553 -8.04802557 -2.14812126
554 -6.68783419 -8.04802557
555 -11.11980052 -6.68783419
556 4.39178305 -11.11980052
557 -0.09812126 4.39178305
558 19.47991241 -0.09812126
559 -2.35821695 19.47991241
560 1.67178305 -2.35821695
561 -1.01802557 1.67178305
562 16.94197443 -1.01802557
563 -0.91812126 16.94197443
564 -4.22802557 -0.91812126
565 4.15178305 -4.22802557
566 3.82197443 4.15178305
567 8.62207012 3.82197443
568 0.99178305 8.62207012
569 2.04197443 0.99178305
570 4.19178305 2.04197443
571 -8.38792988 4.19178305
572 -3.14802557 -8.38792988
573 5.81981672 -3.14802557
574 -2.01802557 5.81981672
575 -1.02821695 -2.01802557
576 -5.16812126 -1.02821695
577 -6.37802557 -5.16812126
578 12.26187874 -6.37802557
579 -1.35802557 12.26187874
580 4.61178305 -1.35802557
581 2.07207012 4.61178305
582 0.85197443 2.07207012
583 3.76187874 0.85197443
584 -2.44812126 3.76187874
585 -0.81812126 -2.44812126
586 -7.63812126 -0.81812126
587 10.66197443 -7.63812126
588 -6.05783419 10.66197443
589 4.30991241 -6.05783419
590 -3.41812126 4.30991241
591 -0.54802557 -3.41812126
592 2.03216581 -0.54802557
593 6.49991241 2.03216581
594 -0.02989621 6.49991241
595 -4.45812126 -0.02989621
596 -2.79821695 -4.45812126
597 -2.28792988 -2.79821695
598 -5.57802557 -2.28792988
599 4.14991241 -5.57802557
600 8.90010379 4.14991241
601 13.75019948 8.90010379
602 -5.70821695 13.75019948
603 3.79991241 -5.70821695
604 -1.18812126 3.79991241
605 4.48226150 -1.18812126
606 8.85029517 4.48226150
607 -16.28745143 8.85029517
608 3.34264426 -16.28745143
609 -10.33932207 3.34264426
610 -1.65960914 -10.33932207
611 1.27245288 -1.65960914
612 4.43226150 1.27245288
613 6.64226150 4.43226150
614 -12.00773850 6.64226150
615 -9.41970483 -12.00773850
616 -14.22970483 -9.41970483
617 -12.07932207 -14.22970483
618 -14.52745143 -12.07932207
619 -6.81764281 -14.52745143
620 -10.04951345 -6.81764281
621 -9.50960914 -10.04951345
622 -6.16970483 -9.50960914
623 19.57226150 -6.16970483
624 1.49226150 19.57226150
625 -3.14745143 1.49226150
626 4.22245288 -3.14745143
627 -1.04773850 4.22245288
628 -6.80951345 -1.04773850
629 -8.01754712 -6.80951345
630 -3.94764281 -8.01754712
631 13.31245288 -3.94764281
632 -4.98970483 13.31245288
633 -2.90970483 -4.98970483
634 -12.45960914 -2.90970483
635 -4.23764281 -12.45960914
636 1.78029517 -4.23764281
637 -2.33941776 1.78029517
638 -9.48773850 -2.33941776
639 2.16039086 -9.48773850
640 -3.41754712 2.16039086
641 2.30235719 -3.41754712
642 -10.87764281 2.30235719
643 -13.56970483 -10.87764281
644 0.84058224 -13.56970483
645 -13.42970483 0.84058224
646 4.35235719 -13.42970483
647 13.36245288 4.35235719
648 10.18039086 13.36245288
649 4.16039086 10.18039086
650 1.06029517 4.16039086
651 -16.81745143 1.06029517
652 -10.58960914 -16.81745143
653 -6.24773850 -10.58960914
> 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()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/7ca0b1292683540.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/8ca0b1292683540.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/freestat/rcomp/tmp/9ca0b1292683540.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/html/freestat/rcomp/tmp/10iyaq1292683540.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/rcomp/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, mysum$coefficients[i,1], 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.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/11lhqe1292683540.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/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,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/12pzpk1292683540.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, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> 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, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/13lr5a1292683540.tab")
> 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,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/1469lg1292683540.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,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/html/freestat/rcomp/tmp/159a141292683540.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,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ 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,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ 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,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ 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="/var/www/html/freestat/rcomp/tmp/16daia1292683540.tab")
+ }
>
> try(system("convert tmp/1tfvw1292683540.ps tmp/1tfvw1292683540.png",intern=TRUE))
character(0)
> try(system("convert tmp/2tfvw1292683540.ps tmp/2tfvw1292683540.png",intern=TRUE))
character(0)
> try(system("convert tmp/346cz1292683540.ps tmp/346cz1292683540.png",intern=TRUE))
character(0)
> try(system("convert tmp/446cz1292683540.ps tmp/446cz1292683540.png",intern=TRUE))
character(0)
> try(system("convert tmp/546cz1292683540.ps tmp/546cz1292683540.png",intern=TRUE))
character(0)
> try(system("convert tmp/6110q1292683540.ps tmp/6110q1292683540.png",intern=TRUE))
character(0)
> try(system("convert tmp/7ca0b1292683540.ps tmp/7ca0b1292683540.png",intern=TRUE))
character(0)
> try(system("convert tmp/8ca0b1292683540.ps tmp/8ca0b1292683540.png",intern=TRUE))
character(0)
> try(system("convert tmp/9ca0b1292683540.ps tmp/9ca0b1292683540.png",intern=TRUE))
character(0)
> try(system("convert tmp/10iyaq1292683540.ps tmp/10iyaq1292683540.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
22.757 3.589 23.614