R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
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.
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(1958
+ ,3
+ ,315.71
+ ,1958
+ ,4
+ ,317.45
+ ,1958
+ ,5
+ ,317.5
+ ,1958
+ ,6
+ ,317.12
+ ,1958
+ ,7
+ ,315.86
+ ,1958
+ ,8
+ ,314.93
+ ,1958
+ ,9
+ ,313.2
+ ,1958
+ ,10
+ ,312.6
+ ,1958
+ ,11
+ ,313.33
+ ,1958
+ ,12
+ ,314.67
+ ,1959
+ ,1
+ ,315.62
+ ,1959
+ ,2
+ ,316.38
+ ,1959
+ ,3
+ ,316.71
+ ,1959
+ ,4
+ ,317.72
+ ,1959
+ ,5
+ ,318.29
+ ,1959
+ ,6
+ ,318.16
+ ,1959
+ ,7
+ ,316.55
+ ,1959
+ ,8
+ ,314.8
+ ,1959
+ ,9
+ ,313.84
+ ,1959
+ ,10
+ ,313.26
+ ,1959
+ ,11
+ ,314.8
+ ,1959
+ ,12
+ ,315.59
+ ,1960
+ ,1
+ ,316.43
+ ,1960
+ ,2
+ ,316.97
+ ,1960
+ ,3
+ ,317.58
+ ,1960
+ ,4
+ ,319.02
+ ,1960
+ ,5
+ ,320.02
+ ,1960
+ ,6
+ ,319.59
+ ,1960
+ ,7
+ ,318.18
+ ,1960
+ ,8
+ ,315.91
+ ,1960
+ ,9
+ ,314.16
+ ,1960
+ ,10
+ ,313.83
+ ,1960
+ ,11
+ ,315
+ ,1960
+ ,12
+ ,316.19
+ ,1961
+ ,1
+ ,316.93
+ ,1961
+ ,2
+ ,317.7
+ ,1961
+ ,3
+ ,318.54
+ ,1961
+ ,4
+ ,319.48
+ ,1961
+ ,5
+ ,320.58
+ ,1961
+ ,6
+ ,319.77
+ ,1961
+ ,7
+ ,318.58
+ ,1961
+ ,8
+ ,316.79
+ ,1961
+ ,9
+ ,314.8
+ ,1961
+ ,10
+ ,315.38
+ ,1961
+ ,11
+ ,316.1
+ ,1961
+ ,12
+ ,317.01
+ ,1962
+ ,1
+ ,317.94
+ ,1962
+ ,2
+ ,318.55
+ ,1962
+ ,3
+ ,319.68
+ ,1962
+ ,4
+ ,320.63
+ ,1962
+ ,5
+ ,321.01
+ ,1962
+ ,6
+ ,320.55
+ ,1962
+ ,7
+ ,319.58
+ ,1962
+ ,8
+ ,317.4
+ ,1962
+ ,9
+ ,316.26
+ ,1962
+ ,10
+ ,315.42
+ ,1962
+ ,11
+ ,316.69
+ ,1962
+ ,12
+ ,317.7
+ ,1963
+ ,1
+ ,318.74
+ ,1963
+ ,2
+ ,319.08
+ ,1963
+ ,3
+ ,319.86
+ ,1963
+ ,4
+ ,321.39
+ ,1963
+ ,5
+ ,322.24
+ ,1963
+ ,6
+ ,321.47
+ ,1963
+ ,7
+ ,319.74
+ ,1963
+ ,8
+ ,317.77
+ ,1963
+ ,9
+ ,316.21
+ ,1963
+ ,10
+ ,315.99
+ ,1963
+ ,11
+ ,317.12
+ ,1963
+ ,12
+ ,318.31
+ ,1964
+ ,1
+ ,319.57
+ ,1964
+ ,2
+ ,320.08
+ ,1964
+ ,3
+ ,320.75
+ ,1964
+ ,4
+ ,321.8
+ ,1964
+ ,5
+ ,322.24
+ ,1964
+ ,6
+ ,321.89
+ ,1964
+ ,7
+ ,320.44
+ ,1964
+ ,8
+ ,318.7
+ ,1964
+ ,9
+ ,316.7
+ ,1964
+ ,10
+ ,316.79
+ ,1964
+ ,11
+ ,317.79
+ ,1964
+ ,12
+ ,318.71
+ ,1965
+ ,1
+ ,319.44
+ ,1965
+ ,2
+ ,320.44
+ ,1965
+ ,3
+ ,320.89
+ ,1965
+ ,4
+ ,322.13
+ ,1965
+ ,5
+ ,322.16
+ ,1965
+ ,6
+ ,321.87
+ ,1965
+ ,7
+ ,321.39
+ ,1965
+ ,8
+ ,318.8
+ ,1965
+ ,9
+ ,317.81
+ ,1965
+ ,10
+ ,317.3
+ ,1965
+ ,11
+ ,318.87
+ ,1965
+ ,12
+ ,319.42
+ ,1966
+ ,1
+ ,320.62
+ ,1966
+ ,2
+ ,321.59
+ ,1966
+ ,3
+ ,322.39
+ ,1966
+ ,4
+ ,323.87
+ ,1966
+ ,5
+ ,324.01
+ ,1966
+ ,6
+ ,323.75
+ ,1966
+ ,7
+ ,322.4
+ ,1966
+ ,8
+ ,320.37
+ ,1966
+ ,9
+ ,318.64
+ ,1966
+ ,10
+ ,318.1
+ ,1966
+ ,11
+ ,319.78
+ ,1966
+ ,12
+ ,321.08
+ ,1967
+ ,1
+ ,322.06
+ ,1967
+ ,2
+ ,322.5
+ ,1967
+ ,3
+ ,323.04
+ ,1967
+ ,4
+ ,324.42
+ ,1967
+ ,5
+ ,325
+ ,1967
+ ,6
+ ,324.09
+ ,1967
+ ,7
+ ,322.55
+ ,1967
+ ,8
+ ,320.92
+ ,1967
+ ,9
+ ,319.31
+ ,1967
+ ,10
+ ,319.31
+ ,1967
+ ,11
+ ,320.72
+ ,1967
+ ,12
+ ,321.96
+ ,1968
+ ,1
+ ,322.57
+ ,1968
+ ,2
+ ,323.15
+ ,1968
+ ,3
+ ,323.89
+ ,1968
+ ,4
+ ,325.02
+ ,1968
+ ,5
+ ,325.57
+ ,1968
+ ,6
+ ,325.36
+ ,1968
+ ,7
+ ,324.14
+ ,1968
+ ,8
+ ,322.03
+ ,1968
+ ,9
+ ,320.41
+ ,1968
+ ,10
+ ,320.25
+ ,1968
+ ,11
+ ,321.31
+ ,1968
+ ,12
+ ,322.84
+ ,1969
+ ,1
+ ,324
+ ,1969
+ ,2
+ ,324.42
+ ,1969
+ ,3
+ ,325.64
+ ,1969
+ ,4
+ ,326.66
+ ,1969
+ ,5
+ ,327.34
+ ,1969
+ ,6
+ ,326.76
+ ,1969
+ ,7
+ ,325.88
+ ,1969
+ ,8
+ ,323.67
+ ,1969
+ ,9
+ ,322.38
+ ,1969
+ ,10
+ ,321.78
+ ,1969
+ ,11
+ ,322.85
+ ,1969
+ ,12
+ ,324.12
+ ,1970
+ ,1
+ ,325.03
+ ,1970
+ ,2
+ ,325.99
+ ,1970
+ ,3
+ ,326.87
+ ,1970
+ ,4
+ ,328.14
+ ,1970
+ ,5
+ ,328.07
+ ,1970
+ ,6
+ ,327.66
+ ,1970
+ ,7
+ ,326.35
+ ,1970
+ ,8
+ ,324.69
+ ,1970
+ ,9
+ ,323.1
+ ,1970
+ ,10
+ ,323.16
+ ,1970
+ ,11
+ ,323.98
+ ,1970
+ ,12
+ ,325.13
+ ,1971
+ ,1
+ ,326.17
+ ,1971
+ ,2
+ ,326.68
+ ,1971
+ ,3
+ ,327.18
+ ,1971
+ ,4
+ ,327.78
+ ,1971
+ ,5
+ ,328.92
+ ,1971
+ ,6
+ ,328.57
+ ,1971
+ ,7
+ ,327.34
+ ,1971
+ ,8
+ ,325.46
+ ,1971
+ ,9
+ ,323.36
+ ,1971
+ ,10
+ ,323.56
+ ,1971
+ ,11
+ ,324.8
+ ,1971
+ ,12
+ ,326.01
+ ,1972
+ ,1
+ ,326.77
+ ,1972
+ ,2
+ ,327.63
+ ,1972
+ ,3
+ ,327.75
+ ,1972
+ ,4
+ ,329.72
+ ,1972
+ ,5
+ ,330.07
+ ,1972
+ ,6
+ ,329.09
+ ,1972
+ ,7
+ ,328.05
+ ,1972
+ ,8
+ ,326.32
+ ,1972
+ ,9
+ ,324.93
+ ,1972
+ ,10
+ ,325.06
+ ,1972
+ ,11
+ ,326.5
+ ,1972
+ ,12
+ ,327.55
+ ,1973
+ ,1
+ ,328.55
+ ,1973
+ ,2
+ ,329.56
+ ,1973
+ ,3
+ ,330.3
+ ,1973
+ ,4
+ ,331.5
+ ,1973
+ ,5
+ ,332.48
+ ,1973
+ ,6
+ ,332.07
+ ,1973
+ ,7
+ ,330.87
+ ,1973
+ ,8
+ ,329.31
+ ,1973
+ ,9
+ ,327.51
+ ,1973
+ ,10
+ ,327.18
+ ,1973
+ ,11
+ ,328.16
+ ,1973
+ ,12
+ ,328.64
+ ,1974
+ ,1
+ ,329.35
+ ,1974
+ ,2
+ ,330.71
+ ,1974
+ ,3
+ ,331.48
+ ,1974
+ ,4
+ ,332.65
+ ,1974
+ ,5
+ ,333.15
+ ,1974
+ ,6
+ ,332.13
+ ,1974
+ ,7
+ ,330.99
+ ,1974
+ ,8
+ ,329.17
+ ,1974
+ ,9
+ ,327.41
+ ,1974
+ ,10
+ ,327.21
+ ,1974
+ ,11
+ ,328.34
+ ,1974
+ ,12
+ ,329.5
+ ,1975
+ ,1
+ ,330.68
+ ,1975
+ ,2
+ ,331.41
+ ,1975
+ ,3
+ ,331.85
+ ,1975
+ ,4
+ ,333.29
+ ,1975
+ ,5
+ ,333.91
+ ,1975
+ ,6
+ ,333.4
+ ,1975
+ ,7
+ ,331.74
+ ,1975
+ ,8
+ ,329.88
+ ,1975
+ ,9
+ ,328.57
+ ,1975
+ ,10
+ ,328.35
+ ,1975
+ ,11
+ ,329.33
+ ,1975
+ ,12
+ ,330.58
+ ,1976
+ ,1
+ ,331.66
+ ,1976
+ ,2
+ ,332.75
+ ,1976
+ ,3
+ ,333.46
+ ,1976
+ ,4
+ ,334.78
+ ,1976
+ ,5
+ ,334.79
+ ,1976
+ ,6
+ ,334.05
+ ,1976
+ ,7
+ ,332.95
+ ,1976
+ ,8
+ ,330.64
+ ,1976
+ ,9
+ ,328.96
+ ,1976
+ ,10
+ ,328.77
+ ,1976
+ ,11
+ ,330.18
+ ,1976
+ ,12
+ ,331.65
+ ,1977
+ ,1
+ ,332.69
+ ,1977
+ ,2
+ ,333.23
+ ,1977
+ ,3
+ ,334.97
+ ,1977
+ ,4
+ ,336.03
+ ,1977
+ ,5
+ ,336.82
+ ,1977
+ ,6
+ ,336.1
+ ,1977
+ ,7
+ ,334.79
+ ,1977
+ ,8
+ ,332.53
+ ,1977
+ ,9
+ ,331.19
+ ,1977
+ ,10
+ ,331.21
+ ,1977
+ ,11
+ ,332.35
+ ,1977
+ ,12
+ ,333.47
+ ,1978
+ ,1
+ ,335.09
+ ,1978
+ ,2
+ ,335.26
+ ,1978
+ ,3
+ ,336.62
+ ,1978
+ ,4
+ ,337.77
+ ,1978
+ ,5
+ ,338
+ ,1978
+ ,6
+ ,337.98
+ ,1978
+ ,7
+ ,336.48
+ ,1978
+ ,8
+ ,334.37
+ ,1978
+ ,9
+ ,332.33
+ ,1978
+ ,10
+ ,332.4
+ ,1978
+ ,11
+ ,333.76
+ ,1978
+ ,12
+ ,334.83
+ ,1979
+ ,1
+ ,336.21
+ ,1979
+ ,2
+ ,336.64
+ ,1979
+ ,3
+ ,338.13
+ ,1979
+ ,4
+ ,338.96
+ ,1979
+ ,5
+ ,339.02
+ ,1979
+ ,6
+ ,339.2
+ ,1979
+ ,7
+ ,337.6
+ ,1979
+ ,8
+ ,335.56
+ ,1979
+ ,9
+ ,333.93
+ ,1979
+ ,10
+ ,334.12
+ ,1979
+ ,11
+ ,335.26
+ ,1979
+ ,12
+ ,336.77
+ ,1980
+ ,1
+ ,337.8
+ ,1980
+ ,2
+ ,338.28
+ ,1980
+ ,3
+ ,340.04
+ ,1980
+ ,4
+ ,340.86
+ ,1980
+ ,5
+ ,341.47
+ ,1980
+ ,6
+ ,341.26
+ ,1980
+ ,7
+ ,339.34
+ ,1980
+ ,8
+ ,337.45
+ ,1980
+ ,9
+ ,336.1
+ ,1980
+ ,10
+ ,336.05
+ ,1980
+ ,11
+ ,337.21
+ ,1980
+ ,12
+ ,338.29
+ ,1981
+ ,1
+ ,339.36
+ ,1981
+ ,2
+ ,340.51
+ ,1981
+ ,3
+ ,341.57
+ ,1981
+ ,4
+ ,342.56
+ ,1981
+ ,5
+ ,343.01
+ ,1981
+ ,6
+ ,342.52
+ ,1981
+ ,7
+ ,340.71
+ ,1981
+ ,8
+ ,338.51
+ ,1981
+ ,9
+ ,336.96
+ ,1981
+ ,10
+ ,337.13
+ ,1981
+ ,11
+ ,338.58
+ ,1981
+ ,12
+ ,339.91
+ ,1982
+ ,1
+ ,340.92
+ ,1982
+ ,2
+ ,341.69
+ ,1982
+ ,3
+ ,342.87
+ ,1982
+ ,4
+ ,343.83
+ ,1982
+ ,5
+ ,344.3
+ ,1982
+ ,6
+ ,343.42
+ ,1982
+ ,7
+ ,341.85
+ ,1982
+ ,8
+ ,339.82
+ ,1982
+ ,9
+ ,337.98
+ ,1982
+ ,10
+ ,338.09
+ ,1982
+ ,11
+ ,339.24
+ ,1982
+ ,12
+ ,340.67
+ ,1983
+ ,1
+ ,341.42
+ ,1983
+ ,2
+ ,342.67
+ ,1983
+ ,3
+ ,343.45
+ ,1983
+ ,4
+ ,345.08
+ ,1983
+ ,5
+ ,345.76
+ ,1983
+ ,6
+ ,345.32
+ ,1983
+ ,7
+ ,343.93
+ ,1983
+ ,8
+ ,342.08
+ ,1983
+ ,9
+ ,340
+ ,1983
+ ,10
+ ,340.12
+ ,1983
+ ,11
+ ,341.35
+ ,1983
+ ,12
+ ,342.89
+ ,1984
+ ,1
+ ,343.87
+ ,1984
+ ,2
+ ,344.59
+ ,1984
+ ,3
+ ,345.29
+ ,1984
+ ,4
+ ,346.59
+ ,1984
+ ,5
+ ,347.36
+ ,1984
+ ,6
+ ,346.8
+ ,1984
+ ,7
+ ,345.37
+ ,1984
+ ,8
+ ,343.06
+ ,1984
+ ,9
+ ,341.24
+ ,1984
+ ,10
+ ,341.54
+ ,1984
+ ,11
+ ,342.9
+ ,1984
+ ,12
+ ,344.36
+ ,1985
+ ,1
+ ,345.08
+ ,1985
+ ,2
+ ,345.89
+ ,1985
+ ,3
+ ,347.49
+ ,1985
+ ,4
+ ,348.02
+ ,1985
+ ,5
+ ,348.75
+ ,1985
+ ,6
+ ,348.19
+ ,1985
+ ,7
+ ,346.49
+ ,1985
+ ,8
+ ,344.7
+ ,1985
+ ,9
+ ,343.04
+ ,1985
+ ,10
+ ,342.92
+ ,1985
+ ,11
+ ,344.22
+ ,1985
+ ,12
+ ,345.61
+ ,1986
+ ,1
+ ,346.42
+ ,1986
+ ,2
+ ,346.95
+ ,1986
+ ,3
+ ,347.88
+ ,1986
+ ,4
+ ,349.57
+ ,1986
+ ,5
+ ,350.35
+ ,1986
+ ,6
+ ,349.7
+ ,1986
+ ,7
+ ,347.78
+ ,1986
+ ,8
+ ,345.89
+ ,1986
+ ,9
+ ,344.88
+ ,1986
+ ,10
+ ,344.34
+ ,1986
+ ,11
+ ,345.67
+ ,1986
+ ,12
+ ,346.89
+ ,1987
+ ,1
+ ,348.2
+ ,1987
+ ,2
+ ,348.55
+ ,1987
+ ,3
+ ,349.56
+ ,1987
+ ,4
+ ,351.12
+ ,1987
+ ,5
+ ,351.84
+ ,1987
+ ,6
+ ,351.45
+ ,1987
+ ,7
+ ,349.77
+ ,1987
+ ,8
+ ,347.62
+ ,1987
+ ,9
+ ,346.37
+ ,1987
+ ,10
+ ,346.48
+ ,1987
+ ,11
+ ,347.8
+ ,1987
+ ,12
+ ,349.03
+ ,1988
+ ,1
+ ,350.23
+ ,1988
+ ,2
+ ,351.58
+ ,1988
+ ,3
+ ,352.22
+ ,1988
+ ,4
+ ,353.53
+ ,1988
+ ,5
+ ,354.14
+ ,1988
+ ,6
+ ,353.64
+ ,1988
+ ,7
+ ,352.53
+ ,1988
+ ,8
+ ,350.42
+ ,1988
+ ,9
+ ,348.84
+ ,1988
+ ,10
+ ,348.94
+ ,1988
+ ,11
+ ,349.99
+ ,1988
+ ,12
+ ,351.29
+ ,1989
+ ,1
+ ,352.72
+ ,1989
+ ,2
+ ,353.1
+ ,1989
+ ,3
+ ,353.64
+ ,1989
+ ,4
+ ,355.43
+ ,1989
+ ,5
+ ,355.7
+ ,1989
+ ,6
+ ,355.11
+ ,1989
+ ,7
+ ,353.79
+ ,1989
+ ,8
+ ,351.42
+ ,1989
+ ,9
+ ,349.83
+ ,1989
+ ,10
+ ,350.1
+ ,1989
+ ,11
+ ,351.26
+ ,1989
+ ,12
+ ,352.66
+ ,1990
+ ,1
+ ,353.63
+ ,1990
+ ,2
+ ,354.72
+ ,1990
+ ,3
+ ,355.49
+ ,1990
+ ,4
+ ,356.1
+ ,1990
+ ,5
+ ,357.08
+ ,1990
+ ,6
+ ,356.11
+ ,1990
+ ,7
+ ,354.67
+ ,1990
+ ,8
+ ,352.67
+ ,1990
+ ,9
+ ,351.05
+ ,1990
+ ,10
+ ,351.36
+ ,1990
+ ,11
+ ,352.81
+ ,1990
+ ,12
+ ,354.21
+ ,1991
+ ,1
+ ,354.87
+ ,1991
+ ,2
+ ,355.67
+ ,1991
+ ,3
+ ,357
+ ,1991
+ ,4
+ ,358.4
+ ,1991
+ ,5
+ ,359
+ ,1991
+ ,6
+ ,357.99
+ ,1991
+ ,7
+ ,355.96
+ ,1991
+ ,8
+ ,353.78
+ ,1991
+ ,9
+ ,352.2
+ ,1991
+ ,10
+ ,352.22
+ ,1991
+ ,11
+ ,353.7
+ ,1991
+ ,12
+ ,354.98
+ ,1992
+ ,1
+ ,356.08
+ ,1992
+ ,2
+ ,356.84
+ ,1992
+ ,3
+ ,357.73
+ ,1992
+ ,4
+ ,358.91
+ ,1992
+ ,5
+ ,359.45
+ ,1992
+ ,6
+ ,359.19
+ ,1992
+ ,7
+ ,356.72
+ ,1992
+ ,8
+ ,354.77
+ ,1992
+ ,9
+ ,352.8
+ ,1992
+ ,10
+ ,353.21
+ ,1992
+ ,11
+ ,354.15
+ ,1992
+ ,12
+ ,355.39
+ ,1993
+ ,1
+ ,356.76
+ ,1993
+ ,2
+ ,357.17
+ ,1993
+ ,3
+ ,358.26
+ ,1993
+ ,4
+ ,359.17
+ ,1993
+ ,5
+ ,360.07
+ ,1993
+ ,6
+ ,359.41
+ ,1993
+ ,7
+ ,357.36
+ ,1993
+ ,8
+ ,355.29
+ ,1993
+ ,9
+ ,353.96
+ ,1993
+ ,10
+ ,354.03
+ ,1993
+ ,11
+ ,355.27
+ ,1993
+ ,12
+ ,356.7
+ ,1994
+ ,1
+ ,358.05
+ ,1994
+ ,2
+ ,358.8
+ ,1994
+ ,3
+ ,359.67
+ ,1994
+ ,4
+ ,361.13
+ ,1994
+ ,5
+ ,361.48
+ ,1994
+ ,6
+ ,360.6
+ ,1994
+ ,7
+ ,359.2
+ ,1994
+ ,8
+ ,357.23
+ ,1994
+ ,9
+ ,355.42
+ ,1994
+ ,10
+ ,355.89
+ ,1994
+ ,11
+ ,357.41
+ ,1994
+ ,12
+ ,358.74
+ ,1995
+ ,1
+ ,359.73
+ ,1995
+ ,2
+ ,360.61
+ ,1995
+ ,3
+ ,361.6
+ ,1995
+ ,4
+ ,363.05
+ ,1995
+ ,5
+ ,363.62
+ ,1995
+ ,6
+ ,363.03
+ ,1995
+ ,7
+ ,361.55
+ ,1995
+ ,8
+ ,358.94
+ ,1995
+ ,9
+ ,357.93
+ ,1995
+ ,10
+ ,357.8
+ ,1995
+ ,11
+ ,359.22
+ ,1995
+ ,12
+ ,360.42
+ ,1996
+ ,1
+ ,361.83
+ ,1996
+ ,2
+ ,362.94
+ ,1996
+ ,3
+ ,363.91
+ ,1996
+ ,4
+ ,364.28
+ ,1996
+ ,5
+ ,364.93
+ ,1996
+ ,6
+ ,364.7
+ ,1996
+ ,7
+ ,363.31
+ ,1996
+ ,8
+ ,361.15
+ ,1996
+ ,9
+ ,359.41
+ ,1996
+ ,10
+ ,359.34
+ ,1996
+ ,11
+ ,360.62
+ ,1996
+ ,12
+ ,361.96
+ ,1997
+ ,1
+ ,362.81
+ ,1997
+ ,2
+ ,363.87
+ ,1997
+ ,3
+ ,364.25
+ ,1997
+ ,4
+ ,366.02
+ ,1997
+ ,5
+ ,366.47
+ ,1997
+ ,6
+ ,365.37
+ ,1997
+ ,7
+ ,364.1
+ ,1997
+ ,8
+ ,361.89
+ ,1997
+ ,9
+ ,360.05
+ ,1997
+ ,10
+ ,360.49
+ ,1997
+ ,11
+ ,362.21
+ ,1997
+ ,12
+ ,364.12
+ ,1998
+ ,1
+ ,365
+ ,1998
+ ,2
+ ,365.82
+ ,1998
+ ,3
+ ,366.95
+ ,1998
+ ,4
+ ,368.42
+ ,1998
+ ,5
+ ,369.33
+ ,1998
+ ,6
+ ,368.78
+ ,1998
+ ,7
+ ,367.59
+ ,1998
+ ,8
+ ,365.81
+ ,1998
+ ,9
+ ,363.83
+ ,1998
+ ,10
+ ,364.18
+ ,1998
+ ,11
+ ,365.36
+ ,1998
+ ,12
+ ,366.88
+ ,1999
+ ,1
+ ,367.97
+ ,1999
+ ,2
+ ,368.83
+ ,1999
+ ,3
+ ,369.46
+ ,1999
+ ,4
+ ,370.77
+ ,1999
+ ,5
+ ,370.66
+ ,1999
+ ,6
+ ,370.1
+ ,1999
+ ,7
+ ,369.1
+ ,1999
+ ,8
+ ,366.7
+ ,1999
+ ,9
+ ,364.61
+ ,1999
+ ,10
+ ,365.17
+ ,1999
+ ,11
+ ,366.51
+ ,1999
+ ,12
+ ,367.86
+ ,2000
+ ,1
+ ,369.07
+ ,2000
+ ,2
+ ,369.32
+ ,2000
+ ,3
+ ,370.38
+ ,2000
+ ,4
+ ,371.63
+ ,2000
+ ,5
+ ,371.32
+ ,2000
+ ,6
+ ,371.51
+ ,2000
+ ,7
+ ,369.69
+ ,2000
+ ,8
+ ,368.18
+ ,2000
+ ,9
+ ,366.87
+ ,2000
+ ,10
+ ,366.94
+ ,2000
+ ,11
+ ,368.27
+ ,2000
+ ,12
+ ,369.62
+ ,2001
+ ,1
+ ,370.47
+ ,2001
+ ,2
+ ,371.44
+ ,2001
+ ,3
+ ,372.39
+ ,2001
+ ,4
+ ,373.32
+ ,2001
+ ,5
+ ,373.77
+ ,2001
+ ,6
+ ,373.13
+ ,2001
+ ,7
+ ,371.51
+ ,2001
+ ,8
+ ,369.59
+ ,2001
+ ,9
+ ,368.12
+ ,2001
+ ,10
+ ,368.38
+ ,2001
+ ,11
+ ,369.64
+ ,2001
+ ,12
+ ,371.11
+ ,2002
+ ,1
+ ,372.38
+ ,2002
+ ,2
+ ,373.08
+ ,2002
+ ,3
+ ,373.87
+ ,2002
+ ,4
+ ,374.93
+ ,2002
+ ,5
+ ,375.58
+ ,2002
+ ,6
+ ,375.44
+ ,2002
+ ,7
+ ,373.91
+ ,2002
+ ,8
+ ,371.77
+ ,2002
+ ,9
+ ,370.72
+ ,2002
+ ,10
+ ,370.5
+ ,2002
+ ,11
+ ,372.19
+ ,2002
+ ,12
+ ,373.71
+ ,2003
+ ,1
+ ,374.92
+ ,2003
+ ,2
+ ,375.63
+ ,2003
+ ,3
+ ,376.51
+ ,2003
+ ,4
+ ,377.75
+ ,2003
+ ,5
+ ,378.54
+ ,2003
+ ,6
+ ,378.21
+ ,2003
+ ,7
+ ,376.65
+ ,2003
+ ,8
+ ,374.28
+ ,2003
+ ,9
+ ,373.12
+ ,2003
+ ,10
+ ,373.1
+ ,2003
+ ,11
+ ,374.67
+ ,2003
+ ,12
+ ,375.97
+ ,2004
+ ,1
+ ,377.03
+ ,2004
+ ,2
+ ,377.87
+ ,2004
+ ,3
+ ,378.88
+ ,2004
+ ,4
+ ,380.42
+ ,2004
+ ,5
+ ,380.62
+ ,2004
+ ,6
+ ,379.66
+ ,2004
+ ,7
+ ,377.48
+ ,2004
+ ,8
+ ,376.07
+ ,2004
+ ,9
+ ,374.1
+ ,2004
+ ,10
+ ,374.47
+ ,2004
+ ,11
+ ,376.15
+ ,2004
+ ,12
+ ,377.51
+ ,2005
+ ,1
+ ,378.43
+ ,2005
+ ,2
+ ,379.7
+ ,2005
+ ,3
+ ,380.91
+ ,2005
+ ,4
+ ,382.2
+ ,2005
+ ,5
+ ,382.45
+ ,2005
+ ,6
+ ,382.14
+ ,2005
+ ,7
+ ,380.6
+ ,2005
+ ,8
+ ,378.6
+ ,2005
+ ,9
+ ,376.72
+ ,2005
+ ,10
+ ,376.98
+ ,2005
+ ,11
+ ,378.29
+ ,2005
+ ,12
+ ,380.07
+ ,2006
+ ,1
+ ,381.36
+ ,2006
+ ,2
+ ,382.19
+ ,2006
+ ,3
+ ,382.65
+ ,2006
+ ,4
+ ,384.65
+ ,2006
+ ,5
+ ,384.94
+ ,2006
+ ,6
+ ,384.01
+ ,2006
+ ,7
+ ,382.15
+ ,2006
+ ,8
+ ,380.33
+ ,2006
+ ,9
+ ,378.81
+ ,2006
+ ,10
+ ,379.06
+ ,2006
+ ,11
+ ,380.17
+ ,2006
+ ,12
+ ,381.85
+ ,2007
+ ,1
+ ,382.88
+ ,2007
+ ,2
+ ,383.77
+ ,2007
+ ,3
+ ,384.42
+ ,2007
+ ,4
+ ,386.36
+ ,2007
+ ,5
+ ,386.53
+ ,2007
+ ,6
+ ,386.01
+ ,2007
+ ,7
+ ,384.45
+ ,2007
+ ,8
+ ,381.96
+ ,2007
+ ,9
+ ,380.81
+ ,2007
+ ,10
+ ,381.09
+ ,2007
+ ,11
+ ,382.37
+ ,2007
+ ,12
+ ,383.84
+ ,2008
+ ,1
+ ,385.42
+ ,2008
+ ,2
+ ,385.72
+ ,2008
+ ,3
+ ,385.96
+ ,2008
+ ,4
+ ,387.18
+ ,2008
+ ,5
+ ,388.5
+ ,2008
+ ,6
+ ,387.88
+ ,2008
+ ,7
+ ,386.38
+ ,2008
+ ,8
+ ,384.15
+ ,2008
+ ,9
+ ,383.07
+ ,2008
+ ,10
+ ,382.98
+ ,2008
+ ,11
+ ,384.11
+ ,2008
+ ,12
+ ,385.54
+ ,2009
+ ,1
+ ,386.92
+ ,2009
+ ,2
+ ,387.41
+ ,2009
+ ,3
+ ,388.77
+ ,2009
+ ,4
+ ,389.46
+ ,2009
+ ,5
+ ,390.18
+ ,2009
+ ,6
+ ,389.43
+ ,2009
+ ,7
+ ,387.74
+ ,2009
+ ,8
+ ,385.91
+ ,2009
+ ,9
+ ,384.77
+ ,2009
+ ,10
+ ,384.38
+ ,2009
+ ,11
+ ,385.99
+ ,2009
+ ,12
+ ,387.26
+ ,2010
+ ,1
+ ,388.45
+ ,2010
+ ,2
+ ,389.7
+ ,2010
+ ,3
+ ,391.08
+ ,2010
+ ,4
+ ,392.46
+ ,2010
+ ,5
+ ,392.96
+ ,2010
+ ,6
+ ,392.03
+ ,2010
+ ,7
+ ,390.13
+ ,2010
+ ,8
+ ,388.15
+ ,2010
+ ,9
+ ,386.8
+ ,2010
+ ,10
+ ,387.18
+ ,2010
+ ,11
+ ,388.59)
+ ,dim=c(3
+ ,633)
+ ,dimnames=list(c('year'
+ ,'month'
+ ,'CO2')
+ ,1:633))
> y <- array(NA,dim=c(3,633),dimnames=list(c('year','month','CO2'),1:633))
> 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 = '3'
> #'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
> 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
CO2 year month
1 315.71 1958 3
2 317.45 1958 4
3 317.50 1958 5
4 317.12 1958 6
5 315.86 1958 7
6 314.93 1958 8
7 313.20 1958 9
8 312.60 1958 10
9 313.33 1958 11
10 314.67 1958 12
11 315.62 1959 1
12 316.38 1959 2
13 316.71 1959 3
14 317.72 1959 4
15 318.29 1959 5
16 318.16 1959 6
17 316.55 1959 7
18 314.80 1959 8
19 313.84 1959 9
20 313.26 1959 10
21 314.80 1959 11
22 315.59 1959 12
23 316.43 1960 1
24 316.97 1960 2
25 317.58 1960 3
26 319.02 1960 4
27 320.02 1960 5
28 319.59 1960 6
29 318.18 1960 7
30 315.91 1960 8
31 314.16 1960 9
32 313.83 1960 10
33 315.00 1960 11
34 316.19 1960 12
35 316.93 1961 1
36 317.70 1961 2
37 318.54 1961 3
38 319.48 1961 4
39 320.58 1961 5
40 319.77 1961 6
41 318.58 1961 7
42 316.79 1961 8
43 314.80 1961 9
44 315.38 1961 10
45 316.10 1961 11
46 317.01 1961 12
47 317.94 1962 1
48 318.55 1962 2
49 319.68 1962 3
50 320.63 1962 4
51 321.01 1962 5
52 320.55 1962 6
53 319.58 1962 7
54 317.40 1962 8
55 316.26 1962 9
56 315.42 1962 10
57 316.69 1962 11
58 317.70 1962 12
59 318.74 1963 1
60 319.08 1963 2
61 319.86 1963 3
62 321.39 1963 4
63 322.24 1963 5
64 321.47 1963 6
65 319.74 1963 7
66 317.77 1963 8
67 316.21 1963 9
68 315.99 1963 10
69 317.12 1963 11
70 318.31 1963 12
71 319.57 1964 1
72 320.08 1964 2
73 320.75 1964 3
74 321.80 1964 4
75 322.24 1964 5
76 321.89 1964 6
77 320.44 1964 7
78 318.70 1964 8
79 316.70 1964 9
80 316.79 1964 10
81 317.79 1964 11
82 318.71 1964 12
83 319.44 1965 1
84 320.44 1965 2
85 320.89 1965 3
86 322.13 1965 4
87 322.16 1965 5
88 321.87 1965 6
89 321.39 1965 7
90 318.80 1965 8
91 317.81 1965 9
92 317.30 1965 10
93 318.87 1965 11
94 319.42 1965 12
95 320.62 1966 1
96 321.59 1966 2
97 322.39 1966 3
98 323.87 1966 4
99 324.01 1966 5
100 323.75 1966 6
101 322.40 1966 7
102 320.37 1966 8
103 318.64 1966 9
104 318.10 1966 10
105 319.78 1966 11
106 321.08 1966 12
107 322.06 1967 1
108 322.50 1967 2
109 323.04 1967 3
110 324.42 1967 4
111 325.00 1967 5
112 324.09 1967 6
113 322.55 1967 7
114 320.92 1967 8
115 319.31 1967 9
116 319.31 1967 10
117 320.72 1967 11
118 321.96 1967 12
119 322.57 1968 1
120 323.15 1968 2
121 323.89 1968 3
122 325.02 1968 4
123 325.57 1968 5
124 325.36 1968 6
125 324.14 1968 7
126 322.03 1968 8
127 320.41 1968 9
128 320.25 1968 10
129 321.31 1968 11
130 322.84 1968 12
131 324.00 1969 1
132 324.42 1969 2
133 325.64 1969 3
134 326.66 1969 4
135 327.34 1969 5
136 326.76 1969 6
137 325.88 1969 7
138 323.67 1969 8
139 322.38 1969 9
140 321.78 1969 10
141 322.85 1969 11
142 324.12 1969 12
143 325.03 1970 1
144 325.99 1970 2
145 326.87 1970 3
146 328.14 1970 4
147 328.07 1970 5
148 327.66 1970 6
149 326.35 1970 7
150 324.69 1970 8
151 323.10 1970 9
152 323.16 1970 10
153 323.98 1970 11
154 325.13 1970 12
155 326.17 1971 1
156 326.68 1971 2
157 327.18 1971 3
158 327.78 1971 4
159 328.92 1971 5
160 328.57 1971 6
161 327.34 1971 7
162 325.46 1971 8
163 323.36 1971 9
164 323.56 1971 10
165 324.80 1971 11
166 326.01 1971 12
167 326.77 1972 1
168 327.63 1972 2
169 327.75 1972 3
170 329.72 1972 4
171 330.07 1972 5
172 329.09 1972 6
173 328.05 1972 7
174 326.32 1972 8
175 324.93 1972 9
176 325.06 1972 10
177 326.50 1972 11
178 327.55 1972 12
179 328.55 1973 1
180 329.56 1973 2
181 330.30 1973 3
182 331.50 1973 4
183 332.48 1973 5
184 332.07 1973 6
185 330.87 1973 7
186 329.31 1973 8
187 327.51 1973 9
188 327.18 1973 10
189 328.16 1973 11
190 328.64 1973 12
191 329.35 1974 1
192 330.71 1974 2
193 331.48 1974 3
194 332.65 1974 4
195 333.15 1974 5
196 332.13 1974 6
197 330.99 1974 7
198 329.17 1974 8
199 327.41 1974 9
200 327.21 1974 10
201 328.34 1974 11
202 329.50 1974 12
203 330.68 1975 1
204 331.41 1975 2
205 331.85 1975 3
206 333.29 1975 4
207 333.91 1975 5
208 333.40 1975 6
209 331.74 1975 7
210 329.88 1975 8
211 328.57 1975 9
212 328.35 1975 10
213 329.33 1975 11
214 330.58 1975 12
215 331.66 1976 1
216 332.75 1976 2
217 333.46 1976 3
218 334.78 1976 4
219 334.79 1976 5
220 334.05 1976 6
221 332.95 1976 7
222 330.64 1976 8
223 328.96 1976 9
224 328.77 1976 10
225 330.18 1976 11
226 331.65 1976 12
227 332.69 1977 1
228 333.23 1977 2
229 334.97 1977 3
230 336.03 1977 4
231 336.82 1977 5
232 336.10 1977 6
233 334.79 1977 7
234 332.53 1977 8
235 331.19 1977 9
236 331.21 1977 10
237 332.35 1977 11
238 333.47 1977 12
239 335.09 1978 1
240 335.26 1978 2
241 336.62 1978 3
242 337.77 1978 4
243 338.00 1978 5
244 337.98 1978 6
245 336.48 1978 7
246 334.37 1978 8
247 332.33 1978 9
248 332.40 1978 10
249 333.76 1978 11
250 334.83 1978 12
251 336.21 1979 1
252 336.64 1979 2
253 338.13 1979 3
254 338.96 1979 4
255 339.02 1979 5
256 339.20 1979 6
257 337.60 1979 7
258 335.56 1979 8
259 333.93 1979 9
260 334.12 1979 10
261 335.26 1979 11
262 336.77 1979 12
263 337.80 1980 1
264 338.28 1980 2
265 340.04 1980 3
266 340.86 1980 4
267 341.47 1980 5
268 341.26 1980 6
269 339.34 1980 7
270 337.45 1980 8
271 336.10 1980 9
272 336.05 1980 10
273 337.21 1980 11
274 338.29 1980 12
275 339.36 1981 1
276 340.51 1981 2
277 341.57 1981 3
278 342.56 1981 4
279 343.01 1981 5
280 342.52 1981 6
281 340.71 1981 7
282 338.51 1981 8
283 336.96 1981 9
284 337.13 1981 10
285 338.58 1981 11
286 339.91 1981 12
287 340.92 1982 1
288 341.69 1982 2
289 342.87 1982 3
290 343.83 1982 4
291 344.30 1982 5
292 343.42 1982 6
293 341.85 1982 7
294 339.82 1982 8
295 337.98 1982 9
296 338.09 1982 10
297 339.24 1982 11
298 340.67 1982 12
299 341.42 1983 1
300 342.67 1983 2
301 343.45 1983 3
302 345.08 1983 4
303 345.76 1983 5
304 345.32 1983 6
305 343.93 1983 7
306 342.08 1983 8
307 340.00 1983 9
308 340.12 1983 10
309 341.35 1983 11
310 342.89 1983 12
311 343.87 1984 1
312 344.59 1984 2
313 345.29 1984 3
314 346.59 1984 4
315 347.36 1984 5
316 346.80 1984 6
317 345.37 1984 7
318 343.06 1984 8
319 341.24 1984 9
320 341.54 1984 10
321 342.90 1984 11
322 344.36 1984 12
323 345.08 1985 1
324 345.89 1985 2
325 347.49 1985 3
326 348.02 1985 4
327 348.75 1985 5
328 348.19 1985 6
329 346.49 1985 7
330 344.70 1985 8
331 343.04 1985 9
332 342.92 1985 10
333 344.22 1985 11
334 345.61 1985 12
335 346.42 1986 1
336 346.95 1986 2
337 347.88 1986 3
338 349.57 1986 4
339 350.35 1986 5
340 349.70 1986 6
341 347.78 1986 7
342 345.89 1986 8
343 344.88 1986 9
344 344.34 1986 10
345 345.67 1986 11
346 346.89 1986 12
347 348.20 1987 1
348 348.55 1987 2
349 349.56 1987 3
350 351.12 1987 4
351 351.84 1987 5
352 351.45 1987 6
353 349.77 1987 7
354 347.62 1987 8
355 346.37 1987 9
356 346.48 1987 10
357 347.80 1987 11
358 349.03 1987 12
359 350.23 1988 1
360 351.58 1988 2
361 352.22 1988 3
362 353.53 1988 4
363 354.14 1988 5
364 353.64 1988 6
365 352.53 1988 7
366 350.42 1988 8
367 348.84 1988 9
368 348.94 1988 10
369 349.99 1988 11
370 351.29 1988 12
371 352.72 1989 1
372 353.10 1989 2
373 353.64 1989 3
374 355.43 1989 4
375 355.70 1989 5
376 355.11 1989 6
377 353.79 1989 7
378 351.42 1989 8
379 349.83 1989 9
380 350.10 1989 10
381 351.26 1989 11
382 352.66 1989 12
383 353.63 1990 1
384 354.72 1990 2
385 355.49 1990 3
386 356.10 1990 4
387 357.08 1990 5
388 356.11 1990 6
389 354.67 1990 7
390 352.67 1990 8
391 351.05 1990 9
392 351.36 1990 10
393 352.81 1990 11
394 354.21 1990 12
395 354.87 1991 1
396 355.67 1991 2
397 357.00 1991 3
398 358.40 1991 4
399 359.00 1991 5
400 357.99 1991 6
401 355.96 1991 7
402 353.78 1991 8
403 352.20 1991 9
404 352.22 1991 10
405 353.70 1991 11
406 354.98 1991 12
407 356.08 1992 1
408 356.84 1992 2
409 357.73 1992 3
410 358.91 1992 4
411 359.45 1992 5
412 359.19 1992 6
413 356.72 1992 7
414 354.77 1992 8
415 352.80 1992 9
416 353.21 1992 10
417 354.15 1992 11
418 355.39 1992 12
419 356.76 1993 1
420 357.17 1993 2
421 358.26 1993 3
422 359.17 1993 4
423 360.07 1993 5
424 359.41 1993 6
425 357.36 1993 7
426 355.29 1993 8
427 353.96 1993 9
428 354.03 1993 10
429 355.27 1993 11
430 356.70 1993 12
431 358.05 1994 1
432 358.80 1994 2
433 359.67 1994 3
434 361.13 1994 4
435 361.48 1994 5
436 360.60 1994 6
437 359.20 1994 7
438 357.23 1994 8
439 355.42 1994 9
440 355.89 1994 10
441 357.41 1994 11
442 358.74 1994 12
443 359.73 1995 1
444 360.61 1995 2
445 361.60 1995 3
446 363.05 1995 4
447 363.62 1995 5
448 363.03 1995 6
449 361.55 1995 7
450 358.94 1995 8
451 357.93 1995 9
452 357.80 1995 10
453 359.22 1995 11
454 360.42 1995 12
455 361.83 1996 1
456 362.94 1996 2
457 363.91 1996 3
458 364.28 1996 4
459 364.93 1996 5
460 364.70 1996 6
461 363.31 1996 7
462 361.15 1996 8
463 359.41 1996 9
464 359.34 1996 10
465 360.62 1996 11
466 361.96 1996 12
467 362.81 1997 1
468 363.87 1997 2
469 364.25 1997 3
470 366.02 1997 4
471 366.47 1997 5
472 365.37 1997 6
473 364.10 1997 7
474 361.89 1997 8
475 360.05 1997 9
476 360.49 1997 10
477 362.21 1997 11
478 364.12 1997 12
479 365.00 1998 1
480 365.82 1998 2
481 366.95 1998 3
482 368.42 1998 4
483 369.33 1998 5
484 368.78 1998 6
485 367.59 1998 7
486 365.81 1998 8
487 363.83 1998 9
488 364.18 1998 10
489 365.36 1998 11
490 366.88 1998 12
491 367.97 1999 1
492 368.83 1999 2
493 369.46 1999 3
494 370.77 1999 4
495 370.66 1999 5
496 370.10 1999 6
497 369.10 1999 7
498 366.70 1999 8
499 364.61 1999 9
500 365.17 1999 10
501 366.51 1999 11
502 367.86 1999 12
503 369.07 2000 1
504 369.32 2000 2
505 370.38 2000 3
506 371.63 2000 4
507 371.32 2000 5
508 371.51 2000 6
509 369.69 2000 7
510 368.18 2000 8
511 366.87 2000 9
512 366.94 2000 10
513 368.27 2000 11
514 369.62 2000 12
515 370.47 2001 1
516 371.44 2001 2
517 372.39 2001 3
518 373.32 2001 4
519 373.77 2001 5
520 373.13 2001 6
521 371.51 2001 7
522 369.59 2001 8
523 368.12 2001 9
524 368.38 2001 10
525 369.64 2001 11
526 371.11 2001 12
527 372.38 2002 1
528 373.08 2002 2
529 373.87 2002 3
530 374.93 2002 4
531 375.58 2002 5
532 375.44 2002 6
533 373.91 2002 7
534 371.77 2002 8
535 370.72 2002 9
536 370.50 2002 10
537 372.19 2002 11
538 373.71 2002 12
539 374.92 2003 1
540 375.63 2003 2
541 376.51 2003 3
542 377.75 2003 4
543 378.54 2003 5
544 378.21 2003 6
545 376.65 2003 7
546 374.28 2003 8
547 373.12 2003 9
548 373.10 2003 10
549 374.67 2003 11
550 375.97 2003 12
551 377.03 2004 1
552 377.87 2004 2
553 378.88 2004 3
554 380.42 2004 4
555 380.62 2004 5
556 379.66 2004 6
557 377.48 2004 7
558 376.07 2004 8
559 374.10 2004 9
560 374.47 2004 10
561 376.15 2004 11
562 377.51 2004 12
563 378.43 2005 1
564 379.70 2005 2
565 380.91 2005 3
566 382.20 2005 4
567 382.45 2005 5
568 382.14 2005 6
569 380.60 2005 7
570 378.60 2005 8
571 376.72 2005 9
572 376.98 2005 10
573 378.29 2005 11
574 380.07 2005 12
575 381.36 2006 1
576 382.19 2006 2
577 382.65 2006 3
578 384.65 2006 4
579 384.94 2006 5
580 384.01 2006 6
581 382.15 2006 7
582 380.33 2006 8
583 378.81 2006 9
584 379.06 2006 10
585 380.17 2006 11
586 381.85 2006 12
587 382.88 2007 1
588 383.77 2007 2
589 384.42 2007 3
590 386.36 2007 4
591 386.53 2007 5
592 386.01 2007 6
593 384.45 2007 7
594 381.96 2007 8
595 380.81 2007 9
596 381.09 2007 10
597 382.37 2007 11
598 383.84 2007 12
599 385.42 2008 1
600 385.72 2008 2
601 385.96 2008 3
602 387.18 2008 4
603 388.50 2008 5
604 387.88 2008 6
605 386.38 2008 7
606 384.15 2008 8
607 383.07 2008 9
608 382.98 2008 10
609 384.11 2008 11
610 385.54 2008 12
611 386.92 2009 1
612 387.41 2009 2
613 388.77 2009 3
614 389.46 2009 4
615 390.18 2009 5
616 389.43 2009 6
617 387.74 2009 7
618 385.91 2009 8
619 384.77 2009 9
620 384.38 2009 10
621 385.99 2009 11
622 387.26 2009 12
623 388.45 2010 1
624 389.70 2010 2
625 391.08 2010 3
626 392.46 2010 4
627 392.96 2010 5
628 392.03 2010 6
629 390.13 2010 7
630 388.15 2010 8
631 386.80 2010 9
632 387.18 2010 10
633 388.59 2010 11
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) year month
-2515.809 1.444 -0.266
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-5.4735 -2.3456 -0.4072 1.8711 8.3789
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.516e+03 1.590e+01 -158.267 < 2e-16 ***
year 1.444e+00 8.009e-03 180.243 < 2e-16 ***
month -2.660e-01 3.544e-02 -7.505 2.11e-13 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 3.068 on 630 degrees of freedom
Multiple R-squared: 0.981, Adjusted R-squared: 0.981
F-statistic: 1.629e+04 on 2 and 630 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,] 8.319892e-02 1.663978e-01 9.168011e-01
[2,] 6.612877e-02 1.322575e-01 9.338712e-01
[3,] 3.373141e-02 6.746282e-02 9.662686e-01
[4,] 1.366116e-02 2.732233e-02 9.863388e-01
[5,] 1.363555e-02 2.727110e-02 9.863644e-01
[6,] 5.427123e-03 1.085425e-02 9.945729e-01
[7,] 2.607672e-03 5.215344e-03 9.973923e-01
[8,] 1.478676e-03 2.957352e-03 9.985213e-01
[9,] 1.856986e-03 3.713972e-03 9.981430e-01
[10,] 2.894822e-03 5.789644e-03 9.971052e-01
[11,] 2.997857e-03 5.995714e-03 9.970021e-01
[12,] 1.489797e-03 2.979593e-03 9.985102e-01
[13,] 8.777695e-04 1.755539e-03 9.991222e-01
[14,] 6.145924e-04 1.229185e-03 9.993854e-01
[15,] 4.060664e-04 8.121328e-04 9.995939e-01
[16,] 2.150607e-04 4.301215e-04 9.997849e-01
[17,] 1.804292e-04 3.608584e-04 9.998196e-01
[18,] 1.289924e-04 2.579848e-04 9.998710e-01
[19,] 6.134902e-05 1.226980e-04 9.999387e-01
[20,] 3.139654e-05 6.279307e-05 9.999686e-01
[21,] 4.240619e-05 8.481238e-05 9.999576e-01
[22,] 1.340952e-04 2.681904e-04 9.998659e-01
[23,] 2.126924e-04 4.253847e-04 9.997873e-01
[24,] 1.459728e-04 2.919456e-04 9.998540e-01
[25,] 9.412330e-05 1.882466e-04 9.999059e-01
[26,] 1.144530e-04 2.289060e-04 9.998855e-01
[27,] 1.129978e-04 2.259957e-04 9.998870e-01
[28,] 6.467873e-05 1.293575e-04 9.999353e-01
[29,] 4.745770e-05 9.491541e-05 9.999525e-01
[30,] 4.981255e-05 9.962511e-05 9.999502e-01
[31,] 2.985362e-05 5.970723e-05 9.999701e-01
[32,] 1.756519e-05 3.513039e-05 9.999824e-01
[33,] 1.560035e-05 3.120070e-05 9.999844e-01
[34,] 3.196262e-05 6.392524e-05 9.999680e-01
[35,] 3.623205e-05 7.246411e-05 9.999638e-01
[36,] 2.572946e-05 5.145891e-05 9.999743e-01
[37,] 1.647119e-05 3.294238e-05 9.999835e-01
[38,] 1.978846e-05 3.957691e-05 9.999802e-01
[39,] 1.429178e-05 2.858356e-05 9.999857e-01
[40,] 8.746454e-06 1.749291e-05 9.999913e-01
[41,] 6.760307e-06 1.352061e-05 9.999932e-01
[42,] 6.536561e-06 1.307312e-05 9.999935e-01
[43,] 4.203516e-06 8.407031e-06 9.999958e-01
[44,] 2.774680e-06 5.549359e-06 9.999972e-01
[45,] 2.984837e-06 5.969673e-06 9.999970e-01
[46,] 4.511531e-06 9.023062e-06 9.999955e-01
[47,] 5.431285e-06 1.086257e-05 9.999946e-01
[48,] 4.450316e-06 8.900633e-06 9.999955e-01
[49,] 3.112676e-06 6.225351e-06 9.999969e-01
[50,] 2.786296e-06 5.572592e-06 9.999972e-01
[51,] 3.023627e-06 6.047255e-06 9.999970e-01
[52,] 1.978515e-06 3.957030e-06 9.999980e-01
[53,] 1.567950e-06 3.135900e-06 9.999984e-01
[54,] 1.601411e-06 3.202821e-06 9.999984e-01
[55,] 1.175239e-06 2.350478e-06 9.999988e-01
[56,] 7.509279e-07 1.501856e-06 9.999992e-01
[57,] 8.449413e-07 1.689883e-06 9.999992e-01
[58,] 1.942225e-06 3.884450e-06 9.999981e-01
[59,] 2.736977e-06 5.473955e-06 9.999973e-01
[60,] 2.039948e-06 4.079895e-06 9.999980e-01
[61,] 1.608950e-06 3.217900e-06 9.999984e-01
[62,] 2.104841e-06 4.209682e-06 9.999979e-01
[63,] 2.389384e-06 4.778768e-06 9.999976e-01
[64,] 1.683974e-06 3.367949e-06 9.999983e-01
[65,] 1.395777e-06 2.791554e-06 9.999986e-01
[66,] 1.339143e-06 2.678286e-06 9.999987e-01
[67,] 9.606451e-07 1.921290e-06 9.999990e-01
[68,] 6.573919e-07 1.314784e-06 9.999993e-01
[69,] 6.495443e-07 1.299089e-06 9.999994e-01
[70,] 9.210925e-07 1.842185e-06 9.999991e-01
[71,] 1.192811e-06 2.385623e-06 9.999988e-01
[72,] 9.429457e-07 1.885891e-06 9.999991e-01
[73,] 7.036069e-07 1.407214e-06 9.999993e-01
[74,] 9.790584e-07 1.958117e-06 9.999990e-01
[75,] 1.021215e-06 2.042431e-06 9.999990e-01
[76,] 7.325560e-07 1.465112e-06 9.999993e-01
[77,] 5.846024e-07 1.169205e-06 9.999994e-01
[78,] 8.430273e-07 1.686055e-06 9.999992e-01
[79,] 6.471877e-07 1.294375e-06 9.999994e-01
[80,] 4.460013e-07 8.920027e-07 9.999996e-01
[81,] 4.027989e-07 8.055979e-07 9.999996e-01
[82,] 4.091837e-07 8.183675e-07 9.999996e-01
[83,] 4.081992e-07 8.163983e-07 9.999996e-01
[84,] 3.722203e-07 7.444406e-07 9.999996e-01
[85,] 2.988123e-07 5.976247e-07 9.999997e-01
[86,] 2.975935e-07 5.951870e-07 9.999997e-01
[87,] 3.106217e-07 6.212435e-07 9.999997e-01
[88,] 2.201697e-07 4.403395e-07 9.999998e-01
[89,] 1.812746e-07 3.625492e-07 9.999998e-01
[90,] 1.820842e-07 3.641684e-07 9.999998e-01
[91,] 1.257387e-07 2.514775e-07 9.999999e-01
[92,] 9.262547e-08 1.852509e-07 9.999999e-01
[93,] 1.488327e-07 2.976655e-07 9.999999e-01
[94,] 3.024335e-07 6.048671e-07 9.999997e-01
[95,] 5.794805e-07 1.158961e-06 9.999994e-01
[96,] 5.902798e-07 1.180560e-06 9.999994e-01
[97,] 4.304363e-07 8.608727e-07 9.999996e-01
[98,] 4.197271e-07 8.394543e-07 9.999996e-01
[99,] 4.343711e-07 8.687421e-07 9.999996e-01
[100,] 3.178780e-07 6.357561e-07 9.999997e-01
[101,] 3.669502e-07 7.339005e-07 9.999996e-01
[102,] 3.012391e-07 6.024781e-07 9.999997e-01
[103,] 2.138384e-07 4.276768e-07 9.999998e-01
[104,] 1.549913e-07 3.099826e-07 9.999998e-01
[105,] 1.996948e-07 3.993896e-07 9.999998e-01
[106,] 4.309344e-07 8.618688e-07 9.999996e-01
[107,] 5.883125e-07 1.176625e-06 9.999994e-01
[108,] 4.839124e-07 9.678249e-07 9.999995e-01
[109,] 3.602405e-07 7.204811e-07 9.999996e-01
[110,] 3.671809e-07 7.343618e-07 9.999996e-01
[111,] 3.263451e-07 6.526901e-07 9.999997e-01
[112,] 2.425140e-07 4.850279e-07 9.999998e-01
[113,] 2.851287e-07 5.702574e-07 9.999997e-01
[114,] 2.590670e-07 5.181340e-07 9.999997e-01
[115,] 1.879671e-07 3.759341e-07 9.999998e-01
[116,] 1.369477e-07 2.738955e-07 9.999999e-01
[117,] 1.545000e-07 3.090000e-07 9.999998e-01
[118,] 2.738897e-07 5.477795e-07 9.999997e-01
[119,] 4.844131e-07 9.688261e-07 9.999995e-01
[120,] 5.075610e-07 1.015122e-06 9.999995e-01
[121,] 3.727112e-07 7.454225e-07 9.999996e-01
[122,] 3.443870e-07 6.887739e-07 9.999997e-01
[123,] 2.989392e-07 5.978783e-07 9.999997e-01
[124,] 2.174441e-07 4.348882e-07 9.999998e-01
[125,] 2.526826e-07 5.053652e-07 9.999997e-01
[126,] 1.956356e-07 3.912712e-07 9.999998e-01
[127,] 1.379172e-07 2.758344e-07 9.999999e-01
[128,] 1.220359e-07 2.440719e-07 9.999999e-01
[129,] 2.002708e-07 4.005416e-07 9.999998e-01
[130,] 6.174525e-07 1.234905e-06 9.999994e-01
[131,] 1.383490e-06 2.766981e-06 9.999986e-01
[132,] 2.044561e-06 4.089123e-06 9.999980e-01
[133,] 1.591390e-06 3.182779e-06 9.999984e-01
[134,] 1.217871e-06 2.435742e-06 9.999988e-01
[135,] 9.616793e-07 1.923359e-06 9.999990e-01
[136,] 7.516195e-07 1.503239e-06 9.999992e-01
[137,] 9.541079e-07 1.908216e-06 9.999990e-01
[138,] 7.616419e-07 1.523284e-06 9.999992e-01
[139,] 5.600165e-07 1.120033e-06 9.999994e-01
[140,] 5.221908e-07 1.044382e-06 9.999995e-01
[141,] 1.033370e-06 2.066739e-06 9.999990e-01
[142,] 2.214862e-06 4.429725e-06 9.999978e-01
[143,] 4.144122e-06 8.288243e-06 9.999959e-01
[144,] 4.566707e-06 9.133413e-06 9.999954e-01
[145,] 3.593105e-06 7.186211e-06 9.999964e-01
[146,] 2.913379e-06 5.826758e-06 9.999971e-01
[147,] 2.250791e-06 4.501583e-06 9.999977e-01
[148,] 1.800125e-06 3.600250e-06 9.999982e-01
[149,] 2.189726e-06 4.379452e-06 9.999978e-01
[150,] 1.756591e-06 3.513182e-06 9.999982e-01
[151,] 1.297376e-06 2.594753e-06 9.999987e-01
[152,] 1.010240e-06 2.020480e-06 9.999990e-01
[153,] 9.735761e-07 1.947152e-06 9.999990e-01
[154,] 1.823251e-06 3.646501e-06 9.999982e-01
[155,] 3.140263e-06 6.280527e-06 9.999969e-01
[156,] 3.381864e-06 6.763729e-06 9.999966e-01
[157,] 2.589448e-06 5.178895e-06 9.999974e-01
[158,] 2.492969e-06 4.985939e-06 9.999975e-01
[159,] 2.089852e-06 4.179704e-06 9.999979e-01
[160,] 1.619136e-06 3.238272e-06 9.999984e-01
[161,] 1.838613e-06 3.677226e-06 9.999982e-01
[162,] 1.602747e-06 3.205494e-06 9.999984e-01
[163,] 1.178307e-06 2.356615e-06 9.999988e-01
[164,] 8.696783e-07 1.739357e-06 9.999991e-01
[165,] 1.229186e-06 2.458373e-06 9.999988e-01
[166,] 2.401209e-06 4.802418e-06 9.999976e-01
[167,] 3.107677e-06 6.215353e-06 9.999969e-01
[168,] 2.976009e-06 5.952018e-06 9.999970e-01
[169,] 2.241711e-06 4.483421e-06 9.999978e-01
[170,] 1.850778e-06 3.701555e-06 9.999981e-01
[171,] 1.424715e-06 2.849430e-06 9.999986e-01
[172,] 1.220523e-06 2.441046e-06 9.999988e-01
[173,] 1.688189e-06 3.376378e-06 9.999983e-01
[174,] 1.278624e-06 2.557248e-06 9.999987e-01
[175,] 9.861526e-07 1.972305e-06 9.999990e-01
[176,] 9.747109e-07 1.949422e-06 9.999990e-01
[177,] 1.944766e-06 3.889533e-06 9.999981e-01
[178,] 8.740602e-06 1.748120e-05 9.999913e-01
[179,] 3.048103e-05 6.096207e-05 9.999695e-01
[180,] 5.536195e-05 1.107239e-04 9.999446e-01
[181,] 5.823272e-05 1.164654e-04 9.999418e-01
[182,] 4.643810e-05 9.287621e-05 9.999536e-01
[183,] 3.684006e-05 7.368011e-05 9.999632e-01
[184,] 3.542933e-05 7.085866e-05 9.999646e-01
[185,] 4.300722e-05 8.601445e-05 9.999570e-01
[186,] 3.639438e-05 7.278876e-05 9.999636e-01
[187,] 2.892881e-05 5.785763e-05 9.999711e-01
[188,] 2.778260e-05 5.556520e-05 9.999722e-01
[189,] 4.576299e-05 9.152599e-05 9.999542e-01
[190,] 1.088482e-04 2.176963e-04 9.998912e-01
[191,] 1.683069e-04 3.366137e-04 9.998317e-01
[192,] 1.858735e-04 3.717469e-04 9.998141e-01
[193,] 1.524718e-04 3.049436e-04 9.998475e-01
[194,] 1.304597e-04 2.609194e-04 9.998695e-01
[195,] 1.098890e-04 2.197779e-04 9.998901e-01
[196,] 8.972058e-05 1.794412e-04 9.999103e-01
[197,] 9.674328e-05 1.934866e-04 9.999033e-01
[198,] 8.083436e-05 1.616687e-04 9.999192e-01
[199,] 6.320232e-05 1.264046e-04 9.999368e-01
[200,] 5.182382e-05 1.036476e-04 9.999482e-01
[201,] 6.581022e-05 1.316204e-04 9.999342e-01
[202,] 1.215233e-04 2.430466e-04 9.998785e-01
[203,] 1.933882e-04 3.867764e-04 9.998066e-01
[204,] 1.892589e-04 3.785178e-04 9.998107e-01
[205,] 1.513240e-04 3.026480e-04 9.998487e-01
[206,] 1.284901e-04 2.569803e-04 9.998715e-01
[207,] 1.077731e-04 2.155463e-04 9.998922e-01
[208,] 8.637151e-05 1.727430e-04 9.999136e-01
[209,] 8.999980e-05 1.799996e-04 9.999100e-01
[210,] 7.740307e-05 1.548061e-04 9.999226e-01
[211,] 6.071793e-05 1.214359e-04 9.999393e-01
[212,] 5.307021e-05 1.061404e-04 9.999469e-01
[213,] 7.501475e-05 1.500295e-04 9.999250e-01
[214,] 1.170506e-04 2.341012e-04 9.998829e-01
[215,] 1.480898e-04 2.961796e-04 9.998519e-01
[216,] 1.447018e-04 2.894036e-04 9.998553e-01
[217,] 1.156118e-04 2.312236e-04 9.998844e-01
[218,] 1.133268e-04 2.266536e-04 9.998867e-01
[219,] 1.080165e-04 2.160330e-04 9.998920e-01
[220,] 8.517918e-05 1.703584e-04 9.999148e-01
[221,] 8.560972e-05 1.712194e-04 9.999144e-01
[222,] 7.597122e-05 1.519424e-04 9.999240e-01
[223,] 6.044584e-05 1.208917e-04 9.999396e-01
[224,] 5.638566e-05 1.127713e-04 9.999436e-01
[225,] 8.036793e-05 1.607359e-04 9.999196e-01
[226,] 1.837621e-04 3.675243e-04 9.998162e-01
[227,] 3.180895e-04 6.361790e-04 9.996819e-01
[228,] 3.645870e-04 7.291739e-04 9.996354e-01
[229,] 2.950617e-04 5.901234e-04 9.997049e-01
[230,] 2.446618e-04 4.893237e-04 9.997553e-01
[231,] 1.977021e-04 3.954042e-04 9.998023e-01
[232,] 1.683786e-04 3.367571e-04 9.998316e-01
[233,] 1.961963e-04 3.923926e-04 9.998038e-01
[234,] 1.569552e-04 3.139105e-04 9.998430e-01
[235,] 1.253717e-04 2.507433e-04 9.998746e-01
[236,] 1.266221e-04 2.532442e-04 9.998734e-01
[237,] 2.072845e-04 4.145689e-04 9.997927e-01
[238,] 4.093690e-04 8.187381e-04 9.995906e-01
[239,] 8.594764e-04 1.718953e-03 9.991405e-01
[240,] 1.074729e-03 2.149457e-03 9.989253e-01
[241,] 9.181648e-04 1.836330e-03 9.990818e-01
[242,] 7.891079e-04 1.578216e-03 9.992109e-01
[243,] 6.561544e-04 1.312309e-03 9.993438e-01
[244,] 5.717655e-04 1.143531e-03 9.994282e-01
[245,] 6.559952e-04 1.311990e-03 9.993440e-01
[246,] 5.426296e-04 1.085259e-03 9.994574e-01
[247,] 4.440185e-04 8.880370e-04 9.995560e-01
[248,] 4.618896e-04 9.237791e-04 9.995381e-01
[249,] 6.660675e-04 1.332135e-03 9.993339e-01
[250,] 1.059443e-03 2.118886e-03 9.989406e-01
[251,] 1.942791e-03 3.885582e-03 9.980572e-01
[252,] 2.224385e-03 4.448771e-03 9.977756e-01
[253,] 1.896739e-03 3.793479e-03 9.981033e-01
[254,] 1.614096e-03 3.228192e-03 9.983859e-01
[255,] 1.339884e-03 2.679768e-03 9.986601e-01
[256,] 1.197858e-03 2.395716e-03 9.988021e-01
[257,] 1.596740e-03 3.193480e-03 9.984033e-01
[258,] 1.335759e-03 2.671518e-03 9.986642e-01
[259,] 1.125244e-03 2.250487e-03 9.988748e-01
[260,] 1.318730e-03 2.637459e-03 9.986813e-01
[261,] 2.160887e-03 4.321775e-03 9.978391e-01
[262,] 4.779587e-03 9.559174e-03 9.952204e-01
[263,] 9.956626e-03 1.991325e-02 9.900434e-01
[264,] 1.187843e-02 2.375687e-02 9.881216e-01
[265,] 1.080706e-02 2.161412e-02 9.891929e-01
[266,] 9.283157e-03 1.856631e-02 9.907168e-01
[267,] 7.978980e-03 1.595796e-02 9.920210e-01
[268,] 7.757985e-03 1.551597e-02 9.922420e-01
[269,] 9.916252e-03 1.983250e-02 9.900837e-01
[270,] 8.563683e-03 1.712737e-02 9.914363e-01
[271,] 7.886857e-03 1.577371e-02 9.921131e-01
[272,] 9.012953e-03 1.802591e-02 9.909870e-01
[273,] 1.417318e-02 2.834635e-02 9.858268e-01
[274,] 2.671585e-02 5.343170e-02 9.732841e-01
[275,] 4.398429e-02 8.796858e-02 9.560157e-01
[276,] 4.944084e-02 9.888167e-02 9.505592e-01
[277,] 4.474469e-02 8.948938e-02 9.552553e-01
[278,] 4.003635e-02 8.007271e-02 9.599636e-01
[279,] 3.541687e-02 7.083375e-02 9.645831e-01
[280,] 3.424463e-02 6.848926e-02 9.657554e-01
[281,] 4.284869e-02 8.569738e-02 9.571513e-01
[282,] 3.814020e-02 7.628041e-02 9.618598e-01
[283,] 3.489812e-02 6.979623e-02 9.651019e-01
[284,] 3.750299e-02 7.500598e-02 9.624970e-01
[285,] 5.079872e-02 1.015974e-01 9.492013e-01
[286,] 7.950061e-02 1.590012e-01 9.204994e-01
[287,] 1.042332e-01 2.084665e-01 8.957668e-01
[288,] 1.096479e-01 2.192959e-01 8.903521e-01
[289,] 1.004829e-01 2.009658e-01 8.995171e-01
[290,] 9.332434e-02 1.866487e-01 9.066757e-01
[291,] 8.517839e-02 1.703568e-01 9.148216e-01
[292,] 7.836864e-02 1.567373e-01 9.216314e-01
[293,] 8.447714e-02 1.689543e-01 9.155229e-01
[294,] 7.904688e-02 1.580938e-01 9.209531e-01
[295,] 7.194678e-02 1.438936e-01 9.280532e-01
[296,] 6.928333e-02 1.385667e-01 9.307167e-01
[297,] 8.640192e-02 1.728038e-01 9.135981e-01
[298,] 1.272639e-01 2.545278e-01 8.727361e-01
[299,] 1.733584e-01 3.467168e-01 8.266416e-01
[300,] 1.939612e-01 3.879223e-01 8.060388e-01
[301,] 1.868718e-01 3.737437e-01 8.131282e-01
[302,] 1.731175e-01 3.462350e-01 8.268825e-01
[303,] 1.594345e-01 3.188690e-01 8.405655e-01
[304,] 1.538792e-01 3.077583e-01 8.461208e-01
[305,] 1.793937e-01 3.587874e-01 8.206063e-01
[306,] 1.660975e-01 3.321949e-01 8.339025e-01
[307,] 1.559777e-01 3.119553e-01 8.440223e-01
[308,] 1.550546e-01 3.101093e-01 8.449454e-01
[309,] 1.839464e-01 3.678929e-01 8.160536e-01
[310,] 2.513633e-01 5.027265e-01 7.486367e-01
[311,] 3.158672e-01 6.317344e-01 6.841328e-01
[312,] 3.412335e-01 6.824669e-01 6.587665e-01
[313,] 3.254237e-01 6.508473e-01 6.745763e-01
[314,] 3.083172e-01 6.166343e-01 6.916828e-01
[315,] 2.894635e-01 5.789270e-01 7.105365e-01
[316,] 2.822436e-01 5.644872e-01 7.177564e-01
[317,] 3.159088e-01 6.318175e-01 6.840912e-01
[318,] 2.986931e-01 5.973863e-01 7.013069e-01
[319,] 2.829638e-01 5.659276e-01 7.170362e-01
[320,] 2.959456e-01 5.918913e-01 7.040544e-01
[321,] 3.324904e-01 6.649808e-01 6.675096e-01
[322,] 4.116354e-01 8.232708e-01 5.883646e-01
[323,] 4.799028e-01 9.598056e-01 5.200972e-01
[324,] 4.950400e-01 9.900801e-01 5.049600e-01
[325,] 4.791587e-01 9.583174e-01 5.208413e-01
[326,] 4.583604e-01 9.167208e-01 5.416396e-01
[327,] 4.372309e-01 8.744617e-01 5.627691e-01
[328,] 4.257030e-01 8.514059e-01 5.742970e-01
[329,] 4.540401e-01 9.080802e-01 5.459599e-01
[330,] 4.360003e-01 8.720006e-01 5.639997e-01
[331,] 4.160632e-01 8.321264e-01 5.839368e-01
[332,] 4.061601e-01 8.123203e-01 5.938399e-01
[333,] 4.466712e-01 8.933424e-01 5.533288e-01
[334,] 5.324813e-01 9.350374e-01 4.675187e-01
[335,] 5.988269e-01 8.023462e-01 4.011731e-01
[336,] 6.070065e-01 7.859870e-01 3.929935e-01
[337,] 5.886137e-01 8.227727e-01 4.113863e-01
[338,] 5.675242e-01 8.649517e-01 4.324758e-01
[339,] 5.469529e-01 9.060941e-01 4.530471e-01
[340,] 5.342296e-01 9.315407e-01 4.657704e-01
[341,] 5.543999e-01 8.912003e-01 4.456001e-01
[342,] 5.352853e-01 9.294294e-01 4.647147e-01
[343,] 5.154400e-01 9.691200e-01 4.845600e-01
[344,] 5.074732e-01 9.850537e-01 4.925268e-01
[345,] 5.480926e-01 9.038148e-01 4.519074e-01
[346,] 6.281626e-01 7.436749e-01 3.718374e-01
[347,] 6.975387e-01 6.049225e-01 3.024613e-01
[348,] 7.156038e-01 5.687924e-01 2.843962e-01
[349,] 7.007317e-01 5.985367e-01 2.992683e-01
[350,] 6.818672e-01 6.362656e-01 3.181328e-01
[351,] 6.627762e-01 6.744476e-01 3.372238e-01
[352,] 6.602336e-01 6.795329e-01 3.397664e-01
[353,] 6.963754e-01 6.072492e-01 3.036246e-01
[354,] 6.785292e-01 6.429416e-01 3.214708e-01
[355,] 6.743201e-01 6.513597e-01 3.256799e-01
[356,] 6.890588e-01 6.218824e-01 3.109412e-01
[357,] 7.523399e-01 4.953203e-01 2.476601e-01
[358,] 8.371191e-01 3.257618e-01 1.628809e-01
[359,] 8.946969e-01 2.106061e-01 1.053031e-01
[360,] 9.219299e-01 1.561403e-01 7.807013e-02
[361,] 9.231925e-01 1.536150e-01 7.680748e-02
[362,] 9.163525e-01 1.672949e-01 8.364747e-02
[363,] 9.107236e-01 1.785529e-01 8.927644e-02
[364,] 9.158474e-01 1.683052e-01 8.415258e-02
[365,] 9.401962e-01 1.196075e-01 5.980377e-02
[366,] 9.354218e-01 1.291563e-01 6.457817e-02
[367,] 9.334919e-01 1.330163e-01 6.650814e-02
[368,] 9.374535e-01 1.250930e-01 6.254650e-02
[369,] 9.614823e-01 7.703542e-02 3.851771e-02
[370,] 9.811545e-01 3.769110e-02 1.884555e-02
[371,] 9.905361e-01 1.892788e-02 9.463939e-03
[372,] 9.936259e-01 1.274817e-02 6.374087e-03
[373,] 9.932991e-01 1.340180e-02 6.700901e-03
[374,] 9.921780e-01 1.564395e-02 7.821973e-03
[375,] 9.911205e-01 1.775906e-02 8.879528e-03
[376,] 9.915560e-01 1.688800e-02 8.444000e-03
[377,] 9.947756e-01 1.044886e-02 5.224431e-03
[378,] 9.939152e-01 1.216959e-02 6.084797e-03
[379,] 9.937299e-01 1.254013e-02 6.270063e-03
[380,] 9.946926e-01 1.061482e-02 5.307412e-03
[381,] 9.965132e-01 6.973591e-03 3.486795e-03
[382,] 9.986610e-01 2.677915e-03 1.338957e-03
[383,] 9.993532e-01 1.293613e-03 6.468066e-04
[384,] 9.995282e-01 9.436265e-04 4.718132e-04
[385,] 9.994795e-01 1.041077e-03 5.205386e-04
[386,] 9.993600e-01 1.280044e-03 6.400221e-04
[387,] 9.992331e-01 1.533783e-03 7.668917e-04
[388,] 9.992987e-01 1.402549e-03 7.012745e-04
[389,] 9.996432e-01 7.135117e-04 3.567559e-04
[390,] 9.995608e-01 8.783693e-04 4.391847e-04
[391,] 9.995023e-01 9.954045e-04 4.977023e-04
[392,] 9.996032e-01 7.936854e-04 3.968427e-04
[393,] 9.998449e-01 3.102083e-04 1.551042e-04
[394,] 9.999682e-01 6.351195e-05 3.175597e-05
[395,] 9.999910e-01 1.802185e-05 9.010923e-06
[396,] 9.999938e-01 1.240227e-05 6.201133e-06
[397,] 9.999925e-01 1.491392e-05 7.456958e-06
[398,] 9.999902e-01 1.954624e-05 9.773122e-06
[399,] 9.999872e-01 2.557560e-05 1.278780e-05
[400,] 9.999868e-01 2.645115e-05 1.322558e-05
[401,] 9.999924e-01 1.525233e-05 7.626164e-06
[402,] 9.999900e-01 1.996865e-05 9.984324e-06
[403,] 9.999877e-01 2.453197e-05 1.226599e-05
[404,] 9.999882e-01 2.351007e-05 1.175503e-05
[405,] 9.999938e-01 1.239640e-05 6.198198e-06
[406,] 9.999981e-01 3.804017e-06 1.902008e-06
[407,] 9.999995e-01 1.039704e-06 5.198522e-07
[408,] 9.999995e-01 9.165617e-07 4.582808e-07
[409,] 9.999994e-01 1.213140e-06 6.065702e-07
[410,] 9.999992e-01 1.500388e-06 7.501939e-07
[411,] 9.999990e-01 2.008474e-06 1.004237e-06
[412,] 9.999987e-01 2.598844e-06 1.299422e-06
[413,] 9.999988e-01 2.316929e-06 1.158464e-06
[414,] 9.999986e-01 2.897329e-06 1.448664e-06
[415,] 9.999981e-01 3.861110e-06 1.930555e-06
[416,] 9.999976e-01 4.750386e-06 2.375193e-06
[417,] 9.999978e-01 4.308690e-06 2.154345e-06
[418,] 9.999989e-01 2.272655e-06 1.136327e-06
[419,] 9.999993e-01 1.428156e-06 7.140782e-07
[420,] 9.999992e-01 1.695575e-06 8.477877e-07
[421,] 9.999989e-01 2.256983e-06 1.128491e-06
[422,] 9.999988e-01 2.476838e-06 1.238419e-06
[423,] 9.999985e-01 2.914438e-06 1.457219e-06
[424,] 9.999980e-01 3.925702e-06 1.962851e-06
[425,] 9.999980e-01 3.994836e-06 1.997418e-06
[426,] 9.999977e-01 4.554491e-06 2.277245e-06
[427,] 9.999970e-01 5.974356e-06 2.987178e-06
[428,] 9.999962e-01 7.618653e-06 3.809326e-06
[429,] 9.999968e-01 6.317223e-06 3.158612e-06
[430,] 9.999980e-01 3.929081e-06 1.964540e-06
[431,] 9.999984e-01 3.174426e-06 1.587213e-06
[432,] 9.999982e-01 3.665669e-06 1.832835e-06
[433,] 9.999976e-01 4.886558e-06 2.443279e-06
[434,] 9.999976e-01 4.867829e-06 2.433915e-06
[435,] 9.999971e-01 5.767372e-06 2.883686e-06
[436,] 9.999962e-01 7.520794e-06 3.760397e-06
[437,] 9.999965e-01 6.956501e-06 3.478251e-06
[438,] 9.999962e-01 7.645147e-06 3.822574e-06
[439,] 9.999950e-01 9.907873e-06 4.953937e-06
[440,] 9.999939e-01 1.225569e-05 6.127844e-06
[441,] 9.999952e-01 9.515157e-06 4.757579e-06
[442,] 9.999976e-01 4.816818e-06 2.408409e-06
[443,] 9.999986e-01 2.806533e-06 1.403267e-06
[444,] 9.999986e-01 2.730728e-06 1.365364e-06
[445,] 9.999982e-01 3.634267e-06 1.817134e-06
[446,] 9.999978e-01 4.352214e-06 2.176107e-06
[447,] 9.999974e-01 5.286813e-06 2.643406e-06
[448,] 9.999966e-01 6.853189e-06 3.426595e-06
[449,] 9.999968e-01 6.445205e-06 3.222603e-06
[450,] 9.999962e-01 7.548070e-06 3.774035e-06
[451,] 9.999951e-01 9.851796e-06 4.925898e-06
[452,] 9.999945e-01 1.099277e-05 5.496383e-06
[453,] 9.999949e-01 1.025948e-05 5.129741e-06
[454,] 9.999967e-01 6.501416e-06 3.250708e-06
[455,] 9.999980e-01 3.946945e-06 1.973472e-06
[456,] 9.999981e-01 3.824670e-06 1.912335e-06
[457,] 9.999974e-01 5.129763e-06 2.564882e-06
[458,] 9.999971e-01 5.728681e-06 2.864341e-06
[459,] 9.999967e-01 6.588794e-06 3.294397e-06
[460,] 9.999956e-01 8.713511e-06 4.356755e-06
[461,] 9.999956e-01 8.849218e-06 4.424609e-06
[462,] 9.999959e-01 8.256526e-06 4.128263e-06
[463,] 9.999949e-01 1.024561e-05 5.122804e-06
[464,] 9.999934e-01 1.324996e-05 6.624981e-06
[465,] 9.999936e-01 1.273937e-05 6.369687e-06
[466,] 9.999954e-01 9.290452e-06 4.645226e-06
[467,] 9.999954e-01 9.197606e-06 4.598803e-06
[468,] 9.999943e-01 1.140156e-05 5.700778e-06
[469,] 9.999933e-01 1.344245e-05 6.721226e-06
[470,] 9.999954e-01 9.237221e-06 4.618611e-06
[471,] 9.999959e-01 8.102422e-06 4.051211e-06
[472,] 9.999948e-01 1.049697e-05 5.248487e-06
[473,] 9.999950e-01 9.950829e-06 4.975415e-06
[474,] 9.999952e-01 9.551777e-06 4.775888e-06
[475,] 9.999943e-01 1.148854e-05 5.744269e-06
[476,] 9.999929e-01 1.422565e-05 7.112824e-06
[477,] 9.999942e-01 1.159106e-05 5.795530e-06
[478,] 9.999974e-01 5.104897e-06 2.552449e-06
[479,] 9.999987e-01 2.536843e-06 1.268421e-06
[480,] 9.999990e-01 1.993278e-06 9.966390e-07
[481,] 9.999988e-01 2.481740e-06 1.240870e-06
[482,] 9.999985e-01 3.086987e-06 1.543493e-06
[483,] 9.999980e-01 4.067245e-06 2.033622e-06
[484,] 9.999976e-01 4.800041e-06 2.400021e-06
[485,] 9.999986e-01 2.872053e-06 1.436027e-06
[486,] 9.999982e-01 3.606711e-06 1.803356e-06
[487,] 9.999977e-01 4.632822e-06 2.316411e-06
[488,] 9.999974e-01 5.162929e-06 2.581465e-06
[489,] 9.999985e-01 3.082953e-06 1.541477e-06
[490,] 9.999992e-01 1.618852e-06 8.094262e-07
[491,] 9.999995e-01 9.706317e-07 4.853159e-07
[492,] 9.999996e-01 8.262334e-07 4.131167e-07
[493,] 9.999994e-01 1.138002e-06 5.690008e-07
[494,] 9.999994e-01 1.128357e-06 5.641785e-07
[495,] 9.999993e-01 1.379611e-06 6.898056e-07
[496,] 9.999991e-01 1.835107e-06 9.175534e-07
[497,] 9.999992e-01 1.608259e-06 8.041294e-07
[498,] 9.999991e-01 1.791504e-06 8.957520e-07
[499,] 9.999989e-01 2.196203e-06 1.098102e-06
[500,] 9.999986e-01 2.884300e-06 1.442150e-06
[501,] 9.999986e-01 2.723491e-06 1.361745e-06
[502,] 9.999987e-01 2.598712e-06 1.299356e-06
[503,] 9.999990e-01 2.034100e-06 1.017050e-06
[504,] 9.999987e-01 2.587052e-06 1.293526e-06
[505,] 9.999983e-01 3.455392e-06 1.727696e-06
[506,] 9.999981e-01 3.742476e-06 1.871238e-06
[507,] 9.999978e-01 4.326580e-06 2.163290e-06
[508,] 9.999971e-01 5.773550e-06 2.886775e-06
[509,] 9.999972e-01 5.517258e-06 2.758629e-06
[510,] 9.999975e-01 4.936384e-06 2.468192e-06
[511,] 9.999970e-01 6.026087e-06 3.013043e-06
[512,] 9.999961e-01 7.819598e-06 3.909799e-06
[513,] 9.999959e-01 8.253293e-06 4.126646e-06
[514,] 9.999967e-01 6.681091e-06 3.340545e-06
[515,] 9.999969e-01 6.210947e-06 3.105474e-06
[516,] 9.999960e-01 8.047109e-06 4.023555e-06
[517,] 9.999951e-01 9.720957e-06 4.860478e-06
[518,] 9.999961e-01 7.729994e-06 3.864997e-06
[519,] 9.999964e-01 7.201485e-06 3.600742e-06
[520,] 9.999953e-01 9.363772e-06 4.681886e-06
[521,] 9.999946e-01 1.082140e-05 5.410699e-06
[522,] 9.999957e-01 8.562029e-06 4.281015e-06
[523,] 9.999955e-01 9.089111e-06 4.544555e-06
[524,] 9.999942e-01 1.153981e-05 5.769903e-06
[525,] 9.999931e-01 1.387655e-05 6.938276e-06
[526,] 9.999934e-01 1.313918e-05 6.569589e-06
[527,] 9.999941e-01 1.175749e-05 5.878743e-06
[528,] 9.999926e-01 1.472356e-05 7.361778e-06
[529,] 9.999913e-01 1.735812e-05 8.679060e-06
[530,] 9.999921e-01 1.580782e-05 7.903911e-06
[531,] 9.999931e-01 1.388906e-05 6.944528e-06
[532,] 9.999909e-01 1.825585e-05 9.127924e-06
[533,] 9.999906e-01 1.885221e-05 9.426103e-06
[534,] 9.999917e-01 1.655980e-05 8.279899e-06
[535,] 9.999906e-01 1.881879e-05 9.409396e-06
[536,] 9.999879e-01 2.421849e-05 1.210925e-05
[537,] 9.999874e-01 2.526773e-05 1.263386e-05
[538,] 9.999915e-01 1.697413e-05 8.487064e-06
[539,] 9.999944e-01 1.125116e-05 5.625579e-06
[540,] 9.999938e-01 1.242974e-05 6.214868e-06
[541,] 9.999919e-01 1.618568e-05 8.092839e-06
[542,] 9.999916e-01 1.686715e-05 8.433574e-06
[543,] 9.999908e-01 1.831131e-05 9.155653e-06
[544,] 9.999877e-01 2.458845e-05 1.229422e-05
[545,] 9.999882e-01 2.352081e-05 1.176041e-05
[546,] 9.999892e-01 2.156275e-05 1.078137e-05
[547,] 9.999869e-01 2.620493e-05 1.310246e-05
[548,] 9.999828e-01 3.444863e-05 1.722431e-05
[549,] 9.999855e-01 2.891053e-05 1.445527e-05
[550,] 9.999908e-01 1.831941e-05 9.159707e-06
[551,] 9.999922e-01 1.565850e-05 7.829248e-06
[552,] 9.999891e-01 2.178317e-05 1.089159e-05
[553,] 9.999858e-01 2.839497e-05 1.419748e-05
[554,] 9.999904e-01 1.916763e-05 9.583816e-06
[555,] 9.999920e-01 1.604691e-05 8.023453e-06
[556,] 9.999889e-01 2.219068e-05 1.109534e-05
[557,] 9.999867e-01 2.652672e-05 1.326336e-05
[558,] 9.999915e-01 1.695838e-05 8.479189e-06
[559,] 9.999903e-01 1.938254e-05 9.691269e-06
[560,] 9.999866e-01 2.680683e-05 1.340341e-05
[561,] 9.999864e-01 2.712982e-05 1.356491e-05
[562,] 9.999895e-01 2.093465e-05 1.046732e-05
[563,] 9.999923e-01 1.548458e-05 7.742292e-06
[564,] 9.999905e-01 1.894459e-05 9.472296e-06
[565,] 9.999863e-01 2.733558e-05 1.366779e-05
[566,] 9.999880e-01 2.406341e-05 1.203171e-05
[567,] 9.999877e-01 2.465617e-05 1.232808e-05
[568,] 9.999821e-01 3.581035e-05 1.790518e-05
[569,] 9.999813e-01 3.733576e-05 1.866788e-05
[570,] 9.999812e-01 3.768060e-05 1.884030e-05
[571,] 9.999747e-01 5.060767e-05 2.530384e-05
[572,] 9.999634e-01 7.328342e-05 3.664171e-05
[573,] 9.999685e-01 6.306695e-05 3.153347e-05
[574,] 9.999824e-01 3.515585e-05 1.757793e-05
[575,] 9.999869e-01 2.622796e-05 1.311398e-05
[576,] 9.999817e-01 3.659901e-05 1.829950e-05
[577,] 9.999724e-01 5.520670e-05 2.760335e-05
[578,] 9.999723e-01 5.534626e-05 2.767313e-05
[579,] 9.999683e-01 6.336103e-05 3.168052e-05
[580,] 9.999515e-01 9.694913e-05 4.847456e-05
[581,] 9.999444e-01 1.112344e-04 5.561722e-05
[582,] 9.999513e-01 9.735089e-05 4.867545e-05
[583,] 9.999374e-01 1.251897e-04 6.259485e-05
[584,] 9.999052e-01 1.896410e-04 9.482048e-05
[585,] 9.999025e-01 1.950512e-04 9.752560e-05
[586,] 9.999286e-01 1.428414e-04 7.142069e-05
[587,] 9.999476e-01 1.047089e-04 5.235447e-05
[588,] 9.999325e-01 1.350460e-04 6.752302e-05
[589,] 9.998961e-01 2.077193e-04 1.038597e-04
[590,] 9.998868e-01 2.263099e-04 1.131550e-04
[591,] 9.998597e-01 2.805345e-04 1.402672e-04
[592,] 9.997693e-01 4.614670e-04 2.307335e-04
[593,] 9.997374e-01 5.252723e-04 2.626362e-04
[594,] 9.996838e-01 6.323764e-04 3.161882e-04
[595,] 9.995704e-01 8.591360e-04 4.295680e-04
[596,] 9.993597e-01 1.280608e-03 6.403038e-04
[597,] 9.990131e-01 1.973835e-03 9.869173e-04
[598,] 9.992033e-01 1.593456e-03 7.967279e-04
[599,] 9.993626e-01 1.274795e-03 6.373975e-04
[600,] 9.991755e-01 1.648912e-03 8.244560e-04
[601,] 9.986063e-01 2.787438e-03 1.393719e-03
[602,] 9.980732e-01 3.853589e-03 1.926795e-03
[603,] 9.974225e-01 5.154938e-03 2.577469e-03
[604,] 9.957215e-01 8.556936e-03 4.278468e-03
[605,] 9.947963e-01 1.040747e-02 5.203737e-03
[606,] 9.948553e-01 1.028949e-02 5.144743e-03
[607,] 9.940746e-01 1.185075e-02 5.925373e-03
[608,] 9.902361e-01 1.952771e-02 9.763855e-03
[609,] 9.848296e-01 3.034082e-02 1.517041e-02
[610,] 9.843346e-01 3.133075e-02 1.566538e-02
[611,] 9.842475e-01 3.150504e-02 1.575252e-02
[612,] 9.768146e-01 4.637085e-02 2.318542e-02
[613,] 9.610907e-01 7.781869e-02 3.890934e-02
[614,] 9.459296e-01 1.081408e-01 5.407038e-02
[615,] 9.411614e-01 1.176771e-01 5.883857e-02
[616,] 9.127655e-01 1.744689e-01 8.723447e-02
[617,] 8.613002e-01 2.773996e-01 1.386998e-01
[618,] 9.181320e-01 1.637359e-01 8.186796e-02
[619,] 9.530776e-01 9.384482e-02 4.692241e-02
[620,] 9.543043e-01 9.139146e-02 4.569573e-02
[621,] 9.010307e-01 1.979385e-01 9.896926e-02
[622,] 8.183321e-01 3.633359e-01 1.816679e-01
> postscript(file="/var/wessaorg/rcomp/tmp/139221321896144.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/wessaorg/rcomp/tmp/2u51d1321896144.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/wessaorg/rcomp/tmp/324jf1321896144.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/wessaorg/rcomp/tmp/4l3sa1321896144.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/wessaorg/rcomp/tmp/5g2r51321896144.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 = 633
Frequency = 1
1 2 3 4 5 6
5.666357339 7.672314205 7.988271071 7.874227937 6.880184802 6.216141668
7 8 9 10 11 12
4.752098534 4.418055399 5.414012265 7.019969131 3.600801989 4.626758855
13 14 15 16 17 18
5.222715721 6.498672587 7.334629452 7.470586318 6.126543184 4.642500050
19 20 21 22 23 24
3.948456915 3.634413781 5.440370647 6.496327512 2.967160371 3.773117237
25 26 27 28 29 30
4.649074103 6.355030968 7.620987834 7.456944700 6.312901565 4.308858431
31 32 33 34 35 36
2.824815297 2.760772163 4.196729028 5.652685894 2.023518753 3.059475618
37 38 39 40 41 42
4.165432484 5.371389350 6.737346216 6.193303081 5.269259947 3.745216813
43 44 45 46 47 48
2.021173679 2.867130544 3.853087410 5.029044276 1.589877134 2.465834000
49 50 51 52 53 54
3.861790866 5.077747731 5.723704597 5.529661463 4.825618329 2.911575194
55 56 57 58 59 60
2.037532060 1.463488926 2.999445792 4.275402657 0.946235516 1.552192382
61 62 63 64 65 66
2.598149247 4.394106113 5.510062979 5.006019845 3.541976710 1.837933576
67 68 69 70 71 72
0.543890442 0.589847307 1.985804173 3.441761039 0.332593897 1.108550763
73 74 75 76 77 78
2.044507629 3.360464495 4.066421360 3.982378226 2.798335092 1.324291958
79 80 81 82 83 84
-0.409751177 -0.053794311 1.212162555 2.398119420 -1.241047721 0.024909145
85 86 87 88 89 90
0.740866011 2.246822876 2.542779742 2.518736608 2.304693473 -0.019349661
91 92 93 94 95 96
-0.743392795 -0.987435929 0.848520936 1.664477802 -1.504689339 -0.268732474
97 98 99 100 101 102
0.797224392 2.543181258 2.949138124 2.955094989 1.871051855 0.107008721
103 104 105 106 107 108
-1.357034413 -1.631077548 0.314879318 1.880836184 -1.508330958 -0.802374092
109 110 111 112 113 114
0.003582774 1.649539639 2.495496505 1.851453371 0.577410237 -0.786632898
115 116 117 118 119 120
-2.130676032 -1.864719166 -0.188762300 1.317194565 -2.441972576 -1.596015710
121 122 123 124 125 126
-0.590058845 0.805898021 1.621854887 1.677811753 0.723768618 -1.120274516
127 128 129 130 131 132
-2.474317650 -2.368360785 -1.042403919 0.753552947 -2.455614195 -1.769657329
133 134 135 136 137 138
-0.283700463 1.002256403 1.948213268 1.634170134 1.020127000 -0.923916134
139 140 141 142 143 144
-1.947959269 -2.282002403 -0.946045537 0.589911328 -2.869255813 -1.643298947
145 146 147 148 149 150
-0.497342081 1.038614784 1.234571650 1.090528516 0.046485381 -1.347557753
151 152 153 154 155 156
-2.671600887 -2.345644021 -1.259687156 0.156269710 -3.172897431 -2.396940566
157 158 159 160 161 162
-1.630983700 -0.765026834 0.640930032 0.556886897 -0.407156237 -2.021199371
163 164 165 166 167 168
-3.855242505 -3.389285640 -1.883328774 -0.407371908 -4.016539050 -2.890582184
169 170 171 172 173 174
-2.504625318 -0.268668453 0.347288413 -0.366754721 -1.140797855 -2.604840990
175 176 177 178 179 180
-3.728884124 -3.332927258 -1.626970392 -0.311013527 -3.680180668 -2.404223802
181 182 183 184 185 186
-1.398266937 0.067689929 1.313646795 1.169603661 0.235560526 -1.058482608
187 188 189 190 191 192
-2.592525742 -2.656568877 -1.410612011 -0.664655145 -4.323822287 -2.697865421
193 194 195 196 197 198
-1.661908555 -0.225951689 0.540005176 -0.214037958 -1.088081092 -2.642124226
199 200 201 202 203 204
-4.136167361 -4.070210495 -2.674253629 -1.248296764 -4.437463905 -3.441507039
205 206 207 208 209 210
-2.735550173 -1.029593308 -0.143636442 -0.387679576 -1.781722711 -3.375765845
211 212 213 214 215 216
-4.419808979 -4.373852113 -3.127895248 -1.611938382 -4.901105523 -3.545148658
217 218 219 220 221 222
-2.569191792 -0.983234926 -0.707278060 -1.181321195 -2.015364329 -4.059407463
223 224 225 226 227 228
-5.473450598 -5.397493732 -3.721536866 -1.985580000 -5.314747142 -4.508790276
229 230 231 232 233 234
-2.502833410 -1.176876545 -0.120919679 -0.574962813 -1.619005947 -3.613049082
235 236 237 238 239 240
-4.687092216 -4.401135350 -2.995178484 -1.609221619 -4.358388760 -3.922431894
241 242 243 244 245 246
-2.296475029 -0.880518163 -0.384561297 -0.138604431 -1.372647566 -3.216690700
247 248 249 250 251 252
-4.990733834 -4.654776969 -3.028820103 -1.692863237 -4.682030379 -3.986073513
253 254 255 256 257 258
-2.230116647 -1.134159781 -0.808202916 -0.362246050 -1.696289184 -3.470332318
259 260 261 262 263 264
-4.834375453 -4.378418587 -2.972461721 -1.196504856 -4.535671997 -3.789715131
265 266 267 268 269 270
-1.763758265 -0.677801400 0.198155466 0.254112332 -1.399930803 -3.023973937
271 272 273 274 275 276
-4.108017071 -3.892060205 -2.466103340 -1.120146474 -4.419313615 -3.003356750
277 278 279 280 281 282
-1.677399884 -0.421443018 0.294513848 0.070470713 -1.473572421 -3.407615555
283 284 285 286 287 288
-4.691658690 -4.255701824 -2.539744958 -0.943788092 -4.302955234 -3.266998368
289 290 291 292 293 294
-1.821041502 -0.595084637 0.140872229 -0.473170905 -1.777214039 -3.541257174
295 296 297 298 299 300
-5.115300308 -4.739343442 -3.323386576 -1.627429711 -5.246596852 -3.730639986
301 302 303 304 305 306
-2.684683121 -0.788726255 0.157230611 -0.016812523 -1.140855658 -2.724898792
307 308 309 310 311 312
-4.538941926 -4.152985061 -2.657028195 -0.851071329 -4.240238471 -3.254281605
313 314 315 316 317 318
-2.288324739 -0.722367873 0.313588992 0.019545858 -1.144497276 -3.188540410
319 320 321 322 323 324
-4.742583545 -4.176626679 -2.550669813 -0.824712948 -4.473880089 -3.397923223
325 326 327 328 329 330
-1.531966357 -0.736009492 0.259947374 -0.034095760 -1.468138895 -2.992182029
331 332 333 334 335 336
-4.386225163 -4.240268297 -2.674311432 -1.018354566 -4.577521707 -3.781564842
337 338 339 340 341 342
-2.585607976 -0.629651110 0.416305756 0.032262621 -1.621780513 -3.245823647
343 344 345 346 347 348
-3.989866782 -4.263909916 -2.667953050 -1.181996184 -4.241163326 -3.625206460
349 350 351 352 353 354
-2.349249594 -0.523292729 0.462664137 0.338621003 -1.075422131 -2.959465266
355 356 357 358 359 360
-3.943508400 -3.567551534 -1.981594668 -0.485637803 -3.654804944 -2.038848078
361 362 363 364 365 366
-1.132891213 0.443065653 1.319022519 1.084979385 0.240936250 -1.603106884
367 368 369 370 371 372
-2.917150018 -2.551193153 -1.235236287 0.330720579 -2.608446563 -1.962489697
373 374 375 376 377 378
-1.156532831 0.899424035 1.435380900 1.111337766 0.057294632 -2.046748502
379 380 381 382 383 384
-3.370791637 -2.834834771 -1.408877905 0.257078960 -3.142088181 -1.786131315
385 386 387 388 389 390
-0.750174449 0.125782416 1.371739282 0.667696148 -0.506346987 -2.240390121
391 392 393 394 395 396
-3.594433255 -3.018476389 -1.302519524 0.363437342 -3.345729799 -2.279772934
397 398 399 400 401 402
-0.683816068 0.982140798 1.848097664 1.104054529 -0.659988605 -2.574031739
403 404 405 406 407 408
-3.888074874 -3.602118008 -1.856161142 -0.310204276 -3.579371418 -2.553414552
409 410 411 412 413 414
-1.397457686 0.048499179 0.854456045 0.860412911 -1.343630223 -3.027673358
415 416 417 418 419 420
-4.731716492 -4.055759626 -2.849802760 -1.343845895 -4.343013036 -3.667056170
421 422 423 424 425 426
-2.311099305 -1.135142439 0.030814427 -0.363228707 -2.147271842 -3.951314976
427 428 429 430 431 432
-5.015358110 -4.679401245 -3.173444379 -1.477487513 -4.496654655 -3.480697789
433 434 435 436 437 438
-2.344740923 -0.618784057 -0.002827192 -0.616870326 -1.750913460 -3.454956594
439 440 441 442 443 444
-4.998999729 -4.263042863 -2.477085997 -0.881129132 -4.260296273 -3.114339407
445 446 447 448 449 450
-1.858382541 -0.142425676 0.693531190 0.369488056 -0.844555079 -3.188598213
451 452 453 454 455 456
-3.932641347 -3.796684481 -2.110727616 -0.644770750 -3.603937891 -2.227981026
457 458 459 460 461 462
-0.992024160 -0.356067294 0.559889572 0.595846437 -0.528196697 -2.422239831
463 464 465 466 467 468
-3.896282966 -3.700326100 -2.154369234 -0.548412368 -4.067579510 -2.741622644
469 470 471 472 473 474
-2.095665778 -0.059708913 0.656247953 -0.177795181 -1.181838315 -3.125881450
475 476 477 478 479 480
-4.699924584 -3.993967718 -2.008010852 0.167946013 -3.321221128 -2.235264262
481 482 483 484 485 486
-0.839307397 0.896649469 2.072606335 1.788563200 0.864520066 -0.649523068
487 488 489 490 491 492
-2.363566202 -1.747609337 -0.301652471 1.484304395 -1.794862747 -0.668905881
493 494 495 496 497 498
0.227050985 1.803007851 1.958964716 1.664921582 0.930878448 -1.203164686
499 500 501 502 503 504
-3.027207821 -2.201250955 -0.595294089 1.020662776 -2.138504365 -1.622547499
505 506 507 508 509 510
-0.296590633 1.219366232 1.175323098 1.631279964 0.077236829 -1.166806305
511 512 513 514 515 516
-2.210849439 -1.874892573 -0.278935708 1.337021158 -2.182145983 -0.946189118
517 518 519 520 521 522
0.269767748 1.465724614 2.181681480 1.807638345 0.453595211 -1.200447923
523 524 525 526 527 528
-2.404491058 -1.878534192 -0.352577326 1.383379540 -1.715787602 -0.749830736
529 530 531 532 533 534
0.306126130 1.632082995 2.548039861 2.673996727 1.409953593 -0.464089542
535 536 537 538 539 540
-1.248132676 -1.202175810 0.753781056 2.539737921 -0.619429220 0.356527646
541 542 543 544 545 546
1.502484511 3.008441377 4.064398243 4.000355108 2.706311974 0.602268840
547 548 549 550 551 552
-0.291774294 -0.045817429 1.790139437 3.356096303 0.046929161 1.152886027
553 554 555 556 557 558
2.428842893 4.234799759 4.700756624 4.006713490 2.092670356 0.948627222
559 560 561 562 563 564
-0.755415913 -0.119459047 1.826497819 3.452454684 0.003287543 1.539244409
565 566 567 568 569 570
3.015201275 4.571158140 5.087115006 5.043071872 3.769028737 2.034985603
571 572 573 574 575 576
0.420942469 0.946899335 2.522856200 4.568813066 1.489645925 2.585602790
577 578 579 580 581 582
3.311559656 5.577516522 6.133473388 5.469430253 3.875387119 2.321343985
583 584 585 586 587 588
1.067300850 1.583257716 2.959214582 4.905171448 1.566004306 2.721961172
589 590 591 592 593 594
3.637918038 5.843874903 6.279831769 6.025788635 4.731745501 2.507702366
595 596 597 598 599 600
1.623659232 2.169616098 3.715572964 5.451529829 2.662362688 3.228319554
601 602 603 604 605 606
3.734276419 5.220233285 6.806190151 6.452147016 5.218103882 3.254060748
607 608 609 610 611 612
2.440017614 2.615974479 4.011931345 5.707888211 2.718721069 3.474677935
613 614 615 616 617 618
5.100634801 6.056591667 7.042548532 6.558505398 5.134462264 3.570419130
619 620 621 622 623 624
2.696375995 2.572332861 4.448289727 5.984246592 2.805079451 4.321036317
625 626 627 628 629 630
5.966993183 7.612950048 8.378906914 7.714863780 6.080820645 4.366777511
631 632 633
3.282734377 3.928691243 5.604648108
> postscript(file="/var/wessaorg/rcomp/tmp/6uk171321896144.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 = 633
Frequency = 1
lag(myerror, k = 1) myerror
0 5.666357339 NA
1 7.672314205 5.666357339
2 7.988271071 7.672314205
3 7.874227937 7.988271071
4 6.880184802 7.874227937
5 6.216141668 6.880184802
6 4.752098534 6.216141668
7 4.418055399 4.752098534
8 5.414012265 4.418055399
9 7.019969131 5.414012265
10 3.600801989 7.019969131
11 4.626758855 3.600801989
12 5.222715721 4.626758855
13 6.498672587 5.222715721
14 7.334629452 6.498672587
15 7.470586318 7.334629452
16 6.126543184 7.470586318
17 4.642500050 6.126543184
18 3.948456915 4.642500050
19 3.634413781 3.948456915
20 5.440370647 3.634413781
21 6.496327512 5.440370647
22 2.967160371 6.496327512
23 3.773117237 2.967160371
24 4.649074103 3.773117237
25 6.355030968 4.649074103
26 7.620987834 6.355030968
27 7.456944700 7.620987834
28 6.312901565 7.456944700
29 4.308858431 6.312901565
30 2.824815297 4.308858431
31 2.760772163 2.824815297
32 4.196729028 2.760772163
33 5.652685894 4.196729028
34 2.023518753 5.652685894
35 3.059475618 2.023518753
36 4.165432484 3.059475618
37 5.371389350 4.165432484
38 6.737346216 5.371389350
39 6.193303081 6.737346216
40 5.269259947 6.193303081
41 3.745216813 5.269259947
42 2.021173679 3.745216813
43 2.867130544 2.021173679
44 3.853087410 2.867130544
45 5.029044276 3.853087410
46 1.589877134 5.029044276
47 2.465834000 1.589877134
48 3.861790866 2.465834000
49 5.077747731 3.861790866
50 5.723704597 5.077747731
51 5.529661463 5.723704597
52 4.825618329 5.529661463
53 2.911575194 4.825618329
54 2.037532060 2.911575194
55 1.463488926 2.037532060
56 2.999445792 1.463488926
57 4.275402657 2.999445792
58 0.946235516 4.275402657
59 1.552192382 0.946235516
60 2.598149247 1.552192382
61 4.394106113 2.598149247
62 5.510062979 4.394106113
63 5.006019845 5.510062979
64 3.541976710 5.006019845
65 1.837933576 3.541976710
66 0.543890442 1.837933576
67 0.589847307 0.543890442
68 1.985804173 0.589847307
69 3.441761039 1.985804173
70 0.332593897 3.441761039
71 1.108550763 0.332593897
72 2.044507629 1.108550763
73 3.360464495 2.044507629
74 4.066421360 3.360464495
75 3.982378226 4.066421360
76 2.798335092 3.982378226
77 1.324291958 2.798335092
78 -0.409751177 1.324291958
79 -0.053794311 -0.409751177
80 1.212162555 -0.053794311
81 2.398119420 1.212162555
82 -1.241047721 2.398119420
83 0.024909145 -1.241047721
84 0.740866011 0.024909145
85 2.246822876 0.740866011
86 2.542779742 2.246822876
87 2.518736608 2.542779742
88 2.304693473 2.518736608
89 -0.019349661 2.304693473
90 -0.743392795 -0.019349661
91 -0.987435929 -0.743392795
92 0.848520936 -0.987435929
93 1.664477802 0.848520936
94 -1.504689339 1.664477802
95 -0.268732474 -1.504689339
96 0.797224392 -0.268732474
97 2.543181258 0.797224392
98 2.949138124 2.543181258
99 2.955094989 2.949138124
100 1.871051855 2.955094989
101 0.107008721 1.871051855
102 -1.357034413 0.107008721
103 -1.631077548 -1.357034413
104 0.314879318 -1.631077548
105 1.880836184 0.314879318
106 -1.508330958 1.880836184
107 -0.802374092 -1.508330958
108 0.003582774 -0.802374092
109 1.649539639 0.003582774
110 2.495496505 1.649539639
111 1.851453371 2.495496505
112 0.577410237 1.851453371
113 -0.786632898 0.577410237
114 -2.130676032 -0.786632898
115 -1.864719166 -2.130676032
116 -0.188762300 -1.864719166
117 1.317194565 -0.188762300
118 -2.441972576 1.317194565
119 -1.596015710 -2.441972576
120 -0.590058845 -1.596015710
121 0.805898021 -0.590058845
122 1.621854887 0.805898021
123 1.677811753 1.621854887
124 0.723768618 1.677811753
125 -1.120274516 0.723768618
126 -2.474317650 -1.120274516
127 -2.368360785 -2.474317650
128 -1.042403919 -2.368360785
129 0.753552947 -1.042403919
130 -2.455614195 0.753552947
131 -1.769657329 -2.455614195
132 -0.283700463 -1.769657329
133 1.002256403 -0.283700463
134 1.948213268 1.002256403
135 1.634170134 1.948213268
136 1.020127000 1.634170134
137 -0.923916134 1.020127000
138 -1.947959269 -0.923916134
139 -2.282002403 -1.947959269
140 -0.946045537 -2.282002403
141 0.589911328 -0.946045537
142 -2.869255813 0.589911328
143 -1.643298947 -2.869255813
144 -0.497342081 -1.643298947
145 1.038614784 -0.497342081
146 1.234571650 1.038614784
147 1.090528516 1.234571650
148 0.046485381 1.090528516
149 -1.347557753 0.046485381
150 -2.671600887 -1.347557753
151 -2.345644021 -2.671600887
152 -1.259687156 -2.345644021
153 0.156269710 -1.259687156
154 -3.172897431 0.156269710
155 -2.396940566 -3.172897431
156 -1.630983700 -2.396940566
157 -0.765026834 -1.630983700
158 0.640930032 -0.765026834
159 0.556886897 0.640930032
160 -0.407156237 0.556886897
161 -2.021199371 -0.407156237
162 -3.855242505 -2.021199371
163 -3.389285640 -3.855242505
164 -1.883328774 -3.389285640
165 -0.407371908 -1.883328774
166 -4.016539050 -0.407371908
167 -2.890582184 -4.016539050
168 -2.504625318 -2.890582184
169 -0.268668453 -2.504625318
170 0.347288413 -0.268668453
171 -0.366754721 0.347288413
172 -1.140797855 -0.366754721
173 -2.604840990 -1.140797855
174 -3.728884124 -2.604840990
175 -3.332927258 -3.728884124
176 -1.626970392 -3.332927258
177 -0.311013527 -1.626970392
178 -3.680180668 -0.311013527
179 -2.404223802 -3.680180668
180 -1.398266937 -2.404223802
181 0.067689929 -1.398266937
182 1.313646795 0.067689929
183 1.169603661 1.313646795
184 0.235560526 1.169603661
185 -1.058482608 0.235560526
186 -2.592525742 -1.058482608
187 -2.656568877 -2.592525742
188 -1.410612011 -2.656568877
189 -0.664655145 -1.410612011
190 -4.323822287 -0.664655145
191 -2.697865421 -4.323822287
192 -1.661908555 -2.697865421
193 -0.225951689 -1.661908555
194 0.540005176 -0.225951689
195 -0.214037958 0.540005176
196 -1.088081092 -0.214037958
197 -2.642124226 -1.088081092
198 -4.136167361 -2.642124226
199 -4.070210495 -4.136167361
200 -2.674253629 -4.070210495
201 -1.248296764 -2.674253629
202 -4.437463905 -1.248296764
203 -3.441507039 -4.437463905
204 -2.735550173 -3.441507039
205 -1.029593308 -2.735550173
206 -0.143636442 -1.029593308
207 -0.387679576 -0.143636442
208 -1.781722711 -0.387679576
209 -3.375765845 -1.781722711
210 -4.419808979 -3.375765845
211 -4.373852113 -4.419808979
212 -3.127895248 -4.373852113
213 -1.611938382 -3.127895248
214 -4.901105523 -1.611938382
215 -3.545148658 -4.901105523
216 -2.569191792 -3.545148658
217 -0.983234926 -2.569191792
218 -0.707278060 -0.983234926
219 -1.181321195 -0.707278060
220 -2.015364329 -1.181321195
221 -4.059407463 -2.015364329
222 -5.473450598 -4.059407463
223 -5.397493732 -5.473450598
224 -3.721536866 -5.397493732
225 -1.985580000 -3.721536866
226 -5.314747142 -1.985580000
227 -4.508790276 -5.314747142
228 -2.502833410 -4.508790276
229 -1.176876545 -2.502833410
230 -0.120919679 -1.176876545
231 -0.574962813 -0.120919679
232 -1.619005947 -0.574962813
233 -3.613049082 -1.619005947
234 -4.687092216 -3.613049082
235 -4.401135350 -4.687092216
236 -2.995178484 -4.401135350
237 -1.609221619 -2.995178484
238 -4.358388760 -1.609221619
239 -3.922431894 -4.358388760
240 -2.296475029 -3.922431894
241 -0.880518163 -2.296475029
242 -0.384561297 -0.880518163
243 -0.138604431 -0.384561297
244 -1.372647566 -0.138604431
245 -3.216690700 -1.372647566
246 -4.990733834 -3.216690700
247 -4.654776969 -4.990733834
248 -3.028820103 -4.654776969
249 -1.692863237 -3.028820103
250 -4.682030379 -1.692863237
251 -3.986073513 -4.682030379
252 -2.230116647 -3.986073513
253 -1.134159781 -2.230116647
254 -0.808202916 -1.134159781
255 -0.362246050 -0.808202916
256 -1.696289184 -0.362246050
257 -3.470332318 -1.696289184
258 -4.834375453 -3.470332318
259 -4.378418587 -4.834375453
260 -2.972461721 -4.378418587
261 -1.196504856 -2.972461721
262 -4.535671997 -1.196504856
263 -3.789715131 -4.535671997
264 -1.763758265 -3.789715131
265 -0.677801400 -1.763758265
266 0.198155466 -0.677801400
267 0.254112332 0.198155466
268 -1.399930803 0.254112332
269 -3.023973937 -1.399930803
270 -4.108017071 -3.023973937
271 -3.892060205 -4.108017071
272 -2.466103340 -3.892060205
273 -1.120146474 -2.466103340
274 -4.419313615 -1.120146474
275 -3.003356750 -4.419313615
276 -1.677399884 -3.003356750
277 -0.421443018 -1.677399884
278 0.294513848 -0.421443018
279 0.070470713 0.294513848
280 -1.473572421 0.070470713
281 -3.407615555 -1.473572421
282 -4.691658690 -3.407615555
283 -4.255701824 -4.691658690
284 -2.539744958 -4.255701824
285 -0.943788092 -2.539744958
286 -4.302955234 -0.943788092
287 -3.266998368 -4.302955234
288 -1.821041502 -3.266998368
289 -0.595084637 -1.821041502
290 0.140872229 -0.595084637
291 -0.473170905 0.140872229
292 -1.777214039 -0.473170905
293 -3.541257174 -1.777214039
294 -5.115300308 -3.541257174
295 -4.739343442 -5.115300308
296 -3.323386576 -4.739343442
297 -1.627429711 -3.323386576
298 -5.246596852 -1.627429711
299 -3.730639986 -5.246596852
300 -2.684683121 -3.730639986
301 -0.788726255 -2.684683121
302 0.157230611 -0.788726255
303 -0.016812523 0.157230611
304 -1.140855658 -0.016812523
305 -2.724898792 -1.140855658
306 -4.538941926 -2.724898792
307 -4.152985061 -4.538941926
308 -2.657028195 -4.152985061
309 -0.851071329 -2.657028195
310 -4.240238471 -0.851071329
311 -3.254281605 -4.240238471
312 -2.288324739 -3.254281605
313 -0.722367873 -2.288324739
314 0.313588992 -0.722367873
315 0.019545858 0.313588992
316 -1.144497276 0.019545858
317 -3.188540410 -1.144497276
318 -4.742583545 -3.188540410
319 -4.176626679 -4.742583545
320 -2.550669813 -4.176626679
321 -0.824712948 -2.550669813
322 -4.473880089 -0.824712948
323 -3.397923223 -4.473880089
324 -1.531966357 -3.397923223
325 -0.736009492 -1.531966357
326 0.259947374 -0.736009492
327 -0.034095760 0.259947374
328 -1.468138895 -0.034095760
329 -2.992182029 -1.468138895
330 -4.386225163 -2.992182029
331 -4.240268297 -4.386225163
332 -2.674311432 -4.240268297
333 -1.018354566 -2.674311432
334 -4.577521707 -1.018354566
335 -3.781564842 -4.577521707
336 -2.585607976 -3.781564842
337 -0.629651110 -2.585607976
338 0.416305756 -0.629651110
339 0.032262621 0.416305756
340 -1.621780513 0.032262621
341 -3.245823647 -1.621780513
342 -3.989866782 -3.245823647
343 -4.263909916 -3.989866782
344 -2.667953050 -4.263909916
345 -1.181996184 -2.667953050
346 -4.241163326 -1.181996184
347 -3.625206460 -4.241163326
348 -2.349249594 -3.625206460
349 -0.523292729 -2.349249594
350 0.462664137 -0.523292729
351 0.338621003 0.462664137
352 -1.075422131 0.338621003
353 -2.959465266 -1.075422131
354 -3.943508400 -2.959465266
355 -3.567551534 -3.943508400
356 -1.981594668 -3.567551534
357 -0.485637803 -1.981594668
358 -3.654804944 -0.485637803
359 -2.038848078 -3.654804944
360 -1.132891213 -2.038848078
361 0.443065653 -1.132891213
362 1.319022519 0.443065653
363 1.084979385 1.319022519
364 0.240936250 1.084979385
365 -1.603106884 0.240936250
366 -2.917150018 -1.603106884
367 -2.551193153 -2.917150018
368 -1.235236287 -2.551193153
369 0.330720579 -1.235236287
370 -2.608446563 0.330720579
371 -1.962489697 -2.608446563
372 -1.156532831 -1.962489697
373 0.899424035 -1.156532831
374 1.435380900 0.899424035
375 1.111337766 1.435380900
376 0.057294632 1.111337766
377 -2.046748502 0.057294632
378 -3.370791637 -2.046748502
379 -2.834834771 -3.370791637
380 -1.408877905 -2.834834771
381 0.257078960 -1.408877905
382 -3.142088181 0.257078960
383 -1.786131315 -3.142088181
384 -0.750174449 -1.786131315
385 0.125782416 -0.750174449
386 1.371739282 0.125782416
387 0.667696148 1.371739282
388 -0.506346987 0.667696148
389 -2.240390121 -0.506346987
390 -3.594433255 -2.240390121
391 -3.018476389 -3.594433255
392 -1.302519524 -3.018476389
393 0.363437342 -1.302519524
394 -3.345729799 0.363437342
395 -2.279772934 -3.345729799
396 -0.683816068 -2.279772934
397 0.982140798 -0.683816068
398 1.848097664 0.982140798
399 1.104054529 1.848097664
400 -0.659988605 1.104054529
401 -2.574031739 -0.659988605
402 -3.888074874 -2.574031739
403 -3.602118008 -3.888074874
404 -1.856161142 -3.602118008
405 -0.310204276 -1.856161142
406 -3.579371418 -0.310204276
407 -2.553414552 -3.579371418
408 -1.397457686 -2.553414552
409 0.048499179 -1.397457686
410 0.854456045 0.048499179
411 0.860412911 0.854456045
412 -1.343630223 0.860412911
413 -3.027673358 -1.343630223
414 -4.731716492 -3.027673358
415 -4.055759626 -4.731716492
416 -2.849802760 -4.055759626
417 -1.343845895 -2.849802760
418 -4.343013036 -1.343845895
419 -3.667056170 -4.343013036
420 -2.311099305 -3.667056170
421 -1.135142439 -2.311099305
422 0.030814427 -1.135142439
423 -0.363228707 0.030814427
424 -2.147271842 -0.363228707
425 -3.951314976 -2.147271842
426 -5.015358110 -3.951314976
427 -4.679401245 -5.015358110
428 -3.173444379 -4.679401245
429 -1.477487513 -3.173444379
430 -4.496654655 -1.477487513
431 -3.480697789 -4.496654655
432 -2.344740923 -3.480697789
433 -0.618784057 -2.344740923
434 -0.002827192 -0.618784057
435 -0.616870326 -0.002827192
436 -1.750913460 -0.616870326
437 -3.454956594 -1.750913460
438 -4.998999729 -3.454956594
439 -4.263042863 -4.998999729
440 -2.477085997 -4.263042863
441 -0.881129132 -2.477085997
442 -4.260296273 -0.881129132
443 -3.114339407 -4.260296273
444 -1.858382541 -3.114339407
445 -0.142425676 -1.858382541
446 0.693531190 -0.142425676
447 0.369488056 0.693531190
448 -0.844555079 0.369488056
449 -3.188598213 -0.844555079
450 -3.932641347 -3.188598213
451 -3.796684481 -3.932641347
452 -2.110727616 -3.796684481
453 -0.644770750 -2.110727616
454 -3.603937891 -0.644770750
455 -2.227981026 -3.603937891
456 -0.992024160 -2.227981026
457 -0.356067294 -0.992024160
458 0.559889572 -0.356067294
459 0.595846437 0.559889572
460 -0.528196697 0.595846437
461 -2.422239831 -0.528196697
462 -3.896282966 -2.422239831
463 -3.700326100 -3.896282966
464 -2.154369234 -3.700326100
465 -0.548412368 -2.154369234
466 -4.067579510 -0.548412368
467 -2.741622644 -4.067579510
468 -2.095665778 -2.741622644
469 -0.059708913 -2.095665778
470 0.656247953 -0.059708913
471 -0.177795181 0.656247953
472 -1.181838315 -0.177795181
473 -3.125881450 -1.181838315
474 -4.699924584 -3.125881450
475 -3.993967718 -4.699924584
476 -2.008010852 -3.993967718
477 0.167946013 -2.008010852
478 -3.321221128 0.167946013
479 -2.235264262 -3.321221128
480 -0.839307397 -2.235264262
481 0.896649469 -0.839307397
482 2.072606335 0.896649469
483 1.788563200 2.072606335
484 0.864520066 1.788563200
485 -0.649523068 0.864520066
486 -2.363566202 -0.649523068
487 -1.747609337 -2.363566202
488 -0.301652471 -1.747609337
489 1.484304395 -0.301652471
490 -1.794862747 1.484304395
491 -0.668905881 -1.794862747
492 0.227050985 -0.668905881
493 1.803007851 0.227050985
494 1.958964716 1.803007851
495 1.664921582 1.958964716
496 0.930878448 1.664921582
497 -1.203164686 0.930878448
498 -3.027207821 -1.203164686
499 -2.201250955 -3.027207821
500 -0.595294089 -2.201250955
501 1.020662776 -0.595294089
502 -2.138504365 1.020662776
503 -1.622547499 -2.138504365
504 -0.296590633 -1.622547499
505 1.219366232 -0.296590633
506 1.175323098 1.219366232
507 1.631279964 1.175323098
508 0.077236829 1.631279964
509 -1.166806305 0.077236829
510 -2.210849439 -1.166806305
511 -1.874892573 -2.210849439
512 -0.278935708 -1.874892573
513 1.337021158 -0.278935708
514 -2.182145983 1.337021158
515 -0.946189118 -2.182145983
516 0.269767748 -0.946189118
517 1.465724614 0.269767748
518 2.181681480 1.465724614
519 1.807638345 2.181681480
520 0.453595211 1.807638345
521 -1.200447923 0.453595211
522 -2.404491058 -1.200447923
523 -1.878534192 -2.404491058
524 -0.352577326 -1.878534192
525 1.383379540 -0.352577326
526 -1.715787602 1.383379540
527 -0.749830736 -1.715787602
528 0.306126130 -0.749830736
529 1.632082995 0.306126130
530 2.548039861 1.632082995
531 2.673996727 2.548039861
532 1.409953593 2.673996727
533 -0.464089542 1.409953593
534 -1.248132676 -0.464089542
535 -1.202175810 -1.248132676
536 0.753781056 -1.202175810
537 2.539737921 0.753781056
538 -0.619429220 2.539737921
539 0.356527646 -0.619429220
540 1.502484511 0.356527646
541 3.008441377 1.502484511
542 4.064398243 3.008441377
543 4.000355108 4.064398243
544 2.706311974 4.000355108
545 0.602268840 2.706311974
546 -0.291774294 0.602268840
547 -0.045817429 -0.291774294
548 1.790139437 -0.045817429
549 3.356096303 1.790139437
550 0.046929161 3.356096303
551 1.152886027 0.046929161
552 2.428842893 1.152886027
553 4.234799759 2.428842893
554 4.700756624 4.234799759
555 4.006713490 4.700756624
556 2.092670356 4.006713490
557 0.948627222 2.092670356
558 -0.755415913 0.948627222
559 -0.119459047 -0.755415913
560 1.826497819 -0.119459047
561 3.452454684 1.826497819
562 0.003287543 3.452454684
563 1.539244409 0.003287543
564 3.015201275 1.539244409
565 4.571158140 3.015201275
566 5.087115006 4.571158140
567 5.043071872 5.087115006
568 3.769028737 5.043071872
569 2.034985603 3.769028737
570 0.420942469 2.034985603
571 0.946899335 0.420942469
572 2.522856200 0.946899335
573 4.568813066 2.522856200
574 1.489645925 4.568813066
575 2.585602790 1.489645925
576 3.311559656 2.585602790
577 5.577516522 3.311559656
578 6.133473388 5.577516522
579 5.469430253 6.133473388
580 3.875387119 5.469430253
581 2.321343985 3.875387119
582 1.067300850 2.321343985
583 1.583257716 1.067300850
584 2.959214582 1.583257716
585 4.905171448 2.959214582
586 1.566004306 4.905171448
587 2.721961172 1.566004306
588 3.637918038 2.721961172
589 5.843874903 3.637918038
590 6.279831769 5.843874903
591 6.025788635 6.279831769
592 4.731745501 6.025788635
593 2.507702366 4.731745501
594 1.623659232 2.507702366
595 2.169616098 1.623659232
596 3.715572964 2.169616098
597 5.451529829 3.715572964
598 2.662362688 5.451529829
599 3.228319554 2.662362688
600 3.734276419 3.228319554
601 5.220233285 3.734276419
602 6.806190151 5.220233285
603 6.452147016 6.806190151
604 5.218103882 6.452147016
605 3.254060748 5.218103882
606 2.440017614 3.254060748
607 2.615974479 2.440017614
608 4.011931345 2.615974479
609 5.707888211 4.011931345
610 2.718721069 5.707888211
611 3.474677935 2.718721069
612 5.100634801 3.474677935
613 6.056591667 5.100634801
614 7.042548532 6.056591667
615 6.558505398 7.042548532
616 5.134462264 6.558505398
617 3.570419130 5.134462264
618 2.696375995 3.570419130
619 2.572332861 2.696375995
620 4.448289727 2.572332861
621 5.984246592 4.448289727
622 2.805079451 5.984246592
623 4.321036317 2.805079451
624 5.966993183 4.321036317
625 7.612950048 5.966993183
626 8.378906914 7.612950048
627 7.714863780 8.378906914
628 6.080820645 7.714863780
629 4.366777511 6.080820645
630 3.282734377 4.366777511
631 3.928691243 3.282734377
632 5.604648108 3.928691243
633 NA 5.604648108
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 7.672314205 5.666357339
[2,] 7.988271071 7.672314205
[3,] 7.874227937 7.988271071
[4,] 6.880184802 7.874227937
[5,] 6.216141668 6.880184802
[6,] 4.752098534 6.216141668
[7,] 4.418055399 4.752098534
[8,] 5.414012265 4.418055399
[9,] 7.019969131 5.414012265
[10,] 3.600801989 7.019969131
[11,] 4.626758855 3.600801989
[12,] 5.222715721 4.626758855
[13,] 6.498672587 5.222715721
[14,] 7.334629452 6.498672587
[15,] 7.470586318 7.334629452
[16,] 6.126543184 7.470586318
[17,] 4.642500050 6.126543184
[18,] 3.948456915 4.642500050
[19,] 3.634413781 3.948456915
[20,] 5.440370647 3.634413781
[21,] 6.496327512 5.440370647
[22,] 2.967160371 6.496327512
[23,] 3.773117237 2.967160371
[24,] 4.649074103 3.773117237
[25,] 6.355030968 4.649074103
[26,] 7.620987834 6.355030968
[27,] 7.456944700 7.620987834
[28,] 6.312901565 7.456944700
[29,] 4.308858431 6.312901565
[30,] 2.824815297 4.308858431
[31,] 2.760772163 2.824815297
[32,] 4.196729028 2.760772163
[33,] 5.652685894 4.196729028
[34,] 2.023518753 5.652685894
[35,] 3.059475618 2.023518753
[36,] 4.165432484 3.059475618
[37,] 5.371389350 4.165432484
[38,] 6.737346216 5.371389350
[39,] 6.193303081 6.737346216
[40,] 5.269259947 6.193303081
[41,] 3.745216813 5.269259947
[42,] 2.021173679 3.745216813
[43,] 2.867130544 2.021173679
[44,] 3.853087410 2.867130544
[45,] 5.029044276 3.853087410
[46,] 1.589877134 5.029044276
[47,] 2.465834000 1.589877134
[48,] 3.861790866 2.465834000
[49,] 5.077747731 3.861790866
[50,] 5.723704597 5.077747731
[51,] 5.529661463 5.723704597
[52,] 4.825618329 5.529661463
[53,] 2.911575194 4.825618329
[54,] 2.037532060 2.911575194
[55,] 1.463488926 2.037532060
[56,] 2.999445792 1.463488926
[57,] 4.275402657 2.999445792
[58,] 0.946235516 4.275402657
[59,] 1.552192382 0.946235516
[60,] 2.598149247 1.552192382
[61,] 4.394106113 2.598149247
[62,] 5.510062979 4.394106113
[63,] 5.006019845 5.510062979
[64,] 3.541976710 5.006019845
[65,] 1.837933576 3.541976710
[66,] 0.543890442 1.837933576
[67,] 0.589847307 0.543890442
[68,] 1.985804173 0.589847307
[69,] 3.441761039 1.985804173
[70,] 0.332593897 3.441761039
[71,] 1.108550763 0.332593897
[72,] 2.044507629 1.108550763
[73,] 3.360464495 2.044507629
[74,] 4.066421360 3.360464495
[75,] 3.982378226 4.066421360
[76,] 2.798335092 3.982378226
[77,] 1.324291958 2.798335092
[78,] -0.409751177 1.324291958
[79,] -0.053794311 -0.409751177
[80,] 1.212162555 -0.053794311
[81,] 2.398119420 1.212162555
[82,] -1.241047721 2.398119420
[83,] 0.024909145 -1.241047721
[84,] 0.740866011 0.024909145
[85,] 2.246822876 0.740866011
[86,] 2.542779742 2.246822876
[87,] 2.518736608 2.542779742
[88,] 2.304693473 2.518736608
[89,] -0.019349661 2.304693473
[90,] -0.743392795 -0.019349661
[91,] -0.987435929 -0.743392795
[92,] 0.848520936 -0.987435929
[93,] 1.664477802 0.848520936
[94,] -1.504689339 1.664477802
[95,] -0.268732474 -1.504689339
[96,] 0.797224392 -0.268732474
[97,] 2.543181258 0.797224392
[98,] 2.949138124 2.543181258
[99,] 2.955094989 2.949138124
[100,] 1.871051855 2.955094989
[101,] 0.107008721 1.871051855
[102,] -1.357034413 0.107008721
[103,] -1.631077548 -1.357034413
[104,] 0.314879318 -1.631077548
[105,] 1.880836184 0.314879318
[106,] -1.508330958 1.880836184
[107,] -0.802374092 -1.508330958
[108,] 0.003582774 -0.802374092
[109,] 1.649539639 0.003582774
[110,] 2.495496505 1.649539639
[111,] 1.851453371 2.495496505
[112,] 0.577410237 1.851453371
[113,] -0.786632898 0.577410237
[114,] -2.130676032 -0.786632898
[115,] -1.864719166 -2.130676032
[116,] -0.188762300 -1.864719166
[117,] 1.317194565 -0.188762300
[118,] -2.441972576 1.317194565
[119,] -1.596015710 -2.441972576
[120,] -0.590058845 -1.596015710
[121,] 0.805898021 -0.590058845
[122,] 1.621854887 0.805898021
[123,] 1.677811753 1.621854887
[124,] 0.723768618 1.677811753
[125,] -1.120274516 0.723768618
[126,] -2.474317650 -1.120274516
[127,] -2.368360785 -2.474317650
[128,] -1.042403919 -2.368360785
[129,] 0.753552947 -1.042403919
[130,] -2.455614195 0.753552947
[131,] -1.769657329 -2.455614195
[132,] -0.283700463 -1.769657329
[133,] 1.002256403 -0.283700463
[134,] 1.948213268 1.002256403
[135,] 1.634170134 1.948213268
[136,] 1.020127000 1.634170134
[137,] -0.923916134 1.020127000
[138,] -1.947959269 -0.923916134
[139,] -2.282002403 -1.947959269
[140,] -0.946045537 -2.282002403
[141,] 0.589911328 -0.946045537
[142,] -2.869255813 0.589911328
[143,] -1.643298947 -2.869255813
[144,] -0.497342081 -1.643298947
[145,] 1.038614784 -0.497342081
[146,] 1.234571650 1.038614784
[147,] 1.090528516 1.234571650
[148,] 0.046485381 1.090528516
[149,] -1.347557753 0.046485381
[150,] -2.671600887 -1.347557753
[151,] -2.345644021 -2.671600887
[152,] -1.259687156 -2.345644021
[153,] 0.156269710 -1.259687156
[154,] -3.172897431 0.156269710
[155,] -2.396940566 -3.172897431
[156,] -1.630983700 -2.396940566
[157,] -0.765026834 -1.630983700
[158,] 0.640930032 -0.765026834
[159,] 0.556886897 0.640930032
[160,] -0.407156237 0.556886897
[161,] -2.021199371 -0.407156237
[162,] -3.855242505 -2.021199371
[163,] -3.389285640 -3.855242505
[164,] -1.883328774 -3.389285640
[165,] -0.407371908 -1.883328774
[166,] -4.016539050 -0.407371908
[167,] -2.890582184 -4.016539050
[168,] -2.504625318 -2.890582184
[169,] -0.268668453 -2.504625318
[170,] 0.347288413 -0.268668453
[171,] -0.366754721 0.347288413
[172,] -1.140797855 -0.366754721
[173,] -2.604840990 -1.140797855
[174,] -3.728884124 -2.604840990
[175,] -3.332927258 -3.728884124
[176,] -1.626970392 -3.332927258
[177,] -0.311013527 -1.626970392
[178,] -3.680180668 -0.311013527
[179,] -2.404223802 -3.680180668
[180,] -1.398266937 -2.404223802
[181,] 0.067689929 -1.398266937
[182,] 1.313646795 0.067689929
[183,] 1.169603661 1.313646795
[184,] 0.235560526 1.169603661
[185,] -1.058482608 0.235560526
[186,] -2.592525742 -1.058482608
[187,] -2.656568877 -2.592525742
[188,] -1.410612011 -2.656568877
[189,] -0.664655145 -1.410612011
[190,] -4.323822287 -0.664655145
[191,] -2.697865421 -4.323822287
[192,] -1.661908555 -2.697865421
[193,] -0.225951689 -1.661908555
[194,] 0.540005176 -0.225951689
[195,] -0.214037958 0.540005176
[196,] -1.088081092 -0.214037958
[197,] -2.642124226 -1.088081092
[198,] -4.136167361 -2.642124226
[199,] -4.070210495 -4.136167361
[200,] -2.674253629 -4.070210495
[201,] -1.248296764 -2.674253629
[202,] -4.437463905 -1.248296764
[203,] -3.441507039 -4.437463905
[204,] -2.735550173 -3.441507039
[205,] -1.029593308 -2.735550173
[206,] -0.143636442 -1.029593308
[207,] -0.387679576 -0.143636442
[208,] -1.781722711 -0.387679576
[209,] -3.375765845 -1.781722711
[210,] -4.419808979 -3.375765845
[211,] -4.373852113 -4.419808979
[212,] -3.127895248 -4.373852113
[213,] -1.611938382 -3.127895248
[214,] -4.901105523 -1.611938382
[215,] -3.545148658 -4.901105523
[216,] -2.569191792 -3.545148658
[217,] -0.983234926 -2.569191792
[218,] -0.707278060 -0.983234926
[219,] -1.181321195 -0.707278060
[220,] -2.015364329 -1.181321195
[221,] -4.059407463 -2.015364329
[222,] -5.473450598 -4.059407463
[223,] -5.397493732 -5.473450598
[224,] -3.721536866 -5.397493732
[225,] -1.985580000 -3.721536866
[226,] -5.314747142 -1.985580000
[227,] -4.508790276 -5.314747142
[228,] -2.502833410 -4.508790276
[229,] -1.176876545 -2.502833410
[230,] -0.120919679 -1.176876545
[231,] -0.574962813 -0.120919679
[232,] -1.619005947 -0.574962813
[233,] -3.613049082 -1.619005947
[234,] -4.687092216 -3.613049082
[235,] -4.401135350 -4.687092216
[236,] -2.995178484 -4.401135350
[237,] -1.609221619 -2.995178484
[238,] -4.358388760 -1.609221619
[239,] -3.922431894 -4.358388760
[240,] -2.296475029 -3.922431894
[241,] -0.880518163 -2.296475029
[242,] -0.384561297 -0.880518163
[243,] -0.138604431 -0.384561297
[244,] -1.372647566 -0.138604431
[245,] -3.216690700 -1.372647566
[246,] -4.990733834 -3.216690700
[247,] -4.654776969 -4.990733834
[248,] -3.028820103 -4.654776969
[249,] -1.692863237 -3.028820103
[250,] -4.682030379 -1.692863237
[251,] -3.986073513 -4.682030379
[252,] -2.230116647 -3.986073513
[253,] -1.134159781 -2.230116647
[254,] -0.808202916 -1.134159781
[255,] -0.362246050 -0.808202916
[256,] -1.696289184 -0.362246050
[257,] -3.470332318 -1.696289184
[258,] -4.834375453 -3.470332318
[259,] -4.378418587 -4.834375453
[260,] -2.972461721 -4.378418587
[261,] -1.196504856 -2.972461721
[262,] -4.535671997 -1.196504856
[263,] -3.789715131 -4.535671997
[264,] -1.763758265 -3.789715131
[265,] -0.677801400 -1.763758265
[266,] 0.198155466 -0.677801400
[267,] 0.254112332 0.198155466
[268,] -1.399930803 0.254112332
[269,] -3.023973937 -1.399930803
[270,] -4.108017071 -3.023973937
[271,] -3.892060205 -4.108017071
[272,] -2.466103340 -3.892060205
[273,] -1.120146474 -2.466103340
[274,] -4.419313615 -1.120146474
[275,] -3.003356750 -4.419313615
[276,] -1.677399884 -3.003356750
[277,] -0.421443018 -1.677399884
[278,] 0.294513848 -0.421443018
[279,] 0.070470713 0.294513848
[280,] -1.473572421 0.070470713
[281,] -3.407615555 -1.473572421
[282,] -4.691658690 -3.407615555
[283,] -4.255701824 -4.691658690
[284,] -2.539744958 -4.255701824
[285,] -0.943788092 -2.539744958
[286,] -4.302955234 -0.943788092
[287,] -3.266998368 -4.302955234
[288,] -1.821041502 -3.266998368
[289,] -0.595084637 -1.821041502
[290,] 0.140872229 -0.595084637
[291,] -0.473170905 0.140872229
[292,] -1.777214039 -0.473170905
[293,] -3.541257174 -1.777214039
[294,] -5.115300308 -3.541257174
[295,] -4.739343442 -5.115300308
[296,] -3.323386576 -4.739343442
[297,] -1.627429711 -3.323386576
[298,] -5.246596852 -1.627429711
[299,] -3.730639986 -5.246596852
[300,] -2.684683121 -3.730639986
[301,] -0.788726255 -2.684683121
[302,] 0.157230611 -0.788726255
[303,] -0.016812523 0.157230611
[304,] -1.140855658 -0.016812523
[305,] -2.724898792 -1.140855658
[306,] -4.538941926 -2.724898792
[307,] -4.152985061 -4.538941926
[308,] -2.657028195 -4.152985061
[309,] -0.851071329 -2.657028195
[310,] -4.240238471 -0.851071329
[311,] -3.254281605 -4.240238471
[312,] -2.288324739 -3.254281605
[313,] -0.722367873 -2.288324739
[314,] 0.313588992 -0.722367873
[315,] 0.019545858 0.313588992
[316,] -1.144497276 0.019545858
[317,] -3.188540410 -1.144497276
[318,] -4.742583545 -3.188540410
[319,] -4.176626679 -4.742583545
[320,] -2.550669813 -4.176626679
[321,] -0.824712948 -2.550669813
[322,] -4.473880089 -0.824712948
[323,] -3.397923223 -4.473880089
[324,] -1.531966357 -3.397923223
[325,] -0.736009492 -1.531966357
[326,] 0.259947374 -0.736009492
[327,] -0.034095760 0.259947374
[328,] -1.468138895 -0.034095760
[329,] -2.992182029 -1.468138895
[330,] -4.386225163 -2.992182029
[331,] -4.240268297 -4.386225163
[332,] -2.674311432 -4.240268297
[333,] -1.018354566 -2.674311432
[334,] -4.577521707 -1.018354566
[335,] -3.781564842 -4.577521707
[336,] -2.585607976 -3.781564842
[337,] -0.629651110 -2.585607976
[338,] 0.416305756 -0.629651110
[339,] 0.032262621 0.416305756
[340,] -1.621780513 0.032262621
[341,] -3.245823647 -1.621780513
[342,] -3.989866782 -3.245823647
[343,] -4.263909916 -3.989866782
[344,] -2.667953050 -4.263909916
[345,] -1.181996184 -2.667953050
[346,] -4.241163326 -1.181996184
[347,] -3.625206460 -4.241163326
[348,] -2.349249594 -3.625206460
[349,] -0.523292729 -2.349249594
[350,] 0.462664137 -0.523292729
[351,] 0.338621003 0.462664137
[352,] -1.075422131 0.338621003
[353,] -2.959465266 -1.075422131
[354,] -3.943508400 -2.959465266
[355,] -3.567551534 -3.943508400
[356,] -1.981594668 -3.567551534
[357,] -0.485637803 -1.981594668
[358,] -3.654804944 -0.485637803
[359,] -2.038848078 -3.654804944
[360,] -1.132891213 -2.038848078
[361,] 0.443065653 -1.132891213
[362,] 1.319022519 0.443065653
[363,] 1.084979385 1.319022519
[364,] 0.240936250 1.084979385
[365,] -1.603106884 0.240936250
[366,] -2.917150018 -1.603106884
[367,] -2.551193153 -2.917150018
[368,] -1.235236287 -2.551193153
[369,] 0.330720579 -1.235236287
[370,] -2.608446563 0.330720579
[371,] -1.962489697 -2.608446563
[372,] -1.156532831 -1.962489697
[373,] 0.899424035 -1.156532831
[374,] 1.435380900 0.899424035
[375,] 1.111337766 1.435380900
[376,] 0.057294632 1.111337766
[377,] -2.046748502 0.057294632
[378,] -3.370791637 -2.046748502
[379,] -2.834834771 -3.370791637
[380,] -1.408877905 -2.834834771
[381,] 0.257078960 -1.408877905
[382,] -3.142088181 0.257078960
[383,] -1.786131315 -3.142088181
[384,] -0.750174449 -1.786131315
[385,] 0.125782416 -0.750174449
[386,] 1.371739282 0.125782416
[387,] 0.667696148 1.371739282
[388,] -0.506346987 0.667696148
[389,] -2.240390121 -0.506346987
[390,] -3.594433255 -2.240390121
[391,] -3.018476389 -3.594433255
[392,] -1.302519524 -3.018476389
[393,] 0.363437342 -1.302519524
[394,] -3.345729799 0.363437342
[395,] -2.279772934 -3.345729799
[396,] -0.683816068 -2.279772934
[397,] 0.982140798 -0.683816068
[398,] 1.848097664 0.982140798
[399,] 1.104054529 1.848097664
[400,] -0.659988605 1.104054529
[401,] -2.574031739 -0.659988605
[402,] -3.888074874 -2.574031739
[403,] -3.602118008 -3.888074874
[404,] -1.856161142 -3.602118008
[405,] -0.310204276 -1.856161142
[406,] -3.579371418 -0.310204276
[407,] -2.553414552 -3.579371418
[408,] -1.397457686 -2.553414552
[409,] 0.048499179 -1.397457686
[410,] 0.854456045 0.048499179
[411,] 0.860412911 0.854456045
[412,] -1.343630223 0.860412911
[413,] -3.027673358 -1.343630223
[414,] -4.731716492 -3.027673358
[415,] -4.055759626 -4.731716492
[416,] -2.849802760 -4.055759626
[417,] -1.343845895 -2.849802760
[418,] -4.343013036 -1.343845895
[419,] -3.667056170 -4.343013036
[420,] -2.311099305 -3.667056170
[421,] -1.135142439 -2.311099305
[422,] 0.030814427 -1.135142439
[423,] -0.363228707 0.030814427
[424,] -2.147271842 -0.363228707
[425,] -3.951314976 -2.147271842
[426,] -5.015358110 -3.951314976
[427,] -4.679401245 -5.015358110
[428,] -3.173444379 -4.679401245
[429,] -1.477487513 -3.173444379
[430,] -4.496654655 -1.477487513
[431,] -3.480697789 -4.496654655
[432,] -2.344740923 -3.480697789
[433,] -0.618784057 -2.344740923
[434,] -0.002827192 -0.618784057
[435,] -0.616870326 -0.002827192
[436,] -1.750913460 -0.616870326
[437,] -3.454956594 -1.750913460
[438,] -4.998999729 -3.454956594
[439,] -4.263042863 -4.998999729
[440,] -2.477085997 -4.263042863
[441,] -0.881129132 -2.477085997
[442,] -4.260296273 -0.881129132
[443,] -3.114339407 -4.260296273
[444,] -1.858382541 -3.114339407
[445,] -0.142425676 -1.858382541
[446,] 0.693531190 -0.142425676
[447,] 0.369488056 0.693531190
[448,] -0.844555079 0.369488056
[449,] -3.188598213 -0.844555079
[450,] -3.932641347 -3.188598213
[451,] -3.796684481 -3.932641347
[452,] -2.110727616 -3.796684481
[453,] -0.644770750 -2.110727616
[454,] -3.603937891 -0.644770750
[455,] -2.227981026 -3.603937891
[456,] -0.992024160 -2.227981026
[457,] -0.356067294 -0.992024160
[458,] 0.559889572 -0.356067294
[459,] 0.595846437 0.559889572
[460,] -0.528196697 0.595846437
[461,] -2.422239831 -0.528196697
[462,] -3.896282966 -2.422239831
[463,] -3.700326100 -3.896282966
[464,] -2.154369234 -3.700326100
[465,] -0.548412368 -2.154369234
[466,] -4.067579510 -0.548412368
[467,] -2.741622644 -4.067579510
[468,] -2.095665778 -2.741622644
[469,] -0.059708913 -2.095665778
[470,] 0.656247953 -0.059708913
[471,] -0.177795181 0.656247953
[472,] -1.181838315 -0.177795181
[473,] -3.125881450 -1.181838315
[474,] -4.699924584 -3.125881450
[475,] -3.993967718 -4.699924584
[476,] -2.008010852 -3.993967718
[477,] 0.167946013 -2.008010852
[478,] -3.321221128 0.167946013
[479,] -2.235264262 -3.321221128
[480,] -0.839307397 -2.235264262
[481,] 0.896649469 -0.839307397
[482,] 2.072606335 0.896649469
[483,] 1.788563200 2.072606335
[484,] 0.864520066 1.788563200
[485,] -0.649523068 0.864520066
[486,] -2.363566202 -0.649523068
[487,] -1.747609337 -2.363566202
[488,] -0.301652471 -1.747609337
[489,] 1.484304395 -0.301652471
[490,] -1.794862747 1.484304395
[491,] -0.668905881 -1.794862747
[492,] 0.227050985 -0.668905881
[493,] 1.803007851 0.227050985
[494,] 1.958964716 1.803007851
[495,] 1.664921582 1.958964716
[496,] 0.930878448 1.664921582
[497,] -1.203164686 0.930878448
[498,] -3.027207821 -1.203164686
[499,] -2.201250955 -3.027207821
[500,] -0.595294089 -2.201250955
[501,] 1.020662776 -0.595294089
[502,] -2.138504365 1.020662776
[503,] -1.622547499 -2.138504365
[504,] -0.296590633 -1.622547499
[505,] 1.219366232 -0.296590633
[506,] 1.175323098 1.219366232
[507,] 1.631279964 1.175323098
[508,] 0.077236829 1.631279964
[509,] -1.166806305 0.077236829
[510,] -2.210849439 -1.166806305
[511,] -1.874892573 -2.210849439
[512,] -0.278935708 -1.874892573
[513,] 1.337021158 -0.278935708
[514,] -2.182145983 1.337021158
[515,] -0.946189118 -2.182145983
[516,] 0.269767748 -0.946189118
[517,] 1.465724614 0.269767748
[518,] 2.181681480 1.465724614
[519,] 1.807638345 2.181681480
[520,] 0.453595211 1.807638345
[521,] -1.200447923 0.453595211
[522,] -2.404491058 -1.200447923
[523,] -1.878534192 -2.404491058
[524,] -0.352577326 -1.878534192
[525,] 1.383379540 -0.352577326
[526,] -1.715787602 1.383379540
[527,] -0.749830736 -1.715787602
[528,] 0.306126130 -0.749830736
[529,] 1.632082995 0.306126130
[530,] 2.548039861 1.632082995
[531,] 2.673996727 2.548039861
[532,] 1.409953593 2.673996727
[533,] -0.464089542 1.409953593
[534,] -1.248132676 -0.464089542
[535,] -1.202175810 -1.248132676
[536,] 0.753781056 -1.202175810
[537,] 2.539737921 0.753781056
[538,] -0.619429220 2.539737921
[539,] 0.356527646 -0.619429220
[540,] 1.502484511 0.356527646
[541,] 3.008441377 1.502484511
[542,] 4.064398243 3.008441377
[543,] 4.000355108 4.064398243
[544,] 2.706311974 4.000355108
[545,] 0.602268840 2.706311974
[546,] -0.291774294 0.602268840
[547,] -0.045817429 -0.291774294
[548,] 1.790139437 -0.045817429
[549,] 3.356096303 1.790139437
[550,] 0.046929161 3.356096303
[551,] 1.152886027 0.046929161
[552,] 2.428842893 1.152886027
[553,] 4.234799759 2.428842893
[554,] 4.700756624 4.234799759
[555,] 4.006713490 4.700756624
[556,] 2.092670356 4.006713490
[557,] 0.948627222 2.092670356
[558,] -0.755415913 0.948627222
[559,] -0.119459047 -0.755415913
[560,] 1.826497819 -0.119459047
[561,] 3.452454684 1.826497819
[562,] 0.003287543 3.452454684
[563,] 1.539244409 0.003287543
[564,] 3.015201275 1.539244409
[565,] 4.571158140 3.015201275
[566,] 5.087115006 4.571158140
[567,] 5.043071872 5.087115006
[568,] 3.769028737 5.043071872
[569,] 2.034985603 3.769028737
[570,] 0.420942469 2.034985603
[571,] 0.946899335 0.420942469
[572,] 2.522856200 0.946899335
[573,] 4.568813066 2.522856200
[574,] 1.489645925 4.568813066
[575,] 2.585602790 1.489645925
[576,] 3.311559656 2.585602790
[577,] 5.577516522 3.311559656
[578,] 6.133473388 5.577516522
[579,] 5.469430253 6.133473388
[580,] 3.875387119 5.469430253
[581,] 2.321343985 3.875387119
[582,] 1.067300850 2.321343985
[583,] 1.583257716 1.067300850
[584,] 2.959214582 1.583257716
[585,] 4.905171448 2.959214582
[586,] 1.566004306 4.905171448
[587,] 2.721961172 1.566004306
[588,] 3.637918038 2.721961172
[589,] 5.843874903 3.637918038
[590,] 6.279831769 5.843874903
[591,] 6.025788635 6.279831769
[592,] 4.731745501 6.025788635
[593,] 2.507702366 4.731745501
[594,] 1.623659232 2.507702366
[595,] 2.169616098 1.623659232
[596,] 3.715572964 2.169616098
[597,] 5.451529829 3.715572964
[598,] 2.662362688 5.451529829
[599,] 3.228319554 2.662362688
[600,] 3.734276419 3.228319554
[601,] 5.220233285 3.734276419
[602,] 6.806190151 5.220233285
[603,] 6.452147016 6.806190151
[604,] 5.218103882 6.452147016
[605,] 3.254060748 5.218103882
[606,] 2.440017614 3.254060748
[607,] 2.615974479 2.440017614
[608,] 4.011931345 2.615974479
[609,] 5.707888211 4.011931345
[610,] 2.718721069 5.707888211
[611,] 3.474677935 2.718721069
[612,] 5.100634801 3.474677935
[613,] 6.056591667 5.100634801
[614,] 7.042548532 6.056591667
[615,] 6.558505398 7.042548532
[616,] 5.134462264 6.558505398
[617,] 3.570419130 5.134462264
[618,] 2.696375995 3.570419130
[619,] 2.572332861 2.696375995
[620,] 4.448289727 2.572332861
[621,] 5.984246592 4.448289727
[622,] 2.805079451 5.984246592
[623,] 4.321036317 2.805079451
[624,] 5.966993183 4.321036317
[625,] 7.612950048 5.966993183
[626,] 8.378906914 7.612950048
[627,] 7.714863780 8.378906914
[628,] 6.080820645 7.714863780
[629,] 4.366777511 6.080820645
[630,] 3.282734377 4.366777511
[631,] 3.928691243 3.282734377
[632,] 5.604648108 3.928691243
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 7.672314205 5.666357339
2 7.988271071 7.672314205
3 7.874227937 7.988271071
4 6.880184802 7.874227937
5 6.216141668 6.880184802
6 4.752098534 6.216141668
7 4.418055399 4.752098534
8 5.414012265 4.418055399
9 7.019969131 5.414012265
10 3.600801989 7.019969131
11 4.626758855 3.600801989
12 5.222715721 4.626758855
13 6.498672587 5.222715721
14 7.334629452 6.498672587
15 7.470586318 7.334629452
16 6.126543184 7.470586318
17 4.642500050 6.126543184
18 3.948456915 4.642500050
19 3.634413781 3.948456915
20 5.440370647 3.634413781
21 6.496327512 5.440370647
22 2.967160371 6.496327512
23 3.773117237 2.967160371
24 4.649074103 3.773117237
25 6.355030968 4.649074103
26 7.620987834 6.355030968
27 7.456944700 7.620987834
28 6.312901565 7.456944700
29 4.308858431 6.312901565
30 2.824815297 4.308858431
31 2.760772163 2.824815297
32 4.196729028 2.760772163
33 5.652685894 4.196729028
34 2.023518753 5.652685894
35 3.059475618 2.023518753
36 4.165432484 3.059475618
37 5.371389350 4.165432484
38 6.737346216 5.371389350
39 6.193303081 6.737346216
40 5.269259947 6.193303081
41 3.745216813 5.269259947
42 2.021173679 3.745216813
43 2.867130544 2.021173679
44 3.853087410 2.867130544
45 5.029044276 3.853087410
46 1.589877134 5.029044276
47 2.465834000 1.589877134
48 3.861790866 2.465834000
49 5.077747731 3.861790866
50 5.723704597 5.077747731
51 5.529661463 5.723704597
52 4.825618329 5.529661463
53 2.911575194 4.825618329
54 2.037532060 2.911575194
55 1.463488926 2.037532060
56 2.999445792 1.463488926
57 4.275402657 2.999445792
58 0.946235516 4.275402657
59 1.552192382 0.946235516
60 2.598149247 1.552192382
61 4.394106113 2.598149247
62 5.510062979 4.394106113
63 5.006019845 5.510062979
64 3.541976710 5.006019845
65 1.837933576 3.541976710
66 0.543890442 1.837933576
67 0.589847307 0.543890442
68 1.985804173 0.589847307
69 3.441761039 1.985804173
70 0.332593897 3.441761039
71 1.108550763 0.332593897
72 2.044507629 1.108550763
73 3.360464495 2.044507629
74 4.066421360 3.360464495
75 3.982378226 4.066421360
76 2.798335092 3.982378226
77 1.324291958 2.798335092
78 -0.409751177 1.324291958
79 -0.053794311 -0.409751177
80 1.212162555 -0.053794311
81 2.398119420 1.212162555
82 -1.241047721 2.398119420
83 0.024909145 -1.241047721
84 0.740866011 0.024909145
85 2.246822876 0.740866011
86 2.542779742 2.246822876
87 2.518736608 2.542779742
88 2.304693473 2.518736608
89 -0.019349661 2.304693473
90 -0.743392795 -0.019349661
91 -0.987435929 -0.743392795
92 0.848520936 -0.987435929
93 1.664477802 0.848520936
94 -1.504689339 1.664477802
95 -0.268732474 -1.504689339
96 0.797224392 -0.268732474
97 2.543181258 0.797224392
98 2.949138124 2.543181258
99 2.955094989 2.949138124
100 1.871051855 2.955094989
101 0.107008721 1.871051855
102 -1.357034413 0.107008721
103 -1.631077548 -1.357034413
104 0.314879318 -1.631077548
105 1.880836184 0.314879318
106 -1.508330958 1.880836184
107 -0.802374092 -1.508330958
108 0.003582774 -0.802374092
109 1.649539639 0.003582774
110 2.495496505 1.649539639
111 1.851453371 2.495496505
112 0.577410237 1.851453371
113 -0.786632898 0.577410237
114 -2.130676032 -0.786632898
115 -1.864719166 -2.130676032
116 -0.188762300 -1.864719166
117 1.317194565 -0.188762300
118 -2.441972576 1.317194565
119 -1.596015710 -2.441972576
120 -0.590058845 -1.596015710
121 0.805898021 -0.590058845
122 1.621854887 0.805898021
123 1.677811753 1.621854887
124 0.723768618 1.677811753
125 -1.120274516 0.723768618
126 -2.474317650 -1.120274516
127 -2.368360785 -2.474317650
128 -1.042403919 -2.368360785
129 0.753552947 -1.042403919
130 -2.455614195 0.753552947
131 -1.769657329 -2.455614195
132 -0.283700463 -1.769657329
133 1.002256403 -0.283700463
134 1.948213268 1.002256403
135 1.634170134 1.948213268
136 1.020127000 1.634170134
137 -0.923916134 1.020127000
138 -1.947959269 -0.923916134
139 -2.282002403 -1.947959269
140 -0.946045537 -2.282002403
141 0.589911328 -0.946045537
142 -2.869255813 0.589911328
143 -1.643298947 -2.869255813
144 -0.497342081 -1.643298947
145 1.038614784 -0.497342081
146 1.234571650 1.038614784
147 1.090528516 1.234571650
148 0.046485381 1.090528516
149 -1.347557753 0.046485381
150 -2.671600887 -1.347557753
151 -2.345644021 -2.671600887
152 -1.259687156 -2.345644021
153 0.156269710 -1.259687156
154 -3.172897431 0.156269710
155 -2.396940566 -3.172897431
156 -1.630983700 -2.396940566
157 -0.765026834 -1.630983700
158 0.640930032 -0.765026834
159 0.556886897 0.640930032
160 -0.407156237 0.556886897
161 -2.021199371 -0.407156237
162 -3.855242505 -2.021199371
163 -3.389285640 -3.855242505
164 -1.883328774 -3.389285640
165 -0.407371908 -1.883328774
166 -4.016539050 -0.407371908
167 -2.890582184 -4.016539050
168 -2.504625318 -2.890582184
169 -0.268668453 -2.504625318
170 0.347288413 -0.268668453
171 -0.366754721 0.347288413
172 -1.140797855 -0.366754721
173 -2.604840990 -1.140797855
174 -3.728884124 -2.604840990
175 -3.332927258 -3.728884124
176 -1.626970392 -3.332927258
177 -0.311013527 -1.626970392
178 -3.680180668 -0.311013527
179 -2.404223802 -3.680180668
180 -1.398266937 -2.404223802
181 0.067689929 -1.398266937
182 1.313646795 0.067689929
183 1.169603661 1.313646795
184 0.235560526 1.169603661
185 -1.058482608 0.235560526
186 -2.592525742 -1.058482608
187 -2.656568877 -2.592525742
188 -1.410612011 -2.656568877
189 -0.664655145 -1.410612011
190 -4.323822287 -0.664655145
191 -2.697865421 -4.323822287
192 -1.661908555 -2.697865421
193 -0.225951689 -1.661908555
194 0.540005176 -0.225951689
195 -0.214037958 0.540005176
196 -1.088081092 -0.214037958
197 -2.642124226 -1.088081092
198 -4.136167361 -2.642124226
199 -4.070210495 -4.136167361
200 -2.674253629 -4.070210495
201 -1.248296764 -2.674253629
202 -4.437463905 -1.248296764
203 -3.441507039 -4.437463905
204 -2.735550173 -3.441507039
205 -1.029593308 -2.735550173
206 -0.143636442 -1.029593308
207 -0.387679576 -0.143636442
208 -1.781722711 -0.387679576
209 -3.375765845 -1.781722711
210 -4.419808979 -3.375765845
211 -4.373852113 -4.419808979
212 -3.127895248 -4.373852113
213 -1.611938382 -3.127895248
214 -4.901105523 -1.611938382
215 -3.545148658 -4.901105523
216 -2.569191792 -3.545148658
217 -0.983234926 -2.569191792
218 -0.707278060 -0.983234926
219 -1.181321195 -0.707278060
220 -2.015364329 -1.181321195
221 -4.059407463 -2.015364329
222 -5.473450598 -4.059407463
223 -5.397493732 -5.473450598
224 -3.721536866 -5.397493732
225 -1.985580000 -3.721536866
226 -5.314747142 -1.985580000
227 -4.508790276 -5.314747142
228 -2.502833410 -4.508790276
229 -1.176876545 -2.502833410
230 -0.120919679 -1.176876545
231 -0.574962813 -0.120919679
232 -1.619005947 -0.574962813
233 -3.613049082 -1.619005947
234 -4.687092216 -3.613049082
235 -4.401135350 -4.687092216
236 -2.995178484 -4.401135350
237 -1.609221619 -2.995178484
238 -4.358388760 -1.609221619
239 -3.922431894 -4.358388760
240 -2.296475029 -3.922431894
241 -0.880518163 -2.296475029
242 -0.384561297 -0.880518163
243 -0.138604431 -0.384561297
244 -1.372647566 -0.138604431
245 -3.216690700 -1.372647566
246 -4.990733834 -3.216690700
247 -4.654776969 -4.990733834
248 -3.028820103 -4.654776969
249 -1.692863237 -3.028820103
250 -4.682030379 -1.692863237
251 -3.986073513 -4.682030379
252 -2.230116647 -3.986073513
253 -1.134159781 -2.230116647
254 -0.808202916 -1.134159781
255 -0.362246050 -0.808202916
256 -1.696289184 -0.362246050
257 -3.470332318 -1.696289184
258 -4.834375453 -3.470332318
259 -4.378418587 -4.834375453
260 -2.972461721 -4.378418587
261 -1.196504856 -2.972461721
262 -4.535671997 -1.196504856
263 -3.789715131 -4.535671997
264 -1.763758265 -3.789715131
265 -0.677801400 -1.763758265
266 0.198155466 -0.677801400
267 0.254112332 0.198155466
268 -1.399930803 0.254112332
269 -3.023973937 -1.399930803
270 -4.108017071 -3.023973937
271 -3.892060205 -4.108017071
272 -2.466103340 -3.892060205
273 -1.120146474 -2.466103340
274 -4.419313615 -1.120146474
275 -3.003356750 -4.419313615
276 -1.677399884 -3.003356750
277 -0.421443018 -1.677399884
278 0.294513848 -0.421443018
279 0.070470713 0.294513848
280 -1.473572421 0.070470713
281 -3.407615555 -1.473572421
282 -4.691658690 -3.407615555
283 -4.255701824 -4.691658690
284 -2.539744958 -4.255701824
285 -0.943788092 -2.539744958
286 -4.302955234 -0.943788092
287 -3.266998368 -4.302955234
288 -1.821041502 -3.266998368
289 -0.595084637 -1.821041502
290 0.140872229 -0.595084637
291 -0.473170905 0.140872229
292 -1.777214039 -0.473170905
293 -3.541257174 -1.777214039
294 -5.115300308 -3.541257174
295 -4.739343442 -5.115300308
296 -3.323386576 -4.739343442
297 -1.627429711 -3.323386576
298 -5.246596852 -1.627429711
299 -3.730639986 -5.246596852
300 -2.684683121 -3.730639986
301 -0.788726255 -2.684683121
302 0.157230611 -0.788726255
303 -0.016812523 0.157230611
304 -1.140855658 -0.016812523
305 -2.724898792 -1.140855658
306 -4.538941926 -2.724898792
307 -4.152985061 -4.538941926
308 -2.657028195 -4.152985061
309 -0.851071329 -2.657028195
310 -4.240238471 -0.851071329
311 -3.254281605 -4.240238471
312 -2.288324739 -3.254281605
313 -0.722367873 -2.288324739
314 0.313588992 -0.722367873
315 0.019545858 0.313588992
316 -1.144497276 0.019545858
317 -3.188540410 -1.144497276
318 -4.742583545 -3.188540410
319 -4.176626679 -4.742583545
320 -2.550669813 -4.176626679
321 -0.824712948 -2.550669813
322 -4.473880089 -0.824712948
323 -3.397923223 -4.473880089
324 -1.531966357 -3.397923223
325 -0.736009492 -1.531966357
326 0.259947374 -0.736009492
327 -0.034095760 0.259947374
328 -1.468138895 -0.034095760
329 -2.992182029 -1.468138895
330 -4.386225163 -2.992182029
331 -4.240268297 -4.386225163
332 -2.674311432 -4.240268297
333 -1.018354566 -2.674311432
334 -4.577521707 -1.018354566
335 -3.781564842 -4.577521707
336 -2.585607976 -3.781564842
337 -0.629651110 -2.585607976
338 0.416305756 -0.629651110
339 0.032262621 0.416305756
340 -1.621780513 0.032262621
341 -3.245823647 -1.621780513
342 -3.989866782 -3.245823647
343 -4.263909916 -3.989866782
344 -2.667953050 -4.263909916
345 -1.181996184 -2.667953050
346 -4.241163326 -1.181996184
347 -3.625206460 -4.241163326
348 -2.349249594 -3.625206460
349 -0.523292729 -2.349249594
350 0.462664137 -0.523292729
351 0.338621003 0.462664137
352 -1.075422131 0.338621003
353 -2.959465266 -1.075422131
354 -3.943508400 -2.959465266
355 -3.567551534 -3.943508400
356 -1.981594668 -3.567551534
357 -0.485637803 -1.981594668
358 -3.654804944 -0.485637803
359 -2.038848078 -3.654804944
360 -1.132891213 -2.038848078
361 0.443065653 -1.132891213
362 1.319022519 0.443065653
363 1.084979385 1.319022519
364 0.240936250 1.084979385
365 -1.603106884 0.240936250
366 -2.917150018 -1.603106884
367 -2.551193153 -2.917150018
368 -1.235236287 -2.551193153
369 0.330720579 -1.235236287
370 -2.608446563 0.330720579
371 -1.962489697 -2.608446563
372 -1.156532831 -1.962489697
373 0.899424035 -1.156532831
374 1.435380900 0.899424035
375 1.111337766 1.435380900
376 0.057294632 1.111337766
377 -2.046748502 0.057294632
378 -3.370791637 -2.046748502
379 -2.834834771 -3.370791637
380 -1.408877905 -2.834834771
381 0.257078960 -1.408877905
382 -3.142088181 0.257078960
383 -1.786131315 -3.142088181
384 -0.750174449 -1.786131315
385 0.125782416 -0.750174449
386 1.371739282 0.125782416
387 0.667696148 1.371739282
388 -0.506346987 0.667696148
389 -2.240390121 -0.506346987
390 -3.594433255 -2.240390121
391 -3.018476389 -3.594433255
392 -1.302519524 -3.018476389
393 0.363437342 -1.302519524
394 -3.345729799 0.363437342
395 -2.279772934 -3.345729799
396 -0.683816068 -2.279772934
397 0.982140798 -0.683816068
398 1.848097664 0.982140798
399 1.104054529 1.848097664
400 -0.659988605 1.104054529
401 -2.574031739 -0.659988605
402 -3.888074874 -2.574031739
403 -3.602118008 -3.888074874
404 -1.856161142 -3.602118008
405 -0.310204276 -1.856161142
406 -3.579371418 -0.310204276
407 -2.553414552 -3.579371418
408 -1.397457686 -2.553414552
409 0.048499179 -1.397457686
410 0.854456045 0.048499179
411 0.860412911 0.854456045
412 -1.343630223 0.860412911
413 -3.027673358 -1.343630223
414 -4.731716492 -3.027673358
415 -4.055759626 -4.731716492
416 -2.849802760 -4.055759626
417 -1.343845895 -2.849802760
418 -4.343013036 -1.343845895
419 -3.667056170 -4.343013036
420 -2.311099305 -3.667056170
421 -1.135142439 -2.311099305
422 0.030814427 -1.135142439
423 -0.363228707 0.030814427
424 -2.147271842 -0.363228707
425 -3.951314976 -2.147271842
426 -5.015358110 -3.951314976
427 -4.679401245 -5.015358110
428 -3.173444379 -4.679401245
429 -1.477487513 -3.173444379
430 -4.496654655 -1.477487513
431 -3.480697789 -4.496654655
432 -2.344740923 -3.480697789
433 -0.618784057 -2.344740923
434 -0.002827192 -0.618784057
435 -0.616870326 -0.002827192
436 -1.750913460 -0.616870326
437 -3.454956594 -1.750913460
438 -4.998999729 -3.454956594
439 -4.263042863 -4.998999729
440 -2.477085997 -4.263042863
441 -0.881129132 -2.477085997
442 -4.260296273 -0.881129132
443 -3.114339407 -4.260296273
444 -1.858382541 -3.114339407
445 -0.142425676 -1.858382541
446 0.693531190 -0.142425676
447 0.369488056 0.693531190
448 -0.844555079 0.369488056
449 -3.188598213 -0.844555079
450 -3.932641347 -3.188598213
451 -3.796684481 -3.932641347
452 -2.110727616 -3.796684481
453 -0.644770750 -2.110727616
454 -3.603937891 -0.644770750
455 -2.227981026 -3.603937891
456 -0.992024160 -2.227981026
457 -0.356067294 -0.992024160
458 0.559889572 -0.356067294
459 0.595846437 0.559889572
460 -0.528196697 0.595846437
461 -2.422239831 -0.528196697
462 -3.896282966 -2.422239831
463 -3.700326100 -3.896282966
464 -2.154369234 -3.700326100
465 -0.548412368 -2.154369234
466 -4.067579510 -0.548412368
467 -2.741622644 -4.067579510
468 -2.095665778 -2.741622644
469 -0.059708913 -2.095665778
470 0.656247953 -0.059708913
471 -0.177795181 0.656247953
472 -1.181838315 -0.177795181
473 -3.125881450 -1.181838315
474 -4.699924584 -3.125881450
475 -3.993967718 -4.699924584
476 -2.008010852 -3.993967718
477 0.167946013 -2.008010852
478 -3.321221128 0.167946013
479 -2.235264262 -3.321221128
480 -0.839307397 -2.235264262
481 0.896649469 -0.839307397
482 2.072606335 0.896649469
483 1.788563200 2.072606335
484 0.864520066 1.788563200
485 -0.649523068 0.864520066
486 -2.363566202 -0.649523068
487 -1.747609337 -2.363566202
488 -0.301652471 -1.747609337
489 1.484304395 -0.301652471
490 -1.794862747 1.484304395
491 -0.668905881 -1.794862747
492 0.227050985 -0.668905881
493 1.803007851 0.227050985
494 1.958964716 1.803007851
495 1.664921582 1.958964716
496 0.930878448 1.664921582
497 -1.203164686 0.930878448
498 -3.027207821 -1.203164686
499 -2.201250955 -3.027207821
500 -0.595294089 -2.201250955
501 1.020662776 -0.595294089
502 -2.138504365 1.020662776
503 -1.622547499 -2.138504365
504 -0.296590633 -1.622547499
505 1.219366232 -0.296590633
506 1.175323098 1.219366232
507 1.631279964 1.175323098
508 0.077236829 1.631279964
509 -1.166806305 0.077236829
510 -2.210849439 -1.166806305
511 -1.874892573 -2.210849439
512 -0.278935708 -1.874892573
513 1.337021158 -0.278935708
514 -2.182145983 1.337021158
515 -0.946189118 -2.182145983
516 0.269767748 -0.946189118
517 1.465724614 0.269767748
518 2.181681480 1.465724614
519 1.807638345 2.181681480
520 0.453595211 1.807638345
521 -1.200447923 0.453595211
522 -2.404491058 -1.200447923
523 -1.878534192 -2.404491058
524 -0.352577326 -1.878534192
525 1.383379540 -0.352577326
526 -1.715787602 1.383379540
527 -0.749830736 -1.715787602
528 0.306126130 -0.749830736
529 1.632082995 0.306126130
530 2.548039861 1.632082995
531 2.673996727 2.548039861
532 1.409953593 2.673996727
533 -0.464089542 1.409953593
534 -1.248132676 -0.464089542
535 -1.202175810 -1.248132676
536 0.753781056 -1.202175810
537 2.539737921 0.753781056
538 -0.619429220 2.539737921
539 0.356527646 -0.619429220
540 1.502484511 0.356527646
541 3.008441377 1.502484511
542 4.064398243 3.008441377
543 4.000355108 4.064398243
544 2.706311974 4.000355108
545 0.602268840 2.706311974
546 -0.291774294 0.602268840
547 -0.045817429 -0.291774294
548 1.790139437 -0.045817429
549 3.356096303 1.790139437
550 0.046929161 3.356096303
551 1.152886027 0.046929161
552 2.428842893 1.152886027
553 4.234799759 2.428842893
554 4.700756624 4.234799759
555 4.006713490 4.700756624
556 2.092670356 4.006713490
557 0.948627222 2.092670356
558 -0.755415913 0.948627222
559 -0.119459047 -0.755415913
560 1.826497819 -0.119459047
561 3.452454684 1.826497819
562 0.003287543 3.452454684
563 1.539244409 0.003287543
564 3.015201275 1.539244409
565 4.571158140 3.015201275
566 5.087115006 4.571158140
567 5.043071872 5.087115006
568 3.769028737 5.043071872
569 2.034985603 3.769028737
570 0.420942469 2.034985603
571 0.946899335 0.420942469
572 2.522856200 0.946899335
573 4.568813066 2.522856200
574 1.489645925 4.568813066
575 2.585602790 1.489645925
576 3.311559656 2.585602790
577 5.577516522 3.311559656
578 6.133473388 5.577516522
579 5.469430253 6.133473388
580 3.875387119 5.469430253
581 2.321343985 3.875387119
582 1.067300850 2.321343985
583 1.583257716 1.067300850
584 2.959214582 1.583257716
585 4.905171448 2.959214582
586 1.566004306 4.905171448
587 2.721961172 1.566004306
588 3.637918038 2.721961172
589 5.843874903 3.637918038
590 6.279831769 5.843874903
591 6.025788635 6.279831769
592 4.731745501 6.025788635
593 2.507702366 4.731745501
594 1.623659232 2.507702366
595 2.169616098 1.623659232
596 3.715572964 2.169616098
597 5.451529829 3.715572964
598 2.662362688 5.451529829
599 3.228319554 2.662362688
600 3.734276419 3.228319554
601 5.220233285 3.734276419
602 6.806190151 5.220233285
603 6.452147016 6.806190151
604 5.218103882 6.452147016
605 3.254060748 5.218103882
606 2.440017614 3.254060748
607 2.615974479 2.440017614
608 4.011931345 2.615974479
609 5.707888211 4.011931345
610 2.718721069 5.707888211
611 3.474677935 2.718721069
612 5.100634801 3.474677935
613 6.056591667 5.100634801
614 7.042548532 6.056591667
615 6.558505398 7.042548532
616 5.134462264 6.558505398
617 3.570419130 5.134462264
618 2.696375995 3.570419130
619 2.572332861 2.696375995
620 4.448289727 2.572332861
621 5.984246592 4.448289727
622 2.805079451 5.984246592
623 4.321036317 2.805079451
624 5.966993183 4.321036317
625 7.612950048 5.966993183
626 8.378906914 7.612950048
627 7.714863780 8.378906914
628 6.080820645 7.714863780
629 4.366777511 6.080820645
630 3.282734377 4.366777511
631 3.928691243 3.282734377
632 5.604648108 3.928691243
> 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/wessaorg/rcomp/tmp/7a4tb1321896144.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/wessaorg/rcomp/tmp/8qi8t1321896144.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/wessaorg/rcomp/tmp/9sbvm1321896144.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/wessaorg/rcomp/tmp/10gxzi1321896144.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/112iyn1321896144.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/wessaorg/rcomp/tmp/12mqx21321896144.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/wessaorg/rcomp/tmp/13tp7o1321896145.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/wessaorg/rcomp/tmp/1499vk1321896145.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/wessaorg/rcomp/tmp/15dow41321896145.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/wessaorg/rcomp/tmp/1674xz1321896145.tab")
+ }
>
> try(system("convert tmp/139221321896144.ps tmp/139221321896144.png",intern=TRUE))
character(0)
> try(system("convert tmp/2u51d1321896144.ps tmp/2u51d1321896144.png",intern=TRUE))
character(0)
> try(system("convert tmp/324jf1321896144.ps tmp/324jf1321896144.png",intern=TRUE))
character(0)
> try(system("convert tmp/4l3sa1321896144.ps tmp/4l3sa1321896144.png",intern=TRUE))
character(0)
> try(system("convert tmp/5g2r51321896144.ps tmp/5g2r51321896144.png",intern=TRUE))
character(0)
> try(system("convert tmp/6uk171321896144.ps tmp/6uk171321896144.png",intern=TRUE))
character(0)
> try(system("convert tmp/7a4tb1321896144.ps tmp/7a4tb1321896144.png",intern=TRUE))
character(0)
> try(system("convert tmp/8qi8t1321896144.ps tmp/8qi8t1321896144.png",intern=TRUE))
character(0)
> try(system("convert tmp/9sbvm1321896144.ps tmp/9sbvm1321896144.png",intern=TRUE))
character(0)
> try(system("convert tmp/10gxzi1321896144.ps tmp/10gxzi1321896144.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
19.831 0.654 20.561