Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

rwasp_One Factor ANOVA.wasp

Title produced by software

One-Way-Between-Groups ANOVA- Free Statistics Software (Calculator)

Date of computation

Mon, 15 Dec 2014 12:12:09 +0000

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Dec/15/t14186456525778k999iipeipq.htm/, Retrieved Fri, 01 Nov 2024 00:06:49 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=268214, Retrieved Fri, 01 Nov 2024 00:06:49 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

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-15 12:12:09] [18673d63f90870b9c004059cd6229007] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

0 21 149 96 18 68
0 22 139 70 31 39
0 22 148 88 39 32
0 18 158 114 46 62
0 23 128 69 31 33
0 12 224 176 67 52
0 20 159 114 35 62
0 22 105 121 52 77
0 21 159 110 77 76
0 19 167 158 37 41
0 22 165 116 32 48
0 15 159 181 36 63
0 20 119 77 38 30
0 19 176 141 69 78
0 18 54 35 21 19
0 15 91 80 26 31
0 20 163 152 54 66
0 21 124 97 36 35
0 21 137 99 42 42
0 15 121 84 23 45
0 16 153 68 34 21
0 23 148 101 112 25
0 21 221 107 35 44
0 18 188 88 47 69
0 25 149 112 47 54
0 9 244 171 37 74
0 30 148 137 109 80
0 20 92 77 24 42
0 23 150 66 20 61
0 16 153 93 22 41
0 16 94 105 23 46
0 19 156 131 32 39
0 25 132 102 30 34
0 18 161 161 92 51
0 23 105 120 43 42
0 21 97 127 55 31
0 10 151 77 16 39
0 14 131 108 49 20
0 22 166 85 71 49
0 26 157 168 43 53
0 23 111 48 29 31
0 23 145 152 56 39
0 24 162 75 46 54
0 24 163 107 19 49
0 18 59 62 23 34
0 23 187 121 59 46
0 15 109 124 30 55
0 19 90 72 61 42
0 16 105 40 7 50
0 25 83 58 38 13
0 23 116 97 32 37
0 17 42 88 16 25
0 19 148 126 19 30
0 21 155 104 22 28
0 18 125 148 48 45
0 27 116 146 23 35
0 21 128 80 26 28
0 13 138 97 33 41
0 8 49 25 9 6
0 29 96 99 24 45
0 28 164 118 34 73
0 23 162 58 48 17
0 21 99 63 18 40
0 19 202 139 43 64
0 19 186 50 33 37
0 20 66 60 28 25
0 18 183 152 71 65
0 19 214 142 26 100
0 17 188 94 67 28
0 19 104 66 34 35
0 25 177 127 80 56
0 19 126 67 29 29
0 22 76 90 16 43
0 23 99 75 59 59
0 14 139 128 32 50
0 16 162 146 43 59
0 24 108 69 38 27
0 20 159 186 29 61
0 12 74 81 36 28
0 24 110 85 32 51
0 22 96 54 35 35
0 12 116 46 21 29
0 22 87 106 29 48
0 20 97 34 12 25
0 10 127 60 37 44
0 23 106 95 37 64
0 17 80 57 47 32
0 22 74 62 51 20
0 24 91 36 32 28
0 18 133 56 21 34
0 21 74 54 13 31
0 20 114 64 14 26
0 20 140 76 -2 58
0 22 95 98 20 23
0 19 98 88 24 21
0 20 121 35 11 21
0 26 126 102 23 33
0 23 98 61 24 16
0 24 95 80 14 20
0 21 110 49 52 37
0 21 70 78 15 35
0 19 102 90 23 33
0 8 86 45 19 27
0 17 130 55 35 41
0 20 96 96 24 40
0 11 102 43 39 35
0 8 100 52 29 28
0 15 94 60 13 32
0 18 52 54 8 22
0 18 98 51 18 44
0 19 118 51 24 27
0 19 99 38 19 17
1 23 48 41 23 12
1 22 50 146 16 45
1 21 150 182 33 37
1 25 154 192 32 37
1 30 109 263 37 108
1 17 68 35 14 10
1 27 194 439 52 68
1 23 158 214 75 72
1 23 159 341 72 143
1 18 67 58 15 9
1 18 147 292 29 55
1 23 39 85 13 17
1 19 100 200 40 37
1 15 111 158 19 27
1 20 138 199 24 37
1 16 101 297 121 58
1 24 131 227 93 66
1 25 101 108 36 21
1 25 114 86 23 19
1 19 165 302 85 78
1 19 114 148 41 35
1 16 111 178 46 48
1 19 75 120 18 27
1 19 82 207 35 43
1 23 121 157 17 30
1 21 32 128 4 25
1 22 150 296 28 69
1 19 117 323 44 72
1 20 71 79 10 23
1 20 165 70 38 13
1 3 154 146 57 61
1 23 126 246 23 43
1 23 149 196 36 51
1 20 145 199 22 67
1 15 120 127 40 36
1 16 109 153 31 44
1 7 132 299 11 45
1 24 172 228 38 34
1 17 169 190 24 36
1 24 114 180 37 72
1 24 156 212 37 39
1 19 172 269 22 43
1 25 68 130 15 25
1 20 89 179 2 56
1 28 167 243 43 80
1 23 113 190 31 40
1 27 115 299 29 73
1 18 78 121 45 34
1 28 118 137 25 72
1 21 87 305 4 42
1 19 173 157 31 61
1 23 2 96 -4 23
1 27 162 183 66 74
1 22 49 52 61 16
1 28 122 238 32 66
1 25 96 40 31 9
1 21 100 226 39 41
1 22 82 190 19 57
1 28 100 214 31 48
1 20 115 145 36 51
1 29 141 119 42 53
1 25 165 222 21 29
1 25 165 222 21 29
1 20 110 159 25 55
1 20 118 165 32 54
1 16 158 249 26 43
1 20 146 125 28 51
1 20 49 122 32 20
1 23 90 186 41 79
1 18 121 148 29 39
1 25 155 274 33 61
1 18 104 172 17 55
1 19 147 84 13 30
1 25 110 168 32 55
1 25 108 102 30 22
1 25 113 106 34 37
1 24 115 2 59 2
1 19 61 139 13 38
1 26 60 95 23 27
1 10 109 130 10 56
1 17 68 72 5 25
1 13 111 141 31 39
1 17 77 113 19 33
1 30 73 206 32 43
1 25 151 268 30 57
1 4 89 175 25 43
1 16 78 77 48 23
1 21 110 125 35 44
1 23 220 255 67 54
1 22 65 111 15 28
1 17 141 132 22 36
1 20 117 211 18 39
1 20 122 92 33 16
1 22 63 76 46 23
1 16 44 171 24 40
1 23 52 83 14 24
1 0 131 266 12 78
1 18 101 186 38 57
1 25 42 50 12 37
1 23 152 117 28 27
1 12 107 219 41 61
1 18 77 246 12 27
1 24 154 279 31 69
1 11 103 148 33 34
1 18 96 137 34 44
1 23 175 181 21 34
1 24 57 98 20 39
1 29 112 226 44 51
1 18 143 234 52 34
1 15 49 138 7 31
1 29 110 85 29 13
1 16 131 66 11 12
1 19 167 236 26 51
1 22 56 106 24 24
1 16 137 135 7 19
1 23 86 122 60 30
1 23 121 218 13 81
1 19 149 199 20 42
1 4 168 112 52 22
1 20 140 278 28 85
1 24 88 94 25 27
1 20 168 113 39 25
1 4 94 84 9 22
1 24 51 86 19 19
1 22 48 62 13 14
1 16 145 222 60 45
1 3 66 167 19 45
1 15 85 82 34 28
1 24 109 207 14 51
1 17 63 184 17 41
1 20 102 83 45 31
1 27 162 183 66 74
1 26 86 89 48 19
1 23 114 225 29 51
1 17 164 237 -2 73
1 20 119 102 51 24
1 22 126 221 2 61
1 19 132 128 24 23
1 24 142 91 40 14
1 19 83 198 20 54
1 23 94 204 19 51
1 15 81 158 16 62
1 27 166 138 20 36
1 26 110 226 40 59
1 22 64 44 27 24
1 22 93 196 25 26
1 18 104 83 49 54
1 15 105 79 39 39
1 22 49 52 61 16
1 27 88 105 19 36
1 10 95 116 67 31
1 20 102 83 45 31
1 17 99 196 30 42
1 23 63 153 8 39
1 19 76 157 19 25
1 13 109 75 52 31
1 27 117 106 22 38
1 23 57 58 17 31
1 16 120 75 33 17
1 25 73 74 34 22
1 2 91 185 22 55
1 26 108 265 30 62
1 20 105 131 25 51
1 23 117 139 38 30
1 22 119 196 26 49
1 24 31 78 13 16

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 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=268214&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=268214&T=0

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

As an alternative you can also use a QR Code:

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

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'George Udny Yule' @ yule.wessa.net

ANOVA Model
LFM ~ Jaar_Bin
means	126.705	-17.254

\begin{tabular}{lllllllll}
\hline
ANOVA Model \tabularnewline
LFM  ~  Jaar_Bin \tabularnewline
means & 126.705 & -17.254 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=268214&T=1

[TABLE]
[ROW][C]ANOVA Model[/C][/ROW]
[ROW][C]LFM  ~  Jaar_Bin[/C][/ROW]
[ROW][C]means[/C][C]126.705[/C][C]-17.254[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=268214&T=1

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

As an alternative you can also use a QR Code:

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

ANOVA Model
LFM ~ Jaar_Bin
means	126.705	-17.254

ANOVA Statistics
	Df	Sum Sq	Mean Sq	F value	Pr(>F)
Jaar_Bin	1	19908.486	19908.486	13.094	0
Residuals	276	419628.391	1520.393

\begin{tabular}{lllllllll}
\hline
ANOVA Statistics \tabularnewline
  & Df & Sum Sq & Mean Sq & F value & Pr(>F) \tabularnewline
Jaar_Bin & 1 & 19908.486 & 19908.486 & 13.094 & 0 \tabularnewline
Residuals & 276 & 419628.391 & 1520.393 &   &   \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=268214&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_Bin[/C][C]1[/C][C]19908.486[/C][C]19908.486[/C][C]13.094[/C][C]0[/C][/ROW]
[ROW][C]Residuals[/C][C]276[/C][C]419628.391[/C][C]1520.393[/C][C] [/C][C] [/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=268214&T=2

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

As an alternative you can also use a QR Code:

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

ANOVA Statistics
	Df	Sum Sq	Mean Sq	F value	Pr(>F)
Jaar_Bin	1	19908.486	19908.486	13.094	0
Residuals	276	419628.391	1520.393

Tukey Honest Significant Difference Comparisons
	diff	lwr	upr	p adj
1-0	-17.254	-26.64	-7.867	0

\begin{tabular}{lllllllll}
\hline
Tukey Honest Significant Difference Comparisons \tabularnewline
  & diff & lwr & upr & p adj \tabularnewline
1-0 & -17.254 & -26.64 & -7.867 & 0 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=268214&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]1-0[/C][C]-17.254[/C][C]-26.64[/C][C]-7.867[/C][C]0[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=268214&T=3

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

As an alternative you can also use a QR Code:

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

Tukey Honest Significant Difference Comparisons
	diff	lwr	upr	p adj
1-0	-17.254	-26.64	-7.867	0

Levenes Test for Homogeneity of Variance
	Df	F value	Pr(>F)
Group	1	0.412	0.521
	276

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

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

As an alternative you can also use a QR Code:

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

Levenes Test for Homogeneity of Variance
	Df	F value	Pr(>F)
Group	1	0.412	0.521
	276

Figure 1

PNG link

Postscript link

PDF link

Figure 2

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = 3 ; par2 = 1 ; par3 = TRUE ;

Parameters (R input):

par1 = 3 ; 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')

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code