R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

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 <- 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)
> par5 = '1'
> par4 = '0'
> par3 = '1'
> par2 = '1'
> par1 = '48'
> #'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!)
> if (par1 == 'Default') {
+ par1 = 10*log10(length(x))
+ } else {
+ par1 <- as.numeric(par1)
+ }
> par2 <- as.numeric(par2)
> par3 <- as.numeric(par3)
> par4 <- as.numeric(par4)
> par5 <- as.numeric(par5)
> if (par2 == 0) {
+ x <- log(x)
+ } else {
+ x <- (x ^ par2 - 1) / par2
+ }
> if (par3 > 0) x <- diff(x,lag=1,difference=par3)
> if (par4 > 0) x <- diff(x,lag=par5,difference=par4)
> postscript(file="/var/www/html/rcomp/tmp/1nl3y1229550185.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> racf <- acf(x,par1,main='Autocorrelation',xlab='lags',ylab='ACF')
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/2wl431229550185.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> rpacf <- pacf(x,par1,main='Partial Autocorrelation',xlab='lags',ylab='PACF')
> dev.off()
null device 
          1 
> (myacf <- c(racf$acf))
 [1]  1.000000000 -0.226173549 -0.090782781  0.088800029 -0.040499387
 [6]  0.003961423 -0.027461766 -0.065853696  0.113433338 -0.049589011
[11]  0.066575109  0.048933266 -0.054633841  0.036698043 -0.012876184
[16] -0.141861510  0.040161365  0.062497322 -0.074077852 -0.016253277
[21]  0.025607406 -0.062514452  0.032942480 -0.020522620 -0.012387728
[26] -0.007923845  0.032805047 -0.019348800 -0.018655368  0.063073859
[31]  0.012189919 -0.091364161  0.039948448 -0.028824469 -0.023309429
[36] -0.058260691 -0.002443793  0.091022970 -0.035600632 -0.004452556
[41] -0.004288880 -0.016345360  0.046622616 -0.072827401 -0.091520354
[46]  0.063494502  0.055498822 -0.067121246  0.021728446
> (mypacf <- c(rpacf$acf))
 [1] -0.2261735493 -0.1495894242  0.0338101319 -0.0248178277  0.0036341619
 [6] -0.0390267450 -0.0841979821  0.0745179925 -0.0164878961  0.0856621746
[11]  0.0687552853 -0.0047902801  0.0229682775 -0.0077084790 -0.1358107190
[16] -0.0434306536  0.0569408562 -0.0436091010 -0.0507421187 -0.0111936626
[21] -0.0937718016 -0.0093778000 -0.0030588176 -0.0096803960 -0.0157642075
[26]  0.0541550042 -0.0192522835 -0.0235322722  0.0774941851  0.0198756167
[31] -0.0608282440  0.0265587241 -0.0559857772 -0.0571832902 -0.0850617966
[36] -0.0628155882  0.0481638312  0.0074776632 -0.0098623924 -0.0558114059
[41] -0.0004744742  0.0305449716 -0.0490064411 -0.0767512133 -0.0001800220
[46]  0.0514459125 -0.0346525160 -0.0043445311
> lengthx <- length(x)
> sqrtn <- sqrt(lengthx)
> 
> #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,'Autocorrelation Function',4,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Time lag k',header=TRUE)
> a<-table.element(a,hyperlink('http://www.xycoon.com/basics.htm','ACF(k)','click here for more information about the Autocorrelation Function'),header=TRUE)
> a<-table.element(a,'T-STAT',header=TRUE)
> a<-table.element(a,'P-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 2:(par1+1)) {
+ a<-table.row.start(a)
+ a<-table.element(a,i-1,header=TRUE)
+ a<-table.element(a,round(myacf[i],6))
+ mytstat <- myacf[i]*sqrtn
+ a<-table.element(a,round(mytstat,4))
+ a<-table.element(a,round(1-pt(abs(mytstat),lengthx),6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/3wh221229550185.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Partial Autocorrelation Function',4,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Time lag k',header=TRUE)
> a<-table.element(a,hyperlink('http://www.xycoon.com/basics.htm','PACF(k)','click here for more information about the Partial Autocorrelation Function'),header=TRUE)
> a<-table.element(a,'T-STAT',header=TRUE)
> a<-table.element(a,'P-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:par1) {
+ a<-table.row.start(a)
+ a<-table.element(a,i,header=TRUE)
+ a<-table.element(a,round(mypacf[i],6))
+ mytstat <- mypacf[i]*sqrtn
+ a<-table.element(a,round(mytstat,4))
+ a<-table.element(a,round(1-pt(abs(mytstat),lengthx),6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/431jm1229550185.tab") 
> 
> system("convert tmp/1nl3y1229550185.ps tmp/1nl3y1229550185.png")
> system("convert tmp/2wl431229550185.ps tmp/2wl431229550185.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.625   0.312   0.754