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(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) > par5 = '1' > par4 = '0' > par3 = '1' > par2 = '2.0' > par1 = '60' > #'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/freestat/rcomp/tmp/14am01229276943.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/freestat/rcomp/tmp/2hpgx1229276943.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.219110031 -0.031683878 -0.055897689 0.055851598 [6] -0.041618038 -0.051256977 -0.082678060 0.124802189 -0.096877893 [11] -0.011145675 0.086212736 -0.015506112 0.023896151 -0.094920808 [16] -0.033474987 0.095879528 0.006142980 -0.004150048 0.018292940 [21] -0.030252512 0.067682612 -0.062630945 0.034870396 0.030721068 [26] 0.021881864 -0.039380067 0.032061725 -0.070848962 0.033242022 [31] -0.016721901 0.042017063 -0.039601514 -0.128569166 0.055966595 [36] 0.007331961 0.008806526 -0.069080597 0.040563855 0.027429259 [41] 0.044775022 -0.036175918 0.026351634 0.052647095 -0.071677118 [46] -0.004129565 0.001594589 0.084479813 -0.012673678 0.017473131 [51] -0.056324844 0.014718330 -0.070604931 0.058779831 -0.008783573 [56] -0.038079890 -0.133559798 0.061009897 -0.012180566 0.080662597 [61] -0.016540982 > (mypacf <- c(rpacf$acf)) [1] -0.219110031 -0.083712032 -0.086492703 0.021230413 -0.033782183 [6] -0.071745982 -0.119639572 0.067292231 -0.076502526 -0.051618941 [11] 0.078858583 -0.011406316 0.027656523 -0.081862678 -0.081724735 [16] 0.047793560 0.044488934 0.030590356 0.016340430 -0.022955351 [21] 0.044237421 -0.026702854 0.037586473 0.035541396 0.063468759 [26] 0.007703857 0.022281616 -0.061631813 -0.023916217 0.022816269 [31] 0.058636577 -0.019762409 -0.157622387 -0.022588851 -0.026898823 [36] 0.007192823 -0.074718815 -0.001809937 0.020166136 0.034495818 [41] 0.003446537 -0.036400138 0.065451584 -0.023075572 0.015515996 [46] -0.008889926 0.052744610 0.030966027 0.084629356 -0.028437858 [51] -0.034625028 -0.063090461 0.062794449 0.044855296 -0.056228009 [56] -0.179872972 -0.021829733 -0.012238460 0.041935329 0.021667983 > lengthx <- length(x) > sqrtn <- sqrt(lengthx) > > #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'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 1:par1) { + 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(mytstat,lengthx),6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/337p31229276944.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-1,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(mytstat,lengthx),6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/43a5u1229276944.tab") > > system("convert tmp/14am01229276943.ps tmp/14am01229276943.png") > system("convert tmp/2hpgx1229276943.ps tmp/2hpgx1229276943.png") > > > proc.time() user system elapsed 0.840 0.426 0.959