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(-0.0554261742134966 + ,-0.0103723682328288 + ,-0.09774824332988 + ,0.0682838841648196 + ,-0.0598771056498073 + ,-0.0115369164551084 + ,0.0737440480162446 + ,0.0652896783949673 + ,0.0212699541543462 + ,-0.00207374705171083 + ,0.110078189837195 + ,-0.0146254989919655 + ,0.114698886599707 + ,0.122494579401249 + ,-0.093917668042089 + ,0.317437918041701 + ,0.0233064257137833 + ,0.0907672748631918 + ,0.137726863879236 + ,0.0413595313391135 + ,0.0478723668655492 + ,0.0876104346350919 + ,0.0384005884923821 + ,0.0522380787977494 + ,0.0585774158784047 + ,-0.0318052895741221 + ,-0.181701157975095 + ,0.312417016566173 + ,-0.161245884995846 + ,-0.0498906070653545 + ,0.0164895011305049 + ,-0.106532442482812 + ,-0.0111252289051036 + ,-0.0112051724184589 + ,-0.0621266808603795 + ,-0.0194747117256815 + ,-0.005215993087672 + ,0.00868532947685446 + ,0.0189643082939415 + ,0.132126779126779 + ,0.114222267764316 + ,-0.0512424762283885 + ,0.111763287181735 + ,-0.0727789247238114 + ,0.099056920984902 + ,-0.0133229962434981 + ,-0.0540901157284234 + ,0.0131003534096119 + ,-0.105083618822970 + ,-0.0532746009105622 + ,-0.152603426068890 + ,0.169880388606489 + ,-0.134937086209412 + ,-0.156403695781353 + ,0.0220625418036537 + ,-0.122447540039379 + ,-0.147301939947001 + ,0.124893101691225 + ,-0.0981932532160595 + ,-0.0338473507167997 + ,0.120285012723111 + ,0.0619113501360601 + ,-0.127397173085151 + ,0.295513148096302 + ,-0.00395731088007977 + ,-0.134137032646724 + ,0.0800913155281188 + ,0.0420046126531497 + ,0.061555085014176 + ,-0.0129403767474939 + ,0.0972860535162994 + ,0.0877108066588468 + ,-0.0356206200854371 + ,0.0505485461316185 + ,-0.131893718779630 + ,0.343673596883514 + ,-0.00930805173243436 + ,-0.244200397859937 + ,-0.054917088296814 + ,0.0337991899193833 + ,-0.0441366211382616 + ,-0.0401696208416182 + ,0.0586113053187203 + ,-0.0362167322426645 + ,0.0214944435265415 + ,-0.0394242009278998 + ,-0.219189486789616 + ,-0.038750207606925 + ,-0.305953588706861 + ,-0.160680743573754 + ,-0.129309183705573 + ,-0.175325132712971 + ,-0.285383831304565 + ,-0.182559850590123 + ,-0.224103473978491 + ,0.0611619478196247 + ,-0.253257239400946 + ,-0.0667794801013696 + ,-0.0203659156623685 + ,0.0390230291850474 + ,-0.362690960075248 + ,0.0333119163149927 + ,-0.0749260168888494 + ,-0.310091538794336 + ,-0.313174314082963 + ,-0.0446943832910636 + ,0.0884651538474875 + ,0.148008967470383 + ,-0.0581099037570537 + ,-0.265335985463967 + ,-0.142815914648150 + ,-0.0634276688751165 + ,0.301971033848160 + ,0.509200438595075 + ,-0.226853203690718 + ,0.132746336473417 + ,0.113478705894465 + ,-0.0714550273266129 + ,0.107207329251342 + ,-0.070724916726179 + ,-0.0351999533754058 + ,0.0397456022551489 + ,0.211380404750537 + ,-0.349347916179288 + ,0.270214046482663 + ,0.281727228579061 + ,0.0814427167227327 + ,0.0585457504519322 + ,0.131338788505694 + ,0.0584680902381685 + ,0.0145094042468520 + ,0.131142061479022 + ,-0.118748120650334 + ,0.161971940232045 + ,0.217140021180596 + ,-0.36075238703227 + ,0.112369257289203 + ,0.177091978202189 + ,0.0560949481137286 + ,0.141966781350567 + ,0.238943521016553 + ,-0.0169617673898806 + ,-0.217034348244890 + ,-0.386510308927318 + ,0.119169209082527 + ,0.239427960434481 + ,0.619747125797162 + ,-0.229810610454884 + ,0.0509120935252796 + ,-0.212011167464921 + ,-0.0887441758042202 + ,0.216923527514844 + ,-0.0205975632055718 + ,0.047788699609204 + ,0.031052755829211 + ,-0.0774958020653822 + ,0.210654616027303 + ,-0.0984341615804984 + ,0.232160736649630 + ,-0.318032671100063 + ,0.0466865439469646 + ,-0.0217807544657625 + ,0.0373144926595483 + ,0.0649938009681715 + ,0.108204156147779 + ,0.0766226157267711 + ,0.109937091688754 + ,0.145533105760753 + ,0.0217820395557200 + ,-0.119359221276209 + ,-0.131759892290998 + ,-0.318234299422278 + ,0.118288319679319 + ,-0.0502150450132071 + ,-0.0419800466493269 + ,-0.0104536303122343) > 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/1kj8c1261513755.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/2kc981261513755.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/3f3c31261513755.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/4582u1261513755.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/593b11261513755.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/6wh8m1261513755.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/754b41261513755.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. -3.865e-01 -8.031e-02 -3.016e-03 1.584e-18 8.904e-02 6.197e-01 > > #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/87eu21261513755.tab") > try(system("convert tmp/1kj8c1261513755.ps tmp/1kj8c1261513755.png",intern=TRUE)) character(0) > try(system("convert tmp/2kc981261513755.ps tmp/2kc981261513755.png",intern=TRUE)) character(0) > try(system("convert tmp/3f3c31261513755.ps tmp/3f3c31261513755.png",intern=TRUE)) character(0) > try(system("convert tmp/4582u1261513755.ps tmp/4582u1261513755.png",intern=TRUE)) character(0) > try(system("convert tmp/593b11261513755.ps tmp/593b11261513755.png",intern=TRUE)) character(0) > try(system("convert tmp/6wh8m1261513755.ps tmp/6wh8m1261513755.png",intern=TRUE)) character(0) > try(system("convert tmp/754b41261513755.ps tmp/754b41261513755.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.528 1.086 2.071