R version 2.6.0 (2007-10-03) 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(0.0640061800725882 + ,0.055543390675083 + ,0.032601938563982 + ,0.0408184280409284 + ,0.0462166153670518 + ,-0.0137516738094333 + ,-0.00063153303605741 + ,-0.00233761316721433 + ,0.0483117836462727 + ,0.0729374241986308 + ,0.0573703673438476 + ,0.0732313798761501 + ,0.0182045786194426 + ,0.0203594752119119 + ,0.0227844996034191 + ,-0.0378014755630382 + ,-0.0646039043977659 + ,-0.0373526591017259 + ,-0.0233447045393399 + ,-0.0482940754041509 + ,0.00826205351002031 + ,-0.0234440266211367 + ,-0.0252275371441901 + ,-0.0430726405517209 + ,-0.125589104079991 + ,-0.112546393698512 + ,-0.099426252925136 + ,0.000565864129830549 + ,0.0438787023103477 + ,0.0542958079481452 + ,0.0126671582088848 + ,0.0115787640677000 + ,0.0285598815647022 + ,0.0626452665189847 + ,0.0379660234532115 + ,0.00888554806573639 + ,0.0313105724572436 + ,0.00713374836619783 + ,0.0321843775013869 + ,-0.00380766793622727 + ,-0.0143558064793364 + ,-0.0633588514749149 + ,-0.0538295596261275 + ,-0.0658057675563223 + ,-0.0597898942402267 + ,-0.0415734047632801 + ,-0.0616191243316618 + ,-0.0630550766637809 + ,-0.0515952964531803 + ,-0.0327611209862616 + ,0.00676817086252592 + ,-0.0325050019320905 + ,-0.0721257329176413 + ,-0.0613610690889093 + ,-0.0588586143055058 + ,0.00260116590509487 + ,0.0317828995629480 + ,-0.0113512499552947 + ,-0.00444756285018907 + ,0.00956039171219682 + ,-0.0324932466098945 + ,-0.0454346987209958 + ,-0.0467157902693218 + ,0.0301878968357837 + ,0.0413001188554499 + ,0.0635324458398158 + ,0.0644519704523407 + ,0.0643288204920626 + ,0.0764837170845318 + ,0.150066683064087 + ,0.0873953587520234 + ,0.0134886624600154) > par2 = '36' > 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/1uc211200413553.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/2v1em1200413553.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/3hwap1200413553.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/4j2c11200413553.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/5dohc1200413553.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/6vwjh1200413553.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. -1.256e-01 -4.195e-02 4.685e-03 9.287e-18 3.868e-02 1.501e-01 > 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/73jz01200413553.tab") > > system("convert tmp/1uc211200413553.ps tmp/1uc211200413553.png") > system("convert tmp/2v1em1200413553.ps tmp/2v1em1200413553.png") > system("convert tmp/3hwap1200413553.ps tmp/3hwap1200413553.png") > system("convert tmp/4j2c11200413553.ps tmp/4j2c11200413553.png") > system("convert tmp/5dohc1200413553.ps tmp/5dohc1200413553.png") > system("convert tmp/6vwjh1200413553.ps tmp/6vwjh1200413553.png") > > > proc.time() user system elapsed 1.577 0.882 1.785