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(-2.02969015 + ,-2.029239482 + ,-2.033801546 + ,-2.036644359 + ,-2.034490446 + ,-2.032001271 + ,-2.030699746 + ,-2.03252381 + ,-2.038987342 + ,-2.040993691 + ,-2.037075472 + ,-2.03754865 + ,-2.04075 + ,-2.040945274 + ,-2.038261272 + ,-2.037953846 + ,-2.0379748 + ,-2.038783784 + ,-2.039284837 + ,-2.042300613 + ,-2.040853659 + ,-2.041885246 + ,-2.041468447 + ,-2.040445725 + ,-2.037787879 + ,-2.038293124 + ,-2.037360024 + ,-2.037595694 + ,-2.038959451 + ,-2.040858586 + ,-2.039697509 + ,-2.036113085 + ,-2.038323477 + ,-2.038952802 + ,-2.040082547 + ,-2.03908957 + ,-2.040747939 + ,-2.040328638 + ,-2.041483324 + ,-2.040980221 + ,-2.040107621 + ,-2.038953488 + ,-2.039244186 + ,-2.040846868 + ,-2.040780347 + ,-2.042795049 + ,-2.043550058 + ,-2.044298044 + ,-2.044189189 + ,-2.044576174 + ,-2.045057143 + ,-2.045278567 + ,-2.047159091 + ,-2.047747002 + ,-2.048223235 + ,-2.047723138 + ,-2.048265422 + ,-2.048072118 + ,-2.04665247 + ,-2.046713961 + ,-2.043882587 + ,-2.044426322 + ,-2.043627451 + ,-2.043628763 + ,-2.044415042 + ,-2.044977728 + ,-2.044315526 + ,-2.044315526 + ,-2.044043382 + ,-2.041325737 + ,-2.041903441 + ,-2.04097561 + ,-2.042195122 + ,-2.044078875 + ,-2.046973612 + ,-2.04809589 + ,-2.046351943 + ,-2.045496996 + ,-2.044705882 + ,-2.044929232 + ,-2.046186072 + ,-2.047930941 + ,-2.047286801 + ,-2.048763557 + ,-2.052723577 + ,-2.051019365 + ,-2.049533762 + ,-2.046950959 + ,-2.045095847 + ,-2.046538052 + ,-2.049521277 + ,-2.052663303 + ,-2.050647558 + ,-2.050437666 + ,-2.051923077 + ,-2.051871686 + ,-2.049962884 + ,-2.049825766 + ,-2.048526316 + ,-2.048903075 + ,-2.050636364 + ,-2.051165803 + ,-2.051640827 + ,-2.050306859 + ,-2.048872297 + ,-2.048329913 + ,-2.04714982 + ,-2.045046178 + ,-2.045815803 + ,-2.045025575 + ,-2.043469388 + ,-2.044602446 + ,-2.043990826 + ,-2.04098369 + ,-2.040779817 + ,-2.040380711 + ,-2.038004552 + ,-2.04 + ,-2.039359233 + ,-2.042455919 + ,-2.043922457 + ,-2.043721805 + ,-2.04455 + ,-2.041140089 + ,-2.039470821 + ,-2.040559389 + ,-2.041977887 + ,-2.041879607 + ,-2.043921569 + ,-2.043486036 + ,-2.044574156 + ,-2.044339853 + ,-2.040839844 + ,-2.042176987 + ,-2.042019932 + ,-2.046234867 + ,-2.048571774 + ,-2.050750845 + ,-2.048840231 + ,-2.048942308 + ,-2.048956835 + ,-2.048815916 + ,-2.049736842 + ,-2.048448111 + ,-2.049229852 + ,-2.046706387 + ,-2.046615055 + ,-2.047952381 + ,-2.050066414 + ,-2.047336802 + ,-2.047961738 + ,-2.046556604 + ,-2.046517436 + ,-2.047561205 + ,-2.046263504 + ,-2.046846298 + ,-2.047559663 + ,-2.04911215 + ,-2.049727018 + ,-2.046061453 + ,-2.042860465 + ,-2.044532191 + ,-2.043863216 + ,-2.044580452 + ,-2.045653775 + ,-2.041850575 + ,-2.043297384 + ,-2.04459193 + ,-2.044661482 + ,-2.043259479 + ,-2.041000912 + ,-2.042477314 + ,-2.043637392 + ,-2.042013575 + ,-2.040843102 + ,-2.040990991 + ,-2.040778578 + ,-2.043898876 + ,-2.043470166 + ,-2.042046168 + ,-2.043227395 + ,-2.044375 + ,-2.047492226 + ,-2.046308204 + ,-2.045305445 + ,-2.046736633 + ,-2.047935982 + ,-2.047819383 + ,-2.049632644 + ,-2.049688323 + ,-2.048467777 + ,-2.047819615 + ,-2.04773523 + ,-2.047263986 + ,-2.047052402 + ,-2.047230467 + ,-2.049417793 + ,-2.049381533 + ,-2.049040887 + ,-2.050130435 + ,-2.048440487 + ,-2.045294627 + ,-2.042272727 + ,-2.041663789 + ,-2.040079776 + ,-2.039377682 + ,-2.036327095 + ,-2.037415254 + ,-2.033007598 + ,-2.030537748 + ,-2.030115789 + ,-2.032371956 + ,-2.032145249 + ,-2.035259306 + ,-2.041916667 + ,-2.044050791 + ,-2.04535343 + ,-2.04698172 + ,-2.052428453 + ,-2.052514499 + ,-2.056527663 + ,-2.059528613 + ,-2.061907895 + ,-2.058198032) > 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/10ohd1257158602.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/2obvl1257158602.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/3s81y1257158602.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/4qpay1257158602.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/5p7e91257158602.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/6zrgn1257158602.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/7sy4w1257158602.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. -2.062 -2.048 -2.044 -2.044 -2.041 -2.029 > > #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/8uwgt1257158602.tab") > > system("convert tmp/10ohd1257158602.ps tmp/10ohd1257158602.png") > system("convert tmp/2obvl1257158602.ps tmp/2obvl1257158602.png") > system("convert tmp/3s81y1257158602.ps tmp/3s81y1257158602.png") > system("convert tmp/4qpay1257158602.ps tmp/4qpay1257158602.png") > system("convert tmp/5p7e91257158602.ps tmp/5p7e91257158602.png") > system("convert tmp/6zrgn1257158602.ps tmp/6zrgn1257158602.png") > system("convert tmp/7sy4w1257158602.ps tmp/7sy4w1257158602.png") > > > proc.time() user system elapsed 1.539 1.082 1.935