R version 2.12.0 (2010-10-15) Copyright (C) 2010 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(321.61 + ,345.85 + ,338.60 + ,345.64 + ,340.71 + ,342.49 + ,342.65 + ,348.68 + ,377.36 + ,418.05 + ,423.13 + ,397.69 + ,390.80 + ,408.29 + ,401.02 + ,409.24 + ,439.28 + ,459.95 + ,449.66 + ,451.14 + ,460.66 + ,460.23 + ,465.69 + ,468.01 + ,486.74 + ,475.89 + ,441.52 + ,443.63 + ,451.62 + ,451.14 + ,450.88 + ,437.56 + ,431.18 + ,412.02 + ,407.14 + ,420.48 + ,418.92 + ,403.57 + ,387.55 + ,390.02 + ,384.37 + ,370.62 + ,367.90 + ,374.91 + ,365.15 + ,361.91 + ,367.03 + ,395.19 + ,409.25 + ,410.49 + ,416.58 + ,392.21 + ,374.29 + ,369.05 + ,352.19 + ,362.85 + ,395.47 + ,388.82 + ,380.39 + ,381.73 + ,377.69 + ,383.04 + ,363.89 + ,363.23 + ,358.37 + ,356.97 + ,366.87 + ,367.57 + ,355.88 + ,348.88 + ,358.77 + ,360.42 + ,361.08 + ,354.57 + ,353.73 + ,344.20 + ,338.34 + ,337.21 + ,340.96 + ,353.29 + ,342.67 + ,345.71 + ,344.17 + ,334.92 + ,334.81 + ,329.05 + ,329.31 + ,330.25 + ,341.89 + ,367.74 + ,371.93 + ,392.79 + ,377.97 + ,354.93 + ,364.40 + ,374.05 + ,383.63 + ,386.56 + ,381.90 + ,384.08 + ,377.29 + ,381.54 + ,385.60 + ,385.47 + ,380.40 + ,391.74 + ,389.57 + ,384.29 + ,379.26 + ,378.44 + ,376.63 + ,382.48 + ,390.89 + ,385.04 + ,387.58 + ,386.19 + ,383.78 + ,383.10 + ,383.25 + ,385.19 + ,387.35 + ,400.49 + ,404.53 + ,396.15 + ,392.79 + ,391.96 + ,385.04 + ,383.58 + ,387.46 + ,382.90 + ,381.04 + ,377.69 + ,368.95 + ,353.87 + ,347.03 + ,351.49 + ,344.23 + ,344.09 + ,340.51 + ,323.90 + ,324.02 + ,323.11 + ,324.36 + ,305.55 + ,288.59 + ,289.15 + ,297.49 + ,295.94 + ,308.29 + ,299.10 + ,292.32 + ,292.87 + ,284.11 + ,288.98 + ,295.93 + ,294.17 + ,291.68 + ,287.07 + ,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 + ,269.82 + ,267.53 + ,272.39 + ,283.42 + ,283.06 + ,276.16 + ,275.85 + ,281.51 + ,295.50 + ,294.06 + ,302.68 + ,314.58 + ,321.18 + ,313.29 + ,310.25 + ,319.14 + ,316.56 + ,319.07 + ,331.92 + ,356.86 + ,358.97 + ,340.55 + ,328.18 + ,355.68 + ,356.35 + ,350.99 + ,359.77 + ,378.95 + ,378.92 + ,389.91 + ,406.11 + ,413.79 + ,404.95 + ,406.67 + ,403.26 + ,383.78 + ,392.48 + ,398.09 + ,400.51 + ,405.28 + ,420.46 + ,439.38 + ,442.08 + ,424.03 + ,423.35 + ,434.32 + ,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.33 + ,598.06 + ,585.78 + ,627.83 + ,629.42 + ,631.17 + ,664.75 + ,654.90 + ,679.37 + ,666.92 + ,655.49 + ,665.30 + ,665.41 + ,712.65 + ,754.60 + ,806.25 + ,803.20 + ,889.60 + ,922.30 + ,968.43 + ,909.70 + ,890.51 + ,889.49 + ,939.77 + ,838.31 + ,829.93 + ,806.62 + ,760.86 + ,822.00 + ,859.19 + ,943.16 + ,924.27 + ,889.50 + ,930.20 + ,945.67 + ,934.23 + ,949.67 + ,996.59 + ,1043.16 + ,1127.04 + ,1126.22 + ,1116.51 + ,1095.41 + ,1113.34 + ,1148.69 + ,1205.43 + ,1232.92 + ,1192.97 + ,1215.81 + ,1270.98 + ,1342.02) > 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/rcomp/tmp/17gut1291755603.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/www/rcomp/tmp/2i7bw1291755603.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/www/rcomp/tmp/3i7bw1291755603.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/www/rcomp/tmp/4ahaz1291755603.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/www/rcomp/tmp/5l8ak1291755603.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/www/rcomp/tmp/6l8ak1291755603.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/www/rcomp/tmp/77r881291755603.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. 256.1 338.5 383.0 457.0 442.9 1342.0 > > #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/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/rcomp/tmp/83ioh1291755603.tab") > try(system("convert tmp/17gut1291755603.ps tmp/17gut1291755603.png",intern=TRUE)) character(0) > try(system("convert tmp/2i7bw1291755603.ps tmp/2i7bw1291755603.png",intern=TRUE)) character(0) > try(system("convert tmp/3i7bw1291755603.ps tmp/3i7bw1291755603.png",intern=TRUE)) character(0) > try(system("convert tmp/4ahaz1291755603.ps tmp/4ahaz1291755603.png",intern=TRUE)) character(0) > try(system("convert tmp/5l8ak1291755603.ps tmp/5l8ak1291755603.png",intern=TRUE)) character(0) > try(system("convert tmp/6l8ak1291755603.ps tmp/6l8ak1291755603.png",intern=TRUE)) character(0) > try(system("convert tmp/77r881291755603.ps tmp/77r881291755603.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.970 0.990 2.968