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.

> x <- c(-71.15
+ ,-45.95
+ ,-25.90
+ ,13.72
+ ,47.48
+ ,61.03
+ ,56.78
+ ,50.00
+ ,44.88
+ ,39.80
+ ,29.18
+ ,24.96
+ ,20.74
+ ,17.60
+ ,15.77
+ ,4.47
+ ,-0.26
+ ,-8.24
+ ,-13.05
+ ,-27.34
+ ,-34.40
+ ,-36.76
+ ,-36.04
+ ,-36.14
+ ,-39.13
+ ,-40.25
+ ,-29.97
+ ,-32.19
+ ,-36.46
+ ,-38.92
+ ,-20.32
+ ,2.62
+ ,6.46
+ ,8.60
+ ,10.56
+ ,11.17
+ ,-1.01
+ ,6.41
+ ,0.57
+ ,-9.90
+ ,-5.80
+ ,-5.89
+ ,2.37
+ ,0.89
+ ,-8.46
+ ,-12.31
+ ,-24.04
+ ,-38.97
+ ,-41.46
+ ,-55.24
+ ,-59.23
+ ,-60.20
+ ,-61.37
+ ,-58.50
+ ,-60.15
+ ,-59.90
+ ,-58.99
+ ,-65.57
+ ,-73.24
+ ,-62.78
+ ,-62.45
+ ,-49.27
+ ,-54.35
+ ,-53.94
+ ,-51.51
+ ,-51.24
+ ,-52.49
+ ,-49.15
+ ,-48.83
+ ,-62.10
+ ,-59.95
+ ,-56.50
+ ,-52.50
+ ,-48.58
+ ,-45.70
+ ,-52.91
+ ,-52.66
+ ,-50.91
+ ,-49.70
+ ,-45.44
+ ,-52.43
+ ,-52.38
+ ,-49.30
+ ,-48.82
+ ,-43.81
+ ,-35.68
+ ,-36.15
+ ,-35.09
+ ,-33.17
+ ,-28.47
+ ,-23.15
+ ,-26.59
+ ,-26.28
+ ,-22.77
+ ,-22.47
+ ,-22.89
+ ,-20.24
+ ,-15.04
+ ,-17.27
+ ,-12.20
+ ,-7.86
+ ,-10.46
+ ,-6.44
+ ,-8.55
+ ,-13.10
+ ,-8.56
+ ,-15.69
+ ,-23.45
+ ,-13.44
+ ,52.76
+ ,46.15
+ ,36.75
+ ,31.01
+ ,29.73
+ ,25.37
+ ,14.72
+ ,-1.56
+ ,-1.09
+ ,0.57
+ ,3.27
+ ,3.09
+ ,-2.19
+ ,-4.98
+ ,-11.71
+ ,-18.84
+ ,-17.77
+ ,-18.36
+ ,-25.92
+ ,-20.11
+ ,-16.45
+ ,-31.49
+ ,-32.77
+ ,-22.89
+ ,-26.12
+ ,-29.49
+ ,-23.30
+ ,-23.89
+ ,-22.33
+ ,-13.41
+ ,-15.06
+ ,-11.44
+ ,-12.46
+ ,-16.72
+ ,-10.79
+ ,-7.02
+ ,7.35
+ ,7.69
+ ,3.41
+ ,11.31
+ ,10.21
+ ,16.98
+ ,14.34
+ ,9.89
+ ,4.57
+ ,-9.86
+ ,-3.47
+ ,-2.25
+ ,5.87
+ ,-5.18
+ ,-9.82
+ ,-1.41
+ ,5.45
+ ,15.89
+ ,15.95
+ ,4.18
+ ,6.18
+ ,0.28
+ ,2.62
+ ,5.38
+ ,10.09
+ ,0.11
+ ,-3.62
+ ,-5.99
+ ,-5.16
+ ,-9.21
+ ,-15.16
+ ,-15.60
+ ,-17.46
+ ,-19.63
+ ,-22.10
+ ,-17.80
+ ,-8.88
+ ,-11.16
+ ,-11.84
+ ,-5.85
+ ,-8.13
+ ,1.56
+ ,3.82
+ ,-2.75
+ ,-6.22
+ ,-15.32
+ ,-16.99
+ ,-17.27
+ ,-6.39
+ ,-4.52
+ ,-7.83
+ ,-4.73
+ ,-1.55
+ ,0.87
+ ,-15.54
+ ,-15.22
+ ,-22.01
+ ,-19.70
+ ,-13.54
+ ,-7.68
+ ,-20.38
+ ,-57.94
+ ,-55.69
+ ,-52.24
+ ,-47.71
+ ,28.95
+ ,106.79
+ ,78.46
+ ,72.85
+ ,76.55
+ ,74.06
+ ,79.28
+ ,64.54
+ ,55.77
+ ,54.80
+ ,47.53
+ ,52.95
+ ,53.35
+ ,62.49
+ ,59.69
+ ,52.57
+ ,47.25
+ ,41.22
+ ,35.33
+ ,31.84
+ ,28.67
+ ,31.78
+ ,22.10
+ ,23.63
+ ,12.87
+ ,-15.19
+ ,-24.44
+ ,-27.99
+ ,-34.17
+ ,-17.49
+ ,-9.34
+ ,-42.68
+ ,-103.84
+ ,-47.33
+ ,-70.93
+ ,-39.14
+ ,53.06
+ ,69.47
+ ,73.18
+ ,61.22
+ ,61.47
+ ,33.71
+ ,20.51
+ ,2.39
+ ,36.10
+ ,25.61
+ ,38.65
+ ,32.34
+ ,29.02
+ ,24.59
+ ,31.68
+ ,23.00
+ ,10.90
+ ,5.32
+ ,6.46
+ ,11.12
+ ,13.76
+ ,15.49
+ ,21.39
+ ,13.88
+ ,24.76
+ ,30.25
+ ,34.48
+ ,43.44
+ ,12.67
+ ,-8.83
+ ,27.49
+ ,50.55
+ ,44.15
+ ,46.92
+ ,43.74
+ ,43.65
+ ,45.74
+ ,52.88
+ ,43.93
+ ,39.63
+ ,29.82
+ ,29.52
+ ,36.70
+ ,47.66
+ ,46.08
+ ,32.49
+ ,22.19
+ ,47.71
+ ,42.53
+ ,36.29
+ ,32.53
+ ,42.31
+ ,34.11
+ ,28.32
+ ,26.70
+ ,30.41
+ ,12.19
+ ,18.25
+ ,32.68
+ ,30.14
+ ,46.76
+ ,45.34
+ ,34.55
+ ,29.65
+ ,24.47
+ ,33.19
+ ,41.52
+ ,36.56
+ ,43.39
+ ,46.10
+ ,36.22
+ ,41.93
+ ,35.91
+ ,39.00
+ ,37.94
+ ,30.76
+ ,24.64
+ ,31.50
+ ,21.46
+ ,17.34
+ ,20.53
+ ,21.86
+ ,11.25
+ ,-6.82
+ ,21.49
+ ,19.69
+ ,6.78
+ ,5.49
+ ,-15.81
+ ,-33.73
+ ,-7.85
+ ,-43.92
+ ,-62.69
+ ,-25.88
+ ,-11.68
+ ,9.62
+ ,19.47
+ ,34.38
+ ,39.64
+ ,21.66
+ ,-0.65
+ ,-6.58
+ ,-35.05
+ ,-23.54
+ ,0.29
+ ,-2.76
+ ,9.33
+ ,3.09
+ ,0.41
+ ,-6.67
+ ,-2.15
+ ,-3.03
+ ,-23.45
+ ,-42.26)
> par2 = '36'
> par1 = '0'
> #'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!)
> par1 <- as.numeric(par1)
> par2 <- as.numeric(par2)
> x <- as.ts(x)
> library(lattice)
> postscript(file="/var/www/html/rcomp/tmp/1umd31256745977.ps",horizontal=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/www/html/rcomp/tmp/2rm7v1256745977.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> hist(x)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/3izlt1256745977.ps",horizontal=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/www/html/rcomp/tmp/4cs1h1256745977.ps",horizontal=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/www/html/rcomp/tmp/5cpb31256745977.ps",horizontal=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/www/html/rcomp/tmp/6lz9q1256745977.ps",horizontal=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/www/html/rcomp/tmp/740u51256745977.ps",horizontal=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. 
-103.8000  -22.8900   -1.4800    0.1305   28.4100  106.8000 
> 
> #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,'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/www/html/rcomp/tmp/89l5u1256745977.tab") 
> 
> system("convert tmp/1umd31256745977.ps tmp/1umd31256745977.png")
> system("convert tmp/2rm7v1256745977.ps tmp/2rm7v1256745977.png")
> system("convert tmp/3izlt1256745977.ps tmp/3izlt1256745977.png")
> system("convert tmp/4cs1h1256745977.ps tmp/4cs1h1256745977.png")
> system("convert tmp/5cpb31256745977.ps tmp/5cpb31256745977.png")
> system("convert tmp/6lz9q1256745977.ps tmp/6lz9q1256745977.png")
> system("convert tmp/740u51256745977.ps tmp/740u51256745977.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.687   1.169   3.145