Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
26 50 0 57 62 1 37 54 0 67 71 1 43 54 1 52 65 1 52 73 0 43 52 1 84 84 1 67 42 1 49 66 1 70 65 1 52 78 1 58 73 0 68 75 0 62 72 0 43 66 1 56 70 0 56 61 1 74 81 0 65 71 1 63 69 1 58 71 0 57 72 1 63 68 1 53 70 1 57 68 1 51 61 0 64 67 1 53 76 0 29 70 0 54 60 0 58 72 1 43 69 1 51 71 1 53 62 1 54 70 0 56 64 1 61 58 1 47 76 0 39 52 1 48 59 1 50 68 1 35 76 1 30 65 1 68 67 0 49 59 1 61 69 1 67 76 0 47 63 1 56 75 1 50 63 1 43 60 1 67 73 1 62 63 1 57 70 1 41 75 0 54 66 1 45 63 0 48 63 1 61 64 1 56 70 0 41 75 0 43 61 1 53 60 0 44 62 1 66 73 0 58 61 1 46 66 1 37 64 0 51 59 0 51 64 0 56 60 0 66 56 1 37 78 0 59 53 1 42 67 0 38 59 1 66 66 0 34 68 0 53 71 1 49 66 0 55 73 0 49 72 0 59 71 1 40 59 0 58 64 1 60 66 1 63 78 0 56 68 0 54 73 0 52 62 1 34 65 1 69 68 1 32 65 0 48 60 1 67 71 0 58 65 1 57 68 1 42 64 1 64 74 1 58 69 1 66 76 0 26 68 1 61 72 1 52 67 1 51 63 0 55 59 0 50 73 0 60 66 0 56 62 0 63 69 0 61 66 1 52 51 1 16 56 1 46 67 1 56 69 1 52 57 0 55 56 1 50 55 1 59 63 0 60 67 1 52 65 0 44 47 0 67 76 1 52 64 1 55 68 1 37 64 1 54 65 1 72 71 1 51 63 1 48 60 1 60 68 0 50 72 1 63 70 1 33 61 1 67 61 1 46 62 1 54 71 1 59 71 0 61 51 1 33 56 1 47 70 1 69 73 1 52 76 1 55 68 0 41 48 0 73 52 1 52 60 0 50 59 0 51 57 1 60 79 0 56 60 1 56 60 1 29 59 0 66 62 1 66 59 1 73 61 1 55 71 0 64 57 0 40 66 0 46 63 0 58 69 1 43 58 0 61 59 1 51 48 0 50 66 1 52 73 0 54 67 1 66 61 0 61 68 0 80 75 1 51 62 0 56 69 1 56 58 1 56 60 1 53 74 1 47 55 1 25 62 0 47 63 1 46 69 0 50 58 0 39 58 0 51 68 1 58 72 0 35 62 1 58 62 0 60 65 0 62 69 0 63 66 0 53 72 1 46 62 1 67 75 1 59 58 1 64 66 0 38 55 0 50 47 1 48 72 0 48 62 0 47 64 0 66 64 0 47 19 1 63 50 1 58 68 0 44 70 0 51 79 1 43 69 0 55 71 1 38 48 1 45 73 0 50 74 1 54 66 1 57 71 1 60 74 0 55 78 0 56 75 0 49 53 1 37 60 1 59 70 1 46 69 1 51 65 0 58 78 0 64 78 0 53 59 1 48 72 1 51 70 0 47 63 0 59 63 0 62 71 1 62 74 1 51 67 0 64 66 0 52 62 0 67 80 1 50 73 1 54 67 1 58 61 1 56 73 0 63 74 1 31 32 1 65 69 1 71 69 0 50 84 0 57 64 1 47 58 0 47 59 1 57 78 1 43 57 0 41 60 1 63 68 0 63 68 1 56 73 1 51 69 0 50 67 1 22 60 0 41 65 1 59 66 0 56 74 1 66 81 0 53 72 0 42 55 1 52 49 1 54 74 0 44 53 1 62 64 1 53 65 0 50 57 1 36 51 0 76 80 0 66 67 1 62 70 1 59 74 0 47 75 1 55 70 0 58 69 0 60 65 1 44 55 0 57 71 0 45 65 1 58 69 1 51 48 1 57 69 0 30 68 1 46 74 1 51 67 1 56 65 1 58 63 0 44 74 1 14 39 0 53 68 0 42 69 1 49 68 0 44 63 1 62 67 0 30 70 0 46 68 1 50 70 1 54 78 1 48 59 0 55 62 0 35 75 0 55 74 1 41 73 0 59 62 1 54 69 1 55 67 1 45 73 0 51 52 1 47 61 0 42 53 1 53 63 0 53 78 0 41 65 0 55 77 0 55 69 0 46 68 0 63 76 1 43 63 1 65 41 1 59 76 0 39 67 0 44 69 0 60 59 0 57 73 0 67 72 1 52 52 1 52 65 1 69 63 1 46 78 0 46 56 1 53 68 0 40 56 1 70 64 1 54 68 0 77 75 1 45 67 0 60 55 0 47 73 0 50 66 0 66 75 0 60 77 0 41 65 1 53 75 0 34 57 0 51 61 1 69 71 1 60 72 1 45 62 1 58 66 0 39 66 1 51 63 1 52 60 0 49 64 0 63 74 0 44 59 0 51 71 1 52 69 0 60 63 0 53 73 0 53 55 0 52 77 0 31 70 0 51 64 1 65 78 1 51 60 1 49 66 0 61 77 0 58 68 1 62 78 0 54 68 1 52 60 1 72 65 1 50 64 1 65 69 1 53 72 0 56 50 0 63 72 0 62 71 0 66 80 0 50 74 1 45 64 0 58 69 0 52 76 1 53 75 0 68 79 0 59 73 1 58 60 0 52 76 1 45 55 1 58 53 0 70 62 1 69 69 0 71 78 1 46 68 0 58 67 1 39 75 1 46 59 1 64 73 0 67 70 1 44 59 1 54 64 0 41 63 1 68 67 1 63 58 1 57 71 1 61 79 0 39 53 1 69 76 0 64 66 0 38 64 1 59 57 0 51 67 1 59 72 1 51 58 0 65 74 0 47 57 1 50 62 1 57 74 1 21 54 1 47 62 0 51 66 1 37 64 1 67 74 1 43 71 1 58 66 0 51 66 0 40 63 1 41 65 0 58 70 1 64 66 1 64 66 0 58 78 1 50 77 0 59 72 1 55 65 0 59 67 0 58 72 1 41 58 1 56 84 1 63 67 1 77 84 0 60 58 0 58 63 1 64 75 0 46 72 0 62 58 1 60 69 1 50 54 1 46 58 1 44 67 1 58 77 1 56 80 1 43 67 1 54 75 0 54 71 1 56 72 0 65 75 0 66 79 1 62 76 1 58 72 1 67 81 1 25 52 1 56 76 1 53 60 1 56 72 0 59 77 1 46 64 1 49 67 1 56 72 0 76 79 1 33 40 1 49 71 1 53 73 1 58 75 0 72 70 1 51 66 1 42 66 0 69 73 0 51 74 1 54 58 1 52 51 1 59 75 1 51 70 0 67 50 1 64 64 0 58 77 1 53 71 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 1 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 1 NA NA 1 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 1 NA NA 1 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 1 NA NA 0 NA NA 1 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 1 NA NA 0 NA NA 1 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 1 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0 NA NA 0
Names of X columns:
AMS.I AMS.E gender
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 <- 'FALSE' par3 <- '2' par2 <- '1' par1 <- '3' 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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation