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(100.21
+ ,100.36
+ ,100.62
+ ,100.78
+ ,100.93
+ ,100.70
+ ,100.00
+ ,100.20
+ ,99.68
+ ,99.56
+ ,100.06
+ ,100.50
+ ,99.30
+ ,99.37
+ ,99.20
+ ,98.11
+ ,97.60
+ ,97.76
+ ,98.06
+ ,98.25
+ ,98.50
+ ,97.39
+ ,98.09
+ ,97.78
+ ,98.12
+ ,97.50
+ ,97.30
+ ,97.64
+ ,96.88
+ ,97.40
+ ,98.27
+ ,97.94
+ ,98.61
+ ,98.72
+ ,98.62
+ ,98.56
+ ,98.06
+ ,97.40
+ ,97.76
+ ,97.05
+ ,97.85
+ ,97.40
+ ,97.27
+ ,97.93
+ ,98.60
+ ,98.70
+ ,98.88
+ ,98.27
+ ,97.85
+ ,97.70
+ ,96.97
+ ,97.72
+ ,97.66
+ ,99.00
+ ,98.86
+ ,99.56
+ ,100.19
+ ,100.37
+ ,100.01
+ ,99.68
+ ,99.78
+ ,99.36
+ ,99.21
+ ,99.26
+ ,99.26
+ ,100.43
+ ,101.50
+ ,102.27
+ ,102.69
+ ,103.47
+ ,104.02
+ ,103.55
+ ,103.77
+ ,104.19
+ ,103.64
+ ,103.68
+ ,105.39
+ ,106.61
+ ,108.12
+ ,109.22
+ ,110.17
+ ,110.31
+ ,111.06
+ ,111.14
+ ,111.39
+ ,112.51
+ ,111.28
+ ,112.22
+ ,113.19
+ ,114.32
+ ,115.34
+ ,116.61
+ ,117.83
+ ,117.70
+ ,118.51
+ ,118.82
+ ,119.49
+ ,119.57
+ ,120.00
+ ,121.96
+ ,121.45
+ ,123.41
+ ,124.44
+ ,126.25
+ ,127.41
+ ,127.63
+ ,129.19
+ ,129.82
+ ,130.45
+ ,132.02
+ ,132.72
+ ,132.96
+ ,135.06
+ ,137.04
+ ,137.83
+ ,139.17
+ ,140.35
+ ,141.01
+ ,141.89
+ ,143.28
+ ,142.90
+ ,143.37
+ ,145.03
+ ,146.05
+ ,147.39
+ ,149.58
+ ,151.02
+ ,153.57
+ ,155.60
+ ,157.18
+ ,158.77
+ ,159.95
+ ,161.34
+ ,161.95
+ ,163.36
+ ,165.00
+ ,166.65
+ ,168.65
+ ,170.29
+ ,172.70
+ ,173.79
+ ,176.45
+ ,177.58
+ ,179.19
+ ,181.01
+ ,184.08
+ ,185.63
+ ,188.51
+ ,190.18
+ ,192.19
+ ,193.47
+ ,196.73
+ ,200.39
+ ,203.24
+ ,205.53
+ ,208.21
+ ,208.88
+ ,212.85
+ ,216.41
+ ,216.23
+ ,219.27
+ ,222.02
+ ,224.89
+ ,230.37
+ ,232.29
+ ,235.53
+ ,236.92
+ ,242.37
+ ,242.75
+ ,244.19
+ ,247.94
+ ,248.80
+ ,250.18
+ ,251.55
+ ,254.40
+ ,255.72
+ ,257.69
+ ,258.37
+ ,258.22
+ ,258.59
+ ,257.45
+ ,257.45
+ ,256.73
+ ,258.82
+ ,257.99
+ ,262.85
+ ,262.58
+ ,261.55
+ ,261.25
+ ,259.78
+ ,256.26
+ ,254.29
+ ,248.50
+ ,241.88
+ ,238.53
+ ,232.24
+ ,232.46
+ ,225.79
+ ,221.63
+ ,219.62
+ ,215.94
+ ,211.81
+ ,205.57
+ ,201.25
+ ,194.70
+ ,187.94
+ ,185.61
+ ,181.15
+ ,186.50
+ ,183.21
+ ,182.61
+ ,187.09
+ ,189.10
+ ,191.25
+ ,190.74
+ ,190.79)
> par1 = '12'
> #'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)
> n <- length(x)
> sx <- sort(x)
> 
> #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,'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(x,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/html/rcomp/tmp/1ptv31261334879.tab") 
> 
> 
> 
> proc.time()
   user  system elapsed 
  0.232   0.029   0.258