R version 2.10.1 (2009-12-14)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

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(0.953159041
+ ,'gp1'
+ ,'s500'
+ ,3.322440087
+ ,'gp1'
+ ,'s500'
+ ,2.46459695
+ ,'gp1'
+ ,'s500'
+ ,3.131808279
+ ,'gp1'
+ ,'s500'
+ ,4.752178649
+ ,'gp1'
+ ,'s500'
+ ,3.376906318
+ ,'gp1'
+ ,'s500'
+ ,2.096949891
+ ,'gp1'
+ ,'s500'
+ ,2.396514161
+ ,'gp1'
+ ,'s500'
+ ,4.234749455
+ ,'gp1'
+ ,'s500'
+ ,3.649237473
+ ,'gp1'
+ ,'s500'
+ ,1.974400871
+ ,'gp1'
+ ,'s500'
+ ,2.532679739
+ ,'gp1'
+ ,'s500'
+ ,3.853485839
+ ,'gp1'
+ ,'s500'
+ ,3.622004357
+ ,'gp1'
+ ,'s500'
+ ,0.672955121
+ ,'gp1'
+ ,'s1000'
+ ,2.867505127
+ ,'gp1'
+ ,'s1000'
+ ,2.074662203
+ ,'gp1'
+ ,'s1000'
+ ,2.619058391
+ ,'gp1'
+ ,'s1000'
+ ,4.642070817
+ ,'gp1'
+ ,'s1000'
+ ,2.297472554
+ ,'gp1'
+ ,'s1000'
+ ,3.064679093
+ ,'gp1'
+ ,'s1000'
+ ,3.399083122
+ ,'gp1'
+ ,'s1000'
+ ,4.822279527
+ ,'gp1'
+ ,'s1000'
+ ,1.821691398
+ ,'gp1'
+ ,'s1000'
+ ,0.475781156
+ ,'gp1'
+ ,'s1000'
+ ,3.133294125
+ ,'gp1'
+ ,'s1000'
+ ,4.341974303
+ ,'gp1'
+ ,'s1000'
+ ,3.411901315
+ ,'gp1'
+ ,'s1000'
+ ,0.885076253
+ ,'gp1'
+ ,'l500'
+ ,3.553921569
+ ,'gp1'
+ ,'l500'
+ ,4.588779956
+ ,'gp1'
+ ,'l500'
+ ,2.791394336
+ ,'gp1'
+ ,'l500'
+ ,5.133442266
+ ,'gp1'
+ ,'l500'
+ ,3.227124183
+ ,'gp1'
+ ,'l500'
+ ,2.369281046
+ ,'gp1'
+ ,'l500'
+ ,2.723311547
+ ,'gp1'
+ ,'l500'
+ ,3.390522876
+ ,'gp1'
+ ,'l500'
+ ,4.08496732
+ ,'gp1'
+ ,'l500'
+ ,1.579520697
+ ,'gp1'
+ ,'l500'
+ ,2.205882353
+ ,'gp1'
+ ,'l500'
+ ,2.696078431
+ ,'gp1'
+ ,'l500'
+ ,3.376906318
+ ,'gp1'
+ ,'l500'
+ ,2.576079744
+ ,'gp1'
+ ,'l1000'
+ ,2.173060683
+ ,'gp1'
+ ,'l1000'
+ ,2.897665581
+ ,'gp1'
+ ,'l1000'
+ ,2.400395102
+ ,'gp1'
+ ,'l1000'
+ ,4.66356014
+ ,'gp1'
+ ,'l1000'
+ ,2.288801424
+ ,'gp1'
+ ,'l1000'
+ ,3.030371577
+ ,'gp1'
+ ,'l1000'
+ ,2.297472554
+ ,'gp1'
+ ,'l1000'
+ ,3.977786826
+ ,'gp1'
+ ,'l1000'
+ ,4.209268307
+ ,'gp1'
+ ,'l1000'
+ ,3.077497286
+ ,'gp1'
+ ,'l1000'
+ ,3.343286283
+ ,'gp1'
+ ,'l1000'
+ ,4.363463626
+ ,'gp1'
+ ,'l1000'
+ ,4.007570274
+ ,'gp1'
+ ,'l1000'
+ ,1.906318083
+ ,'gp2'
+ ,'s500'
+ ,1.879084967
+ ,'gp2'
+ ,'s500'
+ ,3.118191721
+ ,'gp2'
+ ,'s500'
+ ,4.970043573
+ ,'gp2'
+ ,'s500'
+ ,4.793028322
+ ,'gp2'
+ ,'s500'
+ ,3.649237473
+ ,'gp2'
+ ,'s500'
+ ,4.438997821
+ ,'gp2'
+ ,'s500'
+ ,3.322440087
+ ,'gp2'
+ ,'s500'
+ ,3.118191721
+ ,'gp2'
+ ,'s500'
+ ,1.44335512
+ ,'gp2'
+ ,'s500'
+ ,2.696078431
+ ,'gp2'
+ ,'s500'
+ ,2.423747277
+ ,'gp2'
+ ,'s500'
+ ,2.627995643
+ ,'gp2'
+ ,'s500'
+ ,2.46459695
+ ,'gp2'
+ ,'s500'
+ ,1.984557848
+ ,'gp2'
+ ,'s1000'
+ ,2.190402944
+ ,'gp2'
+ ,'s1000'
+ ,2.807561226
+ ,'gp2'
+ ,'s1000'
+ ,2.910483774
+ ,'gp2'
+ ,'s1000'
+ ,2.623205453
+ ,'gp2'
+ ,'s1000'
+ ,2.473157196
+ ,'gp2'
+ ,'s1000'
+ ,2.391723972
+ ,'gp2'
+ ,'s1000'
+ ,6.318238026
+ ,'gp2'
+ ,'s1000'
+ ,2.507464712
+ ,'gp2'
+ ,'s1000'
+ ,1.495958499
+ ,'gp2'
+ ,'s1000'
+ ,1.855998914
+ ,'gp2'
+ ,'s1000'
+ ,4.689196526
+ ,'gp2'
+ ,'s1000'
+ ,2.674478224
+ ,'gp2'
+ ,'s1000'
+ ,3.339139221
+ ,'gp2'
+ ,'s1000'
+ ,1.974400871
+ ,'gp2'
+ ,'l500'
+ ,1.715686275
+ ,'gp2'
+ ,'l500'
+ ,1.252723312
+ ,'gp2'
+ ,'l500'
+ ,3.322440087
+ ,'gp2'
+ ,'l500'
+ ,4.765795207
+ ,'gp2'
+ ,'l500'
+ ,4.452614379
+ ,'gp2'
+ ,'l500'
+ ,2.600762527
+ ,'gp2'
+ ,'l500'
+ ,3.594771242
+ ,'gp2'
+ ,'l500'
+ ,1.811002179
+ ,'gp2'
+ ,'l500'
+ ,3.06372549
+ ,'gp2'
+ ,'l500'
+ ,2.16503268
+ ,'gp2'
+ ,'l500'
+ ,3.690087146
+ ,'gp2'
+ ,'l500'
+ ,2.369281046
+ ,'gp2'
+ ,'l500'
+ ,3.145424837
+ ,'gp2'
+ ,'l500'
+ ,1.727439981
+ ,'gp2'
+ ,'l1000'
+ ,3.15893051
+ ,'gp2'
+ ,'l1000'
+ ,1.868817107
+ ,'gp2'
+ ,'l1000'
+ ,3.304831705
+ ,'gp2'
+ ,'l1000'
+ ,3.082021354
+ ,'gp2'
+ ,'l1000'
+ ,2.614534323
+ ,'gp2'
+ ,'l1000'
+ ,2.078809265
+ ,'gp2'
+ ,'l1000'
+ ,3.95215044
+ ,'gp2'
+ ,'l1000'
+ ,1.637335626
+ ,'gp2'
+ ,'l1000'
+ ,1.230169502
+ ,'gp2'
+ ,'l1000'
+ ,2.065991073
+ ,'gp2'
+ ,'l1000'
+ ,2.108969719
+ ,'gp2'
+ ,'l1000'
+ ,2.353269393
+ ,'gp2'
+ ,'l1000'
+ ,2.297472554
+ ,'gp2'
+ ,'l1000')
+ ,dim=c(3
+ ,112)
+ ,dimnames=list(c('accuracy'
+ ,'group'
+ ,'cond')
+ ,1:112))
>  y <- array(NA,dim=c(3,112),dimnames=list(c('accuracy','group','cond'),1:112))
>  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 = '3'
> par2 = '2'
> par1 = '1'
> ylab = 'Y-Axis Title'
> xlab = 'X-Axis Title'
> main = 'Chart Title'
> 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] "accuracy"
> (V2<-dimnames(y)[[1]][cat2])
[1] "group"
> (V3 <-dimnames(y)[[1]][cat3])
[1] "cond"
> 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$groupgp2"               
[3] "xdf2$condl500"                "xdf2$conds1000"              
[5] "xdf2$conds500"                "xdf2$groupgp2:xdf2$condl500" 
[7] "xdf2$groupgp2:xdf2$conds1000" "xdf2$groupgp2:xdf2$conds500" 
> (names(lmout$coefficients)<-newnames)
[1] "(Intercept)"                  "xdf2$groupgp2"               
[3] "xdf2$condl500"                "xdf2$conds1000"              
[5] "xdf2$conds500"                "xdf2$groupgp2:xdf2$condl500" 
[7] "xdf2$groupgp2:xdf2$conds1000" "xdf2$groupgp2:xdf2$conds500" 
> 
> #Note: the /var/yougetitorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/yougetitorg/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/yougetitorg/rcomp/tmp/150wv1303234350.tab") 
> (aov.xdf<-aov(lmout) )
Call:
   aov(formula = lmout)

Terms:
                xdf2$group xdf2$cond xdf2$group:xdf2$cond Residuals
Sum of Squares     1.60337   0.88202              3.67021 119.23011
Deg. of Freedom          1         3                    3       104

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

Response: xdf2$accuracy
                      Df  Sum Sq Mean Sq F value Pr(>F)
xdf2$group             1   1.603 1.60337  1.3986 0.2397
xdf2$cond              3   0.882 0.29401  0.2565 0.8566
xdf2$group:xdf2$cond   3   3.670 1.22340  1.0671 0.3664
Residuals            104 119.230 1.14644               
> 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/yougetitorg/rcomp/tmp/2chyw1303234350.tab") 
> postscript(file="/var/yougetitorg/rcomp/tmp/3n4x01303234350.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/yougetitorg/rcomp/tmp/4njap1303234350.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/yougetitorg/rcomp/tmp/5jkl71303234350.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/yougetitorg/rcomp/tmp/6fiim1303234350.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/yougetitorg/rcomp/tmp/7w4tx1303234350.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/yougetitorg/rcomp/tmp/8w8yh1303234350.tab") 
+ }
> library(car)
Loading required package: MASS
Loading required package: nnet
Loading required package: survival
Loading required package: splines
> 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/yougetitorg/rcomp/tmp/9wfnh1303234350.tab") 
> 
> try(system("convert tmp/3n4x01303234350.ps tmp/3n4x01303234350.png",intern=TRUE))
character(0)
> try(system("convert tmp/4njap1303234350.ps tmp/4njap1303234350.png",intern=TRUE))
character(0)
> try(system("convert tmp/5jkl71303234350.ps tmp/5jkl71303234350.png",intern=TRUE))
character(0)
> try(system("convert tmp/6fiim1303234350.ps tmp/6fiim1303234350.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.460   0.790   1.812