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(336766 + ,332201 + ,323529 + ,320287 + ,314768 + ,316870 + ,347093 + ,358764 + ,356615 + ,352559 + ,342807 + ,344952 + ,346958 + ,338371 + ,328380 + ,323111 + ,314764 + ,312235 + ,342631 + ,353912 + ,345948 + ,334046 + ,326268 + ,325939 + ,323758 + ,319275 + ,313328 + ,308612 + ,301096 + ,305610 + ,333005 + ,343322 + ,337653 + ,325595 + ,316726 + ,317578 + ,314484 + ,310166 + ,301792 + ,294495 + ,287256 + ,286224 + ,313654 + ,329656 + ,319784 + ,301995 + ,291228 + ,293041 + ,291908 + ,285121 + ,278310 + ,271158 + ,266586 + ,266426 + ,294354 + ,309380 + ,297710 + ,285409 + ,274580 + ,274402 + ,274452 + ,267700 + ,257841 + ,255124 + ,247377 + ,247823 + ,276919 + ,294271 + ,281758 + ,270434 + ,258848 + ,256674 + ,258882 + ,255060 + ,247698 + ,244779 + ,240901 + ,239933 + ,270247 + ,283893 + ,282348 + ,273570 + ,254756 + ,254354 + ,255843 + ,254490 + ,251995 + ,246339 + ,244019 + ,245953 + ,279806 + ,283111 + ,281097 + ,275964 + ,270694 + ,271901 + ,274412 + ,272433 + ,268361 + ,268586 + ,264768 + ,269974 + ,304744 + ,309365 + ,308347 + ,298427 + ,289231 + ,291975 + ,294912 + ,293488 + ,290555 + ,284736 + ,281818 + ,287854 + ,316263 + ,325412 + ,326011 + ,328282 + ,317480 + ,317539 + ,313737 + ,312276 + ,309391 + ,302950 + ,300316 + ,304035 + ,333476 + ,337698 + ,335932 + ,323931 + ,313927 + ,314485 + ,313218 + ,309664 + ,302963 + ,298989 + ,298423 + ,301631 + ,329765 + ,335083 + ,327616 + ,309119 + ,295916 + ,291413 + ,291542 + ,284678 + ,276475 + ,272566 + ,264981 + ,263290 + ,296806 + ,303598 + ,286994 + ,276427 + ,266424 + ,267153 + ,268381 + ,262522 + ,255542 + ,253158 + ,243803 + ,250741 + ,280445 + ,285257 + ,270976 + ,261076) > 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 a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/1senw1228250095.tab") > > > > proc.time() user system elapsed 0.225 0.029 0.251