R version 2.11.1 (2010-05-31) Copyright (C) 2010 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. > y <- c(563668,586111,604378,600991,544686,537034,551531,563250,574761,580112,575093,557560,564478,580523,596594,586570,536214,523597,536535,536322,532638,528222,516141,501866,506174,517945,533590,528379,477580,469357,490243,492622,507561,516922,514258,509846,527070,541657,564591,555362,498662,511038,525919,531673,548854,560576,557274,565742,587625,619916) > x <- c(276444,289742,303725,298305,266795,259497,266148,271037,276239,279681,277509,271115,275902,287224,300713,293860,264221,256167,262572,263276,264291,263903,260376,255603,261076,270976,285257,280445,250741,243803,253158,255542,262522,268381,267153,266424,276427,286994,303598,296806,263290,264981,272566,276475,284678,291542,291413,295916,309119,327616) > par2 = '0' > 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: > par1 <- as.numeric(par1) > par2 <- as.numeric(par2) > x <- as.ts(x) > y <- as.ts(y) > mylm <- lm(y~x) > cbind(mylm$resid) [,1] 1 19834.4356 2 19781.6570 3 14394.0870 4 20175.9184 5 17175.3238 6 21869.1045 7 25114.8319 8 28563.2753 9 31274.2271 10 30802.5117 11 29457.8109 12 22741.3253 13 21561.3187 14 18453.2727 15 11705.3859 16 13274.3751 17 13057.6728 18 14065.3533 19 16168.2306 20 14764.2975 21 9363.2562 22 5603.6227 23 -510.8702 24 -6711.5469 25 -11662.0366 26 -16638.5331 27 -25152.2196 28 -22222.9213 29 -22772.6653 30 -19258.8845 31 -14198.4228 32 -15852.3553 33 -12721.1861 34 -13271.6590 35 -13858.2928 36 -17037.0680 37 -16734.8061 38 -20023.6440 39 -25178.0714 40 -22917.2738 41 -22919.3858 42 -13403.9935 43 -11354.2825 44 -12213.0061 45 -8908.7448 46 -8798.3424 47 -11882.1174 48 -11031.6908 49 -11483.7611 50 -10483.5128 > library(lattice) > postscript(file="/var/www/rcomp/tmp/1tl5r1290624917.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(y,type='l',main='Run Sequence Plot of Y[t]',xlab='time or index',ylab='value') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/2tl5r1290624917.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x,type='l',main='Run Sequence Plot of X[t]',xlab='time or index',ylab='value') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/3lcmu1290624917.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x,y,main='Scatter Plot',xlab='X[t]',ylab='Y[t]') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/4lcmu1290624917.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mylm$resid,type='l',main='Run Sequence Plot of e[t]',xlab='time or index',ylab='value') > grid() > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5lcmu1290624917.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mylm$resid,main='Histogram of e[t]') > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6lcmu1290624917.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > if (par1 > 0) + { + densityplot(~mylm$resid,col='black',main=paste('Density Plot of e[t] bw = ',par1),bw=par1) + } else { + densityplot(~mylm$resid,col='black',main='Density Plot of e[t]') + } > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7wmmx1290624917.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(mylm$resid,main='QQ plot of e[t]') > qqline(mylm$resid) > grid() > dev.off() null device 1 > if (par2 > 0) + { + postscript(file="/var/www/rcomp/tmp/87dl01290624917.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + acf(mylm$resid,lag.max=par2,main='Residual Autocorrelation Function') + grid() + dev.off() + } > summary(x) Min. 1st Qu. Median Mean 3rd Qu. Max. 243800 263400 271800 275800 287200 327600 > > #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Model: Y[t] = c + b X[t] + e[t]',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'c',1,TRUE) > a<-table.element(a,mylm$coeff[[1]]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'b',1,TRUE) > a<-table.element(a,mylm$coeff[[2]]) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/9im2l1290624917.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Descriptive Statistics about e[t]',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'# observations',header=TRUE) > a<-table.element(a,length(mylm$resid)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'minimum',header=TRUE) > a<-table.element(a,min(mylm$resid)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Q1',header=TRUE) > a<-table.element(a,quantile(mylm$resid,0.25)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'median',header=TRUE) > a<-table.element(a,median(mylm$resid)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'mean',header=TRUE) > a<-table.element(a,mean(mylm$resid)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Q3',header=TRUE) > a<-table.element(a,quantile(mylm$resid,0.75)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'maximum',header=TRUE) > a<-table.element(a,max(mylm$resid)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/10weib1290624917.tab") > > try(system("convert tmp/1tl5r1290624917.ps tmp/1tl5r1290624917.png",intern=TRUE)) character(0) > try(system("convert tmp/2tl5r1290624917.ps tmp/2tl5r1290624917.png",intern=TRUE)) character(0) > try(system("convert tmp/3lcmu1290624917.ps tmp/3lcmu1290624917.png",intern=TRUE)) character(0) > try(system("convert tmp/4lcmu1290624917.ps tmp/4lcmu1290624917.png",intern=TRUE)) character(0) > try(system("convert tmp/5lcmu1290624917.ps tmp/5lcmu1290624917.png",intern=TRUE)) character(0) > try(system("convert tmp/6lcmu1290624917.ps tmp/6lcmu1290624917.png",intern=TRUE)) character(0) > try(system("convert tmp/7wmmx1290624917.ps tmp/7wmmx1290624917.png",intern=TRUE)) character(0) > try(system("convert tmp/87dl01290624917.ps tmp/87dl01290624917.png",intern=TRUE)) convert: unable to open image `tmp/87dl01290624917.ps': No such file or directory @ magick/blob.c/OpenBlob/2418. convert: missing an image filename `tmp/87dl01290624917.png' @ wand/convert.c/ConvertImageCommand/2710. character(0) > > > proc.time() user system elapsed 2.240 1.440 3.674