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.325268940652250 + ,0.174803600738832 + ,0.501377605018886 + ,0.0859193298655126 + ,-0.987159360948653 + ,1.60980144119294 + ,0.171696135380085 + ,0.830277801292865 + ,0.620654556971156 + ,-1.01610741303737 + ,1.79309605723783 + ,0.506513108486633 + ,0.267689829297541 + ,0.733670677239202 + ,0.224067410538229 + ,1.47684633432560 + ,-0.851788714449576 + ,0.797595345615805 + ,-0.27480267237158 + ,0.0838834586911507 + ,-0.121505653550161 + ,0.631202199065368 + ,0.0782854866315957 + ,0.764708828909565 + ,-0.175870900848487 + ,0.303658167461428 + ,0.334813565364407 + ,-0.78903961133057 + ,0.548615446736318 + ,1.10803873923536 + ,-0.317891120456626 + ,0.9901703953356 + ,-0.384820110488252 + ,0.368519360304957 + ,0.254786797238417 + ,-0.174145843270876 + ,0.67535792035457 + ,1.44301427700514 + ,0.760082666024789 + ,1.38938642804491 + ,-0.464062381606388 + ,0.425484373422507 + ,0.391704064700779 + ,0.0515364779754748 + ,0.972875806223118 + ,0.55823387433218 + ,1.03687475153784 + ,-1.73799567103410 + ,0.0273213342655299 + ,-1.14728571478772 + ,-0.270507966751069 + ,-0.679168823987675 + ,-0.939028021983555 + ,-0.2782918641341 + ,-0.137673442418869 + ,-0.612554173793871 + ,-1.40158680260588 + ,0.824977600305121 + ,-0.254161313717333 + ,-0.474379671661846 + ,-1.45888578840602 + ,0.158712807341798 + ,0.42684360859133 + ,0.533838048801584 + ,-0.62756770632153 + ,0.62119253835102 + ,-0.596208725794837 + ,0.307750900244233 + ,0.88681798267684 + ,0.892590640076497 + ,0.661243897634714 + ,0.511811081466187 + ,-1.42764204541458 + ,0.203723265798547 + ,0.722114552246169 + ,-1.78484681219131 + ,-1.02745593312050 + ,-0.0995958911573627 + ,-1.42936754519820 + ,-1.57611650808341 + ,0.854828513620954 + ,0.586553281789383 + ,-1.16196477262354 + ,-2.04343799438004 + ,0.405786986779646 + ,-0.372184126642192 + ,-0.104779273295343 + ,1.82824440149189 + ,1.16865912540444 + ,0.356830409107125 + ,0.119057088663047 + ,1.77565157768908 + ,-1.49131443486346 + ,-0.200217727321004 + ,0.0298928204546468 + ,-1.67852086311231 + ,-1.01495199118528 + ,0.255173067016221 + ,-0.672832416453207 + ,0.659128650175414 + ,-0.0965875846564268 + ,-0.646183795439774 + ,-0.0434058058236342 + ,-1.52225523739414 + ,0.383278501531744 + ,2.54996816693744 + ,-1.40741391717925 + ,-2.48687186705108 + ,-1.01554045514841 + ,0.509294766709254 + ,1.24143160200711 + ,-0.185246318721566 + ,0.39688456930075 + ,-0.242001609810631 + ,0.716286944222023 + ,0.0699549556395686 + ,-0.138806717848701 + ,0.615219389201853 + ,-1.77249350687623 + ,0.087103545072493 + ,-0.96891948815133 + ,-0.50960721877958 + ,0.404822542358774 + ,-0.0736712799438494 + ,-0.905700376721695 + ,1.79089536112799 + ,-0.790284604706349 + ,0.221173229032318 + ,-0.140361653105643 + ,1.80511090275527 + ,-1.15473965738908 + ,0.317515548081165 + ,-1.34186790295649 + ,-1.02981734533104 + ,1.48189631421005 + ,0.875905387613113 + ,0.95472158757434 + ,0.361064161739099 + ,1.04329527780004 + ,-0.876172593661597 + ,0.42120464743111 + ,-0.698371027826438 + ,3.54754447079534 + ,0.45117847934413 + ,-3.48538269140986 + ,-0.155096043035555 + ,-0.079357742699155 + ,1.52281022071598 + ,1.43570676791519 + ,-0.615276485856994 + ,0.356078442212051 + ,0.541469521644917 + ,1.50492205932887 + ,-3.7720621054139 + ,0.213632092332079 + ,1.33994771296635 + ,2.67825338773539 + ,0.063316954315315 + ,0.367484253435371 + ,-0.0283276572277870 + ,-0.110324673504478 + ,1.22374472834180 + ,-0.435851235691176 + ,0.934970648285374 + ,0.577832959108354 + ,-1.20771053506042 + ,0.292021166283429 + ,-1.35108168774920 + ,-1.22130727655887 + ,1.08234779980742 + ,0.464910763590748 + ,0.802771961019954 + ,0.525532038555883 + ,-1.47540220075188 + ,0.344476626736514 + ,0.541653951215843 + ,-1.38774551083857 + ,0.365170239831711 + ,1.29208061611421 + ,2.06853524478187 + ,-0.407308651738156 + ,-1.19873665066815 + ,-2.07543702306924 + ,0.386589322014422 + ,0.753911828685413 + ,-0.123281903338415 + ,0.548886717595287 + ,-0.265125891167525 + ,-0.0222764194572018 + ,-2.09084211290852 + ,0.440453670194936 + ,-1.77122200580385 + ,0.284848672920714 + ,0.316954751725844 + ,-0.244506761135233 + ,0.5319161211438 + ,-0.288866519494615 + ,1.01463395174541 + ,1.37270239692662 + ,-0.301359369584999 + ,-1.19603788903969 + ,-1.47897908540647 + ,-0.414176625776648 + ,-3.40106461721486 + ,-0.472479772630227 + ,-1.65041714824595 + ,2.05175685598916 + ,-0.902464685057427 + ,-0.602947860990191 + ,0.818516482836542 + ,-1.59210559699142 + ,-1.76816394085798 + ,1.93284775246144 + ,0.227267023304705 + ,-2.88157270223229 + ,2.21075106222068 + ,0.674128651962353 + ,1.99963127057584 + ,0.0183623571902152 + ,-0.0312894695077324 + ,-0.39966157084491 + ,0.72139904577073 + ,-1.96806757395163 + ,3.21864918928638 + ,-1.20224319686796 + ,-1.32275450875619 + ,-0.268255724957080 + ,0.560341123135318 + ,2.62120611624231 + ,2.04599601516761 + ,1.49983429778472 + ,1.85337545161847 + ,2.29366444684220 + ,0.116117094486825 + ,-0.579730364351328 + ,0.436696611402682 + ,-1.02200667391198 + ,-0.726229843542471 + ,-1.30420266077865 + ,-1.07800520610845 + ,0.604162662833211) > 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/1s75s1261473005.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/2cfdn1261473005.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/39b5d1261473005.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/4h01m1261473005.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/5wd2r1261473005.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/64uvg1261473005.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/7soz21261473005.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.77200 -0.64620 0.11910 0.03784 0.66120 3.54800 > > #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/8xmtx1261473005.tab") > > try(system("convert tmp/1s75s1261473005.ps tmp/1s75s1261473005.png",intern=TRUE)) character(0) > try(system("convert tmp/2cfdn1261473005.ps tmp/2cfdn1261473005.png",intern=TRUE)) character(0) > try(system("convert tmp/39b5d1261473005.ps tmp/39b5d1261473005.png",intern=TRUE)) character(0) > try(system("convert tmp/4h01m1261473005.ps tmp/4h01m1261473005.png",intern=TRUE)) character(0) > try(system("convert tmp/5wd2r1261473005.ps tmp/5wd2r1261473005.png",intern=TRUE)) character(0) > try(system("convert tmp/64uvg1261473005.ps tmp/64uvg1261473005.png",intern=TRUE)) character(0) > try(system("convert tmp/7soz21261473005.ps tmp/7soz21261473005.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.577 1.091 2.599