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

Thu, 01 Nov 2012 08:30:59 -0400

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/2012/Nov/01/t1351773072sbtx925dtqhc410.htm/, Retrieved Sat, 20 Apr 2024 05:13:41 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=185471, Retrieved Sat, 20 Apr 2024 05:13:41 +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

150

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)] [anova] [2012-11-01 12:30:59] [64d4a1c1c0951cc39d80412399cbae13] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

2	1	63	3	30	68
3	1	63	3	30	68
4	1	63	3	30	68
5	1	63	3	30	68
6	1	67	4	31	70
7	1	68	4	32	70
8	1	69	4	32	70
9	1	70	4	34	72
10	1	70	4	34	73
11	1	71	4	34	75
12	1	71	5	34	75
13	1	71	5	36	75
14	1	71	5	36	75
15	1	71	5	36	75
16	1	72	5	36	77
17	1	72	5	37	77
18	1	73	5	37	77
19	1	74	5	37	77
20	1	74	5	37	78
21	1	74	5	37	78
22	1	74	5	39	78
23	1	75	6	39	79
25	1	75	6	39	79
26	1	75	6	39	79
27	1	75	6	39	80
28	1	75	6	39	80
29	1	75	6	39	80
30	1	76	6	39	81
31	1	76	6	39	82
32	1	76	6	39	84
33	1	76	6	41	84
34	1	77	6	41	85
35	1	77	6	41	85
36	1	77	6	41	85
37	1	77	6	41	85
38	1	77	6	41	86
39	1	77	6	41	86
40	2	77	6	41	86
41	2	78	6	41	86
42	2	78	6	42	86
43	2	79	7	43	86
44	2	79	7	43	86
45	2	79	7	43	86
46	2	79	7	43	86
47	2	79	7	43	87
48	2	79	7	43	87
50	2	80	7	44	87
51	2	80	7	44	87
52	2	80	7	44	87
53	2	80	7	44	87
54	2	80	7	44	87
55	2	80	7	44	88
56	2	81	7	44	88
57	2	81	7	44	88
58	2	81	7	44	88
60	2	81	7	44	88
61	2	81	7	45	88
62	2	81	7	46	88
63	2	81	7	46	90
65	2	82	7	46	90
66	2	82	7	46	90
67	2	82	7	46	90
68	2	83	7	46	90
69	2	83	7	46	90
70	2	83	7	46	90
71	2	83	7	46	90
72	2	83	7	46	90
73	2	83	7	48	90
74	2	83	7	48	91
75	2	84	8	48	91
76	2	84	8	48	91
77	2	84	8	48	91
78	2	84	8	48	91
79	2	84	8	48	92
80	2	84	8	48	92
81	2	84	8	48	92
82	2	85	8	48	92
83	2	85	8	49	92
84	2	85	8	50	92
85	3	85	8	50	94
86	3	85	8	50	94
87	3	85	8	50	94
88	3	85	8	50	94
89	3	85	8	50	94
90	3	85	8	50	94
91	3	85	8	50	94
92	3	85	8	50	94
93	3	85	8	50	95
94	3	86	8	50	95
95	3	86	8	50	95
96	3	86	8	51	95
97	3	86	8	52	95
98	3	86	8	52	95
99	3	86	8	52	95
100	3	86	8	52	95
101	3	87	8	52	95
102	3	87	8	52	95
103	3	87	8	52	95
104	3	87	8	52	96
105	3	87	8	52	96
106	3	87	8	52	96
107	3	88	8	52	96
108	3	88	8	52	96
109	3	89	8	54	97
110	3	89	8	54	97
111	3	89	8	54	97
112	3	89	8	54	97
113	3	89	8	54	97
114	3	90	8	54	97
115	3	90	9	54	97
116	3	91	9	54	97
117	3	91	9	55	97
118	3	91	9	55	97
119	3	92	9	55	97
120	3	92	9	55	97
121	3	92	9	55	98
122	3	92	9	55	98
123	3	92	9	55	98
124	3	92	9	55	98
125	3	93	9	55	98
126	3	93	9	56	100
127	3	93	9	56	100
128	3	93	9	56	100
130	3	93	9	56	101
131	3	93	9	56	101
132	3	93	9	56	101
133	3	95	9	56	101
134	3	95	9	57	101
135	3	95	9	57	101
136	3	96	9	57	101
137	3	96	9	57	102
138	3	97	9	57	102
139	3	98	9	58	102
140	3	98	9	58	103
141	3	99	9	58	103
142	3	99	9	59	105
143	3	99	9	59	106
144	3	101	9	59	106
145	3	101	9	59	107
146	3	102	9	59	107
147	3	103	9	59	107
148	3	103	9	59	108
149	3	103	10	60	109
150	3	104	10	60	109
151	3	104	10	60	111
152	3	104	10	60	111
153	3	107	10	62	111
154	3	109	10	62	111
155	3	110	10	62	112
156	3	111	10	66	115
157	3	111	10	66	115

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Herman Ole Andreas Wold' @ wold.wessa.net
R Engine error message	Error in array(list(2, 1, 63, 3, 30, 68, 3, 1, 63, 3, 30, 68, 4, 1, 63, : length of 'dimnames' [1] not equal to array extent Execution halted

\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 & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
R Engine error message & Error in array(list(2, 1, 63, 3, 30, 68, 3, 1, 63, 3, 30, 68, 4, 1, 63,  : 
  length of 'dimnames' [1] not equal to array extent
Execution halted
 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=185471&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]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[ROW][C]R Engine error message[/C][C]Error in array(list(2, 1, 63, 3, 30, 68, 3, 1, 63, 3, 30, 68, 4, 1, 63,  : 
  length of 'dimnames' [1] not equal to array extent
Execution halted
[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=185471&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=185471&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	1 seconds
R Server	'Herman Ole Andreas Wold' @ wold.wessa.net
R Engine error message	Error in array(list(2, 1, 63, 3, 30, 68, 3, 1, 63, 3, 30, 68, 4, 1, 63, : length of 'dimnames' [1] not equal to array extent Execution halted

Parameters (Session):

par1 = 1 ; par2 = 2 ; par3 = FALSE ;

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code