Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
3 24 3 17 3 22 3 15 3 20 3 24 3 15 3 20 3 20 3 17 3 11 3 21 3 28 3 14 3 13 3 12 3 21 3 13 3 19 3 23 3 27 3 25 3 27 3 16 3 20 3 18 3 19 3 17 3 10 3 11 3 16 3 13 3 14 3 12 3 15 3 15 3 14 3 10 3 13 3 21 3 11 3 14 3 20 3 7 3 24 3 16 3 22 3 25 3 5 3 19 3 23 3 13 3 10 3 12 3 21 3 22 3 20 3 17 3 13 3 9 3 22 3 15 3 12 3 25 3 14 3 14 3 17 3 9 3 10 3 15 3 15 3 15 3 14 3 21 3 13 3 20 3 16 3 12 3 20 3 18 3 21 3 23 3 13 3 22 3 14 3 23 3 16 3 14 3 22 3 19 3 23 3 16 3 20 3 8 3 16 3 11 3 16 3 10 3 17 3 16 3 17 3 10 3 15 3 13 3 19 3 14 3 18 3 25 3 10 3 22 3 15 3 18 3 22 3 18 3 15 3 20 3 18 3 6 3 17 3 12 3 12 3 19 3 23 3 26 3 28 3 19 3 16 3 3 3 11 3 15 3 22 3 12 3 21 3 25 3 12 3 14 3 24 3 12 3 13 3 15 3 17 3 12 3 28 3 25 3 14 3 21 3 18 3 23 3 16 3 15 3 5 3 19 3 22 3 19 3 12 3 22 3 18 3 24 3 19 3 4 3 20 3 24 3 26 3 22 3 19 3 9 3 22 3 18 3 16 3 19 3 20 3 21 3 17 3 9 3 26 3 28 3 13 3 16 3 22 3 18 3 21 3 10 3 15 3 15 3 13 3 10 3 23 3 21 3 14 3 17 3 15 3 15 3 17 3 26 3 12 3 14 3 26 3 18 3 17 3 20 3 16 3 19 3 12 3 20 3 19 3 25 3 19 3 15 3 12 2 15 2 21 2 30 2 20 2 14 2 18 2 19 2 25 2 23 2 17 2 21 2 21 2 8 2 29 2 20 2 19 2 22 2 23 2 24 2 12 2 22 2 12 2 22 2 20 2 10 2 23 2 17 2 22 2 24 2 18 2 21 2 20 2 20 2 22 2 19 2 20 2 26 2 23 2 24 2 21 2 21 2 19 2 8 2 17 2 20 2 11 2 8 2 15 2 18 2 18 2 19 2 19 2 30 2 17 2 24 2 20 2 25 2 20 2 27 2 18 2 28 2 21 2 27 2 22 2 28 2 25 2 21 2 22 2 28 2 20 2 29 2 20 2 20 2 23 2 18 2 18 2 19 2 25 2 25 2 25 2 24 2 19 2 26 2 10 2 17 2 13 2 17 2 30 2 4 2 16 2 21 2 22 2 20 2 22 2 23 2 16 2 0 2 18 2 25 2 18 2 18 2 24 2 29 2 15 2 22 2 23 2 24 2 22 2 15 2 17 2 20 2 27 2 26 2 23 2 23 2 15 2 26 2 22 2 18 2 15 2 22 2 27 2 10 2 20 2 17 2 23 2 19 2 13 2 27 2 23 2 16 2 25 2 2 2 26 2 20 2 22 2 24 2 18 2 21 2 NA 2 24 2 19 2 24 2 19 2 17 2 20 2 21 2 21 2 21 2 16 2 27 2 15 2 21 2 18 2 22 2 20 2 17 2 21 2 23 2 18 2 22 2 24 2 27 2 20 2 27 2 20 2 20 2 21 2 26 2 25 2 18 2 21 2 16 2 25 2 20 2 27 2 20 2 18 2 26 2 18 2 16 2 18 2 21 2 18 2 25 2 20 2 23 2 22 2 10 2 18 2 25 2 23 2 22 2 23 2 19 2 14 2 26 2 17 2 15 2 21 2 20 2 22 2 20 2 26 2 26 2 20 2 24 2 20 2 15 2 25 2 20 2 27 2 20 2 17 2 22 2 24 2 22 2 16 2 22 2 23 2 19 2 20 2 15 2 22 2 12 2 15 2 27 2 24 2 18 2 18 2 25 2 12 2 19 2 24 2 22 2 18 2 28 1 21 1 26 1 22 1 22 1 18 1 23 1 12 1 20 1 22 1 21 1 19 1 22 1 15 1 20 1 19 1 18 1 20 1 21 1 15 1 16 1 23 1 21 1 18 1 25 1 9 1 23 1 16 1 16 1 19 1 25 1 25 1 18 1 23 1 21 1 10 1 22 1 26 1 23 1 23 1 24 1 24 1 23 1 15 1 16 1 19 1 18 1 27 1 13 1 28 1 23 1 21 1 19 1 19 1 18 1 19 1 17 1 25 1 19 1 26 1 14 1 28 1 16 1 20 1 24 1 23 1 22 1 21 1 25 1 27 1 23 1 23 1 18 1 18 1 23 1 19 1 15 1 20 1 16 1 25 1 25 1 19 1 19 1 16 1 19 1 19 1 23 1 21 1 22 1 19 1 20 1 3 1 23 1 14 1 23 1 20 1 15 1 13 1 16 1 7 1 24 1 17 1 24 1 24 1 19 1 28 1 23 1 19 1 23 1 25 1 25 1 20 1 16 1 20 1 25 1 25 1 23 1 17 1 20 1 16 1 23 1 12 1 24 1 11 1 14 1 23 1 18 1 29 1 16 1 19 1 16 1 23 1 19 1 4 1 20 1 20 1 4 1 24 1 16 1 3 1 24 1 23 1 17 1 20 1 22 1 19 1 24 1 19 1 27 1 22 1 23 1 22 1 17 1 23 1 23 1 28 1 29 1 21 1 24 1 20 1 7 1 19 1 28 1 26 1 19 1 20 1 23 1 24 1 16 1 19 1 24 1 21 1 16 1 16 1 21 1 NA 1 28 1 16 1 23 1 26 1 29 1 18 1 19 1 19 1 16 1 16 1 16 1 18 1 22 1 14 1 20 1 15 1 22 1 16 1 15 1 11 1 15 1 20 1 21 1 16 1 17 1 15 1 16 1 18 1 25 1 20 1 24 1 28 1 22 1 20 1 27 1 17 1 22 1 23 1 22 1 13 1 19 1 15 1 20 1 24 1 18 1 19 1 15 1 20 1 13 1 23 1 24 1 23 1 19 1 20 1 22 1 25 1 26 1 24 1 27 1 16 1 15 1 25 1 27 1 23 1 21 1 14 1 24 1 16 1 22 1 13 1 17 1 23 1 22 1 23 1 26 1 14 1 24 1 21 1 16 1 11 1 NA 1 19 1 16 1 19 1 16 1 11 1 23 1 27 1 23 1 25 1 24 1 22 1 26 1 19 1 19 1 19 1 20 1 16 1 22 1 21 1 26 1 23 1 21 1 22 1 26 1 27 1 17 1 22 1 19 1 14 1 20 1 26
Names of X columns:
Groep_n NUMERACYTOT
Response Variable (column number)
Factor Variable (column number)
Include Intercept Term ?
1e-08
TRUE
FALSE
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
par3 <- 'FALSE' 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