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 computationFri, 26 Oct 2007 05:14:28 -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/2007/Oct/26/kksisru32kp9sfi1193400558.htm/, Retrieved Sat, 27 Apr 2024 12:13:28 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=1870, Retrieved Sat, 27 Apr 2024 12:13:28 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsmean plot
Estimated Impact627

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] [d06427f3e67cec1f6334fc93f511b0b4] [Current]
- R  D    [Mean Plot] [task 4] [2008-10-28 19:54:33] [c5a66f1c8528a963efc2b82a8519f117]
F    D    [Mean Plot] [Q2: test the hypo...] [2008-10-28 20:26:19] [1e1d8320a8a1170c475bf6e4ce119de6]
F R         [Mean Plot] [task 4 Q2 if trim...] [2008-10-29 13:18:38] [1e1d8320a8a1170c475bf6e4ce119de6]
F   PD        [Mean Plot] [Task 4 Q2] [2008-10-31 00:34:03] [aa5573c1db401b164e448aef050955a1]
F               [Mean Plot] [Task 4 Q2 Blockwi...] [2008-10-31 00:40:24] [aa5573c1db401b164e448aef050955a1]
F   PD      [Mean Plot] [Q2 Aanpassing] [2008-10-30 19:54:10] [aa5573c1db401b164e448aef050955a1]
-   P         [Mean Plot] [Correctie taak 5 Q2] [2008-11-07 17:11:07] [aa5573c1db401b164e448aef050955a1]
F    D    [Mean Plot] [Q3: test the hypo...] [2008-10-28 20:35:16] [1e1d8320a8a1170c475bf6e4ce119de6]
F R         [Mean Plot] [task 4 Q3 if trim...] [2008-10-29 13:21:15] [1e1d8320a8a1170c475bf6e4ce119de6]
F   PD    [Mean Plot] [Mean Plot] [2008-10-29 11:59:16] [87cabf13a90315c7085b765dcebb7412]
F           [Mean Plot] [Hypothesis testin...] [2008-11-03 19:16:00] [631938996a408f8d8cf3d9850ca0cd03]
F   PD    [Mean Plot] [Mean Plot] [2008-10-29 11:59:16] [58bf45a666dc5198906262e8815a9722]
F           [Mean Plot] [] [2008-11-03 19:13:04] [d2d412c7f4d35ffbf5ee5ee89db327d4]
-   P       [Mean Plot] [Q2 mean plot analyse] [2008-11-07 13:09:32] [3754dd41128068acfc463ebbabce5a9c]
-   P       [Mean Plot] [verbetering works...] [2008-11-07 17:08:41] [cf9c64468d04c2c4dd548cc66b4e3677]
-   P       [Mean Plot] [Verbetering works...] [2008-11-07 17:14:40] [cf9c64468d04c2c4dd548cc66b4e3677]
F    D    [Mean Plot] [Herproducering ta...] [2008-10-29 12:00:40] [819b576fab25b35cfda70f80599828ec]
- R         [Mean Plot] [Mean plot Task 4] [2008-10-29 14:01:28] [819b576fab25b35cfda70f80599828ec]
F R         [Mean Plot] [Mean plot task 4] [2008-10-29 14:43:42] [819b576fab25b35cfda70f80599828ec]
F    D      [Mean Plot] [Mean plot deel 1 ...] [2008-10-29 15:16:12] [819b576fab25b35cfda70f80599828ec]
F   PD    [Mean Plot] [Mean plot product...] [2008-10-29 12:03:27] [b635de6fc42b001d22cbe6e730fec936]
F           [Mean Plot] [Taak 1 Q2] [2008-10-29 15:57:22] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
F           [Mean Plot] [Q2 mean plot] [2008-10-30 08:57:32] [fe7291e888d31b8c4db0b24d6c0f75c6]
-   PD    [Mean Plot] [Mean plot clothin...] [2008-10-29 12:22:58] [b635de6fc42b001d22cbe6e730fec936]
F           [Mean Plot] [Taak 1, Q3] [2008-10-29 16:02:43] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
F           [Mean Plot] [Q3 mean plot] [2008-10-30 09:00:03] [fe7291e888d31b8c4db0b24d6c0f75c6]
-    D    [Mean Plot] [Q2 mean plot] [2008-10-29 13:13:51] [7173087adebe3e3a714c80ea2417b3eb]
-   PD    [Mean Plot] [Q2 mean plot] [2008-10-29 13:19:08] [7173087adebe3e3a714c80ea2417b3eb]
F   PD    [Mean Plot] [Q2 en Q3 mean plot] [2008-10-29 13:47:56] [7173087adebe3e3a714c80ea2417b3eb]
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]
F RMPD    [Notched Boxplots] [Notched boxplot a...] [2008-10-29 13:43:18] [b635de6fc42b001d22cbe6e730fec936]
F R         [Notched Boxplots] [Notched Boxplot l...] [2008-10-29 14:23:11] [b635de6fc42b001d22cbe6e730fec936]
F   P         [Notched Boxplots] [Taak 3, Notched b...] [2008-10-29 16:16:39] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
-               [Notched Boxplots] [Hercomputatie fee...] [2008-11-07 17:23:56] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
F   P         [Notched Boxplots] [taak 3 ] [2008-10-30 09:08:28] [fe7291e888d31b8c4db0b24d6c0f75c6]
F   P       [Notched Boxplots] [Taak 2 Notched bo...] [2008-10-29 16:11:29] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
-             [Notched Boxplots] [Hercomputatie fee...] [2008-11-07 17:14:35] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
-               [Notched Boxplots] [Verbetering] [2008-11-08 18:35:34] [79c17183721a40a589db5f9f561947d8]
F   P       [Notched Boxplots] [taak2 notched box...] [2008-10-30 09:06:42] [fe7291e888d31b8c4db0b24d6c0f75c6]
F    D    [Mean Plot] [Q2 Mean plot] [2008-10-29 13:54:33] [8e4e5f204c24e6d05647858dae308d17]
F   P       [Mean Plot] [q3 median plot] [2008-10-29 14:05:11] [8e4e5f204c24e6d05647858dae308d17]
-   P         [Mean Plot] [] [2008-11-07 10:18:13] [f5709eefd05c649ca6dad46019ffd879]
F   PD    [Mean Plot] [seasonality on cl...] [2008-10-29 14:26:50] [495cd80c1a9baafb03c09cd9ab8d8fb5]
- RMPD      [Harrell-Davis Quantiles] [Task 4: Harrell-D...] [2008-10-29 17:06:21] [495cd80c1a9baafb03c09cd9ab8d8fb5]

[Truncated]
Feedback Forum
2008-11-06 14:47:43 [Nathalie Koulouris] [reply] 
De student heeft deze vraag goed opgelost. Bij de Mean plot is het duidelijk visueel zichtbaar dat er seizoenaliteit aanwezig is. In maand 7 nemen we duidelijk een uitzonderlijk hoge stijging waar terwijl we in maand 10 een uitzonderlijk sterke daling waarnemen. Er is dus sprake van een significant verschil en dit wijst op seizoenaliteit.
2008-11-07 13:04:46 [Siem Van Opstal] [reply] 
goede berekeningen en correcte, uitgebreide conclusie
2008-11-09 15:17:25 [Vincent Vanden Poel] [reply] 
Q2: Je argumentatie is correct maar je zou er op moeten wijzen dat, om op deze vraag te antwoorden, het beter is om te kijken naar de periodieke notched box plots. Zo kunnen we nagaan of deze seizoensverschillen ook wel degelijk significant zijn.
2008-11-09 15:18:42 [Vincent Vanden Poel] [reply] 
Q3: Deze interpretatie is niet volledig correct. Aan de grafiek is duidelijk te zien dat de mediaan over de jaren heen gedaald is. Maar om er zeker van te zijn of deze daling structureel is of louter door toeval bepaald moeten we kijken naar de notches. Bij de individuele vergelijking tussen de opeenvolgende periodes merken we dat de inkepingen elkaar overlappen. De daling is met andere woorden niet significant. Wanneer we echter kijken naar het geheel merken we dat er wel min of meer een significante daling is tussen het 1e en het 5e jaar.
2008-11-09 19:29:21 [006ad2c49b6a7c2ad6ab685cfc1dae56] [reply] 
Q2 is ook goed maar het is niet duidelijk wat je bedoelt met de waarden op de y-as die de mean plot kunnen beïnvloeden. Bij Q3 is er inderdaad een daling in de productie van kleding, elke mediaan is lager dan de vorige, maar er zijn geen significante verschillen omdat de inkepingen overlappen. Over het geheel kan je wel spreken van een significante dazling. Het zesde jaar is onvolledig dus dit kan nog niet worden beoordeeld.
2008-11-09 19:33:41 [006ad2c49b6a7c2ad6ab685cfc1dae56] [reply] 
Bij Q2 kan je nog zeggen dat je bij maart moet beginnen, je had ook een grafiek van notched box plots-periodic subseries kunnen toevoegen om alles nog meer te verduielijken. Op deze box plots kan je zien dat er significante verschillen zijn tussen de medianen.
2008-11-09 19:34:16 [006ad2c49b6a7c2ad6ab685cfc1dae56] [reply] 
Bij Q3 kan je nog zeggen dat de dalingen elk afzonderlijk niet significant zijn maar de daling van het geheel van de periodes is wel significant.
2008-11-10 11:30:37 [a7e076854c32462fd499d2de3f6d4e86] [reply] 
De interpretatie is niet volledig. De mediaan is gedaald over de jaren heen.
Om na te gaan als deze daling structureel of toevallig is moeten we kijken naar de notches. De daling is niet significant (allemaal afzonderlijk). Over het geheel kan je wel spreken van een significante daling

2008-11-11 10:39:42 [Jeroen Michel] [reply] 
Belangrijk om weten is dat het hier gaat over rekenkundige gemiddelden en niet over waarden van de data! Bij 7, wat overeenkomt met maand 7, juli, is er een outlier waar te nemen wat wijst op seizoenaliteit. Net zoals beschreven in het document kan dit er op duiden dat mensen meer kledij zullen aankopen. Vanaf oktober is er een vergelijkbare trend, dit is dan wel van toepassing voor de winterkledij.

De grafieken geven een correct beeld weer.
2008-11-11 10:41:20 [Jeroen Michel] [reply] 
De grafieken zijn goed. Het is duidelijk dat we jaar 6 buiten beschouwing mogen laten aangezien deze data onvolledig zijn. Van groot belang is om te stellen dat er niet echt een significant verschil is, maar eerder sprake van toeval. Dit aangezien de alle jaren in het getrouwheidsinterval zitten.
2008-11-11 11:55:41 [Bernard Femont] [reply] 
DIT IS DE LINK MET DE BEREKENINGEN VAN DE STUDENT VAN VORIG JAAR. IN HET DOCUMENT BIJ CONCLUSIE VANONDER Q2 WERD ER EEN NIEUWE LINK AANGEMAAKT EN DAT IS ONZE LINK, GELIEVE DIE TE BESPREKEN
2008-11-11 14:58:17 [Loïque Verhasselt] [reply] 
Q2:Juiste grafiektype en conclusie maar de student vergeet haar conclusie te bevestigen met de notched box plots-periodic subseries. Deze geven aan dat we kunnen betekenis geven aan de grote verschillen door een significant verschil dat niet toe te schrijven is aan het toeval! Er is dus sprake van seizoenaliteit.
2008-11-11 21:42:55 [c233791e22ae82ed03fa45b0d63a2757] [reply] 
Spijtig dat er geen sequential analyse is. De interpretaties zijn wel correct.
2008-11-11 21:45:25 [c233791e22ae82ed03fa45b0d63a2757] [reply] 
q3: De conclusie klopt hier niet. Als je jaar 1 en jaar 5 met elkaar vergelijk zie je dat er wel een significante daling is.
2008-11-11 22:52:35 [5faab2fc6fb120339944528a32d48a04] [reply] 
De student geeft een conlusie maar geen verantwoording van de conclusie. Er is inderdaad een daling maar deze is enkel significant als we de eerste notched boxplot en de laatste bekijken. Bij de tussenliggende notched boxplots overlappen de betrouwbaarheidsintervallen waardoor de daling te wijten kan zijn aan toeval. En bijgevolg dan niet significant is.

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

[TABLE]
[ROW][C]Summary of compuational 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=1870&T=0

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