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 computationWed, 24 Jan 2018 10:32:07 +0100
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2018/Jan/24/t1516786341okg2zpzfoukd5or.htm/, Retrieved Mon, 06 May 2024 10:18:40 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=312276, Retrieved Mon, 06 May 2024 10:18:40 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact46
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Variance Reduction Matrix] [variance reductio...] [2018-01-24 09:32:07] [1f2811b4b3eab344c60f6986921522b4] [Current]
Feedback Forum

Post a new message
Dataseries X:
1.72923686058208
0.122433126801145
0.523768788982079
-0.50175942601324
-1.86386120887019
0.0717422351702872
-0.380578531491989
-0.29851611788274
1.94316215667028
-1.24779901378139
1.35257277275093
2.84406850014719
1.12274878943606
-2.32540657982629
-2.89902363942181
0.246874012745119
-0.731094375584061
1.97769407549524
-0.089308186457892
-1.64931358279907
1.78988776616251
-0.655026882699021
-0.95518008982139
-0.501427837065876
-1.47468882156665
2.4255103110643
-0.388732570521789
-1.41970373897772
-0.435743764210552
0.588337253253952
-1.73087780677345
0.416192838768035
-0.135602934239543
1.5607508968729
-0.619488852927662
0.130061173484127
-0.485714783249833
1.19042117073251
0.504940688912012
-0.266725732456376
1.8987616014054
0.9835731579105
0.221023228988728
-0.176860990484583
-0.162198287814433
0.78911590880292
1.41209201053576
0.518161844628975
2.36462043434148
-1.24446105804486
1.09138861156198
0.422175612267996
-0.997989397904933
0.754650745534752
-1.44375392644115
0.568632265545353
1.35413305240232
-1.34288885715882
-0.0968579661829751
0.242767533935212
0.255744783559028
-0.243232899636027
0.461607264600104
0.0650797112167817
2.25999343044582
1.80462197944889
2.32548979115297
0.682758803307819
-0.335554613197298
-0.997929583796738
0.97507128242145
-1.17268516900173
-0.000337278483281606
-1.71460596672184
-1.34872342854985
-1.6521456119588
-1.48335049731373
0.925016701443022
-0.226525557050897
0.774905929042345
-0.339789522342372
0.89311055438168
-0.864578780222597
-1.95217915873499
-0.11317194501752
1.02470860347022
-1.25650654087286
-0.602130638597027
-0.83233134657693
-0.446577198612737
-4.01053303468812
0.288513333839785
-1.00798474916531
-0.572816189634855
-0.297410570159231
1.21578539483534
-0.727486991378146
-0.843943364912431
-1.20889338911855
-0.649656515850742
3.19984732249364
0.572738529480399
-0.0199646221520179
0.462316869083462
-1.28699501625665
-1.80757939159691
-2.15431670426712
0.707840316820782
1.36410372036211
0.314718394752158
-4.1002515835693
1.26992188862053
1.80144339707427
-0.42970220964945
0.529993653035631
1.18764422512472
-0.398233352310507
-0.421884162222733
-2.6898323919027
1.0140628812162
-0.542058945797581
-0.368690243263354
1.31411524584048
-0.428705333764204
-0.663582558495647
1.37204511801761
-1.27310996987184
0.810157913937622
1.5741305304336
-1.025150115941
-0.998930862337702
-0.0536208441845616
1.66337771592794
-0.178008462116163
-1.09164775518912
1.96797905166753
0.634220496887586
1.38764487483822
0.75570654001001
-2.3214750727902
1.62023993830839
0.423224828517651
-1.7051130142733
0.541981863740972
-0.593027838089475
0.213567450337977
-0.696840242545287
0.0777251570183703
-0.725328497765991
2.1315776363419
-0.33645518676905
-0.979182554450992
-0.242246084610725
-1.82980075345691
-1.14814874978528
0.496126806273696
-0.37803015357093
0.564631359789787
0.268965939275992
0.752661368724124
-2.30203173253409
2.82635859577673
-0.607627712639158
3.77405193258893
1.08426693162876
-2.2508985634727
0.290663952981133
-1.1346060583543
1.41583984485642
0.590041695581542
1.22056539827692
0.0970995178245223
-0.906020145960084
-0.597681226735446
0.437155726246844
0.742487563055592
1.22917286536281
-0.95259160353182
-0.75528682870179




Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R ServerBig Analytics Cloud Computing Center

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time2 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=312276&T=0

[TABLE]
[ROW]
Summary of computational transaction[/C][/ROW] [ROW]Raw Input[/C] view raw input (R code) [/C][/ROW] [ROW]Raw Output[/C]view raw output of R engine [/C][/ROW] [ROW]Computing time[/C]2 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=312276&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=312276&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 Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R ServerBig Analytics Cloud Computing Center







Variance Reduction Matrix
V(Y[t],d=0,D=0)1.67896Range7.8743Trim Var.0.98386
V(Y[t],d=1,D=0)3.50592Range9.78514Trim Var.2.18297
V(Y[t],d=2,D=0)10.59Range18.3475Trim Var.6.4921
V(Y[t],d=3,D=0)35.404Range33.1235Trim Var.21.9711
V(Y[t],d=0,D=1)3.10523Range9.19837Trim Var.1.80129
V(Y[t],d=1,D=1)6.42812Range14.4613Trim Var.4.18954
V(Y[t],d=2,D=1)19.5682Range22.4996Trim Var.13.2642
V(Y[t],d=3,D=1)65.8704Range37.6432Trim Var.44.8441
V(Y[t],d=0,D=2)9.02429Range16.1259Trim Var.5.35095
V(Y[t],d=1,D=2)18.8789Range23.8394Trim Var.12.242
V(Y[t],d=2,D=2)56.8189Range42.6814Trim Var.39.0947
V(Y[t],d=3,D=2)188.932Range69.4048Trim Var.130.032

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 1.67896 & Range & 7.8743 & Trim Var. & 0.98386 \tabularnewline
V(Y[t],d=1,D=0) & 3.50592 & Range & 9.78514 & Trim Var. & 2.18297 \tabularnewline
V(Y[t],d=2,D=0) & 10.59 & Range & 18.3475 & Trim Var. & 6.4921 \tabularnewline
V(Y[t],d=3,D=0) & 35.404 & Range & 33.1235 & Trim Var. & 21.9711 \tabularnewline
V(Y[t],d=0,D=1) & 3.10523 & Range & 9.19837 & Trim Var. & 1.80129 \tabularnewline
V(Y[t],d=1,D=1) & 6.42812 & Range & 14.4613 & Trim Var. & 4.18954 \tabularnewline
V(Y[t],d=2,D=1) & 19.5682 & Range & 22.4996 & Trim Var. & 13.2642 \tabularnewline
V(Y[t],d=3,D=1) & 65.8704 & Range & 37.6432 & Trim Var. & 44.8441 \tabularnewline
V(Y[t],d=0,D=2) & 9.02429 & Range & 16.1259 & Trim Var. & 5.35095 \tabularnewline
V(Y[t],d=1,D=2) & 18.8789 & Range & 23.8394 & Trim Var. & 12.242 \tabularnewline
V(Y[t],d=2,D=2) & 56.8189 & Range & 42.6814 & Trim Var. & 39.0947 \tabularnewline
V(Y[t],d=3,D=2) & 188.932 & Range & 69.4048 & Trim Var. & 130.032 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=312276&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]1.67896[/C][C]Range[/C][C]7.8743[/C][C]Trim Var.[/C][C]0.98386[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]3.50592[/C][C]Range[/C][C]9.78514[/C][C]Trim Var.[/C][C]2.18297[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]10.59[/C][C]Range[/C][C]18.3475[/C][C]Trim Var.[/C][C]6.4921[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]35.404[/C][C]Range[/C][C]33.1235[/C][C]Trim Var.[/C][C]21.9711[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]3.10523[/C][C]Range[/C][C]9.19837[/C][C]Trim Var.[/C][C]1.80129[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]6.42812[/C][C]Range[/C][C]14.4613[/C][C]Trim Var.[/C][C]4.18954[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]19.5682[/C][C]Range[/C][C]22.4996[/C][C]Trim Var.[/C][C]13.2642[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]65.8704[/C][C]Range[/C][C]37.6432[/C][C]Trim Var.[/C][C]44.8441[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]9.02429[/C][C]Range[/C][C]16.1259[/C][C]Trim Var.[/C][C]5.35095[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]18.8789[/C][C]Range[/C][C]23.8394[/C][C]Trim Var.[/C][C]12.242[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]56.8189[/C][C]Range[/C][C]42.6814[/C][C]Trim Var.[/C][C]39.0947[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]188.932[/C][C]Range[/C][C]69.4048[/C][C]Trim Var.[/C][C]130.032[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=312276&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=312276&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)1.67896Range7.8743Trim Var.0.98386
V(Y[t],d=1,D=0)3.50592Range9.78514Trim Var.2.18297
V(Y[t],d=2,D=0)10.59Range18.3475Trim Var.6.4921
V(Y[t],d=3,D=0)35.404Range33.1235Trim Var.21.9711
V(Y[t],d=0,D=1)3.10523Range9.19837Trim Var.1.80129
V(Y[t],d=1,D=1)6.42812Range14.4613Trim Var.4.18954
V(Y[t],d=2,D=1)19.5682Range22.4996Trim Var.13.2642
V(Y[t],d=3,D=1)65.8704Range37.6432Trim Var.44.8441
V(Y[t],d=0,D=2)9.02429Range16.1259Trim Var.5.35095
V(Y[t],d=1,D=2)18.8789Range23.8394Trim Var.12.242
V(Y[t],d=2,D=2)56.8189Range42.6814Trim Var.39.0947
V(Y[t],d=3,D=2)188.932Range69.4048Trim Var.130.032



Parameters (Session):
Parameters (R input):
par1 = 12 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
n.orig <- length(x)
x <- na.omit(x)
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)
if(n.orig!=n) {
a<-table.row.start(a)
a<-table.element(a,'Warning: NAs were removed from the time series! The results shown below will only be correct if the NAs are all located at the start and/or end of the time series.',6,F)
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,signif(var(myx), digits=6))
a<-table.element(a,'Range',header=TRUE)
a<-table.element(a,signif(max(myx)-min(myx), digits=6))
a<-table.element(a,'Trim Var.',header=TRUE)
smyx <- sort(myx)
sn <- length(smyx)
a<-table.element(a,signif(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()