Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1339.794899 'Control' 1378.629534 'Control' 1203.873678 'Control' 1300.960265 'Control' 1320.377582 'Control' 1135.843276 'Routine 1' 1090.409545 'Routine 1' 1135.843276 'Routine 1' 1044.975814 'Routine 1' 984.397506 'Routine 1' 995.3362922 'Routine 2' 874.6894689 'Routine 2' 889.7703218 'Routine 2' 1040.578851 'Routine 2' 1100.902263 'Routine 2' 1036.613358 'Routine 3' 1105.720915 'Routine 3' 1053.890247 'Routine 3' 1105.720915 'Routine 3' 1036.613358 'Routine 3' 794.529785 'Routine 4' 844.1878966 'Routine 4' 769.7007293 'Routine 4' 732.4571456 'Routine 4' 831.7733687 'Routine 4' 1107.279084 'Routine 5' 1090.502128 'Routine 5' 1023.394305 'Routine 5' 1107.279084 'Routine 5' 973.0634371 'Routine 5' 332.7780898 'Routine 6' 371.9284533 'Routine 6' 391.503635 'Routine 6' 548.105089 'Routine 6' 430.6539985 'Routine 6' 1417.39857 'Routine 7' 1019.835556 'Routine 7' 1348.257176 'Routine 7' 1417.39857 'Routine 7' 1348.257176 'Routine 7' 1279.115783 'Routine 8' 1222.684204 'Routine 8' 1203.873678 'Routine 8' 1373.168414 'Routine 8' 1203.873678 'Routine 8' 774.5510794 'Routine 9' 774.5510794 'Routine 9' 840.9411719 'Routine 9' 708.1609869 'Routine 9' 774.5510794 'Routine 9' 730.9233043 'Routine 10' 687.9278158 'Routine 10' 644.9323273 'Routine 10' 752.4210485 'Routine 10' 816.9142813 'Routine 10' 891.7582797 'Routine 11' 1337.63742 'Routine 11' 1374.794015 'Routine 11' 966.0714697 'Routine 11' 928.9148747 'Routine 11' 1174.510905 'Routine 12' 1211.214371 'Routine 12' 1358.028234 'Routine 12' 1247.917837 'Routine 12' 1358.028234 'Routine 12' 917.5866446 'Routine 13' 1101.103973 'Routine 13' 1192.862638 'Routine 13' 789.1245143 'Routine 13' 990.9935761 'Routine 13' 1254.035081 'Routine 14' 1343.609015 'Routine 14' 1307.779442 'Routine 14' 1236.120294 'Routine 14' 1039.057638 'Routine 14' 1169.367744 'Routine 15' 958.498151 'Routine 15' 862.6483359 'Routine 15' 1016.00804 'Routine 15' 939.328188 'Routine 15' 1041.085728 'Routine 16' 1003.228065 'Routine 16' 908.5839077 'Routine 16' 1116.801053 'Routine 16' 776.0820878 'Routine 16'
Names of X columns:
Response Treatment
Response Variable (column number)
Factor Variable (column number)
Include Intercept Term ?
TRUE
TRUE
FALSE
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
par3 <- 'TRUE' par2 <- '2' par1 <- '1' 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]) (V2<-dimnames(y)[[1]][cat2]) names(xdf)<-c('Response', 'Treatment') if(intercept == FALSE) (lmxdf<-lm(Response ~ Treatment - 1, data = xdf) ) else (lmxdf<-lm(Response ~ Treatment, 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, 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='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) 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='mytable1.tab') bitmap(file='anovaplot.png') boxplot(Response ~ Treatment, data=xdf, xlab=V2, ylab=V1) dev.off() if(intercept==TRUE){ 'Tukey Plot' thsd<-TukeyHSD(aov.xdf) bitmap(file='TukeyHSDPlot.png') 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='mytable2.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='mytable2.tab') } library(car) lt.lmxdf<-leveneTest(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')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation