R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

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(117
+ ,116
+ ,166
+ ,180
+ ,202
+ ,290
+ ,298
+ ,441
+ ,388
+ ,260
+ ,175
+ ,105
+ ,137
+ ,142
+ ,176
+ ,231
+ ,240
+ ,316
+ ,363
+ ,537
+ ,487
+ ,324
+ ,185
+ ,133
+ ,169
+ ,157
+ ,206
+ ,244
+ ,243
+ ,393
+ ,405
+ ,579
+ ,525
+ ,373
+ ,198
+ ,148
+ ,201
+ ,177
+ ,222
+ ,275
+ ,290
+ ,402
+ ,534
+ ,614
+ ,578
+ ,419
+ ,203
+ ,173
+ ,229
+ ,192
+ ,294
+ ,310
+ ,365
+ ,509
+ ,537
+ ,655
+ ,643
+ ,444
+ ,259
+ ,229
+ ,276
+ ,245
+ ,324
+ ,323
+ ,349
+ ,480
+ ,530
+ ,676
+ ,670
+ ,476
+ ,281
+ ,240
+ ,259
+ ,237
+ ,400
+ ,367
+ ,497
+ ,593
+ ,696
+ ,969
+ ,878
+ ,581
+ ,373
+ ,232
+ ,358
+ ,318
+ ,410
+ ,480
+ ,604
+ ,713
+ ,844
+ ,1134
+ ,1013
+ ,755
+ ,371
+ ,280
+ ,417
+ ,417
+ ,514
+ ,548
+ ,583
+ ,839
+ ,924
+ ,1179
+ ,1109
+ ,896
+ ,452
+ ,337
+ ,484
+ ,524
+ ,575
+ ,622
+ ,664
+ ,926
+ ,1028
+ ,1361
+ ,1304
+ ,937
+ ,505
+ ,427
+ ,580
+ ,483
+ ,625
+ ,695
+ ,729
+ ,1099
+ ,1090
+ ,1393
+ ,1261
+ ,988
+ ,525
+ ,416
+ ,516
+ ,454
+ ,629
+ ,755
+ ,706
+ ,951
+ ,1099
+ ,1444
+ ,1316
+ ,1066
+ ,585
+ ,430
+ ,669
+ ,598
+ ,714
+ ,835
+ ,912
+ ,1031
+ ,1210
+ ,1581
+ ,1416
+ ,1120
+ ,652
+ ,505
+ ,741
+ ,675
+ ,782
+ ,956
+ ,996
+ ,1259
+ ,1389
+ ,1868
+ ,1609
+ ,1385
+ ,735
+ ,577
+ ,815
+ ,798
+ ,940
+ ,1007
+ ,1094
+ ,1413
+ ,1552
+ ,2038
+ ,1762
+ ,1411
+ ,805
+ ,729
+ ,912
+ ,753
+ ,989
+ ,1137
+ ,1256
+ ,1554
+ ,1629
+ ,2024
+ ,1900
+ ,1563
+ ,905
+ ,766
+ ,952
+ ,915
+ ,1197
+ ,1242
+ ,1197
+ ,1522
+ ,1591
+ ,2128
+ ,1962
+ ,1653
+ ,987
+ ,877
+ ,990
+ ,880
+ ,1258
+ ,1240
+ ,1312
+ ,1713
+ ,1683
+ ,2220
+ ,1996
+ ,1628
+ ,1119
+ ,890
+ ,1118
+ ,1164
+ ,1364
+ ,1412
+ ,1721
+ ,1752
+ ,1794
+ ,2434
+ ,2390
+ ,1929
+ ,1352
+ ,1060
+ ,1435
+ ,1196
+ ,1478
+ ,1648
+ ,1812
+ ,2118
+ ,2211
+ ,2826
+ ,2534
+ ,2290
+ ,1367
+ ,1105
+ ,1463
+ ,1299
+ ,1576
+ ,1850
+ ,1929
+ ,2367
+ ,2508
+ ,3073
+ ,2922
+ ,2377
+ ,1627
+ ,1259
+ ,1547
+ ,1436
+ ,1905
+ ,2079
+ ,1994
+ ,2501
+ ,2569
+ ,3467
+ ,2885
+ ,2211
+ ,1597
+ ,1141
+ ,1533
+ ,1546
+ ,1967
+ ,2171
+ ,2021
+ ,2753
+ ,2626
+ ,3532
+ ,3096
+ ,2639
+ ,1653
+ ,1425
+ ,1802
+ ,1674
+ ,1970
+ ,2092
+ ,2280
+ ,2715
+ ,2971
+ ,3937
+ ,3110
+ ,2662
+ ,1728
+ ,1609
+ ,1922
+ ,1863
+ ,1945
+ ,2365
+ ,2275
+ ,2962
+ ,2930
+ ,4062
+ ,3445
+ ,2943
+ ,1879
+ ,1694
+ ,2147
+ ,1999
+ ,2266
+ ,2562
+ ,2583
+ ,2965
+ ,3142
+ ,4115
+ ,3654
+ ,2992
+ ,2031
+ ,1699
+ ,2313
+ ,1970
+ ,2382
+ ,2830
+ ,2614
+ ,3321
+ ,3418
+ ,4468
+ ,3657
+ ,3250
+ ,2174
+ ,2014
+ ,2118
+ ,2227
+ ,2563
+ ,2817
+ ,2680
+ ,3337
+ ,3559
+ ,4608
+ ,3930
+ ,3133
+ ,2042
+ ,1999
+ ,2679
+ ,2425
+ ,2693
+ ,2760
+ ,2941
+ ,3611
+ ,3779
+ ,4945
+ ,4034
+ ,2906
+ ,2132
+ ,1932
+ ,2268
+ ,2178
+ ,2317
+ ,2552
+ ,2582
+ ,2886
+ ,3283
+ ,4125
+ ,3536
+ ,2568
+ ,1802
+ ,1598
+ ,2013
+ ,1872
+ ,2227
+ ,2497
+ ,2530
+ ,3119
+ ,3411
+ ,4511
+ ,3528
+ ,2833
+ ,1760
+ ,1517
+ ,1968
+ ,1809
+ ,2104
+ ,2391
+ ,2691
+ ,3023
+ ,3188
+ ,4057
+ ,3476)
> par1 = '12'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Dr. Ian E. Holliday
> #To cite this work: Ian E. Holliday, 2009, 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: 
> #Technical description: 
> par1 <- as.numeric(par1)
> n <- length(x)
> sx <- sort(x)
> 
> #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,'Variance Reduction Matrix',6,TRUE)
> a<-table.row.end(a)
> for (bigd in 0:2) {
+ for (smalld in 0:3) {
+ mylabel <- 'V(Y[t],d='
+ mylabel <- paste(mylabel,as.character(smalld),sep='')
+ mylabel <- paste(mylabel,',D=',sep='')
+ mylabel <- paste(mylabel,as.character(bigd),sep='')
+ mylabel <- paste(mylabel,')',sep='')
+ a<-table.row.start(a)
+ a<-table.element(a,mylabel,header=TRUE)
+ myx <- x
+ if (smalld > 0) myx <- diff(myx,lag=1,differences=smalld)
+ if (bigd > 0) myx <- diff(myx,lag=par1,differences=bigd)
+ a<-table.element(a,var(myx))
+ a<-table.element(a,'Range',header=TRUE)
+ a<-table.element(a,max(myx)-min(myx))
+ a<-table.element(a,'Trim Var.',header=TRUE)
+ smyx <- sort(myx)
+ sn <- length(smyx)
+ a<-table.element(a,var(smyx[smyx>quantile(smyx,0.05) & smyx<quantile(smyx,0.95)]))
+ a<-table.row.end(a)
+ }
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/1fdb01323790226.tab") 
> postscript(file="/var/www/rcomp/tmp/2ofc61323790226.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> op <- par(mfrow=c(2,2))
> plot(x,type='l',xlab='time',ylab='value',main='d=0, D=0')
> plot(diff(x,lag=1,differences=1),type='l',xlab='time',ylab='value',main='d=1, D=0')
> plot(diff(x,lag=par1,differences=1),type='l',xlab='time',ylab='value',main='d=0, D=1')
> plot(diff(diff(x,lag=1,differences=1),lag=par1,differences=1),type='l',xlab='time',ylab='value',main='d=1, D=1')
> par(op)
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/2ofc61323790226.ps tmp/2ofc61323790226.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.720   0.108   0.812