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:29:48 +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/t1418841013effmqqgda69ugaj.htm/, Retrieved Thu, 16 May 2024 19:25:17 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=270561, Retrieved Thu, 16 May 2024 19:25:17 +0000
QR Codes:

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

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:29:48] [18673d63f90870b9c004059cd6229007] [Current]
Feedback Forum

Post a new message

Dataset

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

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'Herman Ole Andreas Wold' @ wold.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 & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270561&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]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270561&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270561&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'Herman Ole Andreas Wold' @ wold.wessa.net






ANOVA Model
CH ~ Jaar
means41.250.304

\begin{tabular}{lllllllll}
\hline
ANOVA Model \tabularnewline
CH  ~  Jaar \tabularnewline
means & 41.25 & 0.304 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270561&T=1

[TABLE]
[ROW][C]ANOVA Model[/C][/ROW]
[ROW][C]CH  ~  Jaar[/C][/ROW]
[ROW][C]means[/C][C]41.25[/C][C]0.304[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270561&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270561&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
CH ~ Jaar
means41.250.304






ANOVA Statistics
DfSum SqMean SqF valuePr(>F)
Jaar16.1896.1890.0170.896
Residuals276100700.012364.855 

\begin{tabular}{lllllllll}
\hline
ANOVA Statistics \tabularnewline
  & Df & Sum Sq & Mean Sq & F value & Pr(>F) \tabularnewline
Jaar & 1 & 6.189 & 6.189 & 0.017 & 0.896 \tabularnewline
Residuals & 276 & 100700.012 & 364.855 &   &   \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270561&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]6.189[/C][C]6.189[/C][C]0.017[/C][C]0.896[/C][/ROW]
[ROW][C]Residuals[/C][C]276[/C][C]100700.012[/C][C]364.855[/C][C] [/C][C] [/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270561&T=2

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






Tukey Honest Significant Difference Comparisons
difflwruprp adj
2012-20110.304-4.2944.9020.896

\begin{tabular}{lllllllll}
\hline
Tukey Honest Significant Difference Comparisons \tabularnewline
  & diff & lwr & upr & p adj \tabularnewline
2012-2011 & 0.304 & -4.294 & 4.902 & 0.896 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270561&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]0.304[/C][C]-4.294[/C][C]4.902[/C][C]0.896[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270561&T=3

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






Levenes Test for Homogeneity of Variance
DfF valuePr(>F)
Group12.5950.108
276 

\begin{tabular}{lllllllll}
\hline
Levenes Test for Homogeneity of Variance \tabularnewline
  & Df & F value & Pr(>F) \tabularnewline
Group & 1 & 2.595 & 0.108 \tabularnewline
  & 276 &   &   \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270561&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]2.595[/C][C]0.108[/C][/ROW]
[ROW][C] [/C][C]276[/C][C] [/C][C] [/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270561&T=4

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


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 4 ; par2 = 1 ; par3 = TRUE ;
Parameters (R input):
par1 = 4 ; 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')