Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
21 4 26 9 22 4 22 5 18 4 23 4 12 9 20 8 22 11 21 4 19 4 22 6 15 4 20 8 19 4 18 4 15 11 20 4 21 4 21 6 15 6 16 4 23 8 21 5 18 4 25 9 9 4 30 7 20 10 23 4 16 4 16 7 19 12 25 4 25 7 18 5 23 8 21 5 10 4 14 9 22 7 26 4 23 4 23 4 24 4 24 4 18 7 23 4 15 7 19 4 16 4 25 4 23 4 17 8 19 4 21 4 18 4 27 4 21 7 13 12 8 4 29 4 28 4 23 5 21 15 19 5 19 10 20 9 18 8 19 4 17 5 19 4 25 9 19 4 22 10 23 4 26 7 14 4 28 6 16 7 24 5 20 4 12 4 24 4 22 4 12 4 22 4 20 6 10 10 23 7 17 4 22 4 24 7 18 4 21 8 20 11 20 6 22 14 19 5 20 4 26 8 23 9 24 4 21 4 21 5 19 4 8 5 17 4 20 4 11 7 8 10 15 4 18 5 18 4 19 4 19 4 23 6 22 4 21 8 25 5 30 4 17 17 27 4 23 4 23 8 18 4 18 7 23 4 19 4 15 5 20 7 16 4 24 4 25 7 25 11 19 7 19 4 16 4 19 4 19 4 23 4 21 4 22 6 19 8 20 23 20 4 3 8 23 6 14 4 23 4 20 7 15 4 13 4 16 4 7 4 24 10 17 6 24 5 24 5 19 4 25 4 20 5 28 5 23 5 27 5 18 4 28 6 21 4 19 4 23 4 27 9 22 18 28 6 25 5 21 4 22 11 28 4 20 10 29 6 25 8 25 8 20 6 20 8 16 4 20 4 20 9 23 9 18 5 25 4 18 4 19 15 25 10 25 9 25 7 24 9 19 6 26 4 10 7 17 4 13 7 17 4 30 15 25 4 4 9 16 4 21 4 23 28 22 4 17 4 20 4 20 5 22 4 16 4 23 12 16 5 0 4 18 6 25 6 23 5 12 4 18 4 24 4 11 10 18 7 14 4 23 4 24 7 29 4 18 4 15 12 29 5 16 8 19 6 22 17 16 4 23 5 23 4 19 5 4 5 20 6 24 4 20 4 4 4 24 6 22 8 16 10 3 4 15 5 24 4 17 4 20 4 27 16 23 4 26 7 23 4 17 4 20 14 22 5 19 5 24 5 19 5 23 7 15 19 27 16 26 4 22 4 22 7 18 9 15 5 22 14 27 4 10 16 20 10 17 5 23 6 19 4 13 4 27 4 23 5 16 4 25 4 2 5 26 4 20 4 23 5 22 8 24 15 22 7 17 5 23 8 23 8 28 5 29 4 21 4 24 11 20 5 7 22 19 4 28 4 18 4 26 5 21 4 19 16 20 5 23 6 24 5 16 4 19 4 24 4 21 7 16 4 16 8 21 7 28 6 16 5 23 8 26 8 29 4 18 7 19 4 19 13 16 4 16 4 16 4 18 4 22 7 14 5 20 4 15 5 22 12 24 8 16 4 19 4 24 8 19 5 15 4 11 4 15 7 17 5 20 13 21 4 16 4 17 4 20 6 15 4 21 4 16 4 18 4 25 4 21 5 21 6 16 4 20 4 24 4 28 6 27 9 22 5 20 6 27 13 17 4 22 7 23 5 15 4 22 4 13 4 21 6 18 6 22 8 19 6 15 5 20 9 17 6 21 4 23 9 20 4 18 4 22 4 24 5 24 4 18 4 27 4 19 5 20 5 15 8 20 4 27 4 20 9 20 4 13 4 21 4 23 4 26 4 24 4 25 4 18 4 21 4 23 4 16 4 19 4 20 4 25 5 22 8 20 7 25 4 27 4 20 4 18 5 26 5 26 6 24 12 27 5 16 9 15 12 25 4 27 16 18 4 16 5 18 4 23 4 21 6 21 4 14 4 24 5 18 6 16 5 25 6 22 4 13 4 20 7 17 9 23 5 22 5 23 4 22 4 23 12 10 4 18 6 25 9 26 4 14 5 23 4 22 4 23 4 19 4 14 4 26 11 24 4 21 6 17 4 16 5 15 4 11 4 19 4 21 7 20 9 16 5 19 14 16 4 11 4 22 4 20 5 26 4 26 4 20 9 24 4 20 4 15 10 23 4 25 4 27 6 23 4 20 9 25 5 24 4 22 5 27 14 20 9 17 4 22 4 26 17 19 4 19 5 24 9 22 7 16 4 22 5 23 7 19 10 20 5 16 4 19 8 20 4 15 4 22 6 26 5
Names of X columns:
NUMERACYTOT AMS.A
Response Variable (column number)
Factor Variable (column number)
Include Intercept Term ?
FALSE
TRUE
FALSE
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation