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 computationFri, 01 Jun 2012 05:47:02 -0400
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Jun/01/t1338544037i6w1iowabhunwyb.htm/, Retrieved Thu, 02 May 2024 09:51:23 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=168405, Retrieved Thu, 02 May 2024 09:51:23 +0000
QR Codes:

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

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)] [male1] [2012-06-01 09:41:08] [2b8093ecf59eda89cc76f527bcb38785]
-    D    [One-Way-Between-Groups ANOVA- Free Statistics Software (Calculator)] [female 1] [2012-06-01 09:47:02] [e736236f6db1679287c406e2c94ebdc5] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
58	161	51	159
53	161	54	158
59	157	59	155
166	57	56	163
51	161	52	158
64	168	64	165
52	163	57	160
65	166	66	165
62	168	62	165
61	175	61	171
61	170	61	170
54	171	59	168
50	166	50	165
63	169	61	168
58	166	60	160
39	157	41	153
71	166	71	165
52	164	52	161
68	169	63	170
56	166	54	165
54	164	53	160
63	163	59	159
54	160	55	158
49	161	NA	NA
54	174	56	173
75	162	75	158
56	165	57	163
66	170	65	NA
78	173	75	169
60	162	59	160
64	165	63	163
64	164	62	161
52	158	51	155
62	175	61	171
55	165	54	163
56	163	57	159
50	166	50	161
50	171	NA	NA
50	160	55	150
63	160	64	158
61	165	60	163
53	169	52	175
60	167	55	163
56	170	56	170
53	165	53	165
57	163	59	160
57	162	56	160
56	161	56	161
56	165	57	160
50	169	50	165
52	159	52	153
55	155	NA	154
55	164	55	163
47	163	47	160
45	163	45	160
62	175	63	173
53	164	51	160
52	152	51	150
57	167	55	164
64	166	64	165
59	166	55	163
55	174	57	171
76	167	77	165
62	168	62	163
68	178	68	175
55	165	55	163
52	169	56	NA
47	153	NA	154
45	157	45	153
68	171	68	169
44	157	44	155
62	166	61	163
56	160	53	158
50	148	47	148
53	162	53	160
64	172	62	168
62	167	NA	NA
52	163	53	160
53	165	55	163
54	176	55	176
64	171	66	171
55	160	55	155
55	165	55	165
59	157	55	158
70	173	67	170
57	168	58	165
47	162	47	160
47	150	45	152
55	162	NA	NA
48	163	44	160
59	170	NA	NA
58	169	NA	NA
57	167	56	165
51	163	50	160
54	161	54	160
53	162	52	158
59	172	58	171
59	159	59	155
63	170	62	168
66	166	66	165
53	158	50	155
54	163	NA	NA
60	174	NA	NA
43	154	NA	NA
63	165	59	160
56	162	56	160
60	172	55	168
58	169	54	166
50	158	49	155
59	164	59	165
51	156	51	158
62	164	61	161

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'AstonUniversity' @ aston.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 & 1 seconds \tabularnewline
R Server & 'AstonUniversity' @ aston.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=168405&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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'AstonUniversity' @ aston.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=168405&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=168405&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 time1 seconds
R Server'AstonUniversity' @ aston.wessa.net






ANOVA Model
X1 ~ X2
means5047524743555149.251.66755.555.653.55752.55757.160.162.461.2557.4296159617456.33361.6675468166

\begin{tabular}{lllllllll}
\hline
ANOVA Model \tabularnewline
X1  ~  X2 \tabularnewline
means & 50 & 47 & 52 & 47 & 43 & 55 & 51 & 49.2 & 51.667 & 55.5 & 55.6 & 53.5 & 57 & 52.5 & 57 & 57.1 & 60.1 & 62.4 & 61.25 & 57.429 & 61 & 59 & 61 & 74 & 56.333 & 61.667 & 54 & 68 & 166 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=168405&T=1

[TABLE]
[ROW]ANOVA Model[/C][/ROW]
[ROW]X1  ~  X2[/C][/ROW]
[ROW][C]means[/C][C]50[/C][C]47[/C][C]52[/C][C]47[/C][C]43[/C][C]55[/C][C]51[/C][C]49.2[/C][C]51.667[/C][C]55.5[/C][C]55.6[/C][C]53.5[/C][C]57[/C][C]52.5[/C][C]57[/C][C]57.1[/C][C]60.1[/C][C]62.4[/C][C]61.25[/C][C]57.429[/C][C]61[/C][C]59[/C][C]61[/C][C]74[/C][C]56.333[/C][C]61.667[/C][C]54[/C][C]68[/C][C]166[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=168405&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=168405&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
X1 ~ X2
means5047524743555149.251.66755.555.653.55752.55757.160.162.461.2557.4296159617456.33361.6675468166






ANOVA Statistics
DfSum SqMean SqF valuePr(>F)
X229389251.03613422.45397.8850
Residuals832799.96433.735 

\begin{tabular}{lllllllll}
\hline
ANOVA Statistics \tabularnewline
  & Df & Sum Sq & Mean Sq & F value & Pr(>F) \tabularnewline
X2 & 29 & 389251.036 & 13422.45 & 397.885 & 0 \tabularnewline
Residuals & 83 & 2799.964 & 33.735 &   &   \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=168405&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]X2[/C][C]29[/C][C]389251.036[/C][C]13422.45[/C][C]397.885[/C][C]0[/C][/ROW]
[ROW][C]Residuals[/C][C]83[/C][C]2799.964[/C][C]33.735[/C][C] [/C][C] [/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=168405&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=168405&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)
X229389251.03613422.45397.8850
Residuals832799.96433.735






Must Include Intercept to use Tukey Test 

\begin{tabular}{lllllllll}
\hline
Must Include Intercept to use Tukey Test  \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=168405&T=3

[TABLE]
[ROW][C]Must Include Intercept to use Tukey Test [/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=168405&T=3

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

Must Include Intercept to use Tukey Test






Levenes Test for Homogeneity of Variance
DfF valuePr(>F)
Group280.9720.516
83 

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

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


Figures (Output of Computation)

Input Parameters & R Code

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