Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
41 F' 39 44 F' 44 44 F' 48 45 F' 45 45 F' 45 45 F' 47 47 F' 47 47 F' 50 47 F' 47 49 F' 50 50 F' 50 50 F' 50 50 F' 50 50 F' 51 50 F' 53 51 F' 58 51 F' 52 51 F' 53 51 F' 52 51 F' 51 52 F' 51 52 F' 52 52 F' 53 52 F' 52 52 F' 53 53 F' 54 53 F' 53 53 F' 56 53 F' 53 53 F' 52 54 F' 53 54 F' 56 54 F' 55 54 F' 54 54 F' 58 55 F' 54 55 F' 50 55 F' 60 55 F' 55 55 F' 57 55 F' 59 55 F' 55 55 F' 53 55 F' 54 55 F' 55 55 F' 55 55 F' 59 55 F' 60 56 F' 54 56 F' 56 56 F' 57 56 F' 56 56 F' 57 56 F' 56 57 F' 52 57 F' 56 57 F' 56 57 F' 56 57 F' 55 58 F' 57 58 F' 59 59 F' 59 59 F' 54 59 F' 63 59 F' 60 59 F' 57 59 F' 59 59 F' 63 59 F' 59 60 F' 58 60 F' 61 61 F' 61 61 F' 61 61 F' 63 61 F' 62 61 F' 62 61 F' 62 62 F' 62 62 F' 64 62 F' 62 62 F' 64 62 F' 63 63 F' 68 63 F' 64 63 F' 62 64 F' 64 64 F' 63 64 F' 64 66 F' 65 66 F' 64 66 F' 66 67 F' 70 68 F' 68 68 F' 68 71 F' 71 75 F' 75 75 F' 78 77 F' 76 56 M' 55 58 M' 57 58 M' 54 58 M' 56 61 M' 59 61 M' 61 63 M' 63 64 M' 65 64 M' 65 64 M' 62 65 M' 64 66 M' 65 66 M' 66 66 M' 62 66 M' 65 67 M' 65 67 M' 67 67 M' 66 68 M' 70 68 M' 66 68 M' 69 68 M' 71 68 M' 68 69 M' 69 69 M' 68 70 M' 68 70 M' 66 70 M' 75 70 M' 69 70 M' 69 70 M' 70 71 M' 71 71 M' 71 71 M' 69 71 M' 74 73 M' 69 73 M' 75 73 M' 74 73 M' 69 74 M' 73 75 M' 70 75 M' 76 75 M' 76 75 M' 76 76 M' 76 76 M' 79 76 M' 71 76 M' 80 76 M' 75 77 M' 77 77 M' 76 77 M' 78 78 M' 80 79 M' 79 80 M' 80 80 M' 78 80 M' 83 81 M' 79 82 M' 85 82 M' 81 83 M' 83 83 M' 85 84 M' 83 85 M' 82 86 M' 88 86 M' 84 86 M' 88 86 M' 89 87 M' 88 89 M' 87 90 M' 84 91 M' 90 91 M' 90 93 M' 88 94 M' 96 95 M' 96 98 M' 97 100 M' 101 101 M' 92 101 M' 103 107 M' 102 124 M' 119
Names of X columns:
repwt sex weight
Response : Variable 1
Factor : Variable 2
Factor : Variable 3
Include Intercept Term ?
TRUE
TRUE
FALSE
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
par4 <- 'TRUE' par3 <- '3' par2 <- '2' par1 <- '1' cat1 <- as.numeric(par1) # cat2<- as.numeric(par2) # cat3 <- as.numeric(par3) intercept<-as.logical(par4) x <- t(x) x1<-as.numeric(x[,cat1]) f1<-as.character(x[,cat2]) f2 <- as.character(x[,cat3]) xdf<-data.frame(x1,f1, f2) (V1<-dimnames(y)[[1]][cat1]) (V2<-dimnames(y)[[1]][cat2]) (V3 <-dimnames(y)[[1]][cat3]) names(xdf)<-c('Response', 'Treatment_A', 'Treatment_B') if(intercept == FALSE) (lmxdf<-lm(Response ~ Treatment_A * Treatment_B- 1, data = xdf) ) else (lmxdf<-lm(Response ~ Treatment_A * Treatment_B, 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, lmxdf$call['formula'],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) for(i in 1 : length(rownames(anova.xdf))-1){ a<-table.row.start(a) a<-table.element(a,rownames(anova.xdf)[i] ,,TRUE) a<-table.element(a, anova.xdf$Df[1],,FALSE) a<-table.element(a, round(anova.xdf$'Sum Sq'[i], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'Mean Sq'[i], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'F value'[i], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'Pr(>F)'[i], 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'[i+1],,FALSE) a<-table.element(a, round(anova.xdf$'Sum Sq'[i+1], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'Mean Sq'[i+1], 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_A + Treatment_B, data=xdf, xlab=V2, ylab=V1, main='Boxplots of ANOVA Groups') dev.off() bitmap(file='designplot.png') xdf2 <- xdf # to preserve xdf make copy for function names(xdf2) <- c(V1, V2, V3) plot.design(xdf2, main='Design Plot of Group Means') dev.off() bitmap(file='interactionplot.png') interaction.plot(xdf$Treatment_A, xdf$Treatment_B, xdf$Response, xlab=V2, ylab=V1, trace.label=V3, main='Possible Interactions Between Anova Groups') dev.off() if(intercept==TRUE){ thsd<-TukeyHSD(aov.xdf) names(thsd) <- c(V2, V3, paste(V2, ':', V3, sep='')) bitmap(file='TukeyHSDPlot.png') layout(matrix(c(1,2,3,3), 2,2)) plot(thsd, las=1) dev.off() } if(intercept==TRUE){ ntables<-length(names(thsd)) 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(nt in 1:ntables){ for(i in 1:length(rownames(thsd[[nt]]))){ a<-table.row.start(a) a<-table.element(a,rownames(thsd[[nt]])[i], 1, TRUE) for(j in 1:4){ a<-table.element(a,round(thsd[[nt]][i,j], digits=3), 1, FALSE) } a<-table.row.end(a) } } # end nt a<-table.end(a) table.save(a,file='hsdtable.tab') }#end if hsd tables 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<-levene.test(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