Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_meanplot.wasp
Title produced by softwareMean Plot
Date of computationMon, 26 Apr 2010 21:57:04 +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/2010/Apr/26/t12723190683ranzaz9g8h985l.htm/, Retrieved Thu, 28 Mar 2024 14:36:31 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=74922, Retrieved Thu, 28 Mar 2024 14:36:31 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKDGP2W12
Estimated Impact173
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Histogram] [Datareeks-Inschri...] [2010-02-08 14:03:15] [c656187b38b4f6e6948d94f8dfe6ded2]
- RMP   [Harrell-Davis Quantiles] [Datareeks-Opdrach...] [2010-03-08 17:35:01] [1f3241a8f2363a866734862cbbf73252]
- RMPD      [Mean Plot] [Opgave 6 oefening...] [2010-04-26 21:57:04] [8c87877ca0a068b5d9f0f8fa9cf6c0e7] [Current]
-    D        [Mean Plot] [Opgave 6 oefening...] [2010-04-28 15:35:23] [1f3241a8f2363a866734862cbbf73252]
- RMPD        [(Partial) Autocorrelation Function] [Opgave 6 BIS oefe...] [2010-05-03 18:06:05] [1f3241a8f2363a866734862cbbf73252]
-   P           [(Partial) Autocorrelation Function] [Opgave 6 BIS oefe...] [2010-05-03 18:08:48] [1f3241a8f2363a866734862cbbf73252]
-   P           [(Partial) Autocorrelation Function] [Opgave 6 BIS oefe...] [2010-05-03 18:11:43] [1f3241a8f2363a866734862cbbf73252]
-   P           [(Partial) Autocorrelation Function] [Opgave 6 BIS oefe...] [2010-05-03 18:13:08] [1f3241a8f2363a866734862cbbf73252]
Feedback Forum

Post a new message
Dataseries X:
14538
18730
22485
20036
16971
19028
22759
20516
26195
27786
24090
25447
11509
15572
22518
20520
17789
20205
26835
25826
31934
30019
30111
31566
12738
19814
24776
20424
18688
20418
25778
25100
25859
30651
26551
31124
9367
17382
20995
18205
17328
18157
23691
26736
27165
34506
29506
31956
10698
18479
19785
19684
18730
17970
27044
22405
26482
29096
25591
29743
13807
19169
22782
20366
17537
18004
24319
22679
32034
36438
29383
28029
15548
17704
25316
20764
18089
21705
23843
22397
26105
29462
27071
31514
11514
17409
24561
21382
17351
22780
23970
26887
34777
35022
35338
36845
13971
22228
28806
20506
22414
25814
28352
29965
30212
32609
30364
37702
13253
21790
27192
21725
22205
24693
29133
35953
37863
43129
39690
41086




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 3 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=74922&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]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=74922&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=74922&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 time3 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24



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