Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
12.9 50 12.2 68 12.8 62 7.4 54 6.7 71 12.6 54 14.8 65 13.3 73 11.1 52 8.2 84 11.4 42 6.4 66 10.6 65 12 78 6.3 73 11.3 75 11.9 72 9.3 66 9.6 70 10 61 6.4 81 13.8 71 10.8 69 13.8 71 11.7 72 10.9 68 16.1 70 13.4 68 9.9 61 11.5 67 8.3 76 11.7 70 9 60 9.7 77 10.8 72 10.3 69 10.4 71 12.7 62 9.3 70 11.8 64 5.9 58 11.4 76 13 52 10.8 59 12.3 68 11.3 76 11.8 65 7.9 67 12.7 59 12.3 69 11.6 76 6.7 63 10.9 75 12.1 63 13.3 60 10.1 73 5.7 63 14.3 70 8 75 13.3 66 9.3 63 12.5 63 7.6 64 15.9 70 9.2 75 9.1 61 11.1 60 13 62 14.5 73 12.2 61 12.3 66 11.4 64 8.8 59 14.6 64 12.6 60 13 66 12.6 78 13.2 53 9.9 67 7.7 59 10.5 66 13.4 68 10.9 71 4.3 66 10.3 73 11.8 72 11.2 71 11.4 59 8.6 64 13.2 66 12.6 78 5.6 68 9.9 73 8.8 62 7.7 65 9 68 7.3 65 11.4 60 13.6 71 7.9 65 10.7 68 10.3 64 8.3 74 9.6 69 14.2 76 8.5 68 13.5 72 4.9 67 6.4 63 9.6 59 11.6 73 11.1 66 4.35 62 12.7 69 18.1 66 17.85 51 16.6 56 12.6 67 17.1 69 19.1 57 16.1 56 13.35 55 18.4 63 14.7 67 10.6 65 12.6 47 16.2 76 13.6 64 18.9 68 14.1 64 14.5 65 16.15 71 14.75 63 14.8 60 12.45 68 12.65 72 17.35 70 8.6 61 18.4 61 16.1 62 11.6 71 17.75 71 15.25 51 17.65 56 16.35 70 17.65 73 13.6 76 14.35 59 14.75 68 18.25 48 9.9 52 16 59 18.25 60 16.85 59 14.6 57 13.85 79 18.95 60 15.6 60 14.85 59 11.75 62 18.45 59 15.9 61 17.1 71 16.1 57 19.9 66 10.95 63 18.45 69 15.1 58 15 59 11.35 48 15.95 66 18.1 73 14.6 67 15.4 61 15.4 68 17.6 75 13.35 62 19.1 69 15.35 58 7.6 60 13.4 74 13.9 55 19.1 62 15.25 63 12.9 69 16.1 58 17.35 58 13.15 68 12.15 72 12.6 62 10.35 62 15.4 65 9.6 69 18.2 66 13.6 72 14.85 62 14.75 75 14.1 58 14.9 66 16.25 55 19.25 47 13.6 72 13.6 62 15.65 64 12.75 64 14.6 19 9.85 50 12.65 68 19.2 70 16.6 79 11.2 69 15.25 71 11.9 48 13.2 66 16.35 73 12.4 74 15.85 66 18.15 71 11.15 74 15.65 78 17.75 75 7.65 53 12.35 60 15.6 50 19.3 70 15.2 69 17.1 65 15.6 78 18.4 78 19.05 59 18.55 72 19.1 70 13.1 63 12.85 63 9.5 71 4.5 74 11.85 67 13.6 66 11.7 62 12.4 80 13.35 73 11.4 67 14.9 61 19.9 73 11.2 74 14.6 32 17.6 69 14.05 69 16.1 84 13.35 64 11.85 58 11.95 60 14.75 59 15.15 78 13.2 57 16.85 60 7.85 68 7.7 68 12.6 73 7.85 69 10.95 67 12.35 60 9.95 65 14.9 66 16.65 74 13.4 81 13.95 72 15.7 55 16.85 49 10.95 74 15.35 53 12.2 64 15.1 65 17.75 57 15.2 51 14.6 80 16.65 67 8.1 70
Names of X columns:
TOT AMS.E
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
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