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(6282929 + ,4324047 + ,4108272 + ,-1212617 + ,1485329 + ,1779876 + ,1367203 + ,2519076 + ,912684 + ,1443586 + ,1220017 + ,984885 + ,1457425 + ,-572920 + ,929144 + ,1151176 + ,790090 + ,774497 + ,990576 + ,454195 + ,876607 + ,711969 + ,702380 + ,264449 + ,450033 + ,541063 + ,588864 + ,-37216 + ,783310 + ,467359 + ,688779 + ,608419 + ,696348 + ,597793 + ,821730 + ,377934 + ,651939 + ,697458 + ,700368 + ,225986 + ,348695 + ,373683 + ,501709 + ,413743 + ,379825 + ,336260 + ,636765 + ,481231 + ,469107 + ,211928 + ,563925 + ,511939 + ,521016 + ,543856 + ,329304 + ,423262 + ,509665 + ,455881 + ,367772 + ,406339 + ,493408 + ,232942 + ,416002 + ,337430 + ,361517 + ,360962 + ,235561 + ,408247 + ,450296 + ,418799 + ,247405 + ,378519 + ,326638 + ,328233 + ,386225 + ,283662 + ,370225 + ,269236 + ,365732 + ,420383 + ,345811 + ,431809 + ,418876 + ,297476 + ,416776 + ,357257 + ,458343 + ,388386 + ,358934 + ,407560 + ,392558 + ,373177 + ,428370 + ,369419 + ,358649 + ,376641 + ,467427 + ,364885 + ,436230 + ,329118 + ,317365 + ,286849 + ,376685 + ,407198 + ,377772 + ,271483 + ,153661 + ,513294 + ,324881 + ,264512 + ,420968 + ,129302 + ,191521 + ,268673 + ,353179 + ,354624 + ,363713 + ,456657 + ,211742 + ,338381 + ,418530 + ,351483 + ,372928 + ,485538 + ,279268 + ,219060 + ,325560 + ,325314 + ,322046 + ,325560 + ,325599 + ,377028 + ,325560 + ,323850 + ,325560 + ,331514 + ,325632 + ,325560 + ,325560 + ,325560 + ,322265 + ,325560 + ,325906 + ,325985 + ,346145 + ,325898 + ,325560 + ,325356 + ,325560 + ,325930 + ,318020 + ,326389 + ,325560 + ,302925 + ,325540 + ,325560 + ,325560 + ,326736 + ,340580 + ,325560 + ,325560 + ,325560 + ,325560 + ,331828 + ,323299 + ,325560 + ,325560 + ,387722 + ,325560 + ,325560 + ,325560 + ,324598 + ,325560 + ,328726 + ,325560 + ,325043 + ,325560 + ,325806 + ,325560 + ,325560 + ,387732 + ,349729 + ,332202 + ,305442 + ,329537 + ,327055 + ,356245 + ,328451 + ,307062 + ,325560 + ,331345 + ,325560 + ,331824 + ,325560 + ,325685 + ,325560 + ,404480 + ,325560 + ,325560 + ,318314 + ,325560 + ,325560 + ,325560 + ,311807 + ,337724 + ,326431 + ,327556 + ,325560 + ,356850 + ,325560 + ,325560 + ,325560 + ,322741 + ,310902 + ,324295 + ,325560 + ,326156 + ,326960 + ,325560 + ,333411 + ,297761 + ,325560 + ,325536 + ,325560 + ,325762 + ,327957 + ,325560 + ,325560 + ,318521 + ,325560 + ,319775 + ,325560 + ,325560 + ,332128 + ,325560 + ,325486 + ,325560 + ,325838 + ,325560 + ,325560 + ,325560 + ,331767 + ,325560 + ,324523 + ,339995 + ,325560 + ,325560 + ,319582 + ,325560 + ,325560 + ,307245 + ,325560 + ,317967 + ,331488 + ,335452 + ,325560 + ,334184 + ,313213 + ,325560 + ,325560 + ,325560 + ,325560 + ,348678 + ,328727 + ,325560 + ,325560 + ,325560 + ,387978 + ,325560 + ,336704 + ,325560 + ,325560 + ,325560 + ,322076 + ,325560 + ,334272 + ,338197 + ,325560 + ,321024 + ,322145 + ,325560 + ,325560 + ,323351 + ,325560 + ,327748 + ,325560 + ,325560 + ,328157 + ,325560 + ,311594 + ,325560 + ,335962 + ,372426 + ,325560 + ,319844 + ,355822 + ,325560 + ,325560 + ,325560 + ,325560 + ,324047 + ,311464 + ,325560 + ,325560 + ,353417 + ,325590 + ,325560 + ,328576 + ,326126 + ,325560 + ,325560 + ,369376 + ,325560 + ,332013 + ,325871 + ,342165 + ,324967 + ,314832 + ,325557 + ,325560 + ,325560 + ,325560 + ,325560 + ,325560 + ,325560 + ,322649 + ,325560 + ,325560 + ,325560 + ,325560 + ,325560 + ,325560 + ,325560 + ,324598 + ,325567 + ,325560 + ,324005 + ,325560 + ,325748 + ,323385 + ,315409 + ,325560 + ,325560 + ,325560 + ,325560 + ,325560 + ,312275 + ,325560 + ,325560 + ,325560 + ,320576 + ,325246 + ,332961 + ,323010 + ,325560 + ,325560 + ,345253 + ,325560 + ,325560 + ,325560 + ,325559 + ,325560 + ,325560 + ,319634 + ,319951 + ,325560 + ,325560 + ,325560 + ,325560 + ,325560 + ,325560 + ,318519 + ,343222 + ,317234 + ,325560 + ,325560 + ,314025 + ,320249 + ,325560 + ,325560 + ,325560 + ,349365 + ,289197 + ,325560 + ,329245 + ,240869 + ,327182 + ,322876 + ,323117 + ,306351 + ,335137 + ,308271 + ,301731 + ,382409 + ,279230 + ,298731 + ,243650 + ,532682 + ,319771 + ,171493 + ,347262 + ,343945 + ,311874 + ,302211 + ,316708 + ,333463 + ,344282 + ,319635 + ,301186 + ,300381 + ,318765 + ,286146 + ,306844 + ,307705 + ,312448 + ,299715 + ,373399 + ,299446 + ,325586 + ,291221 + ,261173 + ,255027 + ,-78375 + ,-58143 + ,227033 + ,235098 + ,21267 + ,238675 + ,197687 + ,418341 + ,-297706) > par2 = '0' > 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/1l8fk1291304557.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/html/rcomp/tmp/2l8fk1291304557.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/html/rcomp/tmp/3l8fk1291304557.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/html/rcomp/tmp/4wze51291304557.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/html/rcomp/tmp/5wze51291304557.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/html/rcomp/tmp/6p8v81291304557.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/html/rcomp/tmp/7zhdb1291304557.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. -1213000 325300 325600 399200 350600 6283000 > > #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/86ir41291304557.tab") > > try(system("convert tmp/1l8fk1291304557.ps tmp/1l8fk1291304557.png",intern=TRUE)) character(0) > try(system("convert tmp/2l8fk1291304557.ps tmp/2l8fk1291304557.png",intern=TRUE)) character(0) > try(system("convert tmp/3l8fk1291304557.ps tmp/3l8fk1291304557.png",intern=TRUE)) character(0) > try(system("convert tmp/4wze51291304557.ps tmp/4wze51291304557.png",intern=TRUE)) character(0) > try(system("convert tmp/5wze51291304557.ps tmp/5wze51291304557.png",intern=TRUE)) convert: unable to open image `tmp/5wze51291304557.ps': No such file or directory. convert: missing an image filename `tmp/5wze51291304557.png'. character(0) > try(system("convert tmp/6p8v81291304557.ps tmp/6p8v81291304557.png",intern=TRUE)) convert: unable to open image `tmp/6p8v81291304557.ps': No such file or directory. convert: missing an image filename `tmp/6p8v81291304557.png'. character(0) > try(system("convert tmp/7zhdb1291304557.ps tmp/7zhdb1291304557.png",intern=TRUE)) convert: unable to open image `tmp/7zhdb1291304557.ps': No such file or directory. convert: missing an image filename `tmp/7zhdb1291304557.png'. character(0) > > > proc.time() user system elapsed 1.245 0.803 5.870