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 computationSun, 02 Nov 2008 11:40:47 -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/02/t1225651305o12khgq6flttsdl.htm/, Retrieved Sat, 18 May 2024 18:53:27 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=20695, Retrieved Sat, 18 May 2024 18:53:27 +0000
QR Codes:

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

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] [Hypothesis testing] [2008-11-02 18:40:47] [d6e9f26c3644bfc30f06303d9993b878] [Current]
Feedback Forum
2008-11-09 14:19:02 [Lindsay Heyndrickx] [reply] 
Hier zien we inderdaad een duidelijke stijging op het einde van de meanplot maar dit is niet significant. Als we naar de box plot periodic subseries dan zien we dat de betrouwbaarheidsintervallen elkaar nog steeds overlappen. Hier is dus geen significant verschil.
2008-11-10 01:22:25 [Kristof Augustyns] [reply] 
Het is inderdaad zo dat er een stijging is van de energie vanaf september tot december, maar die is niet significant.
Er is namelijk een overlapping van de mediaan (inkeping) tussen de maand september en october, november, december.
Misschien is het tussen september en december een twijfelgeval en kan het wel een significant verschil zijn.
Ook zie je in de 'sequential blocks' dat er een permanente daling is van de energie in de eerste vier jaren. Vanaf dan blijft het constant.
Enkele uitschieters in jaar drie zijn niet onoverkomelijk.
Het is wel duidelijk dat de energie in de zomermaanden minder is dan in de andere maanden en dit omdat er toch op het nippertje seizonaliteit van kracht is.
Student heeft hier geen volledig antwoord op gegeven, maar wat de persoon heeft gegeven is wel juist, maar weeral vrij gering.
2008-11-10 15:44:24 [Jens Peeters] [reply] 
De stijging is inderdaad niet significant.Tussen september en december overlappen de betrouwbaarheidsintervallen elkaar. Bij september en december kan je inderdaad nog twijfelen zoals kristof al aanhaalde.
2008-11-12 11:26:24 [df2ed12c9b09685cd516719b004050c5] [reply] 
De daling tussen de eerste en de 4de boxplot is significant (inkepingen van de boxplots overlappen elkaar niet). De stijging op het laatste is echter niet significant, want er is wel een overlapping.

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
100.3
97.6
89.1
99.1
94.9
96.5
92.6
80.8
89.5
101.4
95.9
92.3
91.2
88.3
80.7
89.9
87.2
86.9
82.8
72.6
81.3
91.2
87.3
83.4
81.7
80.2
74.1
80.6
79.0
79.3
71.2
78.1
68.2
81.0
106.9
123.7
73.7
69.2
72.5
75.7
73.5
70.4
65.7
68.1
62.4
64.7
77.7
85.9
61.0
57.4
75.1
75.9
71.8
72.3
67.3
71.5
67.6
74.2
77.6
76.4
74.2

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001

\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 & 4 seconds \tabularnewline
R Server & 'Herman Ole Andreas Wold' @ 193.190.124.10:1001 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=20695&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]4 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Herman Ole Andreas Wold' @ 193.190.124.10:1001[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=20695&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=20695&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 time4 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001


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()