Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_variancereduction.wasp
Title produced by softwareVariance Reduction Matrix
Date of computationThu, 22 Nov 2007 07:04:43 -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/2007/Nov/22/t1195739798bjg82to12ifsghq.htm/, Retrieved Fri, 03 May 2024 03:43:59 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=6001, Retrieved Fri, 03 May 2024 03:43:59 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact177
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Variance Reduction Matrix] [VRM P-index GS'n ...] [2007-11-22 14:04:43] [4df98167d5cf79c69ce763f2d4ef5b15] [Current]
Feedback Forum

Post a new message
Dataseries X:
96,9
98,0
97,9
100,9
103,9
103,1
102,5
104,3
102,6
101,7
102,8
105,4
110,9
113,5
116,3
124,0
128,8
133,5
132,6
128,4
127,3
126,7
123,3
123,2
124,4
128,2
128,7
135,7
139,0
145,4
142,4
137,7
137,0
137,1
139,3
139,6
140,4
142,3
148,3
157,7
161,6
161,7
171,8
185,1
176,7
184,4
183,0
180,9
187,0
189,9
193,8
194,5
198,7
204,7
213,2
214,7
211,0
213,2
206,3
210,8




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

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

[TABLE]
[ROW][C]Summary of compuational 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]0 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=6001&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=6001&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 compuational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135







Variance Reduction Matrix
V(Y[t],d=0,D=0)1351.49031355932Range117.8Trim Var.1111.83358490566
V(Y[t],d=1,D=0)16.9725014611338Range21.7Trim Var.9.69071117561684
V(Y[t],d=2,D=0)31.6386085904416Range37.8Trim Var.13.1696945701357
V(Y[t],d=3,D=0)92.6645488721804Range63Trim Var.33.4361019607843
V(Y[t],d=0,D=1)149.337003546099Range38.3Trim Var.117.164651567944
V(Y[t],d=1,D=1)28.6058926919519Range29.8Trim Var.12.3237195121951
V(Y[t],d=2,D=1)61.2552125603865Range45.1000000000001Trim Var.25.2398717948717
V(Y[t],d=3,D=1)175.715131313132Range71.6Trim Var.60.4453576248313
V(Y[t],d=0,D=2)343.840349206349Range58.5Trim Var.278.237651209677
V(Y[t],d=1,D=2)89.7386890756304Range48.3Trim Var.45.9029892473119
V(Y[t],d=2,D=2)195.068065953655Range69.1000000000001Trim Var.92.7349885057473
V(Y[t],d=3,D=2)574.03484848485Range125Trim Var.257.404211822661

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 1351.49031355932 & Range & 117.8 & Trim Var. & 1111.83358490566 \tabularnewline
V(Y[t],d=1,D=0) & 16.9725014611338 & Range & 21.7 & Trim Var. & 9.69071117561684 \tabularnewline
V(Y[t],d=2,D=0) & 31.6386085904416 & Range & 37.8 & Trim Var. & 13.1696945701357 \tabularnewline
V(Y[t],d=3,D=0) & 92.6645488721804 & Range & 63 & Trim Var. & 33.4361019607843 \tabularnewline
V(Y[t],d=0,D=1) & 149.337003546099 & Range & 38.3 & Trim Var. & 117.164651567944 \tabularnewline
V(Y[t],d=1,D=1) & 28.6058926919519 & Range & 29.8 & Trim Var. & 12.3237195121951 \tabularnewline
V(Y[t],d=2,D=1) & 61.2552125603865 & Range & 45.1000000000001 & Trim Var. & 25.2398717948717 \tabularnewline
V(Y[t],d=3,D=1) & 175.715131313132 & Range & 71.6 & Trim Var. & 60.4453576248313 \tabularnewline
V(Y[t],d=0,D=2) & 343.840349206349 & Range & 58.5 & Trim Var. & 278.237651209677 \tabularnewline
V(Y[t],d=1,D=2) & 89.7386890756304 & Range & 48.3 & Trim Var. & 45.9029892473119 \tabularnewline
V(Y[t],d=2,D=2) & 195.068065953655 & Range & 69.1000000000001 & Trim Var. & 92.7349885057473 \tabularnewline
V(Y[t],d=3,D=2) & 574.03484848485 & Range & 125 & Trim Var. & 257.404211822661 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=6001&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]1351.49031355932[/C][C]Range[/C][C]117.8[/C][C]Trim Var.[/C][C]1111.83358490566[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]16.9725014611338[/C][C]Range[/C][C]21.7[/C][C]Trim Var.[/C][C]9.69071117561684[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]31.6386085904416[/C][C]Range[/C][C]37.8[/C][C]Trim Var.[/C][C]13.1696945701357[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]92.6645488721804[/C][C]Range[/C][C]63[/C][C]Trim Var.[/C][C]33.4361019607843[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]149.337003546099[/C][C]Range[/C][C]38.3[/C][C]Trim Var.[/C][C]117.164651567944[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]28.6058926919519[/C][C]Range[/C][C]29.8[/C][C]Trim Var.[/C][C]12.3237195121951[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]61.2552125603865[/C][C]Range[/C][C]45.1000000000001[/C][C]Trim Var.[/C][C]25.2398717948717[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]175.715131313132[/C][C]Range[/C][C]71.6[/C][C]Trim Var.[/C][C]60.4453576248313[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]343.840349206349[/C][C]Range[/C][C]58.5[/C][C]Trim Var.[/C][C]278.237651209677[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]89.7386890756304[/C][C]Range[/C][C]48.3[/C][C]Trim Var.[/C][C]45.9029892473119[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]195.068065953655[/C][C]Range[/C][C]69.1000000000001[/C][C]Trim Var.[/C][C]92.7349885057473[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]574.03484848485[/C][C]Range[/C][C]125[/C][C]Trim Var.[/C][C]257.404211822661[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=6001&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=6001&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)1351.49031355932Range117.8Trim Var.1111.83358490566
V(Y[t],d=1,D=0)16.9725014611338Range21.7Trim Var.9.69071117561684
V(Y[t],d=2,D=0)31.6386085904416Range37.8Trim Var.13.1696945701357
V(Y[t],d=3,D=0)92.6645488721804Range63Trim Var.33.4361019607843
V(Y[t],d=0,D=1)149.337003546099Range38.3Trim Var.117.164651567944
V(Y[t],d=1,D=1)28.6058926919519Range29.8Trim Var.12.3237195121951
V(Y[t],d=2,D=1)61.2552125603865Range45.1000000000001Trim Var.25.2398717948717
V(Y[t],d=3,D=1)175.715131313132Range71.6Trim Var.60.4453576248313
V(Y[t],d=0,D=2)343.840349206349Range58.5Trim Var.278.237651209677
V(Y[t],d=1,D=2)89.7386890756304Range48.3Trim Var.45.9029892473119
V(Y[t],d=2,D=2)195.068065953655Range69.1000000000001Trim Var.92.7349885057473
V(Y[t],d=3,D=2)574.03484848485Range125Trim Var.257.404211822661



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')