Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_meanplot.wasp
Title produced by softwareMean Plot
Date of computationSun, 15 Aug 2010 15:22:42 +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/Aug/15/t1281885743woxao1i1t15ec0d.htm/, Retrieved Sat, 27 Apr 2024 23:39:58 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=78890, Retrieved Sat, 27 Apr 2024 23:39:58 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact127
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Mean Plot] [Mean vs Median pl...] [2010-08-15 15:07:22] [f5ecd041e4b32af12787a4e421b18aaf]
-   P     [Mean Plot] [Mean&Median Plot ...] [2010-08-15 15:22:42] [05b8da000f2ebbd12b039a4b088dd3f2] [Current]
- RM        [(Partial) Autocorrelation Function] [Partial omzet pro...] [2010-08-15 15:56:49] [f5ecd041e4b32af12787a4e421b18aaf]
-   P         [(Partial) Autocorrelation Function] [Partial gedif. om...] [2010-08-16 13:23:29] [f5ecd041e4b32af12787a4e421b18aaf]
- RM        [Variability] [Variability omzet...] [2010-08-15 16:00:24] [f5ecd041e4b32af12787a4e421b18aaf]
- RM        [Standard Deviation-Mean Plot] [SD Mean Plot omze...] [2010-08-15 16:06:10] [f5ecd041e4b32af12787a4e421b18aaf]
- RM        [Classical Decomposition] [Classic Decomposi...] [2010-08-15 16:13:21] [f5ecd041e4b32af12787a4e421b18aaf]
- RM        [Classical Decomposition] [Classical Decompo...] [2010-08-15 16:13:21] [f5ecd041e4b32af12787a4e421b18aaf]
-             [Classical Decomposition] [Classical Decompo...] [2010-08-15 16:22:14] [f5ecd041e4b32af12787a4e421b18aaf]
- RM            [Exponential Smoothing] [Exponential smoot...] [2010-08-15 16:29:17] [f5ecd041e4b32af12787a4e421b18aaf]
Feedback Forum

Post a new message
Dataseries X:
118
117
116
114
112
111
112
114
115
115
116
118
126
131
122
124
119
112
109
108
117
122
127
124
129
141
127
133
114
98
93
101
111
128
126
134
140
158
144
146
138
119
113
120
127
141
144
150
156
174
163
167
160
141
132
144
155
164
162
181
187
209
189
201
193
177
159
158
155
164
163
185
191
217
193
192
184
166
145
146
138
149
145
166




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk

\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 & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=78890&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]'RServer@AstonUniversity' @ vre.aston.ac.uk[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=78890&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=78890&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'RServer@AstonUniversity' @ vre.aston.ac.uk



Parameters (Session):
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()