R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. 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. > y <- c(122.36 + ,123.33 + ,123.04 + ,124.53 + ,125.13 + ,125.85 + ,126.50 + ,126.53 + ,127.07 + ,124.55 + ,124.90 + ,124.32 + ,122.84 + ,123.31 + ,123.31 + ,124.87 + ,124.64 + ,124.73 + ,124.90 + ,124.04 + ,123.28 + ,123.86 + ,122.29 + ,124.09 + ,124.54 + ,125.65 + ,125.70 + ,125.53 + ,125.61 + ,125.55 + ,125.41 + ,127.60 + ,124.68 + ,124.41 + ,126.43 + ,126.38 + ,125.78 + ,124.70 + ,125.07 + ,125.25 + ,126.58 + ,127.13 + ,125.82 + ,123.70 + ,124.39 + ,123.70 + ,124.42 + ,121.05 + ,121.02 + ,123.23 + ,121.32 + ,120.91 + ,120.72 + ,123.31 + ,119.58 + ,119.53 + ,120.59 + ,118.63 + ,118.47 + ,111.81 + ,114.71 + ,117.34 + ,115.77 + ,118.38 + ,117.84 + ,118.83 + ,120.02 + ,116.21 + ,117.08 + ,120.20 + ,119.83 + ,118.92 + ,118.03 + ,117.71 + ,119.55 + ,116.13 + ,115.97 + ,115.99 + ,114.96 + ,116.46 + ,116.55 + ,113.05 + ,117.44 + ,118.84 + ,117.06 + ,117.54 + ,119.31 + ,118.72 + ,121.55 + ,122.61 + ,121.53 + ,123.31 + ,124.07 + ,123.59 + ,122.97 + ,123.22 + ,123.04 + ,122.96 + ,122.81 + ,122.81 + ,122.62 + ,120.82 + ,119.41 + ,121.56 + ,121.59 + ,118.50 + ,118.77 + ,118.86 + ,117.60 + ,119.90 + ,121.83 + ,121.84 + ,122.12 + ,122.12 + ,121.36 + ,119.66 + ,119.32 + ,120.36 + ,117.06 + ,117.48 + ,115.60 + ,113.86 + ,116.92 + ,117.75 + ,117.75 + ,115.31 + ,116.28 + ,115.22 + ,115.65 + ,115.11 + ,118.67 + ,118.04 + ,116.50 + ,119.78 + ,119.95 + ,120.37 + ,119.79 + ,119.43 + ,121.06 + ,121.74 + ,121.09 + ,122.97 + ,120.50 + ,117.18 + ,115.03 + ,113.36 + ,112.59 + ,111.65 + ,111.98 + ,114.87 + ,114.67 + ,114.09 + ,114.77 + ,117.05 + ,117.22 + ,113.18 + ,110.95 + ,112.14 + ,112.72 + ,110.01 + ,110.29 + ,110.74 + ,110.32 + ,105.89 + ,108.97 + ,109.34 + ,106.57 + ,99.49 + ,101.81 + ,104.29 + ,109.73 + ,105.06 + ,107.97 + ,108.13 + ,109.86 + ,108.95 + ,111.20 + ,110.69 + ,106.10 + ,105.68 + ,104.12 + ,104.71 + ,104.30 + ,103.52 + ,107.76 + ,107.80 + ,107.30 + ,108.64 + ,105.03 + ,108.30 + ,107.21 + ,109.27 + ,109.50 + ,111.68 + ,111.80 + ,111.75 + ,106.68 + ,106.37 + ,105.76 + ,109.01 + ,109.01 + ,109.01 + ,109.01 + ,107.69 + ,105.19 + ,105.48 + ,102.22 + ,100.54 + ,105.00 + ,105.44 + ,107.89 + ,108.64 + ,106.70 + ,109.10 + ,105.23 + ,108.41 + ,108.80 + ,110.39 + ,110.22 + ,110.86 + ,108.58 + ,107.70 + ,106.62 + ,109.84 + ,107.16 + ,107.26 + ,108.70 + ,109.85 + ,109.41 + ,112.36 + ,111.03 + ,110.67 + ,109.21 + ,113.58 + ,113.88 + ,114.08 + ,112.33 + ,113.92 + ,114.41 + ,114.57 + ,115.35 + ,113.13 + ,113.29 + ,112.56 + ,113.06 + ,113.46 + ,115.39 + ,116.62 + ,117.04 + ,117.42 + ,115.62 + ,115.16 + ,115.69 + ,112.85 + ,114.05 + ,112.00 + ,113.74 + ,116.26 + ,118.63 + ,116.49 + ,118.23 + ,116.83 + ,118.82 + ,114.36 + ,112.02 + ,113.24 + ,109.75 + ,110.33 + ,112.86 + ,113.04 + ,113.80 + ,110.90 + ,109.96 + ,108.69 + ,108.84 + ,108.47 + ,108.07 + ,107.94 + ,108.11 + ,108.11 + ,106.81 + ,105.58 + ,105.61 + ,106.52 + ,103.86 + ,104.60 + ,104.73 + ,105.12 + ,104.76 + ,103.85 + ,103.83 + ,103.22 + ,101.64 + ,102.13 + ,104.33 + ,104.92 + ,107.78 + ,104.49 + ,102.80 + ,102.86 + ,104.51 + ,104.73 + ,102.58 + ,99.93 + ,101.41 + ,101.05 + ,99.86 + ,101.11 + ,100.89 + ,101.09 + ,98.31 + ,98.08 + ,99.55 + ,99.62 + ,97.37 + ,98.16 + ,97.98 + ,98.15 + ,97.10 + ,97.24 + ,96.70 + ,96.64 + ,100.65 + ,96.75 + ,97.74 + ,97.92 + ,98.34 + ,93.84 + ,97.80 + ,96.20 + ,95.99 + ,95.18 + ,95.95 + ,92.23 + ,91.78 + ,92.97 + ,89.76 + ,92.88 + ,96.23 + ,95.79 + ,93.97 + ,93.90 + ,93.60 + ,93.96 + ,88.69 + ,88.57 + ,85.62 + ,86.25 + ,85.33 + ,83.33 + ,77.78 + ,78.70 + ,72.05 + ,80.75 + ,81.41 + ,82.65 + ,75.85 + ,75.70 + ,78.25 + ,77.41 + ,76.84 + ,74.25 + ,74.95 + ,68.78 + ,73.21 + ,73.26 + ,78.67 + ,75.63 + ,74.99 + ,83.87 + ,79.62 + ,80.13 + ,79.76 + ,78.20 + ,78.05 + ,79.05 + ,73.32 + ,75.17 + ,73.26 + ,73.72 + ,73.57 + ,70.60 + ,71.25 + ,74.22 + ,73.32 + ,73.01 + ,74.21 + ,75.32 + ,71.73 + ,71.94 + ,72.94 + ,72.47 + ,71.94 + ,74.30 + ,74.30) > x <- c(358.59 + ,362.96 + ,362.42 + ,364.97 + ,364.04 + ,361.06 + ,358.48 + ,352.96 + ,359.59 + ,360.39 + ,357.40 + ,362.93 + ,364.55 + ,365.73 + ,364.70 + ,364.65 + ,359.43 + ,362.14 + ,356.97 + ,354.82 + ,353.17 + ,357.06 + ,356.18 + ,355.01 + ,355.65 + ,357.31 + ,357.07 + ,357.91 + ,358.48 + ,358.97 + ,351.77 + ,352.16 + ,359.08 + ,360.35 + ,359.53 + ,359.30 + ,358.41 + ,359.68 + ,355.31 + ,357.08 + ,349.71 + ,354.13 + ,345.49 + ,341.69 + ,344.25 + ,340.17 + ,342.47 + ,344.43 + ,333.23 + ,339.72 + ,342.61 + ,346.36 + ,339.09 + ,339.73 + ,341.12 + ,335.94 + ,333.46 + ,335.66 + ,341.12 + ,342.21 + ,342.62 + ,346.06 + ,344.43 + ,346.65 + ,343.74 + ,335.67 + ,342.75 + ,341.77 + ,345.84 + ,346.52 + ,350.79 + ,345.44 + ,345.87 + ,338.48 + ,337.21 + ,340.81 + ,339.86 + ,342.86 + ,343.33 + ,341.73 + ,351.38 + ,351.13 + ,345.99 + ,347.55 + ,346.02 + ,345.29 + ,347.03 + ,348.01 + ,345.48 + ,349.40 + ,351.05 + ,349.70 + ,350.86 + ,354.45 + ,355.30 + ,357.48 + ,355.24 + ,351.79 + ,355.22 + ,351.02 + ,350.28 + ,350.17 + ,348.16 + ,340.30 + ,343.75 + ,344.71 + ,344.13 + ,342.14 + ,345.04 + ,346.02 + ,346.43 + ,347.07 + ,339.33 + ,339.10 + ,337.19 + ,339.58 + ,327.85 + ,326.81 + ,321.73 + ,320.45 + ,327.69 + ,323.95 + ,320.47 + ,322.13 + ,316.34 + ,314.78 + ,308.90 + ,308.62 + ,314.41 + ,306.88 + ,310.60 + ,321.60 + ,321.50 + ,325.68 + ,324.35 + ,320.01 + ,326.88 + ,332.39 + ,331.48 + ,332.62 + ,324.79 + ,327.12 + ,328.91 + ,328.37 + ,324.83 + ,325.90 + ,326.18 + ,328.94 + ,333.78 + ,328.06 + ,325.87 + ,325.41 + ,318.86 + ,319.13 + ,310.16 + ,311.73 + ,306.54 + ,311.16 + ,311.98 + ,306.72 + ,308.05 + ,300.76 + ,301.90 + ,293.09 + ,292.76 + ,294.58 + ,289.90 + ,296.69 + ,297.21 + ,293.31 + ,296.25 + ,298.60 + ,296.87 + ,301.02 + ,304.73 + ,301.92 + ,295.72 + ,293.18 + ,298.35 + ,297.99 + ,299.85 + ,299.85 + ,304.45 + ,299.45 + ,298.14 + ,298.78 + ,297.02 + ,301.33 + ,294.96 + ,296.69 + ,300.73 + ,301.96 + ,297.38 + ,293.87 + ,285.96 + ,285.41 + ,283.70 + ,284.76 + ,277.11 + ,274.73 + ,274.73 + ,274.73 + ,274.73 + ,274.69 + ,275.42 + ,264.15 + ,276.24 + ,268.88 + ,277.97 + ,280.49 + ,281.09 + ,276.16 + ,272.58 + ,270.94 + ,284.31 + ,283.94 + ,284.18 + ,282.83 + ,283.84 + ,282.71 + ,279.29 + ,280.70 + ,274.47 + ,273.44 + ,275.49 + ,279.46 + ,280.19 + ,288.21 + ,284.80 + ,281.41 + ,283.39 + ,287.97 + ,290.77 + ,290.60 + ,289.67 + ,289.84 + ,298.55 + ,296.07 + ,297.14 + ,295.34 + ,296.25 + ,294.30 + ,296.15 + ,296.49 + ,298.05 + ,301.03 + ,300.52 + ,301.50 + ,296.93 + ,289.84 + ,291.44 + ,286.88 + ,286.74 + ,288.93 + ,292.19 + ,295.39 + ,295.86 + ,293.36 + ,292.86 + ,292.73 + ,296.73 + ,285.02 + ,285.24 + ,288.62 + ,283.36 + ,285.84 + ,291.48 + ,291.41 + ,287.77 + ,284.97 + ,286.05 + ,278.19 + ,281.21 + ,277.92 + ,280.08 + ,269.24 + ,268.48 + ,268.83 + ,269.54 + ,262.37 + ,265.12 + ,265.34 + ,263.32 + ,267.18 + ,260.75 + ,261.78 + ,257.27 + ,255.63 + ,251.39 + ,259.49 + ,261.18 + ,261.65 + ,262.01 + ,265.23 + ,268.10 + ,262.27 + ,263.59 + ,257.85 + ,265.69 + ,271.15 + ,266.69 + ,265.77 + ,262.32 + ,270.48 + ,273.03 + ,269.13 + ,280.65 + ,282.75 + ,281.44 + ,281.99 + ,282.86 + ,287.21 + ,283.11 + ,280.66 + ,282.39 + ,280.83 + ,284.71 + ,279.99 + ,283.50 + ,284.88 + ,288.60 + ,284.80 + ,287.20 + ,286.22 + ,286.54 + ,279.58 + ,283.08 + ,288.88 + ,280.18 + ,284.16 + ,290.57 + ,286.82 + ,273.00 + ,278.69 + ,264.54 + ,271.92 + ,283.60 + ,269.25 + ,263.58 + ,264.16 + ,268.85 + ,269.67 + ,249.41 + ,268.99 + ,268.65 + ,260.16 + ,256.55 + ,251.47 + ,234.93 + ,232.96 + ,215.49 + ,213.68 + ,236.07 + ,235.41 + ,214.77 + ,225.85 + ,224.64 + ,238.26 + ,232.44 + ,222.50 + ,225.28 + ,220.49 + ,216.86 + ,234.70 + ,230.06 + ,238.27 + ,238.56 + ,242.70 + ,249.14 + ,234.89 + ,227.78 + ,234.04 + ,230.70 + ,230.17 + ,218.23 + ,232.20 + ,220.76 + ,215.60 + ,217.69 + ,204.35 + ,191.44 + ,203.84 + ,211.86 + ,210.57 + ,219.57 + ,219.98 + ,226.01 + ,207.04 + ,212.52 + ,217.92 + ,210.45 + ,218.53 + ,223.32 + ,218.76 + ,217.63) > par7 = '0' > par6 = '1' > par5 = '1' > par4 = '1' > par3 = '0' > par2 = '1' > par1 = '1' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: Wessa P., (2008), Cross Correlation Function (v1.0.6) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_cross.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: > par1 <- as.numeric(par1) > par2 <- as.numeric(par2) > par3 <- as.numeric(par3) > par4 <- as.numeric(par4) > par5 <- as.numeric(par5) > par6 <- as.numeric(par6) > par7 <- as.numeric(par7) > if (par1 == 0) { + x <- log(x) + } else { + x <- (x ^ par1 - 1) / par1 + } > if (par5 == 0) { + y <- log(y) + } else { + y <- (y ^ par5 - 1) / par5 + } > if (par2 > 0) x <- diff(x,lag=1,difference=par2) > if (par6 > 0) y <- diff(y,lag=1,difference=par6) > if (par3 > 0) x <- diff(x,lag=par4,difference=par3) > if (par7 > 0) y <- diff(y,lag=par4,difference=par7) > x [1] 4.37 -0.54 2.55 -0.93 -2.98 -2.58 -5.52 6.63 0.80 -2.99 [11] 5.53 1.62 1.18 -1.03 -0.05 -5.22 2.71 -5.17 -2.15 -1.65 [21] 3.89 -0.88 -1.17 0.64 1.66 -0.24 0.84 0.57 0.49 -7.20 [31] 0.39 6.92 1.27 -0.82 -0.23 -0.89 1.27 -4.37 1.77 -7.37 [41] 4.42 -8.64 -3.80 2.56 -4.08 2.30 1.96 -11.20 6.49 2.89 [51] 3.75 -7.27 0.64 1.39 -5.18 -2.48 2.20 5.46 1.09 0.41 [61] 3.44 -1.63 2.22 -2.91 -8.07 7.08 -0.98 4.07 0.68 4.27 [71] -5.35 0.43 -7.39 -1.27 3.60 -0.95 3.00 0.47 -1.60 9.65 [81] -0.25 -5.14 1.56 -1.53 -0.73 1.74 0.98 -2.53 3.92 1.65 [91] -1.35 1.16 3.59 0.85 2.18 -2.24 -3.45 3.43 -4.20 -0.74 [101] -0.11 -2.01 -7.86 3.45 0.96 -0.58 -1.99 2.90 0.98 0.41 [111] 0.64 -7.74 -0.23 -1.91 2.39 -11.73 -1.04 -5.08 -1.28 7.24 [121] -3.74 -3.48 1.66 -5.79 -1.56 -5.88 -0.28 5.79 -7.53 3.72 [131] 11.00 -0.10 4.18 -1.33 -4.34 6.87 5.51 -0.91 1.14 -7.83 [141] 2.33 1.79 -0.54 -3.54 1.07 0.28 2.76 4.84 -5.72 -2.19 [151] -0.46 -6.55 0.27 -8.97 1.57 -5.19 4.62 0.82 -5.26 1.33 [161] -7.29 1.14 -8.81 -0.33 1.82 -4.68 6.79 0.52 -3.90 2.94 [171] 2.35 -1.73 4.15 3.71 -2.81 -6.20 -2.54 5.17 -0.36 1.86 [181] 0.00 4.60 -5.00 -1.31 0.64 -1.76 4.31 -6.37 1.73 4.04 [191] 1.23 -4.58 -3.51 -7.91 -0.55 -1.71 1.06 -7.65 -2.38 0.00 [201] 0.00 0.00 -0.04 0.73 -11.27 12.09 -7.36 9.09 2.52 0.60 [211] -4.93 -3.58 -1.64 13.37 -0.37 0.24 -1.35 1.01 -1.13 -3.42 [221] 1.41 -6.23 -1.03 2.05 3.97 0.73 8.02 -3.41 -3.39 1.98 [231] 4.58 2.80 -0.17 -0.93 0.17 8.71 -2.48 1.07 -1.80 0.91 [241] -1.95 1.85 0.34 1.56 2.98 -0.51 0.98 -4.57 -7.09 1.60 [251] -4.56 -0.14 2.19 3.26 3.20 0.47 -2.50 -0.50 -0.13 4.00 [261] -11.71 0.22 3.38 -5.26 2.48 5.64 -0.07 -3.64 -2.80 1.08 [271] -7.86 3.02 -3.29 2.16 -10.84 -0.76 0.35 0.71 -7.17 2.75 [281] 0.22 -2.02 3.86 -6.43 1.03 -4.51 -1.64 -4.24 8.10 1.69 [291] 0.47 0.36 3.22 2.87 -5.83 1.32 -5.74 7.84 5.46 -4.46 [301] -0.92 -3.45 8.16 2.55 -3.90 11.52 2.10 -1.31 0.55 0.87 [311] 4.35 -4.10 -2.45 1.73 -1.56 3.88 -4.72 3.51 1.38 3.72 [321] -3.80 2.40 -0.98 0.32 -6.96 3.50 5.80 -8.70 3.98 6.41 [331] -3.75 -13.82 5.69 -14.15 7.38 11.68 -14.35 -5.67 0.58 4.69 [341] 0.82 -20.26 19.58 -0.34 -8.49 -3.61 -5.08 -16.54 -1.97 -17.47 [351] -1.81 22.39 -0.66 -20.64 11.08 -1.21 13.62 -5.82 -9.94 2.78 [361] -4.79 -3.63 17.84 -4.64 8.21 0.29 4.14 6.44 -14.25 -7.11 [371] 6.26 -3.34 -0.53 -11.94 13.97 -11.44 -5.16 2.09 -13.34 -12.91 [381] 12.40 8.02 -1.29 9.00 0.41 6.03 -18.97 5.48 5.40 -7.47 [391] 8.08 4.79 -4.56 -1.13 > y [1] 0.97 -0.29 1.49 0.60 0.72 0.65 0.03 0.54 -2.52 0.35 -0.58 -1.48 [13] 0.47 0.00 1.56 -0.23 0.09 0.17 -0.86 -0.76 0.58 -1.57 1.80 0.45 [25] 1.11 0.05 -0.17 0.08 -0.06 -0.14 2.19 -2.92 -0.27 2.02 -0.05 -0.60 [37] -1.08 0.37 0.18 1.33 0.55 -1.31 -2.12 0.69 -0.69 0.72 -3.37 -0.03 [49] 2.21 -1.91 -0.41 -0.19 2.59 -3.73 -0.05 1.06 -1.96 -0.16 -6.66 2.90 [61] 2.63 -1.57 2.61 -0.54 0.99 1.19 -3.81 0.87 3.12 -0.37 -0.91 -0.89 [73] -0.32 1.84 -3.42 -0.16 0.02 -1.03 1.50 0.09 -3.50 4.39 1.40 -1.78 [85] 0.48 1.77 -0.59 2.83 1.06 -1.08 1.78 0.76 -0.48 -0.62 0.25 -0.18 [97] -0.08 -0.15 0.00 -0.19 -1.80 -1.41 2.15 0.03 -3.09 0.27 0.09 -1.26 [109] 2.30 1.93 0.01 0.28 0.00 -0.76 -1.70 -0.34 1.04 -3.30 0.42 -1.88 [121] -1.74 3.06 0.83 0.00 -2.44 0.97 -1.06 0.43 -0.54 3.56 -0.63 -1.54 [133] 3.28 0.17 0.42 -0.58 -0.36 1.63 0.68 -0.65 1.88 -2.47 -3.32 -2.15 [145] -1.67 -0.77 -0.94 0.33 2.89 -0.20 -0.58 0.68 2.28 0.17 -4.04 -2.23 [157] 1.19 0.58 -2.71 0.28 0.45 -0.42 -4.43 3.08 0.37 -2.77 -7.08 2.32 [169] 2.48 5.44 -4.67 2.91 0.16 1.73 -0.91 2.25 -0.51 -4.59 -0.42 -1.56 [181] 0.59 -0.41 -0.78 4.24 0.04 -0.50 1.34 -3.61 3.27 -1.09 2.06 0.23 [193] 2.18 0.12 -0.05 -5.07 -0.31 -0.61 3.25 0.00 0.00 0.00 -1.32 -2.50 [205] 0.29 -3.26 -1.68 4.46 0.44 2.45 0.75 -1.94 2.40 -3.87 3.18 0.39 [217] 1.59 -0.17 0.64 -2.28 -0.88 -1.08 3.22 -2.68 0.10 1.44 1.15 -0.44 [229] 2.95 -1.33 -0.36 -1.46 4.37 0.30 0.20 -1.75 1.59 0.49 0.16 0.78 [241] -2.22 0.16 -0.73 0.50 0.40 1.93 1.23 0.42 0.38 -1.80 -0.46 0.53 [253] -2.84 1.20 -2.05 1.74 2.52 2.37 -2.14 1.74 -1.40 1.99 -4.46 -2.34 [265] 1.22 -3.49 0.58 2.53 0.18 0.76 -2.90 -0.94 -1.27 0.15 -0.37 -0.40 [277] -0.13 0.17 0.00 -1.30 -1.23 0.03 0.91 -2.66 0.74 0.13 0.39 -0.36 [289] -0.91 -0.02 -0.61 -1.58 0.49 2.20 0.59 2.86 -3.29 -1.69 0.06 1.65 [301] 0.22 -2.15 -2.65 1.48 -0.36 -1.19 1.25 -0.22 0.20 -2.78 -0.23 1.47 [313] 0.07 -2.25 0.79 -0.18 0.17 -1.05 0.14 -0.54 -0.06 4.01 -3.90 0.99 [325] 0.18 0.42 -4.50 3.96 -1.60 -0.21 -0.81 0.77 -3.72 -0.45 1.19 -3.21 [337] 3.12 3.35 -0.44 -1.82 -0.07 -0.30 0.36 -5.27 -0.12 -2.95 0.63 -0.92 [349] -2.00 -5.55 0.92 -6.65 8.70 0.66 1.24 -6.80 -0.15 2.55 -0.84 -0.57 [361] -2.59 0.70 -6.17 4.43 0.05 5.41 -3.04 -0.64 8.88 -4.25 0.51 -0.37 [373] -1.56 -0.15 1.00 -5.73 1.85 -1.91 0.46 -0.15 -2.97 0.65 2.97 -0.90 [385] -0.31 1.20 1.11 -3.59 0.21 1.00 -0.47 -0.53 2.36 0.00 > postscript(file="/var/www/html/rcomp/tmp/1572t1260888763.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > (r <- ccf(x,y,main='Cross Correlation Function',ylab='CCF',xlab='Lag (k)')) Autocorrelations of series 'X', by lag -22 -21 -20 -19 -18 -17 -16 -15 -14 -13 -12 0.071 -0.105 -0.015 0.044 -0.027 0.058 -0.078 -0.001 0.050 -0.039 0.104 -11 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 -0.027 -0.037 0.015 0.000 -0.024 -0.002 -0.017 -0.011 -0.035 0.415 0.009 0 1 2 3 4 5 6 7 8 9 10 -0.126 0.020 0.042 -0.049 0.026 -0.039 0.096 -0.052 -0.032 0.098 -0.013 11 12 13 14 15 16 17 18 19 20 21 -0.044 -0.050 -0.042 0.026 0.066 -0.089 0.101 -0.025 -0.055 0.092 -0.070 22 0.028 > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Cross Correlation Function',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'Value',header=TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Box-Cox transformation parameter (lambda) of X series',header=TRUE) > a<-table.element(a,par1) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Degree of non-seasonal differencing (d) of X series',header=TRUE) > a<-table.element(a,par2) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Degree of seasonal differencing (D) of X series',header=TRUE) > a<-table.element(a,par3) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Seasonal Period (s)',header=TRUE) > a<-table.element(a,par4) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Box-Cox transformation parameter (lambda) of Y series',header=TRUE) > a<-table.element(a,par5) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Degree of non-seasonal differencing (d) of Y series',header=TRUE) > a<-table.element(a,par6) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Degree of seasonal differencing (D) of Y series',header=TRUE) > a<-table.element(a,par7) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'k',header=TRUE) > a<-table.element(a,'rho(Y[t],X[t+k])',header=TRUE) > a<-table.row.end(a) > mylength <- length(r$acf) > myhalf <- floor((mylength-1)/2) > for (i in 1:mylength) { + a<-table.row.start(a) + a<-table.element(a,i-myhalf-1,header=TRUE) + a<-table.element(a,r$acf[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/2t12g1260888763.tab") > > try(system("convert tmp/1572t1260888763.ps tmp/1572t1260888763.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.465 0.178 0.610