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_One Factor ANOVA.wasp
Title produced by softwareOne-Way-Between-Groups ANOVA- Free Statistics Software (Calculator)
Date of computationWed, 17 Dec 2014 18:41:07 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Dec/17/t1418841775ot7s8ct85yt6r8u.htm/, Retrieved Thu, 16 May 2024 17:03:11 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=270571, Retrieved Thu, 16 May 2024 17:03:11 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [One-Way-Between-Groups ANOVA- Free Statistics Software (Calculator)] [] [2014-12-17 18:41:07] [18673d63f90870b9c004059cd6229007] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
2011 149 18 68 7.5
2011 139 31 39 6
2011 148 39 32 6.5
2011 158 46 62 1
2011 128 31 33 1
2011 224 67 52 5.5
2011 159 35 62 8.5
2011 105 52 77 6.5
2011 159 77 76 4.5
2011 167 37 41 2
2011 165 32 48 5
2011 159 36 63 0.5
2011 119 38 30 5
2011 176 69 78 5
2011 54 21 19 2.5
2011 91 26 31 5
2011 163 54 66 5.5
2011 124 36 35 3.5
2011 137 42 42 3
2011 121 23 45 4
2011 153 34 21 0.5
2011 148 112 25 6.5
2011 221 35 44 4.5
2011 188 47 69 7.5
2011 149 47 54 5.5
2011 244 37 74 4
2011 148 109 80 7.5
2011 92 24 42 7
2011 150 20 61 4
2011 153 22 41 5.5
2011 94 23 46 2.5
2011 156 32 39 5.5
2011 132 30 34 3.5
2011 161 92 51 2.5
2011 105 43 42 4.5
2011 97 55 31 4.5
2011 151 16 39 4.5
2011 131 49 20 6
2011 166 71 49 2.5
2011 157 43 53 5
2011 111 29 31 0
2011 145 56 39 5
2011 162 46 54 6.5
2011 163 19 49 5
2011 59 23 34 6
2011 187 59 46 4.5
2011 109 30 55 5.5
2011 90 61 42 1
2011 105 7 50 7.5
2011 83 38 13 6
2011 116 32 37 5
2011 42 16 25 1
2011 148 19 30 5
2011 155 22 28 6.5
2011 125 48 45 7
2011 116 23 35 4.5
2011 128 26 28 0
2011 138 33 41 8.5
2011 49 9 6 3.5
2011 96 24 45 7.5
2011 164 34 73 3.5
2011 162 48 17 6
2011 99 18 40 1.5
2011 202 43 64 9
2011 186 33 37 3.5
2011 66 28 25 3.5
2011 183 71 65 4
2011 214 26 100 6.5
2011 188 67 28 7.5
2011 104 34 35 6
2011 177 80 56 5
2011 126 29 29 5.5
2011 76 16 43 3.5
2011 99 59 59 7.5
2011 139 32 50 6.5
2011 162 43 59 6.5
2011 108 38 27 6.5
2011 159 29 61 7
2011 74 36 28 3.5
2011 110 32 51 1.5
2011 96 35 35 4
2011 116 21 29 7.5
2011 87 29 48 4.5
2011 97 12 25 0
2011 127 37 44 3.5
2011 106 37 64 5.5
2011 80 47 32 5
2011 74 51 20 4.5
2011 91 32 28 2.5
2011 133 21 34 7.5
2011 74 13 31 7
2011 114 14 26 0
2011 140 -2 58 4.5
2011 95 20 23 3
2011 98 24 21 1.5
2011 121 11 21 3.5
2011 126 23 33 2.5
2011 98 24 16 5.5
2011 95 14 20 8
2011 110 52 37 1
2011 70 15 35 5
2011 102 23 33 4.5
2011 86 19 27 3
2011 130 35 41 3
2011 96 24 40 8
2011 102 39 35 2.5
2011 100 29 28 7
2011 94 13 32 0
2011 52 8 22 1
2011 98 18 44 3.5
2011 118 24 27 5.5
2011 99 19 17 5.5
2012 48 23 12 0.5
2012 50 16 45 7.5
2012 150 33 37 9
2012 154 32 37 9.5
2012 109 37 108 8.5
2012 68 14 10 7
2012 194 52 68 8
2012 158 75 72 10
2012 159 72 143 7
2012 67 15 9 8.5
2012 147 29 55 9
2012 39 13 17 9.5
2012 100 40 37 4
2012 111 19 27 6
2012 138 24 37 8
2012 101 121 58 5.5
2012 131 93 66 9.5
2012 101 36 21 7.5
2012 114 23 19 7
2012 165 85 78 7.5
2012 114 41 35 8
2012 111 46 48 7
2012 75 18 27 7
2012 82 35 43 6
2012 121 17 30 10
2012 32 4 25 2.5
2012 150 28 69 9
2012 117 44 72 8
2012 71 10 23 6
2012 165 38 13 8.5
2012 154 57 61 6
2012 126 23 43 9
2012 149 36 51 8
2012 145 22 67 9
2012 120 40 36 5.5
2012 109 31 44 7
2012 132 11 45 5.5
2012 172 38 34 9
2012 169 24 36 2
2012 114 37 72 8.5
2012 156 37 39 9
2012 172 22 43 8.5
2012 68 15 25 9
2012 89 2 56 7.5
2012 167 43 80 10
2012 113 31 40 9
2012 115 29 73 7.5
2012 78 45 34 6
2012 118 25 72 10.5
2012 87 4 42 8.5
2012 173 31 61 8
2012 2 -4 23 10
2012 162 66 74 10.5
2012 49 61 16 6.5
2012 122 32 66 9.5
2012 96 31 9 8.5
2012 100 39 41 7.5
2012 82 19 57 5
2012 100 31 48 8
2012 115 36 51 10
2012 141 42 53 7
2012 165 21 29 7.5
2012 165 21 29 7.5
2012 110 25 55 9.5
2012 118 32 54 6
2012 158 26 43 10
2012 146 28 51 7
2012 49 32 20 3
2012 90 41 79 6
2012 121 29 39 7
2012 155 33 61 10
2012 104 17 55 7
2012 147 13 30 3.5
2012 110 32 55 8
2012 108 30 22 10
2012 113 34 37 5.5
2012 115 59 2 6
2012 61 13 38 6.5
2012 60 23 27 6.5
2012 109 10 56 8.5
2012 68 5 25 4
2012 111 31 39 9.5
2012 77 19 33 8
2012 73 32 43 8.5
2012 151 30 57 5.5
2012 89 25 43 7
2012 78 48 23 9
2012 110 35 44 8
2012 220 67 54 10
2012 65 15 28 8
2012 141 22 36 6
2012 117 18 39 8
2012 122 33 16 5
2012 63 46 23 9
2012 44 24 40 4.5
2012 52 14 24 8.5
2012 131 12 78 9.5
2012 101 38 57 8.5
2012 42 12 37 7.5
2012 152 28 27 7.5
2012 107 41 61 5
2012 77 12 27 7
2012 154 31 69 8
2012 103 33 34 5.5
2012 96 34 44 8.5
2012 175 21 34 9.5
2012 57 20 39 7
2012 112 44 51 8
2012 143 52 34 8.5
2012 49 7 31 3.5
2012 110 29 13 6.5
2012 131 11 12 6.5
2012 167 26 51 10.5
2012 56 24 24 8.5
2012 137 7 19 8
2012 86 60 30 10
2012 121 13 81 10
2012 149 20 42 9.5
2012 168 52 22 9
2012 140 28 85 10
2012 88 25 27 7.5
2012 168 39 25 4.5
2012 94 9 22 4.5
2012 51 19 19 0.5
2012 48 13 14 6.5
2012 145 60 45 4.5
2012 66 19 45 5.5
2012 85 34 28 5
2012 109 14 51 6
2012 63 17 41 4
2012 102 45 31 8
2012 162 66 74 10.5
2012 86 48 19 6.5
2012 114 29 51 8
2012 164 -2 73 8.5
2012 119 51 24 5.5
2012 126 2 61 7
2012 132 24 23 5
2012 142 40 14 3.5
2012 83 20 54 5
2012 94 19 51 9
2012 81 16 62 8.5
2012 166 20 36 5
2012 110 40 59 9.5
2012 64 27 24 3
2012 93 25 26 1.5
2012 104 49 54 6
2012 105 39 39 0.5
2012 49 61 16 6.5
2012 88 19 36 7.5
2012 95 67 31 4.5
2012 102 45 31 8
2012 99 30 42 9
2012 63 8 39 7.5
2012 76 19 25 8.5
2012 109 52 31 7
2012 117 22 38 9.5
2012 57 17 31 6.5
2012 120 33 17 9.5
2012 73 34 22 6
2012 91 22 55 8
2012 108 30 62 9.5
2012 105 25 51 8
2012 117 38 30 8
2012 119 26 49 9
2012 31 13 16 5

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'George Udny Yule' @ yule.wessa.net

\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 & 'George Udny Yule' @ yule.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270571&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]'George Udny Yule' @ yule.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270571&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270571&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'George Udny Yule' @ yule.wessa.net






ANOVA Model
Ex ~ Jaar
means4.5582.722

\begin{tabular}{lllllllll}
\hline
ANOVA Model \tabularnewline
Ex  ~  Jaar \tabularnewline
means & 4.558 & 2.722 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270571&T=1

[TABLE]
[ROW][C]ANOVA Model[/C][/ROW]
[ROW][C]Ex  ~  Jaar[/C][/ROW]
[ROW][C]means[/C][C]4.558[/C][C]2.722[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270571&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270571&T=1

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

ANOVA Model
Ex ~ Jaar
means4.5582.722






ANOVA Statistics
DfSum SqMean SqF valuePr(>F)
Jaar1495.547495.547106.6360
Residuals2761282.5974.647 

\begin{tabular}{lllllllll}
\hline
ANOVA Statistics \tabularnewline
  & Df & Sum Sq & Mean Sq & F value & Pr(>F) \tabularnewline
Jaar & 1 & 495.547 & 495.547 & 106.636 & 0 \tabularnewline
Residuals & 276 & 1282.597 & 4.647 &   &   \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270571&T=2

[TABLE]
[ROW][C]ANOVA Statistics[/C][/ROW]
[ROW][C] [/C][C]Df[/C][C]Sum Sq[/C][C]Mean Sq[/C][C]F value[/C][C]Pr(>F)[/C][/ROW]
[ROW][C]Jaar[/C][C]1[/C][C]495.547[/C][C]495.547[/C][C]106.636[/C][C]0[/C][/ROW]
[ROW][C]Residuals[/C][C]276[/C][C]1282.597[/C][C]4.647[/C][C] [/C][C] [/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270571&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270571&T=2

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

ANOVA Statistics
DfSum SqMean SqF valuePr(>F)
Jaar1495.547495.547106.6360
Residuals2761282.5974.647






Tukey Honest Significant Difference Comparisons
difflwruprp adj
2012-20112.7222.2033.2410

\begin{tabular}{lllllllll}
\hline
Tukey Honest Significant Difference Comparisons \tabularnewline
  & diff & lwr & upr & p adj \tabularnewline
2012-2011 & 2.722 & 2.203 & 3.241 & 0 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270571&T=3

[TABLE]
[ROW][C]Tukey Honest Significant Difference Comparisons[/C][/ROW]
[ROW][C] [/C][C]diff[/C][C]lwr[/C][C]upr[/C][C]p adj[/C][/ROW]
[ROW][C]2012-2011[/C][C]2.722[/C][C]2.203[/C][C]3.241[/C][C]0[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270571&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270571&T=3

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Tukey Honest Significant Difference Comparisons
difflwruprp adj
2012-20112.7222.2033.2410






Levenes Test for Homogeneity of Variance
DfF valuePr(>F)
Group10.6150.434
276 

\begin{tabular}{lllllllll}
\hline
Levenes Test for Homogeneity of Variance \tabularnewline
  & Df & F value & Pr(>F) \tabularnewline
Group & 1 & 0.615 & 0.434 \tabularnewline
  & 276 &   &   \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270571&T=4

[TABLE]
[ROW][C]Levenes Test for Homogeneity of Variance[/C][/ROW]
[ROW][C] [/C][C]Df[/C][C]F value[/C][C]Pr(>F)[/C][/ROW]
[ROW][C]Group[/C][C]1[/C][C]0.615[/C][C]0.434[/C][/ROW]
[ROW][C] [/C][C]276[/C][C] [/C][C] [/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270571&T=4

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270571&T=4

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Levenes Test for Homogeneity of Variance
DfF valuePr(>F)
Group10.6150.434
276


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 5 ; par2 = 1 ; par3 = TRUE ;
Parameters (R input):
par1 = 5 ; par2 = 1 ; par3 = TRUE ;
R code (references can be found in the software module):
cat1 <- as.numeric(par1) #
cat2<- as.numeric(par2) #
intercept<-as.logical(par3)
x <- t(x)
x1<-as.numeric(x[,cat1])
f1<-as.character(x[,cat2])
xdf<-data.frame(x1,f1)
(V1<-dimnames(y)[[1]][cat1])
(V2<-dimnames(y)[[1]][cat2])
names(xdf)<-c('Response', 'Treatment')
if(intercept == FALSE) (lmxdf<-lm(Response ~ Treatment - 1, data = xdf) ) else (lmxdf<-lm(Response ~ Treatment, data = xdf) )
(aov.xdf<-aov(lmxdf) )
(anova.xdf<-anova(lmxdf) )
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'ANOVA Model', length(lmxdf$coefficients)+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, paste(V1, ' ~ ', V2), length(lmxdf$coefficients)+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'means',,TRUE)
for(i in 1:length(lmxdf$coefficients)){
a<-table.element(a, round(lmxdf$coefficients[i], digits=3),,FALSE)
}
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'ANOVA Statistics', 5+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, ' ',,TRUE)
a<-table.element(a, 'Df',,FALSE)
a<-table.element(a, 'Sum Sq',,FALSE)
a<-table.element(a, 'Mean Sq',,FALSE)
a<-table.element(a, 'F value',,FALSE)
a<-table.element(a, 'Pr(>F)',,FALSE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, V2,,TRUE)
a<-table.element(a, anova.xdf$Df[1],,FALSE)
a<-table.element(a, round(anova.xdf$'Sum Sq'[1], digits=3),,FALSE)
a<-table.element(a, round(anova.xdf$'Mean Sq'[1], digits=3),,FALSE)
a<-table.element(a, round(anova.xdf$'F value'[1], digits=3),,FALSE)
a<-table.element(a, round(anova.xdf$'Pr(>F)'[1], digits=3),,FALSE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Residuals',,TRUE)
a<-table.element(a, anova.xdf$Df[2],,FALSE)
a<-table.element(a, round(anova.xdf$'Sum Sq'[2], digits=3),,FALSE)
a<-table.element(a, round(anova.xdf$'Mean Sq'[2], digits=3),,FALSE)
a<-table.element(a, ' ',,FALSE)
a<-table.element(a, ' ',,FALSE)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1.tab')
bitmap(file='anovaplot.png')
boxplot(Response ~ Treatment, data=xdf, xlab=V2, ylab=V1)
dev.off()
if(intercept==TRUE){
'Tukey Plot'
thsd<-TukeyHSD(aov.xdf)
bitmap(file='TukeyHSDPlot.png')
plot(thsd)
dev.off()
}
if(intercept==TRUE){
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Tukey Honest Significant Difference Comparisons', 5,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, ' ', 1, TRUE)
for(i in 1:4){
a<-table.element(a,colnames(thsd[[1]])[i], 1, TRUE)
}
a<-table.row.end(a)
for(i in 1:length(rownames(thsd[[1]]))){
a<-table.row.start(a)
a<-table.element(a,rownames(thsd[[1]])[i], 1, TRUE)
for(j in 1:4){
a<-table.element(a,round(thsd[[1]][i,j], digits=3), 1, FALSE)
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')
}
if(intercept==FALSE){
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'TukeyHSD Message', 1,TRUE)
a<-table.row.end(a)
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Must Include Intercept to use Tukey Test ', 1, FALSE)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable2.tab')
}
library(car)
lt.lmxdf<-leveneTest(lmxdf)
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Levenes Test for Homogeneity of Variance', 4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,' ', 1, TRUE)
for (i in 1:3){
a<-table.element(a,names(lt.lmxdf)[i], 1, FALSE)
}
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Group', 1, TRUE)
for (i in 1:3){
a<-table.element(a,round(lt.lmxdf[[i]][1], digits=3), 1, FALSE)
}
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,' ', 1, TRUE)
a<-table.element(a,lt.lmxdf[[1]][2], 1, FALSE)
a<-table.element(a,' ', 1, FALSE)
a<-table.element(a,' ', 1, FALSE)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable3.tab')