R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-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(414.89 + ,444.50 + ,481.29 + ,491.09 + ,419.70 + ,432.88 + ,438.01 + ,412.84 + ,402.91 + ,416.20 + ,411.80 + ,393.60 + ,381.70 + ,388.34 + ,370.89 + ,385.96 + ,394.26 + ,381.37 + ,377.40 + ,377.70 + ,347.71 + ,347.70 + ,340.20 + ,341.00 + ,342.00 + ,319.54 + ,302.79 + ,299.10 + ,313.50 + ,326.80 + ,316.00 + ,316.50 + ,317.20 + ,330.40 + ,323.35 + ,325.85 + ,321.50 + ,321.90 + ,347.48 + ,338.89 + ,345.70 + ,340.44 + ,342.40 + ,342.70 + ,348.34 + ,376.66 + ,417.73 + ,423.51 + ,397.56 + ,390.92 + ,408.26 + ,401.12 + ,408.91 + ,438.35 + ,460.23 + ,449.59 + ,450.52 + ,461.15 + ,460.38 + ,465.35 + ,467.57 + ,486.24 + ,476.58 + ,442.07 + ,443.61 + ,451.55 + ,451.07 + ,451.33 + ,437.63 + ,431.28 + ,413.46 + ,406.78 + ,420.17 + ,419.05 + ,404.01 + ,387.51 + ,390.15 + ,384.40 + ,371.05 + ,367.60 + ,375.04 + ,365.14 + ,361.75 + ,366.88 + ,394.26 + ,409.39 + ,410.11 + ,416.81 + ,393.06 + ,374.24 + ,369.05 + ,352.33 + ,362.53 + ,394.73 + ,389.32 + ,380.74 + ,381.73 + ,376.95 + ,383.64 + ,363.83 + ,363.34 + ,358.38 + ,356.95 + ,366.72 + ,367.69 + ,356.31 + ,348.74 + ,358.69 + ,360.17 + ,361.73 + ,354.45 + ,353.91 + ,344.34 + ,338.62 + ,337.24 + ,340.81 + ,352.72 + ,343.06 + ,345.43 + ,344.38 + ,335.02 + ,334.82 + ,329.01 + ,329.31 + ,330.08 + ,342.15 + ,367.18 + ,371.89 + ,392.19 + ,378.84 + ,355.28 + ,364.18 + ,373.83 + ,383.30 + ,386.88 + ,381.91 + ,384.13 + ,377.27 + ,381.43 + ,385.64 + ,385.49 + ,380.36 + ,391.58 + ,389.77 + ,384.39 + ,379.29 + ,378.55 + ,376.64 + ,382.12 + ,391.03 + ,385.22 + ,387.56 + ,386.23 + ,383.67 + ,383.06 + ,383.14 + ,385.31 + ,387.44 + ,399.45 + ,404.76 + ,396.21 + ,392.85 + ,391.93 + ,385.27 + ,383.47 + ,387.35 + ,383.14 + ,381.07 + ,377.85 + ,369.00 + ,355.11 + ,346.58 + ,351.81 + ,344.47 + ,343.84 + ,340.76 + ,324.10 + ,324.01 + ,322.82 + ,324.87 + ,306.04 + ,288.74 + ,289.10 + ,297.49 + ,295.94 + ,308.29 + ,299.10 + ,292.32 + ,292.87 + ,284.11 + ,288.98 + ,295.93 + ,294.12 + ,291.68 + ,287.08 + ,287.33 + ,285.96 + ,282.62 + ,276.44 + ,261.31 + ,256.08 + ,256.69 + ,264.74 + ,310.72 + ,293.18 + ,283.07 + ,284.32 + ,299.86 + ,286.39 + ,279.69 + ,275.19 + ,285.73 + ,281.59 + ,274.47 + ,273.68 + ,270.00 + ,266.01 + ,271.45 + ,265.49 + ,261.87 + ,263.03 + ,260.48 + ,272.36 + ,270.23 + ,267.53 + ,272.39 + ,283.42 + ,283.06 + ,276.16 + ,275.85 + ,281.51 + ,295.50 + ,294.06 + ,302.68 + ,314.49 + ,321.18 + ,313.29 + ,310.26 + ,319.14 + ,316.56 + ,319.07 + ,331.92 + ,356.86 + ,358.97 + ,340.55 + ,328.18 + ,355.68 + ,356.35 + ,351.02 + ,359.77 + ,378.95 + ,378.92 + ,389.91 + ,406.95 + ,413.79 + ,404.88 + ,406.67 + ,403.26 + ,383.78 + ,392.37 + ,398.09 + ,400.51 + ,405.28 + ,420.46 + ,439.38 + ,442.08 + ,424.03 + ,423.35 + ,433.85 + ,429.23 + ,421.87 + ,430.66 + ,424.48 + ,437.93 + ,456.05 + ,469.90 + ,476.67 + ,510.10 + ,549.86 + ,555.00 + ,557.09 + ,610.65 + ,675.39 + ,596.15 + ,633.71 + ,632.59 + ,598.19 + ,585.78 + ,627.83 + ,629.79 + ,631.17 + ,664.75 + ,654.90 + ,679.37 + ,667.31 + ,655.66 + ,665.38 + ,665.41 + ,712.65 + ,754.60 + ,806.25 + ,803.20 + ,889.60 + ,922.30 + ,968.43 + ,909.71 + ,888.66 + ,889.49 + ,939.77 + ,839.03 + ,829.93 + ,806.62 + ,760.86 + ,816.09 + ,858.69 + ,943.00 + ,924.27 + ,890.20 + ,928.65 + ,945.67 + ,934.23 + ,949.38 + ,996.59 + ,1043.16 + ,1127.04 + ,1134.72 + ,1117.96 + ,1095.41 + ,1113.34 + ,1148.69 + ,1205.43 + ,1232.92 + ,1192.97 + ,1215.81 + ,1270.98 + ,1342.02 + ,1369.89 + ,1390.55 + ,1356.40 + ,1372.73 + ,1424.00 + ,1479.76 + ,1512.60 + ,1528.66 + ,1572.21 + ,1757.21 + ,1770.95 + ,1665.21 + ,1738.11 + ,1641.84 + ,1652.21 + ,1742.14 + ,1673.77 + ,1649.69 + ,1591.19 + ,1598.76 + ,1589.90 + ,1630.31 + ,1744.81 + ,1746.58 + ,1721.64) > par2 = '0' > par1 = '0' > par2 <- '0' > par1 <- '0' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: root > #To cite this work: Wessa P., (2012), Univariate Explorative Data Analysis (v1.0.9) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_edauni.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > # > par1 <- as.numeric(par1) > par2 <- as.numeric(par2) > x <- as.ts(x) > library(lattice) > postscript(file="/var/wessaorg/rcomp/tmp/149ag1355661709.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/2y3e21355661709.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/3njdh1355661709.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/4blnd1355661709.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/538s61355661709.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/6pbc91355661709.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/7x4761355661709.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() + } > summary(x) Min. 1st Qu. Median Mean 3rd Qu. Max. 256.1 340.6 384.4 527.2 467.6 1771.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/8j3dw1355661709.tab") > > try(system("convert tmp/149ag1355661709.ps tmp/149ag1355661709.png",intern=TRUE)) character(0) > try(system("convert tmp/2y3e21355661709.ps tmp/2y3e21355661709.png",intern=TRUE)) character(0) > try(system("convert tmp/3njdh1355661709.ps tmp/3njdh1355661709.png",intern=TRUE)) character(0) > try(system("convert tmp/4blnd1355661709.ps tmp/4blnd1355661709.png",intern=TRUE)) character(0) > try(system("convert tmp/538s61355661709.ps tmp/538s61355661709.png",intern=TRUE)) convert: unable to open image `tmp/538s61355661709.ps': @ error/blob.c/OpenBlob/2587. convert: missing an image filename `tmp/538s61355661709.png' @ error/convert.c/ConvertImageCommand/3011. character(0) attr(,"status") [1] 1 Warning message: running command 'convert tmp/538s61355661709.ps tmp/538s61355661709.png' had status 1 > try(system("convert tmp/6pbc91355661709.ps tmp/6pbc91355661709.png",intern=TRUE)) convert: unable to open image `tmp/6pbc91355661709.ps': @ error/blob.c/OpenBlob/2587. convert: missing an image filename `tmp/6pbc91355661709.png' @ error/convert.c/ConvertImageCommand/3011. character(0) attr(,"status") [1] 1 Warning message: running command 'convert tmp/6pbc91355661709.ps tmp/6pbc91355661709.png' had status 1 > try(system("convert tmp/7x4761355661709.ps tmp/7x4761355661709.png",intern=TRUE)) convert: unable to open image `tmp/7x4761355661709.ps': @ error/blob.c/OpenBlob/2587. convert: missing an image filename `tmp/7x4761355661709.png' @ error/convert.c/ConvertImageCommand/3011. character(0) attr(,"status") [1] 1 Warning message: running command 'convert tmp/7x4761355661709.ps tmp/7x4761355661709.png' had status 1 > > > proc.time() user system elapsed 2.444 0.531 3.031