Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_variancereduction.wasp
Title produced by softwareVariance Reduction Matrix
Date of computationTue, 28 Dec 2010 08:16:14 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Dec/28/t1293524243its3cezpm6kclt5.htm/, Retrieved Sat, 04 May 2024 20:29:59 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=116215, Retrieved Sat, 04 May 2024 20:29:59 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact140
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Variance Reduction Matrix] [] [2010-12-24 15:11:40] [afd301b68d203992295e6972aed62880]
-    D    [Variance Reduction Matrix] [] [2010-12-28 08:16:14] [5a59313293e5c9f616ad36f6edd018c5] [Current]
Feedback Forum

Post a new message
Dataseries X:
547.344		
554.788	
562.325	
560.854	
555.332	
543.599	
536.662	
542.722	
593.530	
610.763	
612.613	
611.324	
594.167	
595.454	
590.865	
589.379	
584.428	
573.100	
567.456	
569.028	
620.735	
628.884	
628.232	
612.117	
595.404	
597.141	
593.408	
590.072	
579.799	
574.205	
572.775	
572.942	
619.567	
625.809	
619.916	
587.625	
565.742	
557.274	
560.576	
548.854	
531.673	
525.919	
511.038	
498.662	
555.362	
564.591	
541.657	
527.070	
509.846	
514.258	
516.922	
507.561	
492.622	
490.243	
469.357	
477.580	
528.379	
533.590	
517.945	
506.174	
501.866	
516.141	
528.222	
532.638	
536.322		
536.535		
523.597		
536.214		
586.570		
596.594		
580.523		
564.478		
557.560		
575.093		
580.112		
574.761		
563.250		
551.531		
537.034		
544.686		
600.991		
604.378		
586.111		
563.668		
548.604		




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 3 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=116215&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=116215&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=116215&T=0

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135







Variance Reduction Matrix
V(Y[t],d=0,D=0)1400.72047372913Range159.527Trim Var.924.987029972252
V(Y[t],d=1,D=0)351.11785313253Range88.991Trim Var.159.615813683080
V(Y[t],d=2,D=0)510.871951271819Range121.994Trim Var.241.500380612633
V(Y[t],d=3,D=0)1175.05501701536Range183.118Trim Var.480.47374166021
V(Y[t],d=0,D=1)1424.90923824772Range141.263Trim Var.1065.64893049135
V(Y[t],d=1,D=1)67.4066052235915Range37.6400000000000Trim Var.39.5151629166666
V(Y[t],d=2,D=1)106.211192356942Range61.2450000000001Trim Var.50.0450972570403
V(Y[t],d=3,D=1)334.266371575362Range107.877000000000Trim Var.166.272645314648
V(Y[t],d=0,D=2)2617.4206922765Range174.115Trim Var.2004.07575525544
V(Y[t],d=1,D=2)167.883231453108Range67.698Trim Var.95.0896401568834
V(Y[t],d=2,D=2)291.515976119813Range94.7980000000003Trim Var.158.721538791727
V(Y[t],d=3,D=2)934.181099072295Range137.188000000001Trim Var.522.529160519982

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 1400.72047372913 & Range & 159.527 & Trim Var. & 924.987029972252 \tabularnewline
V(Y[t],d=1,D=0) & 351.11785313253 & Range & 88.991 & Trim Var. & 159.615813683080 \tabularnewline
V(Y[t],d=2,D=0) & 510.871951271819 & Range & 121.994 & Trim Var. & 241.500380612633 \tabularnewline
V(Y[t],d=3,D=0) & 1175.05501701536 & Range & 183.118 & Trim Var. & 480.47374166021 \tabularnewline
V(Y[t],d=0,D=1) & 1424.90923824772 & Range & 141.263 & Trim Var. & 1065.64893049135 \tabularnewline
V(Y[t],d=1,D=1) & 67.4066052235915 & Range & 37.6400000000000 & Trim Var. & 39.5151629166666 \tabularnewline
V(Y[t],d=2,D=1) & 106.211192356942 & Range & 61.2450000000001 & Trim Var. & 50.0450972570403 \tabularnewline
V(Y[t],d=3,D=1) & 334.266371575362 & Range & 107.877000000000 & Trim Var. & 166.272645314648 \tabularnewline
V(Y[t],d=0,D=2) & 2617.4206922765 & Range & 174.115 & Trim Var. & 2004.07575525544 \tabularnewline
V(Y[t],d=1,D=2) & 167.883231453108 & Range & 67.698 & Trim Var. & 95.0896401568834 \tabularnewline
V(Y[t],d=2,D=2) & 291.515976119813 & Range & 94.7980000000003 & Trim Var. & 158.721538791727 \tabularnewline
V(Y[t],d=3,D=2) & 934.181099072295 & Range & 137.188000000001 & Trim Var. & 522.529160519982 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=116215&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]1400.72047372913[/C][C]Range[/C][C]159.527[/C][C]Trim Var.[/C][C]924.987029972252[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]351.11785313253[/C][C]Range[/C][C]88.991[/C][C]Trim Var.[/C][C]159.615813683080[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]510.871951271819[/C][C]Range[/C][C]121.994[/C][C]Trim Var.[/C][C]241.500380612633[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]1175.05501701536[/C][C]Range[/C][C]183.118[/C][C]Trim Var.[/C][C]480.47374166021[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]1424.90923824772[/C][C]Range[/C][C]141.263[/C][C]Trim Var.[/C][C]1065.64893049135[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]67.4066052235915[/C][C]Range[/C][C]37.6400000000000[/C][C]Trim Var.[/C][C]39.5151629166666[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]106.211192356942[/C][C]Range[/C][C]61.2450000000001[/C][C]Trim Var.[/C][C]50.0450972570403[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]334.266371575362[/C][C]Range[/C][C]107.877000000000[/C][C]Trim Var.[/C][C]166.272645314648[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]2617.4206922765[/C][C]Range[/C][C]174.115[/C][C]Trim Var.[/C][C]2004.07575525544[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]167.883231453108[/C][C]Range[/C][C]67.698[/C][C]Trim Var.[/C][C]95.0896401568834[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]291.515976119813[/C][C]Range[/C][C]94.7980000000003[/C][C]Trim Var.[/C][C]158.721538791727[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]934.181099072295[/C][C]Range[/C][C]137.188000000001[/C][C]Trim Var.[/C][C]522.529160519982[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=116215&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=116215&T=1

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Variance Reduction Matrix
V(Y[t],d=0,D=0)1400.72047372913Range159.527Trim Var.924.987029972252
V(Y[t],d=1,D=0)351.11785313253Range88.991Trim Var.159.615813683080
V(Y[t],d=2,D=0)510.871951271819Range121.994Trim Var.241.500380612633
V(Y[t],d=3,D=0)1175.05501701536Range183.118Trim Var.480.47374166021
V(Y[t],d=0,D=1)1424.90923824772Range141.263Trim Var.1065.64893049135
V(Y[t],d=1,D=1)67.4066052235915Range37.6400000000000Trim Var.39.5151629166666
V(Y[t],d=2,D=1)106.211192356942Range61.2450000000001Trim Var.50.0450972570403
V(Y[t],d=3,D=1)334.266371575362Range107.877000000000Trim Var.166.272645314648
V(Y[t],d=0,D=2)2617.4206922765Range174.115Trim Var.2004.07575525544
V(Y[t],d=1,D=2)167.883231453108Range67.698Trim Var.95.0896401568834
V(Y[t],d=2,D=2)291.515976119813Range94.7980000000003Trim Var.158.721538791727
V(Y[t],d=3,D=2)934.181099072295Range137.188000000001Trim Var.522.529160519982



Parameters (Session):
par1 = 12 ;
Parameters (R input):
par1 = 12 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
n <- length(x)
sx <- sort(x)
load(file='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(myx,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) & smyxa<-table.row.end(a)
}
}
a<-table.end(a)
table.save(a,file='mytable.tab')
bitmap(file='pic0.png')
op <- par(mfrow=c(2,2))
plot(x,type='l',xlab='time',ylab='value',main='d=0, D=0')
plot(diff(x,lag=1,differences=1),type='l',xlab='time',ylab='value',main='d=1, D=0')
plot(diff(x,lag=par1,differences=1),type='l',xlab='time',ylab='value',main='d=0, D=1')
plot(diff(diff(x,lag=1,differences=1),lag=par1,differences=1),type='l',xlab='time',ylab='value',main='d=1, D=1')
par(op)
dev.off()