Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
'F' 2011 12.9 'F' 2011 7.4 'M' 2011 12.2 'F' 2011 12.8 'M' 2011 7.4 'M' 2011 6.7 'M' 2011 12.6 'F' 2011 14.8 'M' 2011 13.3 'M' 2011 11.1 'M' 2011 8.2 'M' 2011 11.4 'M' 2011 6.4 'M' 2011 10.6 'F' 2011 12.0 'F' 2011 6.3 'F' 2011 11.3 'M' 2011 11.9 'F' 2011 9.3 'M' 2011 9.6 'F' 2011 10.0 'M' 2011 6.4 'M' 2011 13.8 'F' 2011 10.8 'M' 2011 13.8 'M' 2011 11.7 'M' 2011 10.9 'M' 2011 16.1 'F' 2011 13.4 'M' 2011 9.9 'F' 2011 11.5 'F' 2011 8.3 'F' 2011 11.7 'M' 2011 6.1 'M' 2011 9.0 'M' 2011 9.7 'M' 2011 10.8 'M' 2011 10.3 'F' 2011 10.4 'M' 2011 12.7 'M' 2011 9.3 'F' 2011 11.8 'M' 2011 5.9 'M' 2011 11.4 'M' 2011 13.0 'M' 2011 10.8 'M' 2011 12.3 'F' 2011 11.3 'M' 2011 11.8 'M' 2011 7.9 'F' 2011 12.7 'M' 2011 12.3 'M' 2011 11.6 'M' 2011 6.7 'M' 2011 10.9 'M' 2011 12.1 'M' 2011 13.3 'M' 2011 10.1 'F' 2011 5.7 'M' 2011 14.3 'F' 2011 8.0 'M' 2011 13.3 'M' 2011 9.3 'F' 2011 12.5 'F' 2011 7.6 'M' 2011 15.9 'F' 2011 9.2 'M' 2011 9.1 'F' 2011 11.1 'M' 2011 13.0 'M' 2011 14.5 'F' 2011 12.2 'F' 2011 12.3 'F' 2011 11.4 'F' 2011 8.8 'M' 2011 14.6 'M' 2011 7.3 'F' 2011 12.6 'F' 2011 13.0 'M' 2011 12.6 'F' 2011 13.2 'F' 2011 9.9 'M' 2011 7.7 'F' 2011 10.5 'F' 2011 13.4 'F' 2011 10.9 'M' 2011 4.3 'F' 2011 10.3 'M' 2011 11.8 'M' 2011 11.2 'F' 2011 11.4 'F' 2011 8.6 'F' 2011 13.2 'M' 2011 12.6 'M' 2011 5.6 'M' 2011 9.9 'F' 2011 8.8 'M' 2011 7.7 'F' 2011 9.0 'M' 2011 7.3 'M' 2011 11.4 'M' 2011 13.6 'M' 2011 7.9 'M' 2011 10.7 'F' 2011 10.3 'M' 2011 8.3 'M' 2011 9.6 'M' 2011 14.2 'F' 2011 8.5 'F' 2011 13.5 'F' 2011 4.9 'F' 2011 6.4 'F' 2011 9.6 'F' 2011 11.6 'M' 2011 11.1 'M' 2012 4.35 'M' 2012 12.7 'M' 2012 18.1 'M' 2012 17.85 'F' 2012 16.6 'M' 2012 12.6 'M' 2012 17.1 'F' 2012 19.1 'M' 2012 16.1 'F' 2012 13.35 'F' 2012 18.4 'M' 2012 14.7 'M' 2012 10.6 'M' 2012 12.6 'M' 2012 16.2 'M' 2012 13.6 'M' 2012 18.9 'M' 2012 14.1 'M' 2012 14.5 'F' 2012 16.15 'M' 2012 14.75 'M' 2012 14.8 'M' 2012 12.45 'M' 2012 12.65 'M' 2012 17.35 'M' 2012 8.6 'F' 2012 18.4 'M' 2012 16.1 'M' 2012 11.6 'M' 2012 17.75 'M' 2012 15.25 'M' 2012 17.65 'F' 2012 15.6 'F' 2012 16.35 'F' 2012 17.65 'M' 2012 13.6 'F' 2012 11.7 'F' 2012 14.35 'F' 2012 14.75 'M' 2012 18.25 'F' 2012 9.9 'M' 2012 16 'M' 2012 18.25 'F' 2012 16.85 'M' 2012 14.6 'M' 2012 13.85 'M' 2012 18.95 'F' 2012 15.6 'F' 2012 14.85 'F' 2012 11.75 'F' 2012 18.45 'M' 2012 15.9 'F' 2012 17.1 'M' 2012 16.1 'F' 2012 19.9 'M' 2012 10.95 'F' 2012 18.45 'M' 2012 15.1 'F' 2012 15 'F' 2012 11.35 'M' 2012 15.95 'F' 2012 18.1 'M' 2012 14.6 'M' 2012 15.4 'M' 2012 15.4 'M' 2012 17.6 'M' 2012 13.35 'F' 2012 19.1 'M' 2012 15.35 'F' 2012 7.6 'F' 2012 13.4 'F' 2012 13.9 'M' 2012 19.1 'F' 2012 15.25 'M' 2012 12.9 'F' 2012 16.1 'F' 2012 17.35 'F' 2012 13.15 'F' 2012 12.15 'M' 2012 12.6 'M' 2012 10.35 'M' 2012 15.4 'M' 2012 9.6 'F' 2012 18.2 'F' 2012 13.6 'M' 2012 14.85 'F' 2012 14.75 'F' 2012 14.1 'F' 2012 14.9 'F' 2012 16.25 'M' 2012 19.25 'M' 2012 13.6 'F' 2012 13.6 'F' 2012 15.65 'M' 2012 12.75 'F' 2012 14.6 'M' 2012 9.85 'M' 2012 12.65 'M' 2012 11.9 'F' 2012 19.2 'M' 2012 16.6 'M' 2012 11.2 'M' 2012 15.25 'F' 2012 11.9 'F' 2012 13.2 'F' 2012 16.35 'M' 2012 12.4 'M' 2012 15.85 'F' 2012 14.35 'M' 2012 18.15 'M' 2012 11.15 'F' 2012 15.65 'F' 2012 17.75 'F' 2012 7.65 'M' 2012 12.35 'M' 2012 15.6 'F' 2012 19.3 'F' 2012 15.2 'F' 2012 17.1 'M' 2012 15.6 'M' 2012 18.4 'F' 2012 19.05 'F' 2012 18.55 'F' 2012 19.1 'M' 2012 13.1 'M' 2012 12.85 'M' 2012 9.5 'M' 2012 4.5 'F' 2012 11.85 'M' 2012 13.6 'M' 2012 11.7 'M' 2012 12.4 'F' 2012 13.35 'F' 2012 11.4 'M' 2012 14.9 'F' 2012 19.9 'M' 2012 17.75 'M' 2012 11.2 'M' 2012 14.6 'F' 2012 17.6 'M' 2012 14.05 'F' 2012 16.1 'M' 2012 13.35 'M' 2012 11.85 'F' 2012 11.95 'M' 2012 14.75 'F' 2012 15.15 'M' 2012 13.2 'F' 2012 16.85 'M' 2012 7.85 'F' 2012 7.7 'F' 2012 12.6 'M' 2012 7.85 'M' 2012 10.95 'F' 2012 12.35 'M' 2012 9.95 'M' 2012 14.9 'F' 2012 16.65 'M' 2012 13.4 'F' 2012 13.95 'F' 2012 15.7 'M' 2012 16.85 'M' 2012 10.95 'F' 2012 15.35 'M' 2012 12.2 'F' 2012 15.1 'F' 2012 17.75 'M' 2012 15.2 'F' 2012 14.6 'F' 2012 16.65 'M' 2012 8.1
Names of X columns:
Geslacht Jaar TotScores
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
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