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_Two Factor ANOVA.wasp
Title produced by softwareTwo-Way ANOVA
Date of computationWed, 01 Dec 2010 15:58:18 +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/2010/Dec/01/t1291218998afexh5thgt31l7a.htm/, Retrieved Sat, 04 May 2024 20:29:46 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=104081, Retrieved Sat, 04 May 2024 20:29:46 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Explorative Data Analysis] [time effect in su...] [2010-11-17 08:55:33] [b98453cac15ba1066b407e146608df68]
- RMPD    [Two-Way ANOVA] [] [2010-12-01 15:58:18] [74602ad29729a97b452f3d2abd95057e] [Current]
-    D      [Two-Way ANOVA] [] [2010-12-01 19:20:05] [58af523ef9b33032fd2497c80088399b]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
162556	1081	807	213118	6282154
29790	309	444	81767	4321023
87550	458	412	153198	4111912
84738	588	428	-26007	223193
54660	302	315	126942	1491348
42634	156	168	157214	1629616
40949	481	263	129352	1398893
45187	353	267	234817	1926517
37704	452	228	60448	983660
16275	109	129	47818	1443586
25830	115	104	245546	1073089
12679	110	122	48020	984885
18014	239	393	-1710	1405225
43556	247	190	32648	227132
24811	505	280	95350	929118
6575	159	63	151352	1071292
7123	109	102	288170	638830
21950	519	265	114337	856956
37597	248	234	37884	992426
17821	373	277	122844	444477
12988	119	73	82340	857217
22330	84	67	79801	711969
13326	102	103	165548	702380
16189	295	290	116384	358589
7146	105	83	134028	297978
15824	64	56	63838	585715
27664	282	236	74996	657954
11920	182	73	31080	209458
8568	37	34	32168	786690
14416	361	139	49857	439798
3369	28	26	87161	688779
11819	85	70	106113	574339
6984	45	40	80570	741409
4519	49	42	102129	597793
2220	22	12	301670	644190
18562	155	211	102313	377934
10327	91	74	88577	640273
5336	81	80	112477	697458
2365	79	83	191778	550608
4069	145	131	79804	207393
8636	855	203	128294	301607
13718	61	56	96448	345783
4525	226	89	93811	501749
6869	105	88	117520	379983
4628	62	39	69159	387475
3689	25	25	101792	377305
4891	217	49	210568	370837
7489	322	149	136996	430866
4901	84	58	121920	469107
2284	33	41	76403	194493
3160	108	90	108094	530670
4150	150	136	134759	518365
7285	115	97	188873	491303
1134	162	63	146216	527021
4658	158	114	156608	233773
2384	97	77	61348	405972
3748	9	6	50350	652925
5371	66	47	87720	446211
1285	107	51	99489	341340
9327	101	85	87419	387699
5565	47	43	94355	493408
1528	38	32	60326	146494
3122	34	25	94670	414462
7561	87	77	82425	364304
2675	79	54	59017	355178
13253	947	251	90829	357760
880	74	15	80791	261216
2053	53	44	100423	397144
1424	94	73	131116	374943
4036	63	85	100269	424898
3045	58	49	27330	202055
5119	49	38	39039	378525
1431	34	35	106885	310768
554	11	9	79285	325738
1975	35	34	118881	394510
1765	20	20	77623	247060
1012	47	29	114768	368078
810	43	11	74015	236761
1280	117	52	69465	312378
666	171	13	117869	339836
1380	26	29	60982	347385
4677	75	66	90131	426280
876	59	33	138971	352850
814	18	15	39625	301881
514	15	15	102725	377516
5692	72	68	64239	357312
3642	86	100	90262	458343
540	14	13	103960	354228
2099	64	45	106611	308636
567	11	14	103345	386212
2001	52	36	95551	393343
2949	41	40	82903	378509
2253	99	68	63593	452469
6533	75	29	126910	364839
1889	45	43	37527	358649
3055	43	30	60247	376641
272	8	9	112995	429112
1414	198	22	70184	330546
2564	22	19	130140	403560
1383	11	9	73221	317892
1261	33	31	76114	307528
975	23	19	90534	235133
3366	80	55	108479	299243
576	18	8	113761	314073
1686	40	28	68696	368186
746	23	29	71561	269661
3192	60	48	59831	125390
2045	20	16	97890	510834
5702	61	47	101481	321896
1932	36	20	72954	249898
936	30	22	67939	408881
3437	47	33	48022	158492
5131	71	44	86111	292154
2397	14	13	74020	289513
1389	9	6	57530	378049
1503	39	35	56364	343466
402	26	8	84990	332743
2239	21	17	88590	442882
2234	16	11	77200	214215
837	69	21	61262	315688
10579	92	92	110309	375195
875	14	12	67000	334280
1585	107	112	93099	355864
1659	29	25	107577	480382
2647	37	17	62920	353058
3294	23	23	75832	217193
0	0	0	60720	315380
94	7	10	60793	314533
422	28	23	57935	318056
0	0	0	60720	315380
34	8	7	60630	314353
1558	63	25	55637	369448
0	0	0	60720	315380
43	3	20	60887	312846
645	5	4	60720	312075
316	9	4	60505	315009
115	13	10	60945	318903
5	2	1	60720	314887
897	5	4	60720	314913
0	0	0	60720	315380
389	14	8	58990	325506
0	0	0	60720	315380
1002	15	11	56750	298568
36	3	4	60894	315834
460	15	15	63346	329784
309	11	9	56535	312878
0	0	0	60720	315380
9	6	7	60835	314987
271	2	2	60720	325249
14	1	0	61016	315877
520	10	7	58650	291650
1766	73	46	60438	305959
0	0	5	60720	315380
458	11	7	58625	297765
20	3	2	60938	315245
0	0	0	60720	315380
0	0	0	60720	315380
98	2	2	61490	315236
405	7	5	60845	336425
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
483	27	7	60830	306268
454	51	24	63261	302187
47	3	1	60720	314882
0	0	0	60720	315380
757	19	18	45689	382712
4655	393	55	60720	341570
0	0	0	60720	315380
0	0	0	60720	315380
36	4	3	61564	312412
0	0	0	60720	315380
203	9	9	61938	309596
0	0	0	60720	315380
126	10	8	60951	315547
400	152	113	60720	313267
71	1	0	60745	316176
0	0	0	60720	315380
0	0	0	60720	315380
972	34	19	71642	359335
531	10	11	71641	330068
2461	57	25	55792	314289
378	52	16	71873	297413
23	5	5	62555	314806
638	14	11	60370	333210
2300	29	23	64873	352108
149	5	6	62041	313332
226	5	5	65745	291787
0	0	0	60720	315380
275	4	7	59500	318745
0	0	0	60720	315380
141	6	7	61630	315366
0	0	0	60720	315380
28	2	3	60890	315688
0	0	0	60720	315380
4980	91	89	113521	409642
0	0	0	60720	315380
0	0	0	60720	315380
472	20	19	80045	269587
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
203	27	12	50804	300962
496	17	12	87390	325479
10	2	5	61656	316155
63	4	2	65688	318574
0	0	0	60720	315380
1136	32	26	48522	343613
265	31	3	60720	306948
0	0	0	60720	315380
0	0	0	60720	315380
267	32	11	57640	330059
474	20	10	61977	288985
534	7	5	62620	304485
0	0	2	60720	315380
15	8	6	60831	315688
397	28	7	60646	317736
0	0	2	60720	315380
1866	29	28	56225	322331
288	4	3	60510	296656
0	0	0	60720	315380
3	2	1	60698	315354
468	21	20	60720	312161
20	2	1	60805	315576
278	26	22	61404	314922
61	14	9	60720	314551
0	0	0	60720	315380
192	4	2	65276	312339
0	0	0	60720	315380
317	9	7	63915	298700
738	10	9	60720	321376
0	0	0	60720	315380
368	17	13	61686	303230
0	0	0	60720	315380
2	1	0	60743	315487
0	0	0	60720	315380
53	6	6	60349	315793
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
94	3	3	61360	312887
0	0	0	60720	315380
24	8	7	59818	315637
2332	4	2	72680	324385
0	0	0	60720	315380
0	0	0	60720	315380
131	11	15	61808	308989
0	0	0	60720	315380
0	0	0	60720	315380
206	9	9	53110	296702
0	0	0	60720	315380
167	2	1	64245	307322
622	73	38	73007	304376
2328	94	57	82732	253588
0	0	0	60720	315380
365	8	7	54820	309560
364	35	26	47705	298466
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
226	12	13	72835	343929
307	15	10	58856	331955
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
188	11	9	77655	381180
0	0	0	60720	315380
138	6	26	69817	331420
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
125	12	19	60798	310201
0	0	0	60720	315380
282	30	12	62452	320016
335	33	23	64175	320398
0	0	0	60720	315380
1324	117	29	67440	291841
176	28	8	68136	310670
0	0	0	60720	315380
0	0	0	60720	315380
249	72	26	56726	313491
0	0	0	60720	315380
333	13	9	70811	331323
0	0	0	60720	315380
601	6	5	60720	319210
30	4	3	62045	318098
0	0	0	60720	315380
249	62	13	54323	292754
0	0	0	60720	315380
165	24	12	62841	325176
453	21	19	81125	365959
0	0	0	60720	315380
53	14	10	59506	302409
382	21	9	59365	340968
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	9	60720	315380
30	4	4	60798	313164
290	2	1	58790	301164
0	0	1	60720	315380
0	0	0	60720	315380
366	53	14	61808	344425
2	9	12	60735	315394
0	0	0	60720	315380
209	13	19	64016	316647
384	22	17	54683	309836
0	0	0	60720	315380
0	0	0	60720	315380
365	83	32	87192	346611
0	0	0	60720	315380
49	8	14	64107	322031
3	4	8	60761	315656
133	14	4	65990	339445
32	1	0	59988	314964
368	17	20	61167	297141
1	6	5	60719	315372
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
22	2	1	60722	312502
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
96	5	4	60379	313729
1	2	1	60727	315388
314	5	4	60720	315371
844	78	20	60925	296139
0	0	0	60720	315380
26	1	1	60896	313880
125	13	10	59734	317698
304	15	12	62969	295580
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
621	48	13	60720	308256
0	0	0	60720	315380
119	6	3	59118	303677
0	0	0	60720	315380
0	0	0	60720	315380
1595	17	10	60720	319369
312	14	3	58598	318690
60	10	7	61124	314049
587	12	10	59595	325699
135	2	1	62065	314210
0	0	0	60720	315380
0	0	0	60720	315380
514	52	15	78780	322378
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
1	4	4	60722	315398
0	0	0	60720	315380
0	0	0	60720	315380
1763	24	28	61600	308336
180	11	9	59635	316386
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
0	0	0	60720	315380
218	21	7	60720	315553
0	0	0	60720	315380
448	40	7	59781	323361
227	9	7	76644	336639
174	1	3	64820	307424
0	0	0	60720	315380
0	0	0	60720	315380
121	24	11	56178	295370
607	11	7	60436	322340
2212	14	10	60720	319864
0	0	0	60720	315380
0	0	0	60720	315380
530	60	18	73433	317291
571	80	14	41477	280398
0	0	0	60720	315380
78	16	12	62700	317330
2489	40	29	67804	238125
131	6	3	59661	327071
923	8	6	58620	309038
72	3	3	60398	314210
572	16	8	58580	307930
397	10	10	62710	322327
450	8	6	59325	292136
622	7	8	60950	263276
694	8	6	68060	367655
3425	12	9	83620	283910
562	13	8	58456	283587
4917	42	26	52811	243650
1442	118	239	121173	438493
529	9	7	63870	296261
2126	138	41	21001	230621
1061	5	3	70415	304252
776	9	8	64230	333505
611	8	6	59190	296919
1526	25	21	69351	278990
592	7	7	64270	276898
1182	13	11	70694	327007
621	16	11	68005	317046
989	11	12	58930	304555
438	11	9	58320	298096
726	3	3	69980	231861
1303	61	57	69863	309422
7419	29	21	63255	286963
1164	17	15	57320	269753
3310	33	32	75230	448243
1920	15	11	79420	165404
965	3	2	73490	204325
3256	66	23	35250	407159
1135	17	20	62285	290476
1270	26	24	69206	275311
661	3	1	65920	246541
1013	2	1	69770	253468
2844	67	74	72683	240897
11528	70	68	-14545	-83265
6526	26	20	55830	-42143
2264	24	20	55174	272713
5109	97	82	67038	215362
3999	30	21	51252	42754
35624	223	244	157278	306275
9252	48	32	79510	253537
15236	90	86	77440	372631
18073	180	69	27284	-7170

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time280 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 280 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=104081&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]280 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=104081&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=104081&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 time280 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 1 ; par2 = 2 ; par3 = 3 ; par4 = FALSE ;
Parameters (R input):
par1 = 1 ; par2 = 2 ; par3 = 3 ; par4 = FALSE ;
R code (references can be found in the software module):
cat1 <- as.numeric(par1) #
cat2<- as.numeric(par2) #
cat3 <- as.numeric(par3)
intercept<-as.logical(par4)
x <- t(x)
x1<-as.numeric(x[,cat1])
f1<-as.character(x[,cat2])
f2 <- as.character(x[,cat3])
xdf<-data.frame(x1,f1, f2)
(V1<-dimnames(y)[[1]][cat1])
(V2<-dimnames(y)[[1]][cat2])
(V3 <-dimnames(y)[[1]][cat3])
names(xdf)<-c('Response', 'Treatment_A', 'Treatment_B')
if(intercept == FALSE) (lmxdf<-lm(Response ~ Treatment_A * Treatment_B- 1, data = xdf) ) else (lmxdf<-lm(Response ~ Treatment_A * Treatment_B, 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, lmxdf$call['formula'],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)
for(i in 1 : length(rownames(anova.xdf))-1){
a<-table.row.start(a)
a<-table.element(a,rownames(anova.xdf)[i] ,,TRUE)
a<-table.element(a, anova.xdf$Df[1],,FALSE)
a<-table.element(a, round(anova.xdf$'Sum Sq'[i], digits=3),,FALSE)
a<-table.element(a, round(anova.xdf$'Mean Sq'[i], digits=3),,FALSE)
a<-table.element(a, round(anova.xdf$'F value'[i], digits=3),,FALSE)
a<-table.element(a, round(anova.xdf$'Pr(>F)'[i], 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'[i+1],,FALSE)
a<-table.element(a, round(anova.xdf$'Sum Sq'[i+1], digits=3),,FALSE)
a<-table.element(a, round(anova.xdf$'Mean Sq'[i+1], 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_A + Treatment_B, data=xdf, xlab=V2, ylab=V1, main='Boxplots of ANOVA Groups')
dev.off()
bitmap(file='designplot.png')
xdf2 <- xdf # to preserve xdf make copy for function
names(xdf2) <- c(V1, V2, V3)
plot.design(xdf2, main='Design Plot of Group Means')
dev.off()
bitmap(file='interactionplot.png')
interaction.plot(xdf$Treatment_A, xdf$Treatment_B, xdf$Response, xlab=V2, ylab=V1, trace.label=V3, main='Possible Interactions Between Anova Groups')
dev.off()
if(intercept==TRUE){
thsd<-TukeyHSD(aov.xdf)
names(thsd) <- c(V2, V3, paste(V2, ':', V3, sep=''))
bitmap(file='TukeyHSDPlot.png')
layout(matrix(c(1,2,3,3), 2,2))
plot(thsd, las=1)
dev.off()
}
if(intercept==TRUE){
ntables<-length(names(thsd))
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(nt in 1:ntables){
for(i in 1:length(rownames(thsd[[nt]]))){
a<-table.row.start(a)
a<-table.element(a,rownames(thsd[[nt]])[i], 1, TRUE)
for(j in 1:4){
a<-table.element(a,round(thsd[[nt]][i,j], digits=3), 1, FALSE)
}
a<-table.row.end(a)
}
} # end nt
a<-table.end(a)
table.save(a,file='hsdtable.tab')
}#end if hsd tables
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<-levene.test(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')