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 computationWed, 29 Oct 2008 07:47:56 -0600
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/Oct/29/t1225288202a87ncce1gb5kpcw.htm/, Retrieved Tue, 14 May 2024 18:42:14 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=19851, Retrieved Tue, 14 May 2024 18:42:14 +0000
QR Codes:

Original text written by user:stéphanie claes kevin engels lindsay heyndrickx katrien bourdiaudhy
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact279

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   PD    [Mean Plot] [Q2 en Q3 mean plot] [2008-10-29 13:47:56] [95d95b0e883740fcbc85e18ec42dcafb] [Current]
F           [Mean Plot] [q2] [2008-11-03 17:14:50] [e43247bc0ab243a5af99ac7f55ba0b41]
F           [Mean Plot] [q3] [2008-11-03 17:16:25] [e43247bc0ab243a5af99ac7f55ba0b41]
-   P       [Mean Plot] [verbetering Q3] [2008-11-07 16:05:08] [2bd2ad6af3eef3a703e9ec23e39bd695]
-   P       [Mean Plot] [Q2 TAAK 3A] [2008-11-07 16:05:29] [c993f605b206b366f754f7f8c1fcc291]
-   P       [Mean Plot] [Q3 TAAK 3A] [2008-11-07 17:17:11] [c993f605b206b366f754f7f8c1fcc291]
Feedback Forum
2008-11-07 16:05:54 [Mehmet Yilmaz] [reply] 
De student heeft de blockwidth veranderd van 12 naar 60. Hierdoor veranderd ze de perioden. Dit is fout. De blockwidth dient op 12 te blijven. De reproductie is dus verkeerd.

Hieronder vind je de juiste berekening:
http://www.freestatistics.org/blog/index.php?v=date/2008/Nov/07/t1226073928lgg8f29h1y02tac.htm
2008-11-07 16:27:02 [Kevin Engels] [reply] 
Zoals de student van vorig jaar gedaan heeft, is het beter om de blockwidth op 12 te laten staan en niet op 60 te zetten zoals wij gedaan hebben. Via deze link vind je dus de juist boxplot terug:
http://www.freestatistics.org/blog/index.php?v=date/2008/Nov/07/t1226074065ysaj29xr27ieqmz.htm

Hierop vertegenwoordigt het eerste punt het gemiddelde van alle gegevens van de maand maart in 2002, 2003,... het volgende punt dat gemiddelde van van april 2002, 2003,... Het geeft dus een periodiek gemiddelde weer.

We hadden gelijk daar we zeiden dat de seizonaliteit wel degelijk een effect heeft op de productie. Op de link van hierboven vinden we ook de notched boxplot - periodic subseries terug. De notches zijn hier zeer groot en overlappen elkaar soms. We zien enkele significante verschillen: bijvoorbeeld tussen de maanden augustus en september (punt 6 en 7) wat duidt op seizonaliteit.
2008-11-07 17:26:51 [Kevin Engels] [reply] 
Wat Q3 betreft, hebben we hier de foutieve plot gebruikt. Ipv de mean plot kan men beter zoals de student van vorig jaar de 'notched boxplot - sequential blocks' gebruiken om aan te tonen of de kledingproductie gedaald is de voorbije jaren. We gaan deze dus jaar per jaar bekijken.

Net zoals bij Q1 moet men dan de blockwidth op 12 zetten en niet op 60, we krijgen dan deze plots:
http://www.freestatistics.org/blog/index.php?v=date/2008/Nov/07/t1226078502ur1y0fwtm9r8m26.htm

De mediaan daalt duidelijk. Is het significant lager? De eerste 4 notches overlappen elkaar, maar de 5e is een twijfelgeval. Als blok 5 eruit zou vallen, is het niet meer significant.
2008-11-10 11:55:08 [a7e076854c32462fd499d2de3f6d4e86] [reply] 
De blockwidth werd veranderd van 12 naar 60 wat foutief is, hierdoor worden de perioden veranderd. De blockwidth moet op 12 blijven staan. De reproductie is dus verkeerd.
Zie bovenstaande link voor de correcte berekening.

2008-11-11 13:05:49 [Bonifer Spillemaeckers] [reply] 
Deze vraag werd foutief opgelost. Men had hier de blockwidth op 12 moeten laten staan en niet moeten veranderen naar 60, want hierdoor verandert men de perioden.

Het mean plot berekent inderdaad het periodiek gemiddelde van een bepaalde maand over x aantal jaren (bv gemiddelde maart 2000, 20001,...). We kunnen duidelijk vaststellen dat het gemiddelde van de maand augustus anders is dan het gemiddelde van de maand september, maar aan de hand van deze grafiek kunnen we nog niet spreken van een beduidend verschil. We gaan hier ook kijken naar de grafiek van de notched boxplots. Hier kunnen we een groot verschil bemerken tussen de boxplot van de maand maart en deze van de maanden augustus en september. Hier kunnen we dus wel spreken van een beduidend verschil. Hier is dus sprake van seizoenaliteit.

2008-11-11 13:12:19 [Bonifer Spillemaeckers] [reply] 
Q3

De student maakt hier gebruik van het mean plot, maar om het verloop van de dataset over de jaren heen te bekijken, nemen we de grafiek van de sequential blocks erbij. We zien dat de mediaan hier telkens daalt. We kunnen ook zien dat de inkepingen van de boxplots (=betrouwbaarheidsintervallen) elkaar overlappen. We kunnen dit dus toeschrijven aan het toeval. We kunnen dus stellen dat de mediaan steeds daalt, maar dit is niet significant. Enkel over het 1ste jaar tov het 5de jaar kan er twijfel bestaan.
2008-11-11 15:03:11 [Bart Haemels] [reply] 
Je Had de Blockwidth op 12 moeten laten staan. Want elk punt is een gemiddele van die maand over X aantal jaar. AAngezien een jaar 12 maanden heeft, is een blockwidth van 12 dus vanzelfsprekend.

In volgende link:
http://www.freestatistics.org/blog/index.php?v=date/2008/Nov/03/t1225738511w50rcldv9qcznl6.htm
Kun je het resultaat zien zoals het moest zijn. Je kan 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 15:08:19 [Bart Haemels] [reply] 
Q3: Je had beter nu de sequential blocks gebruikt. Hieruit kun je alfeiden dat er een daling is van de mediaan, maar dit is volledig te wijten aan het toeval want alle notches van de 5 jaar overlappen mekaar, zelfs die van het 1ste & het 5de jaar.
2008-11-12 11:25:22 [Nilay Erdogdu] [reply] 
De student heeft de blockwidth veranderd van 12 naar 60. Hierdoor veranderd ze de perioden. Dit is fout. De blockwidth dient op 12 te blijven. De reproductie is dus verkeerd.

Hieronder vind je de juiste berekening:
http://www.freestatistics.org/blog/index.php?v=date/2008/Nov/07/t1226073928lgg8f29h1y02tac.htm

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=19851&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 = 60 ;
Parameters (R input):
par1 = 60 ;
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()