Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software ModuleIan.Hollidayrwasp_Two Factor ANOVA -V2.wasp
Title produced by softwareVariability
Date of computationWed, 02 Jun 2010 09:17:37 +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/Jun/02/t1275470319jisa9ss1c86m8sh.htm/, Retrieved Fri, 29 Mar 2024 12:07:50 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=77018, Retrieved Fri, 29 Mar 2024 12:07:50 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact137
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Two-Way ANOVA] [two-way anova wit...] [2010-05-26 17:02:24] [98fd0e87c3eb04e0cc2efde01dbafab6]
- R PD  [Two-Way ANOVA] [ANOVA with good l...] [2010-05-28 23:09:47] [98fd0e87c3eb04e0cc2efde01dbafab6]
- R       [Variability] [ANOVA with better...] [2010-05-29 09:47:12] [98fd0e87c3eb04e0cc2efde01dbafab6]
- R         [Variability] [ANOVA with better...] [2010-05-29 09:54:40] [98fd0e87c3eb04e0cc2efde01dbafab6]
- R  D          [Variability] [Question 1 ii)] [2010-06-02 09:17:37] [d0320cec1021a7468159c5ccef08dea7] [Current]
Feedback Forum

Post a new message
Dataseries X:
'GOOD'	'HIGH'	25
'GOOD'	'HIGH'	0
'GOOD'	'HIGH'	-16
'GOOD'	'HIGH'	5
'GOOD'	'HIGH'	11
'GOOD'	'HIGH'	-6
'GOOD'	'HIGH'	-2
'GOOD'	'HIGH'	-13
'GOOD'	'HIGH'	14
'GOOD'	'HIGH'	4
'GOOD'	'HIGH'	-22
'GOOD'	'HIGH'	19
'GOOD'	'HIGH'	6
'GOOD'	'HIGH'	-6
'GOOD'	'LOW'	-25
'GOOD'	'LOW'	-23
'GOOD'	'LOW'	-28
'GOOD'	'LOW'	-22
'GOOD'	'LOW'	-22
'GOOD'	'LOW'	-10
'GOOD'	'LOW'	-20
'GOOD'	'LOW'	-24
'GOOD'	'LOW'	-24
'GOOD'	'LOW'	-22
'GOOD'	'LOW'	-23
'GOOD'	'LOW'	-19
'GOOD'	'LOW'	-2
'GOOD'	'LOW'	-8
'GOOD'	'LOW'	-17
'GOOD'	'LOW'	-30
'SCIENTIFIC'	'HIGH'	-19
'SCIENTIFIC'	'HIGH'	-24
'SCIENTIFIC'	'HIGH'	-4
'SCIENTIFIC'	'HIGH'	-24
'SCIENTIFIC'	'HIGH'	0
'SCIENTIFIC'	'HIGH'	-4
'SCIENTIFIC'	'HIGH'	5
'SCIENTIFIC'	'HIGH'	-1
'SCIENTIFIC'	'HIGH'	-9
'SCIENTIFIC'	'HIGH'	-5
'SCIENTIFIC'	'HIGH'	-6
'SCIENTIFIC'	'HIGH'	4
'SCIENTIFIC'	'HIGH'	-13
'SCIENTIFIC'	'HIGH'	-1
'SCIENTIFIC'	'HIGH'	-3
'SCIENTIFIC'	'HIGH'	-11
'SCIENTIFIC'	'HIGH'	-6
'SCIENTIFIC'	'HIGH'	-4
'SCIENTIFIC'	'LOW'	6
'SCIENTIFIC'	'LOW'	-5
'SCIENTIFIC'	'LOW'	14
'SCIENTIFIC'	'LOW'	-11
'SCIENTIFIC'	'LOW'	14
'SCIENTIFIC'	'LOW'	-5
'SCIENTIFIC'	'LOW'	-22
'SCIENTIFIC'	'LOW'	7
'SCIENTIFIC'	'LOW'	14
'SCIENTIFIC'	'LOW'	15
'SCIENTIFIC'	'LOW'	-6
'SCIENTIFIC'	'LOW'	9
'SCIENTIFIC'	'LOW'	-5
'NONE'	'HIGH'	-26
'NONE'	'HIGH'	-1
'NONE'	'HIGH'	22
'NONE'	'HIGH'	3
'NONE'	'HIGH'	-26
'NONE'	'HIGH'	4
'NONE'	'HIGH'	-21
'NONE'	'HIGH'	-19
'NONE'	'HIGH'	-12
'NONE'	'HIGH'	9
'NONE'	'HIGH'	-9
'NONE'	'HIGH'	-27
'NONE'	'HIGH'	-10
'NONE'	'HIGH'	-37
'NONE'	'HIGH'	0
'NONE'	'HIGH'	-10
'NONE'	'LOW'	-12
'NONE'	'LOW'	-4
'NONE'	'LOW'	13
'NONE'	'LOW'	-27
'NONE'	'LOW'	-7
'NONE'	'LOW'	-20
'NONE'	'LOW'	-4
'NONE'	'LOW'	-10
'NONE'	'LOW'	-3
'NONE'	'LOW'	-11
'NONE'	'LOW'	2
'NONE'	'LOW'	-9
'NONE'	'LOW'	20
'NONE'	'LOW'	9
'NONE'	'LOW'	-8
'NONE'	'LOW'	8
'NONE'	'LOW'	-6
'NONE'	'LOW'	6




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\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 & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=77018&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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=77018&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=77018&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 Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135







ANOVA Model
R ~ Exp * Inst
names(Intercept)ExpLOWInstNONEInstSCIENTIFICExpLOW:InstNONEExpLOW:InstSCIENTIFIC
means1.357-21.295-11.357-8.30227.79530.162

\begin{tabular}{lllllllll}
\hline
ANOVA Model \tabularnewline
R ~ Exp * Inst \tabularnewline
names & (Intercept) & ExpLOW & InstNONE & InstSCIENTIFIC & ExpLOW:InstNONE & ExpLOW:InstSCIENTIFIC \tabularnewline
means & 1.357 & -21.295 & -11.357 & -8.302 & 27.795 & 30.162 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=77018&T=1

[TABLE]
[ROW][C]ANOVA Model[/C][/ROW]
[ROW][C]R ~ Exp * Inst[/C][/ROW]
[ROW][C]names[/C][C](Intercept)[/C][C]ExpLOW[/C][C]InstNONE[/C][C]InstSCIENTIFIC[/C][C]ExpLOW:InstNONE[/C][C]ExpLOW:InstSCIENTIFIC[/C][/ROW]
[ROW][C]means[/C][C]1.357[/C][C]-21.295[/C][C]-11.357[/C][C]-8.302[/C][C]27.795[/C][C]30.162[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=77018&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=77018&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
R ~ Exp * Inst
names(Intercept)ExpLOWInstNONEInstSCIENTIFICExpLOW:InstNONEExpLOW:InstSCIENTIFIC
means1.357-21.295-11.357-8.30227.79530.162







ANOVA Statistics
DfSum SqMean SqF valuePr(>F)
1
Exp1100.196100.1960.7430.391
Inst1657.743328.8712.440.093
Exp:Inst14278.9742139.48715.8750
Residuals8911994.519134.77

\begin{tabular}{lllllllll}
\hline
ANOVA Statistics \tabularnewline
  & Df & Sum Sq & Mean Sq & F value & Pr(>F) \tabularnewline
 & 1 &  &  &  &  \tabularnewline
Exp & 1 & 100.196 & 100.196 & 0.743 & 0.391 \tabularnewline
Inst & 1 & 657.743 & 328.871 & 2.44 & 0.093 \tabularnewline
Exp:Inst & 1 & 4278.974 & 2139.487 & 15.875 & 0 \tabularnewline
Residuals & 89 & 11994.519 & 134.77 &   &   \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=77018&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][/C][C]1[/C][C][/C][C][/C][C][/C][C][/C][/ROW]
[ROW][C]Exp[/C][C]1[/C][C]100.196[/C][C]100.196[/C][C]0.743[/C][C]0.391[/C][/ROW]
[ROW][C]Inst[/C][C]1[/C][C]657.743[/C][C]328.871[/C][C]2.44[/C][C]0.093[/C][/ROW]
[ROW][C]Exp:Inst[/C][C]1[/C][C]4278.974[/C][C]2139.487[/C][C]15.875[/C][C]0[/C][/ROW]
[ROW][C]Residuals[/C][C]89[/C][C]11994.519[/C][C]134.77[/C][C] [/C][C] [/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=77018&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=77018&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
DfSum SqMean SqF valuePr(>F)
1
Exp1100.196100.1960.7430.391
Inst1657.743328.8712.440.093
Exp:Inst14278.9742139.48715.8750
Residuals8911994.519134.77







Tukey Honest Significant Difference Comparisons
difflwruprp adj
LOW-HIGH-2.054-6.7882.6790.391
NONE-GOOD3.433-3.49810.3640.468
SCIENTIFIC-GOOD6.54-0.54713.6270.077
SCIENTIFIC-NONE3.107-3.7659.9790.53
LOW:GOOD-HIGH:GOOD-21.295-33.669-8.920
HIGH:NONE-HIGH:GOOD-11.357-23.7321.0180.091
LOW:NONE-HIGH:GOOD-4.857-16.9077.1920.848
HIGH:SCIENTIFIC-HIGH:GOOD-8.302-20.3513.7480.347
LOW:SCIENTIFIC-HIGH:GOOD0.566-12.45813.591
HIGH:NONE-LOW:GOOD9.937-2.01821.8930.16
LOW:NONE-LOW:GOOD16.4374.81928.0560.001
HIGH:SCIENTIFIC-LOW:GOOD12.9931.37524.6110.019
LOW:SCIENTIFIC-LOW:GOOD21.8619.23534.4870
LOW:NONE-HIGH:NONE6.5-5.11818.1180.581
HIGH:SCIENTIFIC-HIGH:NONE3.056-8.56314.6740.972
LOW:SCIENTIFIC-HIGH:NONE11.923-0.70324.5490.075
HIGH:SCIENTIFIC-LOW:NONE-3.444-14.7167.8270.948
LOW:SCIENTIFIC-LOW:NONE5.423-6.88417.7310.793
LOW:SCIENTIFIC-HIGH:SCIENTIFIC8.868-3.4421.1750.298

\begin{tabular}{lllllllll}
\hline
Tukey Honest Significant Difference Comparisons \tabularnewline
  & diff & lwr & upr & p adj \tabularnewline
LOW-HIGH & -2.054 & -6.788 & 2.679 & 0.391 \tabularnewline
NONE-GOOD & 3.433 & -3.498 & 10.364 & 0.468 \tabularnewline
SCIENTIFIC-GOOD & 6.54 & -0.547 & 13.627 & 0.077 \tabularnewline
SCIENTIFIC-NONE & 3.107 & -3.765 & 9.979 & 0.53 \tabularnewline
LOW:GOOD-HIGH:GOOD & -21.295 & -33.669 & -8.92 & 0 \tabularnewline
HIGH:NONE-HIGH:GOOD & -11.357 & -23.732 & 1.018 & 0.091 \tabularnewline
LOW:NONE-HIGH:GOOD & -4.857 & -16.907 & 7.192 & 0.848 \tabularnewline
HIGH:SCIENTIFIC-HIGH:GOOD & -8.302 & -20.351 & 3.748 & 0.347 \tabularnewline
LOW:SCIENTIFIC-HIGH:GOOD & 0.566 & -12.458 & 13.59 & 1 \tabularnewline
HIGH:NONE-LOW:GOOD & 9.937 & -2.018 & 21.893 & 0.16 \tabularnewline
LOW:NONE-LOW:GOOD & 16.437 & 4.819 & 28.056 & 0.001 \tabularnewline
HIGH:SCIENTIFIC-LOW:GOOD & 12.993 & 1.375 & 24.611 & 0.019 \tabularnewline
LOW:SCIENTIFIC-LOW:GOOD & 21.861 & 9.235 & 34.487 & 0 \tabularnewline
LOW:NONE-HIGH:NONE & 6.5 & -5.118 & 18.118 & 0.581 \tabularnewline
HIGH:SCIENTIFIC-HIGH:NONE & 3.056 & -8.563 & 14.674 & 0.972 \tabularnewline
LOW:SCIENTIFIC-HIGH:NONE & 11.923 & -0.703 & 24.549 & 0.075 \tabularnewline
HIGH:SCIENTIFIC-LOW:NONE & -3.444 & -14.716 & 7.827 & 0.948 \tabularnewline
LOW:SCIENTIFIC-LOW:NONE & 5.423 & -6.884 & 17.731 & 0.793 \tabularnewline
LOW:SCIENTIFIC-HIGH:SCIENTIFIC & 8.868 & -3.44 & 21.175 & 0.298 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=77018&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]LOW-HIGH[/C][C]-2.054[/C][C]-6.788[/C][C]2.679[/C][C]0.391[/C][/ROW]
[ROW][C]NONE-GOOD[/C][C]3.433[/C][C]-3.498[/C][C]10.364[/C][C]0.468[/C][/ROW]
[ROW][C]SCIENTIFIC-GOOD[/C][C]6.54[/C][C]-0.547[/C][C]13.627[/C][C]0.077[/C][/ROW]
[ROW][C]SCIENTIFIC-NONE[/C][C]3.107[/C][C]-3.765[/C][C]9.979[/C][C]0.53[/C][/ROW]
[ROW][C]LOW:GOOD-HIGH:GOOD[/C][C]-21.295[/C][C]-33.669[/C][C]-8.92[/C][C]0[/C][/ROW]
[ROW][C]HIGH:NONE-HIGH:GOOD[/C][C]-11.357[/C][C]-23.732[/C][C]1.018[/C][C]0.091[/C][/ROW]
[ROW][C]LOW:NONE-HIGH:GOOD[/C][C]-4.857[/C][C]-16.907[/C][C]7.192[/C][C]0.848[/C][/ROW]
[ROW][C]HIGH:SCIENTIFIC-HIGH:GOOD[/C][C]-8.302[/C][C]-20.351[/C][C]3.748[/C][C]0.347[/C][/ROW]
[ROW][C]LOW:SCIENTIFIC-HIGH:GOOD[/C][C]0.566[/C][C]-12.458[/C][C]13.59[/C][C]1[/C][/ROW]
[ROW][C]HIGH:NONE-LOW:GOOD[/C][C]9.937[/C][C]-2.018[/C][C]21.893[/C][C]0.16[/C][/ROW]
[ROW][C]LOW:NONE-LOW:GOOD[/C][C]16.437[/C][C]4.819[/C][C]28.056[/C][C]0.001[/C][/ROW]
[ROW][C]HIGH:SCIENTIFIC-LOW:GOOD[/C][C]12.993[/C][C]1.375[/C][C]24.611[/C][C]0.019[/C][/ROW]
[ROW][C]LOW:SCIENTIFIC-LOW:GOOD[/C][C]21.861[/C][C]9.235[/C][C]34.487[/C][C]0[/C][/ROW]
[ROW][C]LOW:NONE-HIGH:NONE[/C][C]6.5[/C][C]-5.118[/C][C]18.118[/C][C]0.581[/C][/ROW]
[ROW][C]HIGH:SCIENTIFIC-HIGH:NONE[/C][C]3.056[/C][C]-8.563[/C][C]14.674[/C][C]0.972[/C][/ROW]
[ROW][C]LOW:SCIENTIFIC-HIGH:NONE[/C][C]11.923[/C][C]-0.703[/C][C]24.549[/C][C]0.075[/C][/ROW]
[ROW][C]HIGH:SCIENTIFIC-LOW:NONE[/C][C]-3.444[/C][C]-14.716[/C][C]7.827[/C][C]0.948[/C][/ROW]
[ROW][C]LOW:SCIENTIFIC-LOW:NONE[/C][C]5.423[/C][C]-6.884[/C][C]17.731[/C][C]0.793[/C][/ROW]
[ROW][C]LOW:SCIENTIFIC-HIGH:SCIENTIFIC[/C][C]8.868[/C][C]-3.44[/C][C]21.175[/C][C]0.298[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=77018&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=77018&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
difflwruprp adj
LOW-HIGH-2.054-6.7882.6790.391
NONE-GOOD3.433-3.49810.3640.468
SCIENTIFIC-GOOD6.54-0.54713.6270.077
SCIENTIFIC-NONE3.107-3.7659.9790.53
LOW:GOOD-HIGH:GOOD-21.295-33.669-8.920
HIGH:NONE-HIGH:GOOD-11.357-23.7321.0180.091
LOW:NONE-HIGH:GOOD-4.857-16.9077.1920.848
HIGH:SCIENTIFIC-HIGH:GOOD-8.302-20.3513.7480.347
LOW:SCIENTIFIC-HIGH:GOOD0.566-12.45813.591
HIGH:NONE-LOW:GOOD9.937-2.01821.8930.16
LOW:NONE-LOW:GOOD16.4374.81928.0560.001
HIGH:SCIENTIFIC-LOW:GOOD12.9931.37524.6110.019
LOW:SCIENTIFIC-LOW:GOOD21.8619.23534.4870
LOW:NONE-HIGH:NONE6.5-5.11818.1180.581
HIGH:SCIENTIFIC-HIGH:NONE3.056-8.56314.6740.972
LOW:SCIENTIFIC-HIGH:NONE11.923-0.70324.5490.075
HIGH:SCIENTIFIC-LOW:NONE-3.444-14.7167.8270.948
LOW:SCIENTIFIC-LOW:NONE5.423-6.88417.7310.793
LOW:SCIENTIFIC-HIGH:SCIENTIFIC8.868-3.4421.1750.298







Levenes Test for Homogeneity of Variance
DfF valuePr(>F)
Group52.3090.051
89

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=77018&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
DfF valuePr(>F)
Group52.3090.051
89



Parameters (Session):
par1 = 3 ; par2 = 2 ; par3 = 1 ; par4 = TRUE ;
Parameters (R input):
par1 = 3 ; par2 = 2 ; par3 = 1 ; par4 = TRUE ;
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])
mynames<- c(V1, V2, V3)
xdf2<-xdf
names(xdf2)<-mynames
names(xdf)<-c('R', 'A', 'B')
mynames <- c(V1, V2, V3)
if(intercept == FALSE)eval (substitute(lmout<-lm(xdf$R ~ xdf$A * xdf$B- 1, data = xdf), list(xdf=quote(xdf2),R=mynames[1],A=mynames[2],B=mynames[3]) ))else eval(substitute(lmout<-lm(xdf$R ~ xdf$A * xdf$B, data = xdf), list(xdf=quote(xdf2),R=mynames[1],A=mynames[2],B=mynames[3]) ))

oldnames<-names(lmout$coeff)
newnames<-gsub('xdf2\\$', '', oldnames)
(names(lmout$coeff)<-newnames)
(names(lmout$coefficients)<-newnames)

load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'ANOVA Model', length(lmout$coefficients)+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
callstr<-gsub('xdf2\\$', '',as.character(lmout$call$formula))
callstr<-paste(callstr[2], callstr[1], callstr[3])
a<-table.element(a,callstr ,length(lmout$coefficients)+1,TRUE)
a<-table.row.end(a)

a<-table.row.start(a)
a<-table.element(a, 'names',,TRUE)
for(i in 1:length(lmout$coefficients)){
a<-table.element(a, names(lmout$coefficients[i]),,FALSE)
}
a<-table.row.end(a)


a<-table.row.start(a)
a<-table.element(a, 'means',,TRUE)
for(i in 1:length(lmout$coefficients)){
a<-table.element(a, round(lmout$coefficients[i], digits=3),,FALSE)
}
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')

(aov.xdf<-aov(lmout) )
(anova.xdf<-anova(lmout) )
rownames(anova.xdf)<-gsub('xdf2\\$','',rownames(anova.xdf))

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(R ~ A + B, data=xdf, xlab=V2, ylab=V1, main='Boxplots of ANOVA Groups',cex.axis=0.7 )
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$A, xdf$B, xdf$R, 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')
par(mai=c(1,1.5,1,1))
layout(matrix(c(1,2,1,2,3,3,3,3), 2,4))
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(lmout)
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')