R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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(6392.3,8686.4,9244.7,8182.7,7451.4,7988.8,8243.5,8843,9092.7,8246.7,9311.7,8341.2,7116.7,9635.7,9815.4,8611.3,8297.8,8715.1,8919.9,10085.8,9511.7,8991.3,10311.2,8895.4,7449.8,10084,9859.4,9100.1,8920.8,8502.7,8599.6,10394.4,9290.4,8742.2,10217.3,8639,8139.6,10779.1,10427.7,10349.1,10036.4,9492.1,10638.8,12054.5,10324.7,11817.3,11008.9,9996.6,9419.5,11958.8,12594.6,11890.6,10871.7,11835.7,11542.2,13093.7,11180.2,12035.7,12112,10875.2,9897.3,11672.1,12385.7,11405.6,9830.9,11025.1,10853.8,12252.6,11839.4,11669.1,11601.4,11178.4,9516.4,12102.8,12989,11610.2,10205.5,11356.2,11307.1,12648.6,11947.2,11714.1,12192.5,11268.8,9097.4,12639.8,13040.1,11687.3,11191.7,11391.9,11793.1,13933.2,12778.1,11810.3,13698.4,11956.6,10723.8,13938.9,13979.8,13807.4,12973.9,12509.8,12934.1,14908.3,13772.1,13012.6,14049.9,11816.5,11593.2,14466.2,13615.9,14733.9,13880.7,13527.5,13584,16170.2,13260.6,14741.9,15486.5,13154.5,12621.2)
> 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/1ew5h1229349845.tab") 
> 
> 
> 
> proc.time()
   user  system elapsed 
  0.219   0.031   0.257