FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_meanplot.wasp
Title produced by softwareMean Plot
Date of computationThu, 06 Nov 2008 06:33:18 -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/Nov/06/t1225978529jhan4g5oi280mo1.htm/, Retrieved Sun, 19 May 2024 07:20:30 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=22144, Retrieved Sun, 19 May 2024 07:20:30 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F       [Mean Plot] [EDA (part 1) Task...] [2008-11-06 13:33:18] [51254d789fff0741e6503951f574c682] [Current]
Feedback Forum
2008-11-11 09:34:59 [Lana Van Wesemael] [reply] 
Q2
Met het tweede gedeelte van de conclusie van de student ben ik het niet eens. Het is niet het geval dat er (bij de notched box plots)elke maand een significant verschil is. Bij sommige maanden overlappen de notches elkaar bvb tussen de 7e en de 8e maand is er geen sprake van overlapping, het verschil kan dus aan het toeval te wijten zijn. Wel is het zo dat er tussen de grootste schommelingen geen overlapping van de notches voorkomt bvb tussen de 6e en de 7e maand. Hier is dus wel sprake van een significant verschil. Dit duidt aan seizoenaliteit aanwezig is.

Q3
Hier ben ik het niet helemaal eens met de student. Het klopt inderdaad dat de notches van de box plots elkaar overlappen. Maar in dit geval moeten we kijken naar elke box plot en het basisjaar. Hier kunnen we zien dat de mediaan telkens lager ligt dan het basisjar (wat gelijk is aan 100). Of dit dan een significante daling is kunnen we zien aan de notches. Het eerste jaar is er nog sprake van overlapping van de notches met het basisjaar. Maar de volgende jaren liggen de notches significant lager ten opzichte van het basisjaar.
2008-11-12 09:28:53 [Bert Moons] [reply] 
Q2: Er zou duidelijk gemaakt moeten worden dan de getallen op de x-as maanden voorstellen. En dat elk punt in de grafiek het gemiddelde is van de productiecijfers van die maanden. Ook zou er uitgelegd moeten worden waarom er sprake is van een significante daling(noches overlappen elkaar weinig).
De conclusie is naar mijn mening wel ok.
2008-11-12 09:29:31 [Bert Moons] [reply] 
Q3: Er is inderdaad steeds overlapping van de inkepingen. De inkeping van de laatste “noched box plot” en de eerste overlappen elkaar echter (bijna) niet wat dit tot een potentieel twijfelgeval maakt.
2008-11-12 10:00:13 [Stef Vermeiren] [reply] 
Ik ben het niet volledig eens met de student. het is niet zo dat er voor alle maanden een significant verschil is. bijvoorbeeld tussen maand 6 en 7 is er een significant verschil waarbij we spreken van seizoenaliteit.
2008-11-12 10:04:06 [Stef Vermeiren] [reply] 
q3:De inkepingen overlappen inderdaag. De waarden zijn dus niet significant ook al verschillen de medianen wel.

Jaar 1 en 5 kan je misschien wel vergelijken omdat de notches niet meer overlappen. Maar hier is verder onderzoek en analyse voor nodig.

In het totale opzicht is er een dalende trend met overlappende notches, dus is het meer aan het toeval te wijten.
2008-11-12 10:13:38 [407693b66d7f2e0b350979005057872d] [reply] 
Dit is niet volledig correct geantwoord. Op de boxplots kunnen we duidelijker zien of er sprake is van seizonaliteit dan op de mean plot maar dar kunnen we wel grote verschillen detecteren. We zien dat in de zevende maand de mediaan enorm veel hoger ligt dan de vijf voorgaande maanden(inkepingen overlappen elkaar niet dus mediaan zevende maand is niet toevallig) dus m.a.w. hier is er sprake van seizonaliteit
2008-11-12 10:14:07 [407693b66d7f2e0b350979005057872d] [reply] 
De medianen dalen jaar op jaar maar de inkepingen overlappen elkaar wel dus dit betekent dat het toeval kan zijn of de medianen lager liggen . Jaar 6 moeten we hier geen rekening mee houden omdat er hier te weinig metingenvoor zijn

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
109,20
88,60
94,30
98,30
86,40
80,60
104,10
108,20
93,40
71,90
94,10
94,90
96,40
91,10
84,40
86,40
88,00
75,10
109,70
103,00
82,10
68,00
96,40
94,30
90,00
88,00
76,10
82,50
81,40
66,50
97,20
94,10
80,70
70,50
87,80
89,50
99,60
84,20
75,10
92,00
80,80
73,10
99,80
90,00
83,10
72,40
78,80
87,30
91,00
80,10
73,60
86,40
74,50
71,20
92,40
81,50
85,30
69,90
84,20
90,70
100,30

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=22144&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=22144&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=22144&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 = 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()