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.

> y <- c(3.2,1.9,0,0.6,0.2,0.9,2.4,4.7,9.4,12.5,15.8,18.2,16.8,17.3,19.3,17.9,20.2,18.7,20.1,18.2,18.4,18.2,18.9,19.9,21.3,20,19.5,19.6,20.9,21,19.9,19.6,20.9,21.7,22.9,21.5,21.3,23.5,21.6,24.5,22.2,23.5,20.9,20.7,18.1,17.1,14.8,13.8,15.2,16,17.6,15,15,16.3,19.4,21.3,20.5,21.1,21.6,22.6)
> x <- c(1.74,1.75,1.83,2.09,2.12,2.29,2.4,2.82,3.18,4,4.8,5.28,5.37,5.27,5.33,5.23,5.08,5.11,5.1,4.97,5,5.2,4.9,4.82,5.04,4.82,4.77,4.79,4.58,4.59,4.57,4.35,4.27,4.39,3.97,3.84,3.73,3.58,3.45,3.44,3.25,3.25,3.02,2.87,2.92,2.95,2.75,2.7,2.75,2.72,2.71,2.76,2.68,2.78,2.86,2.75,2.87,2.91,2.79,2.77)
> 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: 
> 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   -7.4076630
2   -8.7393864
3  -10.8931730
4  -11.1179795
5  -11.6131495
6  -11.4524461
7  -10.3014027
8   -9.3337825
9   -5.7758223
10  -5.2771352
11  -4.5150014
12  -3.6377211
13  -5.3232311
14  -4.5059978
15  -2.6963378
16  -3.7791045
17  -1.0032546
18  -2.5984246
19  -1.1667012
20  -2.6542980
21  -2.5494680
22  -3.3839345
23  -1.7322347
24  -0.4784481
25   0.2236387
26  -0.3784481
27  -0.7198314
28  -0.6832781
29   1.2829118
30   1.3511885
31   0.3146351
32   0.7125484
33   2.2663350
34   2.6856550
35   5.2180348
36   4.2304381
37   4.3793947
38   7.0552446
39   5.5676479
40   8.4993712
41   6.8021144
42   8.1021144
43   6.2317510
44   6.5076009
45   3.7489843
46   2.6538143
47   0.9882808
48   0.1468975
49   1.3882808
50   2.2834508
51   3.9151741
52   1.1565575
53   1.4103441
54   2.3931108
55   5.2393242
56   7.4882808
57   6.3076009
58   6.7807076
59   7.6613875
60   8.7248342
> library(lattice)
> postscript(file="/var/www/html/rcomp/tmp/1cb7o1258131726.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/html/rcomp/tmp/2jrhz1258131726.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/html/rcomp/tmp/340pu1258131726.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/html/rcomp/tmp/4xf5j1258131726.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/html/rcomp/tmp/54r131258131726.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/html/rcomp/tmp/69fcu1258131726.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/html/rcomp/tmp/7egzi1258131726.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/html/rcomp/tmp/8oa1l1258131726.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()
+ }
null device 
          1 
> summary(x)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  1.740   2.768   3.445   3.682   4.805   5.370 
> 
> #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,'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/html/rcomp/tmp/9zolo1258131726.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/html/rcomp/tmp/106cwk1258131726.tab") 
> system("convert tmp/1cb7o1258131726.ps tmp/1cb7o1258131726.png")
> system("convert tmp/2jrhz1258131726.ps tmp/2jrhz1258131726.png")
> system("convert tmp/340pu1258131726.ps tmp/340pu1258131726.png")
> system("convert tmp/4xf5j1258131726.ps tmp/4xf5j1258131726.png")
> system("convert tmp/54r131258131726.ps tmp/54r131258131726.png")
> system("convert tmp/69fcu1258131726.ps tmp/69fcu1258131726.png")
> system("convert tmp/7egzi1258131726.ps tmp/7egzi1258131726.png")
> system("convert tmp/8oa1l1258131726.ps tmp/8oa1l1258131726.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.627   1.189   1.918