Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
33024 1 'MVG' 32526 1 'MVG' 31455 1 'MVG' 31524 1 'MVG' 31856 1 'MVG' 32696 1 'MVG' 32584 1 'MVG' 33498 1 'MVG' 34175 1 'MVG' 34172 1 'MVG' 34379 1 'MVG' 34988 1 'MVG' 36158 1 'MVG' 37411 1 'MVG' 38015 1 'MVG' 37577 1 'MVG' 36354 1 'MVG' 36030 1 'MVG' 35636 1 'MVG' 35669 1 'MVG' 34635 1 'MVG' 35496 1 'MVG' 36376 1 'MVG' 37635 1 'MVG' 38875 1 'MVG' 38372 1 'MVG' 38897 1 'MVG' 38018 1 'MVG' 37325 1 'MVG' 36893 1 'MVG' 36117 1 'MVG' 37599 1 'MVG' 39037 1 'MVG' 40809 1 'MVG' 42508 1 'MVG' 44021 1 'MVG' 44088 1 'MVG' 44510 1 'MVG' 45786 1 'MVG' 47349 1 'MVG' 48696 1 'MVG' 50598 1 'MVG' 50066 1 'MVG' 49367 1 'MVG' 48784 1 'MVG' 47841 1 'MVG' 48300 1 'MVG' 47518 1 'MVG' 46504 1 'MVG' 45147 1 'MVG' 44404 1 'MVG' 43455 1 'MVG' 42299 1 'MVG' 42105 1 'MVG' 40152 1 'MVG' 39519 1 'MVG' 39633 1 'MVG' 39376 1 'MVG' 38850 1 'MVG' 39657 1 'MVG' 34804 1 'MVG' 34372 1 'MVG' 32678 1 'MVG' 28420 1 'MVG' 25420 1 'MVG' 27683 1 'MVG' 29904 1 'MVG' 30546 1 'MVG' 29142 1 'MVG' 27724 1 'MVG' 27069 1 'MVG' 26665 1 'MVG' 26004 1 'MVG' 25767 1 'MVG' 24915 1 'MVG' 23689 1 'MVG' 20915 1 'MVG' 19414 1 'MVG' 17824 1 'MVG' 16348 1 'MVG' 15571 1 'MVG' 13929 1 'MVG' 12480 1 'MVG' 10837 1 'MVG' 9473 1 'MVG' 8051 1 'MVG' 5278 1 'MVG' 3008 1 'MVG' 2404 1 'MVG' 2298 1 'MVG' 2260 1 'MVG' 1938 1 'MVG' 1371 1 'MVG' 1009 1 'MVG' 686 1 'MVG' 493 1 'MVG' 285 1 'MVG' 192 1 'MVG' 129 1 'MVG' 60 1 'MVG' 54 1 'MVG' 26 1 'MVG' 11 1 'MVG' 3 1 'MVG' 0 1 'MVG' 2 1 'MVG' 1 1 'MVG' 0 1 'MVG' 0 1 'MVG' 0 1 'MVG' 0 1 'MVG' 31086 0 'VVG' 30839 0 'VVG' 30051 0 'VVG' 29976 0 'VVG' 30463 0 'VVG' 31422 0 'VVG' 31588 0 'VVG' 31900 0 'VVG' 32878 0 'VVG' 33010 0 'VVG' 32954 0 'VVG' 33076 0 'VVG' 35057 0 'VVG' 35906 0 'VVG' 36100 0 'VVG' 35824 0 'VVG' 34579 0 'VVG' 34484 0 'VVG' 33920 0 'VVG' 34059 0 'VVG' 33812 0 'VVG' 34594 0 'VVG' 36083 0 'VVG' 36563 0 'VVG' 37416 0 'VVG' 37953 0 'VVG' 37517 0 'VVG' 37467 0 'VVG' 36963 0 'VVG' 36019 0 'VVG' 35232 0 'VVG' 36857 0 'VVG' 37978 0 'VVG' 40160 0 'VVG' 42165 0 'VVG' 43069 0 'VVG' 43021 0 'VVG' 43376 0 'VVG' 43978 0 'VVG' 45911 0 'VVG' 47107 0 'VVG' 49168 0 'VVG' 48390 0 'VVG' 47678 0 'VVG' 47822 0 'VVG' 46695 0 'VVG' 47185 0 'VVG' 45684 0 'VVG' 44884 0 'VVG' 44256 0 'VVG' 43637 0 'VVG' 42368 0 'VVG' 40892 0 'VVG' 40616 0 'VVG' 39026 0 'VVG' 38921 0 'VVG' 38512 0 'VVG' 38884 0 'VVG' 38406 0 'VVG' 38804 0 'VVG' 34871 0 'VVG' 34660 0 'VVG' 33104 0 'VVG' 28952 0 'VVG' 26488 0 'VVG' 29418 0 'VVG' 32315 0 'VVG' 32885 0 'VVG' 31565 0 'VVG' 30782 0 'VVG' 30442 0 'VVG' 30851 0 'VVG' 30432 0 'VVG' 31260 0 'VVG' 30737 0 'VVG' 30129 0 'VVG' 27672 0 'VVG' 26469 0 'VVG' 24895 0 'VVG' 24427 0 'VVG' 23252 0 'VVG' 21815 0 'VVG' 20837 0 'VVG' 18537 0 'VVG' 17237 0 'VVG' 15476 0 'VVG' 10709 0 'VVG' 6776 0 'VVG' 5810 0 'VVG' 5765 0 'VVG' 5775 0 'VVG' 5589 0 'VVG' 4687 0 'VVG' 3630 0 'VVG' 2552 0 'VVG' 1928 0 'VVG' 1323 0 'VVG' 1005 0 'VVG' 678 0 'VVG' 397 0 'VVG' 286 0 'VVG' 166 0 'VVG' 80 0 'VVG' 53 0 'VVG' 32 0 'VVG' 11 0 'VVG' 6 0 'VVG' 4 0 'VVG' 2 0 'VVG' 0 0 'VVG' 1 0 'VVG' 19828 1 'MWG' 19967 1 'MWG' 19814 1 'MWG' 20053 1 'MWG' 20719 1 'MWG' 21174 1 'MWG' 20648 1 'MWG' 20659 1 'MWG' 20733 1 'MWG' 21069 1 'MWG' 20566 1 'MWG' 20839 1 'MWG' 21615 1 'MWG' 22739 1 'MWG' 23222 1 'MWG' 23031 1 'MWG' 23014 1 'MWG' 22868 1 'MWG' 22182 1 'MWG' 22177 1 'MWG' 21216 1 'MWG' 21031 1 'MWG' 20968 1 'MWG' 21049 1 'MWG' 21033 1 'MWG' 21078 1 'MWG' 20702 1 'MWG' 20309 1 'MWG' 20449 1 'MWG' 20737 1 'MWG' 20849 1 'MWG' 21966 1 'MWG' 23100 1 'MWG' 23975 1 'MWG' 24350 1 'MWG' 24020 1 'MWG' 24005 1 'MWG' 23602 1 'MWG' 24120 1 'MWG' 24847 1 'MWG' 25702 1 'MWG' 26312 1 'MWG' 25891 1 'MWG' 25172 1 'MWG' 25698 1 'MWG' 25833 1 'MWG' 25658 1 'MWG' 25269 1 'MWG' 24846 1 'MWG' 24390 1 'MWG' 23954 1 'MWG' 23828 1 'MWG' 23507 1 'MWG' 23144 1 'MWG' 22302 1 'MWG' 23028 1 'MWG' 22741 1 'MWG' 23129 1 'MWG' 22911 1 'MWG' 22071 1 'MWG' 16466 1 'MWG' 16370 1 'MWG' 15049 1 'MWG' 13174 1 'MWG' 12231 1 'MWG' 13620 1 'MWG' 14317 1 'MWG' 14039 1 'MWG' 13526 1 'MWG' 12826 1 'MWG' 12360 1 'MWG' 12592 1 'MWG' 12381 1 'MWG' 12554 1 'MWG' 12338 1 'MWG' 11768 1 'MWG' 10687 1 'MWG' 9964 1 'MWG' 9338 1 'MWG' 8697 1 'MWG' 8068 1 'MWG' 7295 1 'MWG' 6372 1 'MWG' 5649 1 'MWG' 4926 1 'MWG' 4199 1 'MWG' 2568 1 'MWG' 1461 1 'MWG' 1173 1 'MWG' 1084 1 'MWG' 978 1 'MWG' 947 1 'MWG' 679 1 'MWG' 457 1 'MWG' 262 1 'MWG' 218 1 'MWG' 132 1 'MWG' 70 1 'MWG' 44 1 'MWG' 24 1 'MWG' 20 1 'MWG' 4 1 'MWG' 4 1 'MWG' 1 1 'MWG' 0 1 'MWG' 0 1 'MWG' 0 1 'MWG' 0 1 'MWG' 0 1 'MWG' 0 1 'MWG' 0 1 'MWG' 18932 0 'VWG' 18927 0 'VWG' 19124 0 'VWG' 19066 0 'VWG' 19971 0 'VWG' 20165 0 'VWG' 19705 0 'VWG' 19718 0 'VWG' 19938 0 'VWG' 20039 0 'VWG' 19721 0 'VWG' 19777 0 'VWG' 20505 0 'VWG' 21763 0 'VWG' 22404 0 'VWG' 22038 0 'VWG' 22038 0 'VWG' 21874 0 'VWG' 21269 0 'VWG' 21127 0 'VWG' 20609 0 'VWG' 20565 0 'VWG' 19791 0 'VWG' 20672 0 'VWG' 20938 0 'VWG' 20675 0 'VWG' 19992 0 'VWG' 19801 0 'VWG' 20050 0 'VWG' 20427 0 'VWG' 20815 0 'VWG' 21666 0 'VWG' 22720 0 'VWG' 23650 0 'VWG' 24244 0 'VWG' 23669 0 'VWG' 23881 0 'VWG' 23857 0 'VWG' 23999 0 'VWG' 24780 0 'VWG' 25426 0 'VWG' 26229 0 'VWG' 25973 0 'VWG' 25375 0 'VWG' 25966 0 'VWG' 25391 0 'VWG' 26046 0 'VWG' 25572 0 'VWG' 24900 0 'VWG' 24744 0 'VWG' 24526 0 'VWG' 24274 0 'VWG' 23774 0 'VWG' 23414 0 'VWG' 23002 0 'VWG' 23137 0 'VWG' 22947 0 'VWG' 23733 0 'VWG' 23234 0 'VWG' 22969 0 'VWG' 17708 0 'VWG' 17377 0 'VWG' 16273 0 'VWG' 14342 0 'VWG' 13522 0 'VWG' 15210 0 'VWG' 16493 0 'VWG' 16701 0 'VWG' 15662 0 'VWG' 15526 0 'VWG' 15413 0 'VWG' 15805 0 'VWG' 15802 0 'VWG' 16753 0 'VWG' 16906 0 'VWG' 16891 0 'VWG' 15703 0 'VWG' 15429 0 'VWG' 14762 0 'VWG' 14426 0 'VWG' 14250 0 'VWG' 13267 0 'VWG' 12397 0 'VWG' 11586 0 'VWG' 10888 0 'VWG' 9841 0 'VWG' 6443 0 'VWG' 4019 0 'VWG' 3449 0 'VWG' 3179 0 'VWG' 3341 0 'VWG' 3325 0 'VWG' 2478 0 'VWG' 1982 0 'VWG' 1405 0 'VWG' 1059 0 'VWG' 740 0 'VWG' 533 0 'VWG' 366 0 'VWG' 224 0 'VWG' 147 0 'VWG' 75 0 'VWG' 54 0 'VWG' 23 0 'VWG' 16 0 'VWG' 6 0 'VWG' 7 0 'VWG' 2 0 'VWG' 0 0 'VWG' 0 0 'VWG' 0 0 'VWG'
Names of X columns:
Aantallen Geslacht Gewest
Response : Variable 1
Factor : Variable 2
Factor : Variable 3
Include Intercept Term ?
FALSE
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