FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_meanplot.wasp
Title produced by softwareMean Plot
Date of computationThu, 18 Dec 2008 09:03:30 -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/18/t1229616294zeua68nugvu9hap.htm/, Retrieved Sat, 11 May 2024 06:16:36 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=34867, Retrieved Sat, 11 May 2024 06:16:36 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact224

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [(Partial) Autocorrelation Function] [(Partial) Autocor...] [2008-12-11 14:56:15] [87cabf13a90315c7085b765dcebb7412]
-   P   [(Partial) Autocorrelation Function] [(Partial) Autocor...] [2008-12-11 15:06:01] [87cabf13a90315c7085b765dcebb7412]
-   PD    [(Partial) Autocorrelation Function] [(Partial) Autocor...] [2008-12-11 15:23:01] [87cabf13a90315c7085b765dcebb7412]
-   P       [(Partial) Autocorrelation Function] [(Partial) Autocor...] [2008-12-12 09:41:51] [87cabf13a90315c7085b765dcebb7412]
-   PD        [(Partial) Autocorrelation Function] [autocorrelation 2...] [2008-12-18 15:13:35] [631938996a408f8d8cf3d9850ca0cd03]
- RM D          [Central Tendency] [central tendency ...] [2008-12-18 15:43:24] [631938996a408f8d8cf3d9850ca0cd03]
- RM                [Mean Plot] [mean plot tijdree...] [2008-12-18 16:03:30] [4e8974eee929007194de34cbeefcb780] [Current]
-    D                [Mean Plot] [mean plot tijdree...] [2008-12-18 16:10:04] [631938996a408f8d8cf3d9850ca0cd03]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
12460
12870
14122
15366
15366
15046
10356
10194
11027
11449
12107
12862
12587
13951
15210
17511
17247
16900
12237
12188
12801
13276
14171
14629
16242
17239
18516
19862
19542
18515
13182
13148
13853
14341
14895
15317
15649
16373
17543
18911
19401
18295
12665
12246
12969
13563
15714
16094
15823
15460
16827
18279
18789
18041
12129
11708
12387
12982
13664
14073
14222
14715
15697
16908
17412
16232
10281
10258
11013
11779
12333
12684
12654
13239
14729
16146
16324
13797
9399
9626
10290
10931
11382
11704
11716
12470
13662
15097
15821
13083
9081
9479
10295
10999
11643
12055
12431
12940
14168
16031
16445
14245
9820
10389
11298
12156
12759
13414
13683
14070
15176
16694
17694
14294
9255
9868
10878
11685
12496
12488

Tables (Output of Computation)





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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=34867&T=0

As an alternative you can also use a QR Code:  
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 time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 60 ; par2 = 1 ; par3 = 1 ; par4 = 0 ; par5 = 12 ;
Parameters (R input):
par1 = 12 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
(n <- length(x))
(np <- floor(n / par1))
arr <- array(NA,dim=c(par1,np+1))
darr <- array(NA,dim=c(par1,np+1))
ari <- array(0,dim=par1)
dx <- diff(x)
j <- 0
for (i in 1:n)
{
j = j + 1
ari[j] = ari[j] + 1
arr[j,ari[j]] <- x[i]
darr[j,ari[j]] <- dx[i]
if (j == par1) j = 0
}
ari
arr
darr
arr.mean <- array(NA,dim=par1)
arr.median <- array(NA,dim=par1)
arr.midrange <- array(NA,dim=par1)
for (j in 1:par1)
{
arr.mean[j] <- mean(arr[j,],na.rm=TRUE)
arr.median[j] <- median(arr[j,],na.rm=TRUE)
arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2
}
overall.mean <- mean(x)
overall.median <- median(x)
overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2
bitmap(file='plot1.png')
plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index')
mtext(paste('#blocks = ',np))
abline(overall.mean,0)
dev.off()
bitmap(file='plot2.png')
plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index')
mtext(paste('#blocks = ',np))
abline(overall.median,0)
dev.off()
bitmap(file='plot3.png')
plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index')
mtext(paste('#blocks = ',np))
abline(overall.midrange,0)
dev.off()
bitmap(file='plot4.png')
z <- data.frame(t(arr))
names(z) <- c(1:par1)
(boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries'))
dev.off()
bitmap(file='plot4b.png')
z <- data.frame(t(darr))
names(z) <- c(1:par1)
(boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries'))
dev.off()
bitmap(file='plot5.png')
z <- data.frame(arr)
names(z) <- c(1:np)
(boxplot(z,notch=TRUE,col='grey',xlab='Block Index',ylab='Value',main='Notched Box Plots - Sequential Blocks'))
dev.off()
bitmap(file='plot6.png')
z <- data.frame(cbind(arr.mean,arr.median,arr.midrange))
names(z) <- list('mean','median','midrange')
(boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots'))
dev.off()