R version 2.7.0 (2008-04-22)
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(9568.3
+ ,9920.3
+ ,11353.5
+ ,9247.5
+ ,10114.2
+ ,10763.1
+ ,8456.1
+ ,8071.6
+ ,10328
+ ,10551.4
+ ,10186.1
+ ,8821.6
+ ,9841.3
+ ,10233.6
+ ,10794.6
+ ,10289.3
+ ,10513.4
+ ,10607.6
+ ,9707.4
+ ,8103.5
+ ,10982.6
+ ,11836.9
+ ,10517.5
+ ,9810.5
+ ,10374.8
+ ,10855.3
+ ,11671.3
+ ,11901.2
+ ,10846.4
+ ,11917.5
+ ,11362.8
+ ,9314.5
+ ,12605.9
+ ,12815.1
+ ,11254.5
+ ,11111.8
+ ,11282.9
+ ,11554.5
+ ,12935.6
+ ,12146.3
+ ,11615.3
+ ,13214.8
+ ,11735.5
+ ,9522.3
+ ,12694.8
+ ,12317.6
+ ,11450
+ ,11380.9
+ ,10604.6
+ ,10972.2
+ ,13331.5
+ ,11733.1
+ ,11284.7
+ ,13295.8
+ ,11881.4
+ ,10374.2
+ ,13828
+ ,13490.5
+ ,13092.2
+ ,13184.4
+ ,12398.4
+ ,13882.3
+ ,15861.5
+ ,13286.1
+ ,15634.9
+ ,14211
+ ,13646.8
+ ,12224.6
+ ,15916.4
+ ,16535.9
+ ,15796
+ ,14418.6
+ ,15044.5
+ ,14944.2
+ ,16754.8
+ ,14254
+ ,15454.9
+ ,15644.8
+ ,14568.3
+ ,12520.2
+ ,14803
+ ,15873.2
+ ,14755.3
+ ,12875.1
+ ,14291.1
+ ,14205.3
+ ,15859.4
+ ,15258.9
+ ,15498.6
+ ,15106.5
+ ,15023.6
+ ,12083
+ ,15761.3
+ ,16943
+ ,15070.3
+ ,13659.6
+ ,14768.9
+ ,14725.1
+ ,15998.1
+ ,15370.6
+ ,14956.9
+ ,15469.7
+ ,15101.8
+ ,11703.7
+ ,16283.6
+ ,16726.5
+ ,14968.9
+ ,14861
+ ,14583.3
+ ,15305.8
+ ,17903.9
+ ,16379.4
+ ,15420.3
+ ,17870.5
+ ,15912.8
+ ,13866.5
+ ,17823.2
+ ,17872
+ ,17420.4
+ ,16704.4
+ ,15991.2
+ ,16583.6
+ ,19123.5
+ ,17838.7
+ ,17209.4
+ ,18586.5
+ ,16258.1
+ ,15141.6
+ ,19202.1
+ ,17746.5
+ ,19090.1
+ ,18040.3
+ ,17515.5
+ ,17751.8
+ ,21072.4
+ ,17170
+ ,19439.5
+ ,19795.4
+ ,17574.9
+ ,16165.4
+ ,19464.6
+ ,19932.1
+ ,19961.2
+ ,17343.4
+ ,18924.2
+ ,18574.1
+ ,21350.6
+ ,18594.6
+ ,19823.1
+ ,20844.4
+ ,19640.2
+ ,17735.4
+ ,19813.6
+ ,22238.5
+ ,20682.2
+ ,17818.6
+ ,21872.1
+ ,22117
+ ,21865.9
+ ,23451.3
+ ,20953.7
+ ,22497.3)
> 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/127as1228661351.tab") 
> 
> 
> 
> proc.time()
   user  system elapsed 
  0.891   0.020   0.905