Wessa.net

Send output to:

Data X:
888 63031 51 256 546 66751 25 160 186 7176 17 70 1405 78306 66 360 2047 137944 85 721 3626 261308 130 938 845 69266 36 287 663 83529 33 154 1181 73226 33 311 1836 178519 65 617 855 66476 35 262 1245 98606 46 385 993 50001 69 369 1685 91093 61 558 742 73884 25 220 868 72961 41 315 949 69388 34 229 332 15629 21 88 1602 71693 54 494 525 19920 17 155 629 39403 38 234 1279 99933 51 361 767 56088 28 280 1156 62006 32 331 1120 81665 51 378 622 65223 12 227 1203 88794 98 396 745 90642 53 179 1535 203699 61 509 1234 99340 25 504 757 56695 28 225 1087 108143 23 366 1105 58313 60 341 592 29101 40 171 1305 113060 29 437 0 0 0 0 706 65773 34 313 1188 67047 34 366 1111 41953 22 232 1094 109835 35 389 1087 86584 35 349 748 59588 26 316 404 40064 12 140 1076 70227 45 419 673 60437 29 226 517 47000 36 161 354 40295 13 103 1011 103397 54 356 891 78982 40 293 1067 60206 28 414 518 39887 21 156 697 49791 36 189 1095 129283 44 442 928 104816 44 321 1008 101395 33 367 951 72824 30 309 779 76018 27 235 439 33891 12 137 580 63694 39 198 614 28266 24 220 500 35093 22 149 824 35252 35 306 541 36977 20 178 475 42406 18 145 434 56353 13 144 818 58817 23 270 1173 76053 43 301 1720 70872 49 501 549 42372 20 153 157 19144 12 40 1594 114177 73 500 622 53544 26 199 656 51379 33 242 920 40756 39 265 847 46956 22 298 497 17799 20 141 864 71154 30 234 994 58305 38 336 443 27454 16 124 615 34323 34 241 525 44761 33 127 899 113862 27 327 556 35027 16 175 896 62396 22 331 515 29613 25 176 894 65559 28 281 1336 112138 112 292 505 27883 29 137 472 40181 25 155 639 53398 22 194 795 56435 23 300 1244 77283 46 370 559 71738 30 187 584 48503 31 212 440 25214 23 185 1319 119424 65 449 765 79201 27 234 222 19349 11 67 965 78760 54 316 821 54133 35 336 317 21623 16 116 425 25497 11 141 711 69535 38 236 364 30709 22 98 427 37043 23 97 463 24716 13 152 545 54865 16 132 369 27246 19 97 0 0 0 0 596 38814 13 165 479 27646 32 153 713 65373 18 226 639 43021 34 182 477 43116 40 172 38 3058 4 1 0 0 0 0 593 96347 25 196 691 48626 31 263 1056 73073 50 304 495 45266 22 183 778 43410 19 292 875 83842 34 257 490 39296 21 141 713 38490 25 192 485 39841 19 129 285 19764 12 75 934 59975 20 301 554 64589 16 204 753 63339 14 257 256 11796 9 79 80 7627 8 25 618 68998 30 217 41 6836 3 11 520 33365 16 217 42 5118 3 6 347 20898 13 115 0 0 0 0 441 42690 18 167 281 14507 11 75 81 7131 4 27 61 4194 11 14 314 21416 10 96 327 30591 13 95 554 42419 10 228
Names of X columns:
Total number of PageviewsTotal Time spent in RFC in secondsNumber of LoginsTotal number of Course Compendium Views
Response : Variable 1


Factor : Variable 2


Factor : Variable 3


Include Intercept Term ?


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')

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
R Engine error message	Error in array(list(888, 63031, 51, 256, 546, 66751, 25, 160, 186, 7176, : length of 'dimnames' [1] not equal to array extent Execution halted