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(61.10986116
+ ,54.12031648
+ ,54.95749446
+ ,51.99594378
+ ,56.24546553
+ ,64.77452146
+ ,54.33699821
+ ,53.72119053
+ ,48.83342072
+ ,55.93379725
+ ,64.65918024
+ ,61.60053636
+ ,63.05659263
+ ,67.15166067
+ ,60.02845855
+ ,69.52301201
+ ,65.01551552
+ ,60.87993666
+ ,70.72557348
+ ,73.09285805
+ ,73.29160052
+ ,72.66843666
+ ,68.83636849
+ ,76.67414394
+ ,70.41408315
+ ,77.9401298
+ ,73.20981757
+ ,77.41384989
+ ,79.4630442
+ ,75.42057899
+ ,82.68120243
+ ,94.03299383
+ ,84.0698862
+ ,85.67449188
+ ,81.89949935
+ ,92.9720738
+ ,88.85764153
+ ,80.49481348
+ ,91.71796544
+ ,82.26348598
+ ,82.94116291
+ ,90.99841989
+ ,91.31788491
+ ,96.95308859
+ ,90.04825345
+ ,100.9221214
+ ,99.89956639
+ ,93.01828481
+ ,95.97742386
+ ,103.5617168
+ ,110.7469177
+ ,108.3696381
+ ,107.0012411
+ ,110.7524482
+ ,98.56007743
+ ,100.5808127
+ ,101.9883776
+ ,101.8791562
+ ,108.9640668
+ ,104.5917605
+ ,119.5561823
+ ,119.3482511
+ ,105.3661433
+ ,121.4417904
+ ,117.747784
+ ,116.6637445
+ ,119.6190875
+ ,115.2413454
+ ,106.7405737
+ ,113.7508379
+ ,124.9325395
+ ,124.9165153
+ ,134.0975866
+ ,132.0328883
+ ,124.4484784
+ ,115.0077499
+ ,122.6887566
+ ,125.0399372
+ ,119.2917291
+ ,131.6979609
+ ,132.3570185
+ ,136.9068086
+ ,147.3595766
+ ,135.5406186
+ ,137.5528819
+ ,128.6595377
+ ,135.701906
+ ,128.428847
+ ,134.1136873
+ ,133.3455986
+ ,139.8282064
+ ,152.4304094
+ ,148.7181493
+ ,142.8992652
+ ,152.7563365
+ ,139.7011827
+ ,146.5597201
+ ,146.5053758
+ ,145.5923464
+ ,140.0328201
+ ,154.0235601
+ ,154.9190928
+ ,158.101761
+ ,154.2821531
+ ,163.6434212
+ ,163.7318741
+ ,161.4377475
+ ,156.7145855
+ ,153.561049
+ ,153.8774688
+ ,169.2235043
+ ,156.3795439
+ ,165.490775
+ ,154.0411596
+ ,154.9816657
+ ,169.586015
+ ,160.6091489
+ ,160.6353333
+ ,163.6734875
+ ,166.1373986
+ ,175.3861465
+ ,181.8040035
+ ,172.5480095
+ ,175.5379656
+ ,177.4873292
+ ,178.8761947
+ ,178.375176
+ ,176.7772914
+ ,181.8905564
+ ,180.3690141
+ ,180.2465823
+ ,189.4586351
+ ,190.2788307
+ ,169.6704531
+ ,186.9818197
+ ,179.4104796
+ ,182.1682369
+ ,188.2755288
+ ,176.4747358
+ ,189.2416564
+ ,192.4300356
+ ,203.944427
+ ,203.3403548
+ ,188.1282826
+ ,196.354948
+ ,197.9087125
+ ,201.9367353
+ ,193.445003
+ ,198.8122889
+ ,196.4723325
+ ,204.5398687
+ ,205.0704692
+ ,204.3924263
+ ,207.9987893
+ ,200.9225201
+ ,214.2768736
+ ,208.0597487
+ ,208.5093111
+ ,201.22446
+ ,203.3434561
+ ,212.5706014
+ ,210.5406441
+ ,210.887768
+ ,209.6129575
+ ,209.313722
+ ,215.7808384
+ ,226.4811578
+ ,220.9966323
+ ,205.9677281
+ ,210.4273103
+ ,216.9371659
+ ,226.5364144
+ ,226.6246165
+ ,227.5501626
+ ,230.5976031
+ ,218.4966592
+ ,225.660519
+ ,219.7215854
+ ,224.2551636
+ ,223.4610044
+ ,241.6074588
+ ,220.6620768
+ ,231.2874518
+ ,239.7293117
+ ,235.6812993
+ ,239.2962581
+ ,230.904901
+ ,236.5510765
+ ,248.5657307
+ ,226.3445746
+ ,242.5403835
+ ,242.0843806
+ ,245.795533
+ ,244.5513355
+ ,244.957443
+ ,242.9284176
+ ,246.9975039
+ ,243.5141125
+ ,242.3741637
+ ,249.0753006
+ ,252.3467369
+ ,248.3848335
+ ,256.0723248
+ ,256.1781746
+ ,250.3571524
+ ,258.7296195
+ ,256.3099685
+ ,256.0184129
+ ,245.31261
+ ,248.7418796
+ ,263.5634464
+ ,268.8057802
+ ,259.432041
+ ,280.6897114
+ ,263.5020284
+ ,262.910557
+ ,254.8534255
+ ,262.914486
+ ,262.4803197
+ ,265.5528724
+ ,268.3350219
+ ,274.2269561
+ ,285.0879902
+ ,287.4846281
+ ,276.9574892
+ ,270.2328838
+ ,273.343679
+ ,282.3265576
+ ,271.3683696
+ ,276.7731473
+ ,284.984373
+ ,285.2336093
+ ,286.4188091
+ ,283.9944259
+ ,285.9733102
+ ,284.4558291
+ ,281.7885504
+ ,287.3979095
+ ,290.8197966
+ ,281.4066935
+ ,293.8918414
+ ,295.4459854
+ ,296.4138982
+ ,299.309
+ ,299.1331287
+ ,297.5508534
+ ,299.870089
+ ,298.1574167
+ ,293.7846644
+ ,289.1280973
+ ,302.7972545
+ ,311.6582427
+ ,306.8750536
+ ,300.9725842
+ ,309.210378
+ ,314.0951678
+ ,302.0571523
+ ,307.1139977
+ ,306.9729935
+ ,307.7220354
+ ,306.3991609
+ ,310.2296239
+ ,327.1883072
+ ,323.0723765
+ ,319.3126762
+ ,320.7490193
+ ,316.7516682
+ ,314.156943
+ ,318.148116
+ ,313.1329884
+ ,324.4617348
+ ,323.299863
+ ,323.5194599
+ ,325.766307
+ ,315.8600699
+ ,314.4862665
+ ,318.7478783
+ ,330.3323151
+ ,324.4215348
+ ,319.5462879
+ ,333.9837493
+ ,329.2936641
+ ,332.9469371
+ ,332.2435291
+ ,331.1389228
+ ,329.5469119
+ ,340.7498356
+ ,335.7682144
+ ,336.8801173
+ ,333.1690483
+ ,343.3672947
+ ,339.7428617
+ ,347.8462662
+ ,338.0019601
+ ,339.0321535
+ ,349.7254835
+ ,339.9088441
+ ,341.8686896
+ ,348.8204737
+ ,351.4017628
+ ,354.3571353
+ ,354.7671573
+ ,354.6799216
+ ,346.3211461
+ ,352.8674728
+ ,353.7544304
+ ,352.3393992
+ ,360.1987502
+ ,354.4125537
+ ,353.6133866
+ ,366.9320672
+ ,364.0824985
+ ,363.1392851
+ ,363.203798
+ ,358.7129274
+ ,357.8204599
+ ,366.2509726
+ ,365.6025863
+ ,370.8290372
+ ,366.5541817
+ ,374.8778657
+ ,374.9878806
+ ,375.9899812
+ ,377.0300208
+ ,376.3060598
+ ,381.8212896
+ ,380.7436544
+ ,370.7164768
+ ,368.0023819
+ ,376.9856375
+ ,387.8136742
+ ,383.166558
+ ,386.1462067
+ ,388.3909748
+ ,383.3338225
+ ,379.9889392
+ ,383.332495
+ ,394.2459788
+ ,391.6334594
+ ,390.6626585
+ ,385.1222833
+ ,397.8269152
+ ,394.4597945
+ ,395.180818
+ ,393.6250717
+ ,396.190475
+ ,395.3262294
+ ,392.0564838
+ ,402.4440354
+ ,393.7467618
+ ,394.3378194
+ ,410.4003906
+ ,398.2528185
+ ,407.9129263
+ ,408.5478406
+ ,403.4111597
+ ,398.4169705
+ ,406.8649612
+ ,406.8042436
+ ,405.1187658
+ ,417.866504
+ ,416.6827187
+ ,412.3299343
+ ,409.3072168
+ ,413.2115411
+ ,416.0710076
+ ,408.932371
+ ,404.8629156
+ ,413.4135646
+ ,406.9434141
+ ,426.4176986
+ ,420.0718548
+ ,430.8056774
+ ,419.9182
+ ,431.2393615
+ ,429.5679042
+ ,420.0486899
+ ,428.0971546
+ ,430.9662678
+ ,425.1897311
+ ,434.9097032
+ ,443.4206608
+ ,435.6845397
+ ,445.0445079
+ ,433.9489262
+ ,423.4262595
+ ,432.0529621
+ ,446.8513244
+ ,432.7851253
+ ,433.2442466
+ ,447.2221488
+ ,446.0185703
+ ,456.7722076
+ ,441.2170262
+ ,444.4915125
+ ,439.5461461
+ ,443.1393853
+ ,444.2187665
+ ,443.1337644
+ ,447.4276883
+ ,453.5096647
+ ,458.1323521
+ ,453.2685384
+ ,461.58333
+ ,454.3458605
+ ,460.8233618
+ ,456.7880646
+ ,451.2040437
+ ,450.7924225
+ ,450.3340228)
> x <- c(1000
+ ,1000
+ ,1000
+ ,1000
+ ,1000
+ ,1000
+ ,1000
+ ,1000
+ ,1000
+ ,1000
+ ,1100
+ ,1100
+ ,1100
+ ,1100
+ ,1100
+ ,1100
+ ,1100
+ ,1100
+ ,1100
+ ,1100
+ ,1200
+ ,1200
+ ,1200
+ ,1200
+ ,1200
+ ,1200
+ ,1200
+ ,1200
+ ,1200
+ ,1200
+ ,1300
+ ,1300
+ ,1300
+ ,1300
+ ,1300
+ ,1300
+ ,1300
+ ,1300
+ ,1300
+ ,1300
+ ,1400
+ ,1400
+ ,1400
+ ,1400
+ ,1400
+ ,1400
+ ,1400
+ ,1400
+ ,1400
+ ,1400
+ ,1500
+ ,1500
+ ,1500
+ ,1500
+ ,1500
+ ,1500
+ ,1500
+ ,1500
+ ,1500
+ ,1500
+ ,1600
+ ,1600
+ ,1600
+ ,1600
+ ,1600
+ ,1600
+ ,1600
+ ,1600
+ ,1600
+ ,1600
+ ,1700
+ ,1700
+ ,1700
+ ,1700
+ ,1700
+ ,1700
+ ,1700
+ ,1700
+ ,1700
+ ,1700
+ ,1800
+ ,1800
+ ,1800
+ ,1800
+ ,1800
+ ,1800
+ ,1800
+ ,1800
+ ,1800
+ ,1800
+ ,1900
+ ,1900
+ ,1900
+ ,1900
+ ,1900
+ ,1900
+ ,1900
+ ,1900
+ ,1900
+ ,1900
+ ,2000
+ ,2000
+ ,2000
+ ,2000
+ ,2000
+ ,2000
+ ,2000
+ ,2000
+ ,2000
+ ,2000
+ ,2100
+ ,2100
+ ,2100
+ ,2100
+ ,2100
+ ,2100
+ ,2100
+ ,2100
+ ,2100
+ ,2100
+ ,2200
+ ,2200
+ ,2200
+ ,2200
+ ,2200
+ ,2200
+ ,2200
+ ,2200
+ ,2200
+ ,2200
+ ,2300
+ ,2300
+ ,2300
+ ,2300
+ ,2300
+ ,2300
+ ,2300
+ ,2300
+ ,2300
+ ,2300
+ ,2400
+ ,2400
+ ,2400
+ ,2400
+ ,2400
+ ,2400
+ ,2400
+ ,2400
+ ,2400
+ ,2400
+ ,2500
+ ,2500
+ ,2500
+ ,2500
+ ,2500
+ ,2500
+ ,2500
+ ,2500
+ ,2500
+ ,2500
+ ,2600
+ ,2600
+ ,2600
+ ,2600
+ ,2600
+ ,2600
+ ,2600
+ ,2600
+ ,2600
+ ,2600
+ ,2700
+ ,2700
+ ,2700
+ ,2700
+ ,2700
+ ,2700
+ ,2700
+ ,2700
+ ,2700
+ ,2700
+ ,2800
+ ,2800
+ ,2800
+ ,2800
+ ,2800
+ ,2800
+ ,2800
+ ,2800
+ ,2800
+ ,2800
+ ,2900
+ ,2900
+ ,2900
+ ,2900
+ ,2900
+ ,2900
+ ,2900
+ ,2900
+ ,2900
+ ,2900
+ ,3000
+ ,3000
+ ,3000
+ ,3000
+ ,3000
+ ,3000
+ ,3000
+ ,3000
+ ,3000
+ ,3000
+ ,3100
+ ,3100
+ ,3100
+ ,3100
+ ,3100
+ ,3100
+ ,3100
+ ,3100
+ ,3100
+ ,3100
+ ,3200
+ ,3200
+ ,3200
+ ,3200
+ ,3200
+ ,3200
+ ,3200
+ ,3200
+ ,3200
+ ,3200
+ ,3300
+ ,3300
+ ,3300
+ ,3300
+ ,3300
+ ,3300
+ ,3300
+ ,3300
+ ,3300
+ ,3300
+ ,3400
+ ,3400
+ ,3400
+ ,3400
+ ,3400
+ ,3400
+ ,3400
+ ,3400
+ ,3400
+ ,3400
+ ,3500
+ ,3500
+ ,3500
+ ,3500
+ ,3500
+ ,3500
+ ,3500
+ ,3500
+ ,3500
+ ,3500
+ ,3600
+ ,3600
+ ,3600
+ ,3600
+ ,3600
+ ,3600
+ ,3600
+ ,3600
+ ,3600
+ ,3600
+ ,3700
+ ,3700
+ ,3700
+ ,3700
+ ,3700
+ ,3700
+ ,3700
+ ,3700
+ ,3700
+ ,3700
+ ,3800
+ ,3800
+ ,3800
+ ,3800
+ ,3800
+ ,3800
+ ,3800
+ ,3800
+ ,3800
+ ,3800
+ ,3900
+ ,3900
+ ,3900
+ ,3900
+ ,3900
+ ,3900
+ ,3900
+ ,3900
+ ,3900
+ ,3900
+ ,4000
+ ,4000
+ ,4000
+ ,4000
+ ,4000
+ ,4000
+ ,4000
+ ,4000
+ ,4000
+ ,4000
+ ,4100
+ ,4100
+ ,4100
+ ,4100
+ ,4100
+ ,4100
+ ,4100
+ ,4100
+ ,4100
+ ,4100
+ ,4200
+ ,4200
+ ,4200
+ ,4200
+ ,4200
+ ,4200
+ ,4200
+ ,4200
+ ,4200
+ ,4200
+ ,4300
+ ,4300
+ ,4300
+ ,4300
+ ,4300
+ ,4300
+ ,4300
+ ,4300
+ ,4300
+ ,4300
+ ,4400
+ ,4400
+ ,4400
+ ,4400
+ ,4400
+ ,4400
+ ,4400
+ ,4400
+ ,4400
+ ,4400
+ ,4500
+ ,4500
+ ,4500
+ ,4500
+ ,4500
+ ,4500
+ ,4500
+ ,4500
+ ,4500
+ ,4500
+ ,4600
+ ,4600
+ ,4600
+ ,4600
+ ,4600
+ ,4600
+ ,4600
+ ,4600
+ ,4600
+ ,4600
+ ,4700
+ ,4700
+ ,4700
+ ,4700
+ ,4700
+ ,4700
+ ,4700
+ ,4700
+ ,4700
+ ,4700
+ ,4800
+ ,4800
+ ,4800
+ ,4800
+ ,4800
+ ,4800
+ ,4800
+ ,4800
+ ,4800
+ ,4800
+ ,4900
+ ,4900
+ ,4900
+ ,4900
+ ,4900
+ ,4900
+ ,4900
+ ,4900
+ ,4900
+ ,4900
+ ,5000
+ ,5000
+ ,5000
+ ,5000
+ ,5000
+ ,5000
+ ,5000
+ ,5000
+ ,5000
+ ,5000)
> 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)
> library(lattice)
> z <- as.data.frame(cbind(x,y))
> m <- lm(y~x)
> summary(m)

Call:
lm(formula = y ~ x)

Residuals:
     Min       1Q   Median       3Q      Max 
-15.4885  -2.9879  -0.1656   3.0983  15.6577 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)    
(Intercept) -44.476036   0.664269  -66.95   <2e-16 ***
x             0.099841   0.000206  484.71   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 4.935 on 408 degrees of freedom
Multiple R-squared: 0.9983,	Adjusted R-squared: 0.9983 
F-statistic: 2.349e+05 on 1 and 408 DF,  p-value: < 2.2e-16 

> postscript(file="/var/www/rcomp/tmp/1sr001289938742.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(z,main='Scatterplot, lowess, and regression line')
> lines(lowess(z),col='red')
> abline(m)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/2sr001289938742.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> m2 <- lm(m$fitted.values ~ x)
> summary(m2)

Call:
lm(formula = m$fitted.values ~ x)

Residuals:
       Min         1Q     Median         3Q        Max 
-3.852e-14 -6.174e-15  2.424e-15  8.223e-15  1.103e-13 

Coefficients:
              Estimate Std. Error    t value Pr(>|t|)    
(Intercept) -4.448e+01  1.694e-15 -2.626e+16   <2e-16 ***
x            9.984e-02  5.252e-19  1.901e+17   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 1.258e-14 on 408 degrees of freedom
Multiple R-squared:     1,	Adjusted R-squared:     1 
F-statistic: 3.614e+34 on 1 and 408 DF,  p-value: < 2.2e-16 

> z2 <- as.data.frame(cbind(x,m$fitted.values))
> names(z2) <- list('x','Fitted')
> plot(z2,main='Scatterplot, lowess, and regression line')
> lines(lowess(z2),col='red')
> abline(m2)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/3k0il1289938742.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> m3 <- lm(m$residuals ~ x)
> summary(m3)

Call:
lm(formula = m$residuals ~ x)

Residuals:
     Min       1Q   Median       3Q      Max 
-15.4885  -2.9879  -0.1656   3.0983  15.6577 

Coefficients:
              Estimate Std. Error   t value Pr(>|t|)
(Intercept) -2.647e-17  6.643e-01 -3.99e-17        1
x            2.064e-20  2.060e-04  1.00e-16        1

Residual standard error: 4.935 on 408 degrees of freedom
Multiple R-squared: 1.887e-34,	Adjusted R-squared: -0.002451 
F-statistic: 7.7e-32 on 1 and 408 DF,  p-value: 1 

> z3 <- as.data.frame(cbind(x,m$residuals))
> names(z3) <- list('x','Residuals')
> plot(z3,main='Scatterplot, lowess, and regression line')
> lines(lowess(z3),col='red')
> abline(m3)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/4v9zo1289938742.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> m4 <- lm(m$fitted.values ~ m$residuals)
> summary(m4)

Call:
lm(formula = m$fitted.values ~ m$residuals)

Residuals:
       Min         1Q     Median         3Q        Max 
-1.997e+02 -9.984e+01  8.857e-15  9.984e+01  1.997e+02 

Coefficients:
              Estimate Std. Error   t value Pr(>|t|)    
(Intercept)  2.550e+02  5.848e+00     43.61   <2e-16 ***
m$residuals -9.606e-17  1.188e+00 -8.09e-17        1    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 118.4 on 408 degrees of freedom
Multiple R-squared: 5.125e-32,	Adjusted R-squared: -0.002451 
F-statistic: 2.091e-29 on 1 and 408 DF,  p-value: 1 

> z4 <- as.data.frame(cbind(m$residuals,m$fitted.values))
> names(z4) <- list('Residuals','Fitted')
> plot(z4,main='Scatterplot, lowess, and regression line')
> lines(lowess(z4),col='red')
> abline(m4)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5v9zo1289938742.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> myr <- as.ts(m$residuals)
> z5 <- as.data.frame(cbind(lag(myr,1),myr))
> names(z5) <- list('Lagged Residuals','Residuals')
> plot(z5,main='Lag plot')
> m5 <- lm(z5)
> summary(m5)

Call:
lm(formula = z5)

Residuals:
     Min       1Q   Median       3Q      Max 
-15.6126  -3.0411  -0.1319   3.0951  15.8234 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.01434    0.24382  -0.059    0.953
Residuals    0.02703    0.04952   0.546    0.585

Residual standard error: 4.931 on 407 degrees of freedom
  (2 observations deleted due to missingness)
Multiple R-squared: 0.0007316,	Adjusted R-squared: -0.001724 
F-statistic: 0.298 on 1 and 407 DF,  p-value: 0.5855 

> abline(m5)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6gajm1289938743.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> hist(m$residuals,main='Residual Histogram',xlab='Residuals')
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/7gajm1289938743.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> if (par1 > 0)
+ {
+ densityplot(~m$residuals,col='black',main=paste('Density Plot   bw = ',par1),bw=par1)
+ } else {
+ densityplot(~m$residuals,col='black',main='Density Plot')
+ }
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/8r1io1289938743.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> acf(m$residuals,main='Residual Autocorrelation Function')
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/9r1io1289938743.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> qqnorm(x)
> qqline(x)
> grid()
> dev.off()
null device 
          1 
> 
> #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,'Simple Linear Regression',5,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Statistics',1,TRUE)
> a<-table.element(a,'Estimate',1,TRUE)
> a<-table.element(a,'S.D.',1,TRUE)
> a<-table.element(a,'T-STAT (H0: coeff=0)',1,TRUE)
> a<-table.element(a,'P-value (two-sided)',1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'constant term',header=TRUE)
> a<-table.element(a,m$coefficients[[1]])
> sd <- sqrt(vcov(m)[1,1])
> a<-table.element(a,sd)
> tstat <- m$coefficients[[1]]/sd
> a<-table.element(a,tstat)
> pval <- 2*(1-pt(abs(tstat),length(x)-2))
> a<-table.element(a,pval)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'slope',header=TRUE)
> a<-table.element(a,m$coefficients[[2]])
> sd <- sqrt(vcov(m)[2,2])
> a<-table.element(a,sd)
> tstat <- m$coefficients[[2]]/sd
> a<-table.element(a,tstat)
> pval <- 2*(1-pt(abs(tstat),length(x)-2))
> a<-table.element(a,pval)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/10g3fi1289938743.tab") 
> 
> try(system("convert tmp/1sr001289938742.ps tmp/1sr001289938742.png",intern=TRUE))
character(0)
> try(system("convert tmp/2sr001289938742.ps tmp/2sr001289938742.png",intern=TRUE))
character(0)
> try(system("convert tmp/3k0il1289938742.ps tmp/3k0il1289938742.png",intern=TRUE))
character(0)
> try(system("convert tmp/4v9zo1289938742.ps tmp/4v9zo1289938742.png",intern=TRUE))
character(0)
> try(system("convert tmp/5v9zo1289938742.ps tmp/5v9zo1289938742.png",intern=TRUE))
character(0)
> try(system("convert tmp/6gajm1289938743.ps tmp/6gajm1289938743.png",intern=TRUE))
character(0)
> try(system("convert tmp/7gajm1289938743.ps tmp/7gajm1289938743.png",intern=TRUE))
character(0)
> try(system("convert tmp/8r1io1289938743.ps tmp/8r1io1289938743.png",intern=TRUE))
character(0)
> try(system("convert tmp/9r1io1289938743.ps tmp/9r1io1289938743.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.260   1.860   5.209