FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

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, 03 Nov 2008 11:53:26 -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/03/t1225738511w50rcldv9qcznl6.htm/, Retrieved Sun, 19 May 2024 09:17:50 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=20974, Retrieved Sun, 19 May 2024 09:17:50 +0000
QR Codes:

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

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] [workshop 3] [2007-10-26 12:14:28] [e9ffc5de6f8a7be62f22b142b5b6b1a8]
F    D    [Mean Plot] [Q2] [2008-11-03 18:53:26] [07b7cf1321bc38017c2c7efcf91ca696] [Current]
Feedback Forum
2008-11-09 11:57:46 [Kristof Augustyns] [reply] 
Het is inderdaad juist wat er gezegd wordt.
Er zijn twee pieken om al de grafieken.
Die pieken zijn er in de zomermaanden (maand 7 en 8).
In juli en augustus ligt gemiddelde, mediaan en midrange een pak hoger dan in de andere maanden.
Dit is dus duidelijk een significant verschil door dat dit te maken heeft met seizonaliteit.
Men zal blijkbaar ieder jaar in die twee maanden iets speciaals doen om zo de kleding productie een duwtje in de rug te geven.
2008-11-09 13:55:53 [Tim Loyens] [reply] 
Op de mean plot kunnen we duidelijk zien dat er een verschil is tussen, bijvoorbeeld, maanden 6 en 7. Is dit verschil ook significant? Hiervoor moeten we kijken naar de notched boxplot. Net zoals al gedacht is het duidelijk dat in deze maanden de gemiddelden significant verschillen. Hier weer merkbaar aan de inkepingen met betrouwbaarheid 95% die ver van elkaar liggen.
2008-11-09 14:57:51 [Natascha Meeus] [reply] 
De conclussie is niet correct. Het gaat hier over een significant verschil tussen de maand augustus en september. Maand 1 is hier maart. Er is sprake van seizonaliteit.
2008-11-10 12:21:25 [Glenn De Maeyer] [reply] 
Elk punt op het mean plot toont ons het gemiddelde per maand over de duur van tijdreeks. Zo is het eerste puntje het gemiddelde van alle maanden maart.
Er is een grote stijging waar te nemen in maand 7 (aangezien de eerste maand maart is, is de 7e maand dus september). We nemen ook een sterke daling waar in de 10e maand.
We kunnen ook eens kijken naar de notched box plots periodic subseries. Hier zijn we ook dat de maanden 7 en 10 respectievelijk hoger en lager liggen dan de rest. Maar om na te gaan of de verschillen significant zijn dienen we te kijken naar de betrouwbaarheidsintervallen (hier de notches).
De notches van maand 7 bijvoorbeeld overlappen enkel met die van maanden 1 en 8. Dus het verschil in mediaan bij maand 7 t.o.v. maanden 1 en 8 is niet significant. T.o.v. de andere maanden wel.
Als we dan ook eens kijken naar maand 10 dan merken we dat er hier een significant verschil is. Significant omdat de notches van de maand 10 met geen enkele andere maand overlappen.
2008-11-11 15:51:21 [Bart Haemels] [reply] 
Ik trek een correcte conclusie. Enkel heb ik de reeste maand als januari genomen, wat fout was aangezien dit de maand Maart is.
Ik kon ook op de grafiek van de notched box plots zien dat de mediaan van 7 &8 significant honger ligt dan de 5 voorgaande en de 2 maanden achteraf.
2008-11-11 22:26:24 [Liese Tormans] [reply] 
De student heeft de juiste techniek gebruikt ook is zijn conclusie goed maar een beetje onvolledig.

Aan de hand van deze grafiek gaan we kijken of we kunnen spreken van seizoenaliteit.
Als we naar de grafiek zien, zien we een uitschieter rond de maand 7 en maand 10. Zoals de student al gezegd heeft kunnen we spreken van seizoenaliteit.
De student had misschien nog iets duidelijker kunnen zijn, door te onderzoeken of het verschil toevallig of significant is. Dit kunnen we aflezen aan de notched boxplot. We kunnen zien dat in maand 7 en maand 10 de notches sterk afwijken van een denkbeeldige horizontale door het midden. De andere notches liggen dus op ongeveer dezelfde lijn terwijl maand 7 en 10 er sterk van afwijken. We kunnen dus spreken van een significant verschil.

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=20974&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 ; par2 = ; par3 = ; par4 = ; par5 = ; par6 = ; par7 = ; par8 = ; par9 = ; par10 = ; par11 = ; par12 = ; par13 = ; par14 = ; par15 = ; par16 = ; par17 = ; par18 = ; par19 = ; par20 = ;
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))
ari <- array(0,dim=par1)
j <- 0
for (i in 1:n)
{
j = j + 1
ari[j] = ari[j] + 1
arr[j,ari[j]] <- x[i]
if (j == par1) j = 0
}
ari
arr
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='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()