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