Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_variancereduction.wasp
Title produced by softwareVariance Reduction Matrix
Date of computationTue, 01 Feb 2011 11:45:45 +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/2011/Feb/01/t1296560666ctyygzd8kk0gb7z.htm/, Retrieved Wed, 08 May 2024 20:34:00 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=117991, Retrieved Wed, 08 May 2024 20:34:00 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact235
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Data Series] [] [2008-12-08 19:22:39] [d2d412c7f4d35ffbf5ee5ee89db327d4]
- RMP   [Variance Reduction Matrix] [ws9 tim damen vrm] [2010-12-05 12:31:01] [74be16979710d4c4e7c6647856088456]
- R  D      [Variance Reduction Matrix] [PaperTimDamen] [2011-02-01 11:45:45] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
349
336
331
327
323
322
385
405
412
411
410
415
414
411
408
410
411
416
479
498
502
498
499
506
510
509
502
495
490
490
553
570
573
572
575
580
580
574
563
556
546
545
605
628
631
626
614
606
602
589
574
558
552
546
607
636
631
623
618
605
619
596
570
546
528
506
555
568
564
553
541
542
540
521
505
491
482
478
523
531
532
540
525
533
531
508
495
482
470
466
515
518
516
511
500
498
494
476
458
443
430
424
476
481
470
460
451
450
444
429
421
400
389
384
432
446
431
423
416
416
413
399
386
374
365
365
418
428
424
421
417
423
423
419
406
398
390
391
444
460
455
456
452
459
461
451
443
439
430
436
488
506
502
501
501
515
521
520
512
509
505
511
570
592
594
586
586
592
594
594
586
586
572
572
563
563
555
555
554
554
601
601
622
622
617
617
606
606
595
595
599
599
600
600
592
592
575
575
567
567
555
555
555
555
608
608
631
631
629
629
624
624
610
610
616
616
621
621
604
604
584
584
574
574
555
555
545
545
599
599
620
620
608
608
590
590
579
579
580
580
579
579
572
572
560
560
551
551
537
537
541
541
588
588
607
607
599
599
578
578
563
563
566
566
561
561
554
554
540
540
526
526
512
512
505
505
554
554
584
584
569
569
540
540
522
522
526
526
527
527
516
516
503
503
489
489
479
479
475
475
524
524
552
552
532
532
511
511
492
492
492
492
493
493
481
481
462
462
457
457
442
442
439
439
488
488
521
521
501
501
485
485
464
464
460
460
467
467
460
460
448
448
443
443
436
436
431
431
484
484
510
510
513
513
503
503
471
471
471
471
476
476
475
475
470
470
461
461
455
455
456
456
517
517
525
525
523
523
519
519
509
509
512
512
519
519
517
517
510
510
509
509
501
501
507
507
569
569
580
580
578
578
565
565
547
547
555
555
562
561
555
544
537
543
594
611
613
611
594
595
591
589
584
573
567
569
621
629
628
612
595
597
593
590
580
574
573
573
620
626
620
588
566
557
561
549
532
526
511
499
555
565
542
527
510
514
517
508
493
490
469
478
528
534
518
506
502
516
528
533
536
537
524
536
587
597
581
564
558
575
580
575
563
552
537
545
601
604
586
564
549
551




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Gwilym Jenkins' @ www.wessa.org

\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 & 0 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ www.wessa.org \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=117991&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]0 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ www.wessa.org[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=117991&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=117991&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 time0 seconds
R Server'Gwilym Jenkins' @ www.wessa.org







Variance Reduction Matrix
V(Y[t],d=0,D=0)4193.35843558028Range314Trim Var.2810.1509688724
V(Y[t],d=1,D=0)273.761012510911Range95Trim Var.93.6619736214842
V(Y[t],d=2,D=0)442.602332957723Range133Trim Var.169.587746493726
V(Y[t],d=3,D=0)1242.30908679473Range206Trim Var.530.869565217391
V(Y[t],d=0,D=1)1957.92953636047Range177Trim Var.1455.07568627451
V(Y[t],d=1,D=1)205.172517710362Range132Trim Var.33.7591706975421
V(Y[t],d=2,D=1)384.934365762304Range120Trim Var.149.649292567897
V(Y[t],d=3,D=1)1272.24367986188Range222Trim Var.492.074679214866
V(Y[t],d=0,D=2)3817.17997675653Range250Trim Var.2882.37340896514
V(Y[t],d=1,D=2)735.511097223626Range219Trim Var.102.828824022122
V(Y[t],d=2,D=2)1417.65115141446Range220Trim Var.483.404767145505
V(Y[t],d=3,D=2)4733.52292361884Range438Trim Var.1710.73776747832

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 4193.35843558028 & Range & 314 & Trim Var. & 2810.1509688724 \tabularnewline
V(Y[t],d=1,D=0) & 273.761012510911 & Range & 95 & Trim Var. & 93.6619736214842 \tabularnewline
V(Y[t],d=2,D=0) & 442.602332957723 & Range & 133 & Trim Var. & 169.587746493726 \tabularnewline
V(Y[t],d=3,D=0) & 1242.30908679473 & Range & 206 & Trim Var. & 530.869565217391 \tabularnewline
V(Y[t],d=0,D=1) & 1957.92953636047 & Range & 177 & Trim Var. & 1455.07568627451 \tabularnewline
V(Y[t],d=1,D=1) & 205.172517710362 & Range & 132 & Trim Var. & 33.7591706975421 \tabularnewline
V(Y[t],d=2,D=1) & 384.934365762304 & Range & 120 & Trim Var. & 149.649292567897 \tabularnewline
V(Y[t],d=3,D=1) & 1272.24367986188 & Range & 222 & Trim Var. & 492.074679214866 \tabularnewline
V(Y[t],d=0,D=2) & 3817.17997675653 & Range & 250 & Trim Var. & 2882.37340896514 \tabularnewline
V(Y[t],d=1,D=2) & 735.511097223626 & Range & 219 & Trim Var. & 102.828824022122 \tabularnewline
V(Y[t],d=2,D=2) & 1417.65115141446 & Range & 220 & Trim Var. & 483.404767145505 \tabularnewline
V(Y[t],d=3,D=2) & 4733.52292361884 & Range & 438 & Trim Var. & 1710.73776747832 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=117991&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]4193.35843558028[/C][C]Range[/C][C]314[/C][C]Trim Var.[/C][C]2810.1509688724[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]273.761012510911[/C][C]Range[/C][C]95[/C][C]Trim Var.[/C][C]93.6619736214842[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]442.602332957723[/C][C]Range[/C][C]133[/C][C]Trim Var.[/C][C]169.587746493726[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]1242.30908679473[/C][C]Range[/C][C]206[/C][C]Trim Var.[/C][C]530.869565217391[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]1957.92953636047[/C][C]Range[/C][C]177[/C][C]Trim Var.[/C][C]1455.07568627451[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]205.172517710362[/C][C]Range[/C][C]132[/C][C]Trim Var.[/C][C]33.7591706975421[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]384.934365762304[/C][C]Range[/C][C]120[/C][C]Trim Var.[/C][C]149.649292567897[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]1272.24367986188[/C][C]Range[/C][C]222[/C][C]Trim Var.[/C][C]492.074679214866[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]3817.17997675653[/C][C]Range[/C][C]250[/C][C]Trim Var.[/C][C]2882.37340896514[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]735.511097223626[/C][C]Range[/C][C]219[/C][C]Trim Var.[/C][C]102.828824022122[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]1417.65115141446[/C][C]Range[/C][C]220[/C][C]Trim Var.[/C][C]483.404767145505[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]4733.52292361884[/C][C]Range[/C][C]438[/C][C]Trim Var.[/C][C]1710.73776747832[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=117991&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=117991&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)4193.35843558028Range314Trim Var.2810.1509688724
V(Y[t],d=1,D=0)273.761012510911Range95Trim Var.93.6619736214842
V(Y[t],d=2,D=0)442.602332957723Range133Trim Var.169.587746493726
V(Y[t],d=3,D=0)1242.30908679473Range206Trim Var.530.869565217391
V(Y[t],d=0,D=1)1957.92953636047Range177Trim Var.1455.07568627451
V(Y[t],d=1,D=1)205.172517710362Range132Trim Var.33.7591706975421
V(Y[t],d=2,D=1)384.934365762304Range120Trim Var.149.649292567897
V(Y[t],d=3,D=1)1272.24367986188Range222Trim Var.492.074679214866
V(Y[t],d=0,D=2)3817.17997675653Range250Trim Var.2882.37340896514
V(Y[t],d=1,D=2)735.511097223626Range219Trim Var.102.828824022122
V(Y[t],d=2,D=2)1417.65115141446Range220Trim Var.483.404767145505
V(Y[t],d=3,D=2)4733.52292361884Range438Trim Var.1710.73776747832



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