R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-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(1,44,1,30,1,46,1,58,1,36,1,52,1,44,1,55,1,41,1,48,1,39,1,39,1,57,1,54,1,60,1,59,1,44,1,56,1,36,1,32,1,60,1,48,1,39,1,37,1,31,1,48,1,55,1,43,1,51,1,48,1,54,1,57,1,30,1,41,1,52,1,46,1,41,1,55,1,39,1,43,1,41,1,46,1,44,1,49,1,62,1,55,1,50,1,45,1,46,2,44,2,34,2,56,2,59,2,50,2,50,2,48,2,46,2,59,2,54,2,44,2,50,2,43,2,39,2,42,2,52,2,41,2,46,2,57,2,43,2,52,2,57,2,59,2,56,2,41,2,50,2,50,2,41,2,48,2,48,2,37,2,55,2,44,2,30,2,44,2,54,2,55,2,59,2,48,2,59,2,37,2,39,2,48,2,36,2,56,2,54,2,52,2,66,2,55,2,60,2,52,2,36,2,32,2,34,2,36,2,56,2,50,2,39,2,52,2,44,2,46,2,44,2,41,2,56,2,54,2,50,2,52,2,55,2,62,2,52,2,57,2,34,2,30,2,39,2,34,2,37,2,52,2,58,2,48,2,46,2,46,2,52,2,50,3,37,3,52,3,54,3,41,3,60,3,56,3,43,3,41,3,55,3,46,3,39,3,48,3,39,3,54,3,48,3,62,3,50,3,56,3,66,3,66,3,48,3,52,3,59,3,50,3,50,3,66,3,58,3,43),dim=c(2,160),dimnames=list(c('MVRBIQ0','MC30VRB'),1:160)) > y <- array(NA,dim=c(2,160),dimnames=list(c('MVRBIQ0','MC30VRB'),1:160)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'FALSE' > par2 = '2' > par1 = '1' > ylab = 'Y Variable Name' > xlab = 'X Variable Name' > main = 'Title Goes Here' > 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]) [1] "MVRBIQ0" > (V2<-dimnames(y)[[1]][cat2]) [1] "MC30VRB" > names(xdf)<-c('Response', 'Treatment') > if(intercept == FALSE) (lmxdf<-lm(Response ~ Treatment - 1, data = xdf) ) else (lmxdf<-lm(Response ~ Treatment, data = xdf) ) Call: lm(formula = Response ~ Treatment - 1, data = xdf) Coefficients: Treatment30 Treatment31 Treatment32 Treatment34 Treatment36 Treatment37 1.500 1.000 1.500 2.000 1.600 2.000 Treatment39 Treatment41 Treatment42 Treatment43 Treatment44 Treatment45 1.800 1.800 2.000 2.000 1.600 1.000 Treatment46 Treatment48 Treatment49 Treatment50 Treatment51 Treatment52 1.700 1.923 1.000 2.167 1.000 2.000 Treatment54 Treatment55 Treatment56 Treatment57 Treatment58 Treatment59 2.000 1.667 2.125 1.600 2.000 2.000 Treatment60 Treatment62 Treatment66 1.750 2.000 2.750 > (aov.xdf<-aov(lmxdf) ) Call: aov(formula = lmxdf) Terms: Treatment Residuals Sum of Squares 570.4353 62.5647 Deg. of Freedom 27 133 Residual standard error: 0.6858656 Estimated effects are balanced > (anova.xdf<-anova(lmxdf) ) Analysis of Variance Table Response: Response Df Sum Sq Mean Sq F value Pr(>F) Treatment 27 570.44 21.1272 44.912 < 2.2e-16 *** Residuals 133 62.56 0.4704 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/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="/var/wessaorg/rcomp/tmp/1mqk71322507477.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="/var/wessaorg/rcomp/tmp/264171322507477.tab") > postscript(file="/var/wessaorg/rcomp/tmp/3lgbf1322507477.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > boxplot(Response ~ Treatment, data=xdf, xlab=V2, ylab=V1) > dev.off() null device 1 > if(intercept==TRUE){ + thsd<-TukeyHSD(aov.xdf) + postscript(file="/var/wessaorg/rcomp/tmp/4rmrq1322507477.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + 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="/var/wessaorg/rcomp/tmp/5q5aw1322507477.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="/var/wessaorg/rcomp/tmp/64yru1322507477.tab") + } > library(car) Loading required package: MASS Loading required package: nnet Loading required package: survival Loading required package: splines > lt.lmxdf<-levene.test(lmxdf) 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,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="/var/wessaorg/rcomp/tmp/7mz7n1322507477.tab") > > try(system("convert tmp/3lgbf1322507477.ps tmp/3lgbf1322507477.png",intern=TRUE)) character(0) > try(system("convert tmp/4rmrq1322507477.ps tmp/4rmrq1322507477.png",intern=TRUE)) convert: unable to open image `tmp/4rmrq1322507477.ps': No such file or directory @ blob.c/OpenBlob/2480. convert: missing an image filename `tmp/4rmrq1322507477.png' @ convert.c/ConvertImageCommand/2838. character(0) Warning message: running command 'convert tmp/4rmrq1322507477.ps tmp/4rmrq1322507477.png' had status 1 > > > proc.time() user system elapsed 0.699 0.095 2.104