## Free Statistics

of Irreproducible Research!

Author's title
Author*Unverified author*
R Software Modulerwasp_meanplot.wasp
Title produced by softwareMean Plot
Date of computationThu, 27 Dec 2012 14:36:32 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Dec/27/t1356637021qieaonepus4hus4.htm/, Retrieved Tue, 07 Feb 2023 19:15:49 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=204783, Retrieved Tue, 07 Feb 2023 19:15:49 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact109
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Mean Plot] [Mean - Median Plo...] [2012-03-09 13:42:04] [cb8210d8f054d565da1456a9f4ea534d]
- R  D    [Mean Plot] [Mean - Median Plo...] [2012-12-27 19:36:32] [97ff841fcf87514e420f2e9629cfd808] [Current]
Feedback Forum

Post a new message
Dataseries X:
0,75
0,75
0,77
0,78
0,79
1,01
1,16
1,14
1,12
1,1
1,1
1,1
1,1
1,09
1,09
1,1
1,1
1,17
1,15
1,04
0,94
0,88
0,85
0,85
0,85
0,84
0,83
0,8
0,78
1,02
1,19
1,1
0,96
0,87
0,83
0,82
0,81
0,78
0,79
0,8
0,79
0,97
1,01
0,92
0,87
0,84
0,81
0,81
0,83
0,83
0,85
0,88
0,89
1,21
1,32
1,33
1,23
1,16
1,12
1,06
1,08
1,09
1,03
1,04
1,05
1,19
1,14
1,05
0,95
0,87
0,86
0,85

 Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 4 seconds R Server 'Gertrude Mary Cox' @ cox.wessa.net

\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 & 4 seconds \tabularnewline
R Server & 'Gertrude Mary Cox' @ cox.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=204783&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]4 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gertrude Mary Cox' @ cox.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=204783&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=204783&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 4 seconds R Server 'Gertrude Mary Cox' @ cox.wessa.net

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