## Free Statistics

of Irreproducible Research!

Author's title
Author*The author of this computation has been verified*
R Software Modulerwasp_meanplot.wasp
Title produced by softwareMean Plot
Date of computationMon, 16 Nov 2009 12:37:47 -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/2009/Nov/16/t12584004554t1nxj0x2579i9q.htm/, Retrieved Fri, 02 Dec 2022 10:23:51 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=57354, Retrieved Fri, 02 Dec 2022 10:23:51 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsws6 mean en median plot
Estimated Impact156
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Mean Plot] [3/11/2009] [2009-11-02 22:07:54] [b98453cac15ba1066b407e146608df68]
-   PD  [Mean Plot] [ws 6 mean plot] [2009-11-10 21:19:13] [616e2df490b611f6cb7080068870ecbd]
-   PD      [Mean Plot] [ws6 mean en media...] [2009-11-16 19:37:47] [88e98f4c87ea17c4967db8279bda8533] [Current]
Feedback Forum

Post a new message
Dataseries X:
4.8
4.8
4.7
4.7
4.7
4.6
5.0
5.4
5.5
5.6
5.6
5.8
6.0
6.1
6.1
6.0
6.0
6.1
6.5
7.1
7.4
7.4
7.5
7.6
7.8
7.8
7.7
7.6
7.5
7.3
7.6
8.0
8.8
7.9
7.8
7.7
7.8
7.7
7.5
7.3
7.1
7.0
7.3
7.8
7.9
7.9
7.8
7.8
7.9
7.8
7.6
7.4
7.2
6.9
7.1
7.5
7.6
7.4
7.3
7.2
7.3
7.2
7.1
7.0
6.9
6.8
7.2
7.6
7.7
7.6
7.5
7.5
7.6
7.6
7.6
7.5
7.3
7.2
7.4
8.0
8.2
8.0
7.7
7.7
7.8
7.8
7.7
7.5
7.3
7.1
7.1
7.2
6.8
6.6
6.4
6.4
6.5
6.3
5.9
5.5
5.2
4.9
5.4
5.8
5.7
5.6
5.5
5.4
5.4
5.4
5.5
5.8
5.7
5.4
5.6
5.8
6.2
6.8
6.7
6.7
6.4
6.3
6.3
6.4
6.3
6.0
6.3
6.3
6.6
7.5
7.8
7.9
7.8
7.6
7.5
7.6
7.5
7.3
7.6
7.5
7.6
7.9
7.9
8.1
8.2
8.0
7.5
6.8
6.5
6.6
7.6
8.0
8.1
7.7
7.5
7.6
7.8
7.8
7.8
7.5
7.5
7.1
7.5
7.5
7.6
7.7
7.7
7.9
8.1
8.2
8.2
8.2
7.9
7.3
6.9
6.6
6.7
6.9
7.0
7.1

 Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 2 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=57354&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=57354&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=57354&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 Input view raw input (R code) Raw Output view raw output of R engine Computing time 2 seconds R Server 'Gwilym Jenkins' @ 72.249.127.135

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 <- 0for (i in 1:n){j = j + 1ari[j] = ari[j] + 1arr[j,ari[j]] <- x[i]darr[j,ari[j]] <- dx[i]if (j == par1) j = 0}ariarrdarrarr.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)) / 2bitmap(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()