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 computationSun, 14 Dec 2008 08:44:05 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Dec/14/t1229269479wqjdvxaun6ha7im.htm/, Retrieved Wed, 15 May 2024 11:34:42 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=33434, Retrieved Wed, 15 May 2024 11:34:42 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact152
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Variance Reduction Matrix] [Paper: Variance R...] [2008-12-14 15:44:05] [8da7502cfecb272886bc60b3f290b8b8] [Current]
Feedback Forum

Post a new message
Dataseries X:
17,3
15,4
16,9
20,8
16,4
11,3
17,5
16,6
17,5
19,5
18,8
20,2
19,2
14,4
24,5
25,7
27,1
21
18,6
20
21,8
20,4
18
21,5
19,1
19,7
26
26,3
24,6
22,4
32
24
30
24,1
26,3
29,8
21,9
22,8
29,2
27,5
27,4
31
26,1
22,2
34
26,9
31,9
34,2
31,2
28,5
37,1
36
34,8
32,1
37,2
36,3
39,5
37,1
35,6
36,2
35,9
32,5
39,2
39,4
42,8
34,5
43,7
46,3
40,8
48,4
43,2
48,1
42,8




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'George Udny Yule' @ 72.249.76.132

\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 & 1 seconds \tabularnewline
R Server & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=33434&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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=33434&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=33434&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 time1 seconds
R Server'George Udny Yule' @ 72.249.76.132







Variance Reduction Matrix
V(Y[t],d=0,D=0)82.428302891933Range37.1Trim Var.59.0439903846154
V(Y[t],d=1,D=0)21.4027992957746Range20.1Trim Var.13.6638070436508
V(Y[t],d=2,D=0)64.4948169014085Range36.4Trim Var.41.4731694828469
V(Y[t],d=3,D=0)213.689242236025Range66.2Trim Var.131.228358011634
V(Y[t],d=0,D=1)16.9993661202186Range20Trim Var.10.0360957910015
V(Y[t],d=1,D=1)26.9180225988701Range26.5Trim Var.16.4094199860238
V(Y[t],d=2,D=1)80.4492285213326Range40.1Trim Var.50.9047750362845
V(Y[t],d=3,D=1)264.580012099214Range71.3Trim Var.173.444886877828
V(Y[t],d=0,D=2)46.7299659863946Range36.3Trim Var.22.3211074197121
V(Y[t],d=1,D=2)71.9615425531915Range51Trim Var.27.1292973286876
V(Y[t],d=2,D=2)199.397243293247Range72.3Trim Var.82.5900609756097
V(Y[t],d=3,D=2)625.162550724638Range130.6Trim Var.262.402846153846

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 82.428302891933 & Range & 37.1 & Trim Var. & 59.0439903846154 \tabularnewline
V(Y[t],d=1,D=0) & 21.4027992957746 & Range & 20.1 & Trim Var. & 13.6638070436508 \tabularnewline
V(Y[t],d=2,D=0) & 64.4948169014085 & Range & 36.4 & Trim Var. & 41.4731694828469 \tabularnewline
V(Y[t],d=3,D=0) & 213.689242236025 & Range & 66.2 & Trim Var. & 131.228358011634 \tabularnewline
V(Y[t],d=0,D=1) & 16.9993661202186 & Range & 20 & Trim Var. & 10.0360957910015 \tabularnewline
V(Y[t],d=1,D=1) & 26.9180225988701 & Range & 26.5 & Trim Var. & 16.4094199860238 \tabularnewline
V(Y[t],d=2,D=1) & 80.4492285213326 & Range & 40.1 & Trim Var. & 50.9047750362845 \tabularnewline
V(Y[t],d=3,D=1) & 264.580012099214 & Range & 71.3 & Trim Var. & 173.444886877828 \tabularnewline
V(Y[t],d=0,D=2) & 46.7299659863946 & Range & 36.3 & Trim Var. & 22.3211074197121 \tabularnewline
V(Y[t],d=1,D=2) & 71.9615425531915 & Range & 51 & Trim Var. & 27.1292973286876 \tabularnewline
V(Y[t],d=2,D=2) & 199.397243293247 & Range & 72.3 & Trim Var. & 82.5900609756097 \tabularnewline
V(Y[t],d=3,D=2) & 625.162550724638 & Range & 130.6 & Trim Var. & 262.402846153846 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=33434&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]82.428302891933[/C][C]Range[/C][C]37.1[/C][C]Trim Var.[/C][C]59.0439903846154[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]21.4027992957746[/C][C]Range[/C][C]20.1[/C][C]Trim Var.[/C][C]13.6638070436508[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]64.4948169014085[/C][C]Range[/C][C]36.4[/C][C]Trim Var.[/C][C]41.4731694828469[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]213.689242236025[/C][C]Range[/C][C]66.2[/C][C]Trim Var.[/C][C]131.228358011634[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]16.9993661202186[/C][C]Range[/C][C]20[/C][C]Trim Var.[/C][C]10.0360957910015[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]26.9180225988701[/C][C]Range[/C][C]26.5[/C][C]Trim Var.[/C][C]16.4094199860238[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]80.4492285213326[/C][C]Range[/C][C]40.1[/C][C]Trim Var.[/C][C]50.9047750362845[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]264.580012099214[/C][C]Range[/C][C]71.3[/C][C]Trim Var.[/C][C]173.444886877828[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]46.7299659863946[/C][C]Range[/C][C]36.3[/C][C]Trim Var.[/C][C]22.3211074197121[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]71.9615425531915[/C][C]Range[/C][C]51[/C][C]Trim Var.[/C][C]27.1292973286876[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]199.397243293247[/C][C]Range[/C][C]72.3[/C][C]Trim Var.[/C][C]82.5900609756097[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]625.162550724638[/C][C]Range[/C][C]130.6[/C][C]Trim Var.[/C][C]262.402846153846[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=33434&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=33434&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)82.428302891933Range37.1Trim Var.59.0439903846154
V(Y[t],d=1,D=0)21.4027992957746Range20.1Trim Var.13.6638070436508
V(Y[t],d=2,D=0)64.4948169014085Range36.4Trim Var.41.4731694828469
V(Y[t],d=3,D=0)213.689242236025Range66.2Trim Var.131.228358011634
V(Y[t],d=0,D=1)16.9993661202186Range20Trim Var.10.0360957910015
V(Y[t],d=1,D=1)26.9180225988701Range26.5Trim Var.16.4094199860238
V(Y[t],d=2,D=1)80.4492285213326Range40.1Trim Var.50.9047750362845
V(Y[t],d=3,D=1)264.580012099214Range71.3Trim Var.173.444886877828
V(Y[t],d=0,D=2)46.7299659863946Range36.3Trim Var.22.3211074197121
V(Y[t],d=1,D=2)71.9615425531915Range51Trim Var.27.1292973286876
V(Y[t],d=2,D=2)199.397243293247Range72.3Trim Var.82.5900609756097
V(Y[t],d=3,D=2)625.162550724638Range130.6Trim Var.262.402846153846



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(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) & smyxa<-table.row.end(a)
}
}
a<-table.end(a)
table.save(a,file='mytable.tab')