Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
40 1 27 5 26 49 35 45 1 36 4 25 45 34 38 1 25 4 17 54 13 28 1 27 3 37 36 35 NA 2 25 3 35 36 28 38 2 44 3 15 53 32 39 1 50 4 27 46 35 37 1 41 4 36 42 36 30 1 48 5 25 41 27 30 2 43 4 30 45 29 30 2 47 2 27 47 27 26 2 41 3 33 42 28 29 1 44 2 29 45 29 31 2 47 5 30 40 28 27 2 40 3 25 45 30 25 2 46 3 23 40 25 39 1 28 3 26 42 15 35 1 56 3 24 45 33 27 2 49 4 35 47 31 40 2 25 4 39 31 37 34 2 41 4 23 46 37 32 2 26 3 32 34 34 34 1 50 5 29 43 32 38 1 47 4 26 45 21 21 1 52 2 21 42 25 33 2 37 5 35 51 32 27 2 41 3 23 44 28 35 1 45 4 21 47 22 33 2 26 4 28 47 25 36 1 NA 3 30 41 26 NA 1 52 4 21 44 34 37 1 46 2 29 51 34 37 1 58 3 28 46 36 37 1 54 5 19 47 36 32 1 29 3 26 46 26 25 2 50 3 33 38 26 31 1 43 2 34 50 34 33 2 30 3 33 48 33 18 2 47 2 40 36 31 42 1 45 3 24 51 33 26 2 48 1 35 35 22 26 2 48 3 35 49 29 32 2 26 4 32 38 24 31 1 46 5 20 47 37 29 2 NA 3 35 36 32 35 2 50 3 35 47 23 44 1 25 4 21 46 29 35 1 47 2 33 43 35 30 2 47 2 40 53 20 32 1 41 3 22 55 28 24 2 45 2 35 39 26 34 2 41 4 20 55 36 27 2 45 5 28 41 26 31 2 40 3 46 33 33 38 1 29 4 18 52 25 41 2 34 5 22 42 29 40 1 45 5 20 56 32 25 2 52 3 25 46 35 19 2 41 4 31 33 24 33 2 48 3 21 51 31 27 2 45 3 23 46 29 45 1 54 2 26 46 27 27 2 25 3 34 50 29 30 2 26 4 31 46 29 42 1 28 4 23 51 27 21 2 50 4 31 48 34 32 2 48 4 26 44 32 31 2 51 3 36 38 31 36 2 53 3 28 42 31 34 1 37 3 34 39 31 11 1 56 2 25 45 16 35 1 43 3 33 31 25 39 1 34 3 46 29 27 32 1 42 3 24 48 32 28 2 32 3 32 38 28 45 2 31 5 33 55 25 18 1 46 3 42 32 25 35 2 30 5 17 51 36 35 2 47 4 36 53 36 36 2 33 4 40 47 36 34 1 25 4 30 45 27 34 1 25 5 19 33 29 38 2 21 4 33 49 32 28 2 36 5 35 46 29 23 2 50 3 23 42 31 37 2 48 3 15 56 34 29 2 48 2 38 35 27 28 1 25 3 37 40 28 30 1 48 4 23 44 32 24 2 49 5 41 46 33 36 1 27 5 34 46 29 40 1 28 3 38 39 32 37 2 43 2 45 35 35 27 2 48 3 27 48 33 25 2 48 4 46 42 27 22 1 25 1 26 39 16 21 2 49 4 44 39 32 28 1 26 3 36 41 26 34 1 51 3 20 52 32 32 2 25 4 44 45 38 23 1 29 3 27 42 24 29 1 29 4 27 44 26 35 1 43 2 41 33 19 31 2 46 3 30 42 37 36 1 44 3 33 46 25 32 1 25 3 37 45 24 35 1 51 2 30 40 23 45 1 42 5 20 48 28 29 2 53 5 44 32 38 41 1 25 4 20 53 28 36 2 49 2 33 39 28 37 1 51 3 31 45 26 25 2 20 3 23 36 21 36 2 44 3 33 38 35 34 2 38 4 33 49 31 33 1 46 5 32 46 34 32 2 42 4 25 43 30 40 1 29 NA 22 37 30 27 2 46 4 16 48 24 24 2 49 2 36 45 27 26 2 51 3 35 32 26 NA 1 38 3 25 46 30 13 1 41 1 27 20 15 22 2 47 3 32 42 28 29 2 44 3 36 45 34 30 2 47 3 51 29 29 24 2 46 3 30 51 26 NA 1 44 4 20 55 31 26 2 28 3 29 50 28 37 2 47 4 26 44 33 36 2 28 4 20 41 32 38 1 41 5 40 40 33 34 2 45 4 29 47 31 35 2 46 4 32 42 37 32 1 46 4 33 40 27 44 2 22 3 32 51 19 40 2 33 3 34 43 27 24 1 41 4 24 45 31 36 2 47 5 25 41 38 20 1 25 3 41 41 22 28 2 42 3 39 37 35 18 2 47 3 21 46 35 23 2 50 3 38 38 30 28 1 55 5 28 39 41 30 1 21 3 37 45 25 30 1 NA 3 26 46 28 43 1 52 3 30 39 45 20 2 49 4 25 21 21 37 2 46 4 38 31 33 24 1 NA 4 31 35 25 33 2 45 3 31 49 29 43 2 52 3 27 40 31 27 1 NA 3 21 45 29 22 2 40 4 26 46 31 28 2 49 4 37 45 31 18 1 38 5 28 34 25 38 1 32 5 29 41 27 23 2 46 4 33 43 26 38 2 32 3 41 45 26 21 2 41 3 19 48 23 25 2 43 3 37 43 27 NA 1 44 4 36 45 24 30 1 47 5 27 45 35 25 2 28 3 33 34 24 17 1 52 1 29 40 32 26 1 27 2 42 40 24 39 2 45 5 27 55 24 27 1 27 4 47 44 38 33 1 25 4 17 44 36 47 1 28 4 34 48 24 37 1 25 3 32 51 18 34 1 52 4 25 49 34 24 1 44 3 27 33 23 25 2 43 3 37 43 35 20 2 47 4 34 44 22 34 2 52 4 27 44 34 22 2 40 2 37 41 28 39 1 42 3 32 45 34 33 1 45 5 26 44 32 35 1 45 2 29 44 24 26 1 50 5 28 40 34 32 1 49 3 19 48 33 22 1 52 2 46 49 33 39 2 48 3 31 46 29 35 2 51 3 42 49 38 21 2 49 4 33 55 24 27 2 31 4 39 51 25 31 2 43 3 27 46 37 20 2 31 3 35 37 33 28 2 28 4 23 43 30 26 2 43 4 32 41 22 36 2 31 3 22 45 28 16 2 51 3 17 39 24 34 2 58 4 35 38 33 30 2 25 5 34 41 37
Names of X columns:
Intrinsieke_waarden geslacht leeftijd opleiding Neuroticisme Extraversie Openheid
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