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 <- c(315.71
+ ,317.45
+ ,317.5
+ ,317.12
+ ,315.86
+ ,314.93
+ ,313.2
+ ,312.6
+ ,313.33
+ ,314.67
+ ,315.62
+ ,316.38
+ ,316.71
+ ,317.72
+ ,318.29
+ ,318.16
+ ,316.55
+ ,314.8
+ ,313.84
+ ,313.26
+ ,314.8
+ ,315.59
+ ,316.43
+ ,316.97
+ ,317.58
+ ,319.02
+ ,320.02
+ ,319.59
+ ,318.18
+ ,315.91
+ ,314.16
+ ,313.83
+ ,315
+ ,316.19
+ ,316.93
+ ,317.7
+ ,318.54
+ ,319.48
+ ,320.58
+ ,319.77
+ ,318.58
+ ,316.79
+ ,314.8
+ ,315.38
+ ,316.1
+ ,317.01
+ ,317.94
+ ,318.55
+ ,319.68
+ ,320.63
+ ,321.01
+ ,320.55
+ ,319.58
+ ,317.4
+ ,316.26
+ ,315.42
+ ,316.69
+ ,317.7
+ ,318.74
+ ,319.08
+ ,319.86
+ ,321.39
+ ,322.24
+ ,321.47
+ ,319.74
+ ,317.77
+ ,316.21
+ ,315.99
+ ,317.12
+ ,318.31
+ ,319.57
+ ,320.08
+ ,320.75
+ ,321.8
+ ,322.24
+ ,321.89
+ ,320.44
+ ,318.7
+ ,316.7
+ ,316.79
+ ,317.79
+ ,318.71
+ ,319.44
+ ,320.44
+ ,320.89
+ ,322.13
+ ,322.16
+ ,321.87
+ ,321.39
+ ,318.8
+ ,317.81
+ ,317.3
+ ,318.87
+ ,319.42
+ ,320.62
+ ,321.59
+ ,322.39
+ ,323.87
+ ,324.01
+ ,323.75
+ ,322.4
+ ,320.37
+ ,318.64
+ ,318.1
+ ,319.78
+ ,321.08
+ ,322.06
+ ,322.5
+ ,323.04
+ ,324.42
+ ,325
+ ,324.09
+ ,322.55
+ ,320.92
+ ,319.31
+ ,319.31
+ ,320.72
+ ,321.96
+ ,322.57
+ ,323.15
+ ,323.89
+ ,325.02
+ ,325.57
+ ,325.36
+ ,324.14
+ ,322.03
+ ,320.41
+ ,320.25
+ ,321.31
+ ,322.84
+ ,324
+ ,324.42
+ ,325.64
+ ,326.66
+ ,327.34
+ ,326.76
+ ,325.88
+ ,323.67
+ ,322.38
+ ,321.78
+ ,322.85
+ ,324.12
+ ,325.03
+ ,325.99
+ ,326.87
+ ,328.14
+ ,328.07
+ ,327.66
+ ,326.35
+ ,324.69
+ ,323.1
+ ,323.16
+ ,323.98
+ ,325.13
+ ,326.17
+ ,326.68
+ ,327.18
+ ,327.78
+ ,328.92
+ ,328.57
+ ,327.34
+ ,325.46
+ ,323.36
+ ,323.56
+ ,324.8
+ ,326.01
+ ,326.77
+ ,327.63
+ ,327.75
+ ,329.72
+ ,330.07
+ ,329.09
+ ,328.05
+ ,326.32
+ ,324.93
+ ,325.06
+ ,326.5
+ ,327.55
+ ,328.55
+ ,329.56
+ ,330.3
+ ,331.5
+ ,332.48
+ ,332.07
+ ,330.87
+ ,329.31
+ ,327.51
+ ,327.18
+ ,328.16
+ ,328.64
+ ,329.35
+ ,330.71
+ ,331.48
+ ,332.65
+ ,333.15
+ ,332.13
+ ,330.99
+ ,329.17
+ ,327.41
+ ,327.21
+ ,328.34
+ ,329.5
+ ,330.68
+ ,331.41
+ ,331.85
+ ,333.29
+ ,333.91
+ ,333.4
+ ,331.74
+ ,329.88
+ ,328.57
+ ,328.35
+ ,329.33
+ ,330.58
+ ,331.66
+ ,332.75
+ ,333.46
+ ,334.78
+ ,334.79
+ ,334.05
+ ,332.95
+ ,330.64
+ ,328.96
+ ,328.77
+ ,330.18
+ ,331.65
+ ,332.69
+ ,333.23
+ ,334.97
+ ,336.03
+ ,336.82
+ ,336.1
+ ,334.79
+ ,332.53
+ ,331.19
+ ,331.21
+ ,332.35
+ ,333.47
+ ,335.09
+ ,335.26
+ ,336.62
+ ,337.77
+ ,338
+ ,337.98
+ ,336.48
+ ,334.37
+ ,332.33
+ ,332.4
+ ,333.76
+ ,334.83
+ ,336.21
+ ,336.64
+ ,338.13
+ ,338.96
+ ,339.02
+ ,339.2
+ ,337.6
+ ,335.56
+ ,333.93
+ ,334.12
+ ,335.26
+ ,336.77
+ ,337.8
+ ,338.28
+ ,340.04
+ ,340.86
+ ,341.47
+ ,341.26
+ ,339.34
+ ,337.45
+ ,336.1
+ ,336.05
+ ,337.21
+ ,338.29
+ ,339.36
+ ,340.51
+ ,341.57
+ ,342.56
+ ,343.01
+ ,342.52
+ ,340.71
+ ,338.51
+ ,336.96
+ ,337.13
+ ,338.58
+ ,339.91
+ ,340.92
+ ,341.69
+ ,342.87
+ ,343.83
+ ,344.3
+ ,343.42
+ ,341.85
+ ,339.82
+ ,337.98
+ ,338.09
+ ,339.24
+ ,340.67
+ ,341.42
+ ,342.67
+ ,343.45
+ ,345.08
+ ,345.76
+ ,345.32
+ ,343.93
+ ,342.08
+ ,340
+ ,340.12
+ ,341.35
+ ,342.89
+ ,343.87
+ ,344.59
+ ,345.29
+ ,346.59
+ ,347.36
+ ,346.8
+ ,345.37
+ ,343.06
+ ,341.24
+ ,341.54
+ ,342.9
+ ,344.36
+ ,345.08
+ ,345.89
+ ,347.49
+ ,348.02
+ ,348.75
+ ,348.19
+ ,346.49
+ ,344.7
+ ,343.04
+ ,342.92
+ ,344.22
+ ,345.61
+ ,346.42
+ ,346.95
+ ,347.88
+ ,349.57
+ ,350.35
+ ,349.7
+ ,347.78
+ ,345.89
+ ,344.88
+ ,344.34
+ ,345.67
+ ,346.89
+ ,348.2
+ ,348.55
+ ,349.56
+ ,351.12
+ ,351.84
+ ,351.45
+ ,349.77
+ ,347.62
+ ,346.37
+ ,346.48
+ ,347.8
+ ,349.03
+ ,350.23
+ ,351.58
+ ,352.22
+ ,353.53
+ ,354.14
+ ,353.64
+ ,352.53
+ ,350.42
+ ,348.84
+ ,348.94
+ ,349.99
+ ,351.29
+ ,352.72
+ ,353.1
+ ,353.64
+ ,355.43
+ ,355.7
+ ,355.11
+ ,353.79
+ ,351.42
+ ,349.83
+ ,350.1
+ ,351.26
+ ,352.66
+ ,353.63
+ ,354.72
+ ,355.49
+ ,356.1
+ ,357.08
+ ,356.11
+ ,354.67
+ ,352.67
+ ,351.05
+ ,351.36
+ ,352.81
+ ,354.21
+ ,354.87
+ ,355.67
+ ,357
+ ,358.4
+ ,359
+ ,357.99
+ ,355.96
+ ,353.78
+ ,352.2
+ ,352.22
+ ,353.7
+ ,354.98
+ ,356.08
+ ,356.84
+ ,357.73
+ ,358.91
+ ,359.45
+ ,359.19
+ ,356.72
+ ,354.77
+ ,352.8
+ ,353.21
+ ,354.15
+ ,355.39
+ ,356.76
+ ,357.17
+ ,358.26
+ ,359.17
+ ,360.07
+ ,359.41
+ ,357.36
+ ,355.29
+ ,353.96
+ ,354.03
+ ,355.27
+ ,356.7
+ ,358.05
+ ,358.8
+ ,359.67
+ ,361.13
+ ,361.48
+ ,360.6
+ ,359.2
+ ,357.23
+ ,355.42
+ ,355.89
+ ,357.41
+ ,358.74
+ ,359.73
+ ,360.61
+ ,361.6
+ ,363.05
+ ,363.62
+ ,363.03
+ ,361.55
+ ,358.94
+ ,357.93
+ ,357.8
+ ,359.22
+ ,360.42
+ ,361.83
+ ,362.94
+ ,363.91
+ ,364.28
+ ,364.93
+ ,364.7
+ ,363.31
+ ,361.15
+ ,359.41
+ ,359.34
+ ,360.62
+ ,361.96
+ ,362.81
+ ,363.87
+ ,364.25
+ ,366.02
+ ,366.47
+ ,365.37
+ ,364.1
+ ,361.89
+ ,360.05
+ ,360.49
+ ,362.21
+ ,364.12
+ ,365
+ ,365.82
+ ,366.95
+ ,368.42
+ ,369.33
+ ,368.78
+ ,367.59
+ ,365.81
+ ,363.83
+ ,364.18
+ ,365.36
+ ,366.88
+ ,367.97
+ ,368.83
+ ,369.46
+ ,370.77
+ ,370.66
+ ,370.1
+ ,369.1
+ ,366.7
+ ,364.61
+ ,365.17
+ ,366.51
+ ,367.86
+ ,369.07
+ ,369.32
+ ,370.38
+ ,371.63
+ ,371.32
+ ,371.51
+ ,369.69
+ ,368.18
+ ,366.87
+ ,366.94
+ ,368.27
+ ,369.62
+ ,370.47
+ ,371.44
+ ,372.39
+ ,373.32
+ ,373.77
+ ,373.13
+ ,371.51
+ ,369.59
+ ,368.12
+ ,368.38
+ ,369.64
+ ,371.11
+ ,372.38
+ ,373.08
+ ,373.87
+ ,374.93
+ ,375.58
+ ,375.44
+ ,373.91
+ ,371.77
+ ,370.72
+ ,370.5
+ ,372.19
+ ,373.71
+ ,374.92
+ ,375.63
+ ,376.51
+ ,377.75
+ ,378.54
+ ,378.21
+ ,376.65
+ ,374.28
+ ,373.12
+ ,373.1
+ ,374.67
+ ,375.97
+ ,377.03
+ ,377.87
+ ,378.88
+ ,380.42
+ ,380.62
+ ,379.66
+ ,377.48
+ ,376.07
+ ,374.1
+ ,374.47
+ ,376.15
+ ,377.51
+ ,378.43
+ ,379.7
+ ,380.91
+ ,382.2
+ ,382.45
+ ,382.14
+ ,380.6
+ ,378.6
+ ,376.72
+ ,376.98
+ ,378.29
+ ,380.07
+ ,381.36
+ ,382.19
+ ,382.65
+ ,384.65
+ ,384.94
+ ,384.01
+ ,382.15
+ ,380.33
+ ,378.81
+ ,379.06
+ ,380.17
+ ,381.85
+ ,382.88
+ ,383.77
+ ,384.42
+ ,386.36
+ ,386.53
+ ,386.01
+ ,384.45
+ ,381.96
+ ,380.81
+ ,381.09
+ ,382.37
+ ,383.84
+ ,385.42
+ ,385.72
+ ,385.96
+ ,387.18
+ ,388.5
+ ,387.88
+ ,386.38
+ ,384.15
+ ,383.07
+ ,382.98
+ ,384.11
+ ,385.54
+ ,386.92
+ ,387.41
+ ,388.77
+ ,389.46
+ ,390.18
+ ,389.43
+ ,387.74
+ ,385.91
+ ,384.77
+ ,384.38
+ ,385.99
+ ,387.26
+ ,388.45
+ ,389.7
+ ,391.08
+ ,392.46
+ ,392.96
+ ,392.03
+ ,390.13
+ ,388.15
+ ,386.8
+ ,387.18
+ ,388.59)
> par2 = '36'
> par1 = '0'
> par1 <- as.numeric(par1)
> par2 <- as.numeric(par2)
> x <- as.ts(x)
> library(lattice)
> postscript(file="/var/wessaorg/rcomp/tmp/1br281321999313.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(x,type='l',main='Run Sequence Plot',xlab='time or index',ylab='value')
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/2d2sg1321999313.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> hist(x)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/3fy391321999313.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> if (par1 > 0)
+ {
+ densityplot(~x,col='black',main=paste('Density Plot   bw = ',par1),bw=par1)
+ } else {
+ densityplot(~x,col='black',main='Density Plot')
+ }
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/4bak91321999313.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> qqnorm(x)
> qqline(x)
> grid()
> dev.off()
null device 
          1 
> if (par2 > 0)
+ {
+ postscript(file="/var/wessaorg/rcomp/tmp/5ux2c1321999313.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
+ dum <- cbind(lag(x,k=1),x)
+ dum
+ dum1 <- dum[2:length(x),]
+ dum1
+ z <- as.data.frame(dum1)
+ z
+ plot(z,main='Lag plot (k=1), lowess, and regression line')
+ lines(lowess(z))
+ abline(lm(z))
+ dev.off()
+ if (par2 > 1) {
+ postscript(file="/var/wessaorg/rcomp/tmp/64cxv1321999313.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
+ dum <- cbind(lag(x,k=par2),x)
+ dum
+ dum1 <- dum[(par2+1):length(x),]
+ dum1
+ z <- as.data.frame(dum1)
+ z
+ mylagtitle <- 'Lag plot (k='
+ mylagtitle <- paste(mylagtitle,par2,sep='')
+ mylagtitle <- paste(mylagtitle,'), and lowess',sep='')
+ plot(z,main=mylagtitle)
+ lines(lowess(z))
+ dev.off()
+ }
+ postscript(file="/var/wessaorg/rcomp/tmp/7kk2u1321999313.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
+ acf(x,lag.max=par2,main='Autocorrelation Function')
+ grid()
+ dev.off()
+ }
null device 
          1 
> summary(x)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  312.6   326.7   344.4   346.7   364.3   393.0 
> 
> #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,'Descriptive Statistics',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'# observations',header=TRUE)
> a<-table.element(a,length(x))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'minimum',header=TRUE)
> a<-table.element(a,min(x))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Q1',header=TRUE)
> a<-table.element(a,quantile(x,0.25))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'median',header=TRUE)
> a<-table.element(a,median(x))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'mean',header=TRUE)
> a<-table.element(a,mean(x))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Q3',header=TRUE)
> a<-table.element(a,quantile(x,0.75))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'maximum',header=TRUE)
> a<-table.element(a,max(x))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/84hnf1321999313.tab") 
> 
> try(system("convert tmp/1br281321999313.ps tmp/1br281321999313.png",intern=TRUE))
character(0)
> try(system("convert tmp/2d2sg1321999313.ps tmp/2d2sg1321999313.png",intern=TRUE))
character(0)
> try(system("convert tmp/3fy391321999313.ps tmp/3fy391321999313.png",intern=TRUE))
character(0)
> try(system("convert tmp/4bak91321999313.ps tmp/4bak91321999313.png",intern=TRUE))
character(0)
> try(system("convert tmp/5ux2c1321999313.ps tmp/5ux2c1321999313.png",intern=TRUE))
character(0)
> try(system("convert tmp/64cxv1321999313.ps tmp/64cxv1321999313.png",intern=TRUE))
character(0)
> try(system("convert tmp/7kk2u1321999313.ps tmp/7kk2u1321999313.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.394   0.330   2.731