x <- c(5.81,5.76,5.99,6.12,6.03,6.25,5.80,5.67,5.89,5.91,5.86,6.07,6.27,6.68,6.77,6.71,6.62,6.50,5.89,6.05,6.43,6.47,6.62,6.77,6.70,6.95,6.73,7.07,7.28,7.32,6.76,6.93,6.99,7.16,7.28,7.08,7.34,7.87,6.28,6.30,6.36,6.28,5.89,6.04,5.96,6.10,6.26,6.02,6.25,6.41,6.22,6.57,6.18,6.26,6.10,6.02,6.06,6.35,6.21,6.48,6.74,6.53,6.80,6.75,6.56,6.66,6.18,6.40,6.43,6.54,6.44,6.64,6.82,6.97,7.00,6.91,6.74,6.98,6.37,6.56,6.63,6.87,6.68,6.75,6.84,7.15,7.09,6.97,7.15)
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/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/1rbf31291658954.tab")