R version 2.6.1 (2007-11-26) Copyright (C) 2007 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(-64.7052505627806 + ,-3456.92146791498 + ,1221.53818880044 + ,-3425.92181819955 + ,476.517449571596 + ,771.237252580258 + ,1486.03059810756 + ,-1420.53604746731 + ,-1231.65983588785 + ,-1072.35237321869 + ,-1028.79999310919 + ,1689.70160953415 + ,104.781462266857 + ,3056.00468005053 + ,2058.12890000085 + ,-2120.15571952206 + ,857.951846139995 + ,-7653.62307012801 + ,-3625.47032737070 + ,-1066.78246083915 + ,1480.43801029269 + ,383.807285529813 + ,-183.460775788262 + ,-412.66366558563 + ,-3386.49397717813 + ,-868.872254604358 + ,-4519.64765815746 + ,3043.34922909291 + ,4526.15618595789 + ,2298.41482668304 + ,1992.23133104446 + ,1442.95792211978 + ,-893.403718302956 + ,3167.99873285149 + ,231.403932376355 + ,-534.568799433291 + ,-3856.42358359645 + ,-418.841481409176 + ,-1376.94550043225 + ,1524.94083007199 + ,-6240.6385168242 + ,512.555201154925 + ,1722.16374385519 + ,756.960403479329 + ,2999.69977354755 + ,-2634.01843507769 + ,-1578.85998252178 + ,78.6364326059933 + ,-376.690663222133 + ,992.978747098893 + ,994.51938653588 + ,-10.5193017588915 + ,-921.573813091965 + ,606.066588722763 + ,1003.79725596858 + ,1201.6546911613 + ,1805.44593151825 + ,596.116762365994 + ,-1360.77445520603 + ,3383.24617420951 + ,-5515.05833999521 + ,-2431.17906069741 + ,4733.10028012447 + ,4010.28330224708 + ,-2898.76736676388 + ,76.4283006144755 + ,-2221.26939354851 + ,2373.68897140019 + ,-1244.72571223942 + ,537.169792121625 + ,2246.81073824300 + ,864.588365965909 + ,611.893000896821 + ,-1922.54556007542 + ,-2682.31131104567 + ,86.6180766216501 + ,2273.08966203739 + ,3100.18749454311 + ,-3642.00381742456 + ,-661.164658428561 + ,-728.23512638401 + ,-3035.53289896783 + ,585.916426040515 + ,-2542.63265568714) > par2 = '0' > par1 = '0' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: Wessa P., (2007), Univariate Explorative Data Analysis (v1.0.5) 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 > #Technical description: Write here your technical program description > par1 <- as.numeric(par1) > par2 <- as.numeric(par2) > x <- as.ts(x) > library(lattice) > postscript(file="/var/www/html/rcomp/tmp/10ukn1197551978.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/2m8jy1197551978.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/31ob11197551979.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/4px9l1197551979.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(x) > grid() > dev.off() null device 1 > if (par2 > 0) + { + postscript(file="/var/www/html/rcomp/tmp/57fn31197551979.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=paste('Lag plot, lowess, and regression line')) + lines(lowess(z)) + abline(lm(z)) + dev.off() + postscript(file="/var/www/html/rcomp/tmp/6rq2t1197551979.ps",horizontal=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. -7654.00 -1388.00 82.63 -182.10 1452.00 4733.00 > 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/7pdrw1197551979.tab") > > system("convert tmp/10ukn1197551978.ps tmp/10ukn1197551978.png") > system("convert tmp/2m8jy1197551978.ps tmp/2m8jy1197551978.png") > system("convert tmp/31ob11197551979.ps tmp/31ob11197551979.png") > system("convert tmp/4px9l1197551979.ps tmp/4px9l1197551979.png") > system("convert tmp/57fn31197551979.ps tmp/57fn31197551979.png") convert: unable to open image `tmp/57fn31197551979.ps': No such file or directory. convert: missing an image filename `tmp/57fn31197551979.png'. > system("convert tmp/6rq2t1197551979.ps tmp/6rq2t1197551979.png") convert: unable to open image `tmp/6rq2t1197551979.ps': No such file or directory. convert: missing an image filename `tmp/6rq2t1197551979.png'. > > > proc.time() user system elapsed 3.529 1.766 3.792