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 computationThu, 15 Jan 2015 21:27:20 +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/2015/Jan/15/t1421357319j00rjnue63zpv1s.htm/, Retrieved Thu, 31 Oct 2024 22:49:28 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=273006, Retrieved Thu, 31 Oct 2024 22:49:28 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact64
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Variance Reduction Matrix] [Berekening 1 time...] [2015-01-15 21:27:20] [8d9a76f7eb5407bb245b7faf6d068dcd] [Current]
Feedback Forum

Post a new message
Dataseries X:
1775
2197
2920
4240
5415
6136
6719
6234
7152
3646
2165
2803
1615
2350
3350
3536
5834
6767
5993
7276
5641
3477
2247
2466
1567
2237
2598
3729
5715
5776
5852
6878
5488
3583
2054
2282
1552
2261
2446
3519
5161
5085
5711
6057
5224
3363
1899
2115
1491
2061
2419
3430
4778
4862
6176
5664
5529
3418
1941
2402
1579
2146
2462
3695
4831
5134
6250
5760
6249
2917
1741
2359
1511
2059
2635
2867
4403
5720
4502
5749
5627
2846
1762
2429
1169
2154
2249
2687
4359
5382
4459
6398
4596
3024
1887
2070
1351
2218
2461
3028
4784
4975
4607
6249
4809
3157
1910
2228
1594
2467
2222
3607
4685
4962
5770
5480
5000
3228
1993
2288
1588
2105
2191
3591
4668
4885
5822
5599
5340
3082
2010
2301




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gertrude Mary Cox' @ cox.wessa.net

\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 & 'Gertrude Mary Cox' @ cox.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=273006&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]'Gertrude Mary Cox' @ cox.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=273006&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=273006&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'Gertrude Mary Cox' @ cox.wessa.net







Variance Reduction Matrix
V(Y[t],d=0,D=0)2781150Range6107Trim Var.2241210
V(Y[t],d=1,D=0)1315990Range5804Trim Var.841009
V(Y[t],d=2,D=0)2106150Range7286Trim Var.1343060
V(Y[t],d=3,D=0)6182220Range13052Trim Var.3791230
V(Y[t],d=0,D=1)172027Range2911Trim Var.73169.5
V(Y[t],d=1,D=1)439812Range4321Trim Var.169031
V(Y[t],d=2,D=1)1484010Range8392Trim Var.535149
V(Y[t],d=3,D=1)5200960Range15662Trim Var.2012440
V(Y[t],d=0,D=2)329958Range3527Trim Var.146145
V(Y[t],d=1,D=2)796680Range4934Trim Var.376653
V(Y[t],d=2,D=2)2615530Range9055Trim Var.1216430
V(Y[t],d=3,D=2)9053510Range17759Trim Var.4263350

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 2781150 & Range & 6107 & Trim Var. & 2241210 \tabularnewline
V(Y[t],d=1,D=0) & 1315990 & Range & 5804 & Trim Var. & 841009 \tabularnewline
V(Y[t],d=2,D=0) & 2106150 & Range & 7286 & Trim Var. & 1343060 \tabularnewline
V(Y[t],d=3,D=0) & 6182220 & Range & 13052 & Trim Var. & 3791230 \tabularnewline
V(Y[t],d=0,D=1) & 172027 & Range & 2911 & Trim Var. & 73169.5 \tabularnewline
V(Y[t],d=1,D=1) & 439812 & Range & 4321 & Trim Var. & 169031 \tabularnewline
V(Y[t],d=2,D=1) & 1484010 & Range & 8392 & Trim Var. & 535149 \tabularnewline
V(Y[t],d=3,D=1) & 5200960 & Range & 15662 & Trim Var. & 2012440 \tabularnewline
V(Y[t],d=0,D=2) & 329958 & Range & 3527 & Trim Var. & 146145 \tabularnewline
V(Y[t],d=1,D=2) & 796680 & Range & 4934 & Trim Var. & 376653 \tabularnewline
V(Y[t],d=2,D=2) & 2615530 & Range & 9055 & Trim Var. & 1216430 \tabularnewline
V(Y[t],d=3,D=2) & 9053510 & Range & 17759 & Trim Var. & 4263350 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=273006&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]2781150[/C][C]Range[/C][C]6107[/C][C]Trim Var.[/C][C]2241210[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]1315990[/C][C]Range[/C][C]5804[/C][C]Trim Var.[/C][C]841009[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]2106150[/C][C]Range[/C][C]7286[/C][C]Trim Var.[/C][C]1343060[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]6182220[/C][C]Range[/C][C]13052[/C][C]Trim Var.[/C][C]3791230[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]172027[/C][C]Range[/C][C]2911[/C][C]Trim Var.[/C][C]73169.5[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]439812[/C][C]Range[/C][C]4321[/C][C]Trim Var.[/C][C]169031[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]1484010[/C][C]Range[/C][C]8392[/C][C]Trim Var.[/C][C]535149[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]5200960[/C][C]Range[/C][C]15662[/C][C]Trim Var.[/C][C]2012440[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]329958[/C][C]Range[/C][C]3527[/C][C]Trim Var.[/C][C]146145[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]796680[/C][C]Range[/C][C]4934[/C][C]Trim Var.[/C][C]376653[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]2615530[/C][C]Range[/C][C]9055[/C][C]Trim Var.[/C][C]1216430[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]9053510[/C][C]Range[/C][C]17759[/C][C]Trim Var.[/C][C]4263350[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=273006&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=273006&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)2781150Range6107Trim Var.2241210
V(Y[t],d=1,D=0)1315990Range5804Trim Var.841009
V(Y[t],d=2,D=0)2106150Range7286Trim Var.1343060
V(Y[t],d=3,D=0)6182220Range13052Trim Var.3791230
V(Y[t],d=0,D=1)172027Range2911Trim Var.73169.5
V(Y[t],d=1,D=1)439812Range4321Trim Var.169031
V(Y[t],d=2,D=1)1484010Range8392Trim Var.535149
V(Y[t],d=3,D=1)5200960Range15662Trim Var.2012440
V(Y[t],d=0,D=2)329958Range3527Trim Var.146145
V(Y[t],d=1,D=2)796680Range4934Trim Var.376653
V(Y[t],d=2,D=2)2615530Range9055Trim Var.1216430
V(Y[t],d=3,D=2)9053510Range17759Trim Var.4263350



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,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()