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 computationThu, 30 Oct 2008 13:54:10 -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/30/t1225396718rz0uorqv9gr4fpn.htm/, Retrieved Tue, 14 May 2024 01:25:18 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=20165, Retrieved Tue, 14 May 2024 01:25:18 +0000
QR Codes:

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

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: test the hypo...] [2008-10-28 20:26:19] [1e1d8320a8a1170c475bf6e4ce119de6]
F   PD      [Mean Plot] [Q2 Aanpassing] [2008-10-30 19:54:10] [8a1195ff8db4df756ce44b463a631c76] [Current]
-   P         [Mean Plot] [Correctie taak 5 Q2] [2008-11-07 17:11:07] [aa5573c1db401b164e448aef050955a1]
Feedback Forum
2008-11-08 10:19:57 [Stéphanie Claes] [reply] 
De student heeft de blockwidth aangepast, waardoor je de seizonaliteit niet goed meer kan zien (in dit geval valt dat nog mee). Je moet dit op 12 laten staan zodat de calculator de waarden van de overeenstemende maanden kan vergelijken. Op de X-as staan weldegelijk de maanden.
De maand bij 1 is maart (voor alle jaren) en daar wordt een gemiddelde voor berekend = periodiek gemiddelde. We zien dat het gemiddelde van augustus helemaal anders is dan dat van september (voor mean, median en midrange).
Als we naar de boxplot (periodic subseries) gaan kijken dan zien we ook hier dat het verschil tussen augustus en september (dit kan verklaard worden door het feit dat de zomer bijna voorbij is en mensen geen zomerkledij meer gaan kopen) significant is => seizonaliteit.

http://www.freestatistics.org/blog/index.php?v=date/2008/Nov/08/t122613876548txjhg0ufd7a35.htm => geen blockwidth aanpassing
  2008-11-08 14:33:35 [Stéphanie Claes] [reply] 
Even snel een correctie, de student heeft gelijk wanneer hij erop wijst dat we op de x as de groepen van data vinden en dus niet de maanden.
2008-11-08 10:35:06 [Stéphanie Claes] [reply] 
Q3: Opnieuw de opmerking dat je dus ook voor deze vraag de blockwidth niet moest aanpassen. We gaan dus jaar per jaar vergelijken (kijken naar de notched boxplot, sequential blocks). Bij kledingproductie gaat de mediaan dalen tegenover het basisjaar. Bij 2 ligt het toevallig lager dan bij 1.
Het 5de jaar ten opzichte van het 1ste jaar is bijna significant (dit kan je op deze blog niet zien omdat de blockwidth werd aangepast). Dus we kunnen stellen dat de daling toch wel 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 computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 3 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=20165&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]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=20165&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=20165&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 time3 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 15 ;
Parameters (R input):
par1 = 15 ;
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()