R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-pc-linux-gnu (32-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> x <- array(list('M','G1',547,'M','G1',538,'M','G1',437,'M','G1',534,'M','G1',529,'M','G1',539,'M','G2',590,'M','G2',489,'M','G2',349,'M','G2',458,'M','G2',593,'M','G2',600,'M','G3',586,'M','G3',499,'M','G3',586,'M','G3',538,'M','G3',528,'M','G3',487,'F','G1',278,'F','G1',389,'F','G1',484,'F','G1',530,'F','G1',540,'F','G1',489,'F','G2',599,'F','G2',454,'F','G2',511,'F','G2',489,'F','G2',379,'F','G2',547,'F','G3',432,'F','G3',357,'F','G3',688,'F','G3',538,'F','G3',290,'F','G3',375),dim=c(3,36),dimnames=list(c('Geslacht','Groep','Liters'),1:36))
>  y <- array(NA,dim=c(3,36),dimnames=list(c('Geslacht','Groep','Liters'),1:36))
>  for (i in 1:dim(x)[1])
+  {
+  	for (j in 1:dim(x)[2])
+  	{
+  		y[i,j] <- as.numeric(x[i,j])
+  	}
+  }
There were 50 or more warnings (use warnings() to see the first 50)
> par4 = 'TRUE'
> par3 = '1'
> par2 = '2'
> par1 = '3'
> ylab = 'Y Variable Name'
> xlab = 'X Variable Name'
> main = 'Title Goes Here'
> 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])
[1] "Liters"
> (V2<-dimnames(y)[[1]][cat2])
[1] "Groep"
> (V3 <-dimnames(y)[[1]][cat3])
[1] "Geslacht"
> 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)
[1] "(Intercept)"                 "xdf2$GroepG2"               
[3] "xdf2$GroepG3"                "xdf2$GeslachtM"             
[5] "xdf2$GroepG2:xdf2$GeslachtM" "xdf2$GroepG3:xdf2$GeslachtM"
> (names(lmout$coefficients)<-newnames)
[1] "(Intercept)"                 "xdf2$GroepG2"               
[3] "xdf2$GroepG3"                "xdf2$GeslachtM"             
[5] "xdf2$GroepG2:xdf2$GeslachtM" "xdf2$GroepG3:xdf2$GeslachtM"
> 
> #Note: the /var/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/fisher/rcomp/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, signif(lmout$coefficients[i], digits=5),,FALSE)
+ }
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/fisher/rcomp/tmp/1lnl61355348654.tab") 
> (aov.xdf<-aov(lmout) )
Call:
   aov(formula = lmout)

Terms:
                xdf2$Groep xdf2$Geslacht xdf2$Groep:xdf2$Geslacht Residuals
Sum of Squares     2188.67      31093.44                  8684.22 251719.67
Deg. of Freedom          2             1                        2        30

Residual standard error: 91.60052 
Estimated effects may be unbalanced
> (anova.xdf<-anova(lmout) )
Analysis of Variance Table

Response: xdf2$Liters
                         Df Sum Sq Mean Sq F value  Pr(>F)  
xdf2$Groep                2   2189  1094.3  0.1304 0.87822  
xdf2$Geslacht             1  31093 31093.4  3.7057 0.06375 .
xdf2$Groep:xdf2$Geslacht  2   8684  4342.1  0.5175 0.60124  
Residuals                30 251720  8390.7                  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 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, signif(anova.xdf$'Sum Sq'[i], digits=5),,FALSE)
+ a<-table.element(a, signif(anova.xdf$'Mean Sq'[i], digits=5),,FALSE)
+ a<-table.element(a, signif(anova.xdf$'F value'[i], digits=5),,FALSE)
+ a<-table.element(a, signif(anova.xdf$'Pr(>F)'[i], digits=5),,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, signif(anova.xdf$'Sum Sq'[i+1], digits=5),,FALSE)
> a<-table.element(a, signif(anova.xdf$'Mean Sq'[i+1], digits=5),,FALSE)
> a<-table.element(a, ' ',,FALSE)
> a<-table.element(a, ' ',,FALSE)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/fisher/rcomp/tmp/2z6d01355348654.tab") 
> postscript(file="/var/fisher/rcomp/tmp/334211355348654.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> boxplot(R ~ A + B, data=xdf, xlab=V2, ylab=V1, main='Boxplots of ANOVA Groups')
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/4rgav1355348654.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> 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()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/5zg341355348654.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> interaction.plot(xdf$A, xdf$B, xdf$R, xlab=V2, ylab=V1, trace.label=V3, main='Possible Interactions Between Anova Groups')
> dev.off()
null device 
          1 
> if(intercept==TRUE){
+ thsd<-TukeyHSD(aov.xdf)
+ names(thsd) <- c(V2, V3, paste(V2, ':', V3, sep=''))
+ postscript(file="/var/fisher/rcomp/tmp/6ana01355348654.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
+ layout(matrix(c(0,0,1,2,1,2,0,0,3,3,3,3), 2,6))
+ plot(thsd, las=1)
+ dev.off()
+ }
null device 
          1 
> 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,signif(thsd[[nt]][i,j], digits=5), 1, FALSE)
+ }
+ a<-table.row.end(a)
+ }
+ } # end nt
+ a<-table.end(a)
+ table.save(a,file="/var/fisher/rcomp/tmp/77nfv1355348654.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="/var/fisher/rcomp/tmp/8a78u1355348654.tab") 
+ }
> library(car)
Loading required package: MASS
Loading required package: nnet
> lt.lmxdf<-levene.test(lmout)
Warning message:
'levene.test' is deprecated.
Use 'leveneTest' instead.
See help("Deprecated") and help("car-deprecated"). 
> 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,signif(lt.lmxdf[[i]][1], digits=5), 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="/var/fisher/rcomp/tmp/9b4cm1355348654.tab") 
> 
> try(system("convert tmp/334211355348654.ps tmp/334211355348654.png",intern=TRUE))
character(0)
> try(system("convert tmp/4rgav1355348654.ps tmp/4rgav1355348654.png",intern=TRUE))
character(0)
> try(system("convert tmp/5zg341355348654.ps tmp/5zg341355348654.png",intern=TRUE))
character(0)
> try(system("convert tmp/6ana01355348654.ps tmp/6ana01355348654.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.453   0.731   3.173