Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1418 210907 56 396 81 3 79 30 115 94 112285 24188 146283 144 145 11 869 120982 56 297 55 4 58 28 109 103 84786 18273 98364 103 101 15 1530 176508 54 559 50 12 60 38 146 93 83123 14130 86146 98 98 19 3201 385534 92 1562 63 0 121 25 96 91 119182 33251 195663 150 144 23 1583 149061 44 656 66 5 43 26 100 93 116174 27101 95757 84 84 16 1439 165446 33 511 57 0 69 25 93 60 57635 16373 85584 80 79 21 1764 237213 84 655 74 0 78 38 140 123 66198 19716 143983 130 127 24 1373 133131 55 525 52 7 44 30 99 90 57793 9028 59238 60 60 15 4041 324799 154 1436 108 0 158 47 181 168 97668 29498 151511 140 133 17 1706 230964 53 612 43 4 102 30 116 115 133824 27563 136368 151 150 19 2152 236785 119 865 75 3 77 31 116 71 101481 18293 112642 91 91 19 1036 135473 41 385 32 0 82 23 88 66 99645 22530 94728 138 132 25 1929 215147 58 639 85 0 101 36 135 117 99052 35082 121527 124 124 19 2242 344297 75 963 86 1 80 30 108 108 67654 16116 127766 119 118 28 1220 153935 33 398 56 5 50 25 89 84 65553 15849 98958 73 70 24 2515 174724 92 966 135 0 123 34 129 120 69112 26569 85646 123 119 26 2147 174415 100 801 63 0 73 31 118 114 82753 24785 98579 90 89 15 2352 225548 112 892 81 5 81 31 118 94 85323 17569 130767 116 112 21 1638 223632 73 513 52 0 105 33 125 120 72654 23825 131741 113 108 26 1222 124817 40 469 44 0 47 25 95 81 30727 7869 53907 56 52 16 1677 210767 60 643 39 3 94 35 135 133 117478 37791 146761 119 116 16 1579 170266 62 535 73 4 44 42 154 122 74007 9605 82036 129 123 20 2452 294424 77 992 59 2 107 33 127 124 101494 34461 171975 175 162 24 2662 325107 99 937 64 0 84 36 136 126 79215 24919 159676 96 92 10 186 7176 17 70 1 0 0 0 0 0 1423 603 1929 0 0 19 865 106408 30 260 32 1 33 14 46 37 31081 12558 58391 41 41 25 1793 96560 76 503 129 0 42 17 54 38 22996 7784 31580 47 47 22 2527 265769 146 927 37 2 96 32 124 120 83122 28522 136815 126 120 15 1324 149112 56 537 65 6 56 35 128 95 60578 14459 69107 80 79 21 2702 175824 107 910 107 0 57 20 80 77 39992 14526 50495 70 65 22 1383 152871 58 532 74 5 59 28 97 90 79892 22240 108016 73 70 27 1179 111665 34 345 54 4 39 28 104 80 49810 11802 46341 57 55 26 4308 362301 119 1635 715 2 76 34 125 110 100708 11912 79336 68 67 26 1831 183167 66 557 66 0 91 39 149 138 82875 18220 93176 127 127 22 1438 168809 66 452 32 0 76 28 118 100 72260 21884 127969 102 99 20 496 24188 24 218 20 0 8 4 12 7 5950 2694 15049 7 7 22 2253 329267 259 764 71 8 79 39 144 140 115762 15808 155135 148 141 21 2352 218946 41 866 112 1 76 29 108 96 80670 25239 102996 112 109 22 2144 244052 68 574 66 5 101 44 166 164 143558 29801 160604 137 133 20 4691 341570 168 1276 190 1 94 21 80 78 117105 18450 158051 135 123 21 1112 103597 43 379 66 1 27 16 60 49 23789 7132 44547 26 26 20 1973 256462 105 798 56 0 123 35 127 124 105195 35940 174141 181 166 25 2474 235800 94 921 127 8 105 23 84 62 149193 46230 184301 190 179 18 1226 196553 57 503 50 2 41 29 111 99 95260 30546 129847 107 108 22 1389 174184 53 382 52 0 72 25 98 70 55183 19746 117286 94 90 25 1496 143246 103 464 42 5 67 27 105 104 106671 15977 71180 116 114 21 2269 187559 121 717 76 8 75 36 135 116 73511 22583 109377 106 103 20 1833 187681 62 690 67 2 114 28 107 91 92945 17274 85298 143 142 20 893 73566 32 385 39 6 22 23 88 67 22618 3007 23824 26 25 18 1403 167488 45 619 77 2 69 28 104 72 83737 21113 82981 113 113 8 1425 143756 46 479 57 0 105 34 132 120 69094 17401 73815 120 118 22 1840 243199 75 752 34 3 88 28 108 105 95536 23567 132190 134 129 26 1502 182999 88 430 39 6 73 34 129 104 225920 13065 128754 54 51 18 1420 152299 53 537 63 0 62 33 122 98 61370 14587 67808 78 76 20 2970 346485 90 1000 106 0 118 38 147 111 106117 24021 131722 121 118 24 1644 193339 78 465 47 2 100 35 87 71 84651 20537 106175 145 141 17 1654 122774 45 711 162 0 24 24 90 69 15986 4527 25157 27 27 20 1054 130585 46 299 57 5 67 29 109 107 95364 30495 76669 91 91 23 937 112611 41 248 36 0 46 20 78 73 26706 7117 57283 48 48 20 3004 286468 144 1162 263 1 57 29 111 107 89691 17719 105805 68 63 22 2547 148446 91 905 63 1 135 37 141 129 126846 33473 72413 150 144 20 1626 182079 63 512 63 2 124 33 124 118 102860 21115 96971 181 168 19 1468 140344 53 472 77 6 33 25 93 73 51715 7236 71299 65 64 15 2445 220516 62 905 79 1 98 32 124 119 55801 13790 77494 97 97 20 1964 243060 63 786 110 4 58 29 112 104 111813 32902 120336 121 117 22 1381 162765 32 489 56 2 68 28 108 107 120293 25131 93913 99 100 13 1659 232138 62 617 43 0 131 31 117 90 161647 35947 181248 188 187 20 2888 265318 117 925 111 10 110 52 199 197 115929 29848 146123 138 127 17 1290 85574 34 351 71 0 37 21 78 36 24266 6943 32036 40 37 14 2845 310839 92 1144 62 9 130 24 91 85 162901 42705 186646 254 245 22 1982 225060 93 669 56 7 93 41 158 139 109825 31808 102255 87 87 24 1904 232317 54 707 74 0 118 33 126 106 129838 26675 168237 178 177 22 1391 144966 144 458 60 0 39 32 122 50 37510 8435 64219 51 49 23 1559 164709 109 572 53 0 81 31 115 63 87771 36867 115338 176 177 17 2146 220801 75 720 105 1 51 18 72 63 44418 12607 84845 66 60 23 874 99466 50 273 32 0 28 23 91 69 192565 22609 153197 56 55 25 1590 92661 61 508 133 1 40 17 45 41 35232 5892 29877 39 39 16 1590 133328 55 506 79 0 56 20 78 56 40909 17014 63506 66 64 18 1210 61361 77 451 51 0 27 12 39 25 13294 5394 22445 27 26 20 1281 100750 72 407 67 0 83 30 119 93 140867 6440 68370 58 58 18 1105 102010 53 370 66 3 28 13 50 44 44332 5818 42071 26 26 24 1272 101523 42 316 76 0 59 22 88 87 61056 18647 50517 77 76 23 1944 243511 71 603 65 0 133 42 155 110 101338 20556 103950 130 129 24 391 22938 10 154 9 0 12 1 0 0 1168 238 5841 11 11 23 1605 152474 65 577 45 0 106 32 123 83 65567 22392 84396 101 101 23 1988 99923 66 617 115 0 44 25 99 80 32334 12237 35753 36 36 13 1386 132487 41 411 97 0 71 36 136 98 40735 8388 55515 120 89 20 2395 317394 86 975 53 1 116 31 117 82 91413 22120 209056 195 193 18 387 21054 16 146 2 0 4 0 0 0 855 338 6622 4 4 21 1742 209641 42 705 52 5 62 24 88 60 97068 11727 115814 89 84 17 620 22648 19 184 44 0 12 13 39 28 44339 3704 11609 24 23 20 449 31414 19 200 22 0 18 8 25 9 14116 3988 13155 39 39 19 800 46698 45 274 35 0 14 13 52 33 10288 3030 18274 14 14 18 1684 131698 65 502 74 0 60 19 75 59 65622 13520 72875 78 78 19 2699 244749 95 964 144 2 98 33 124 115 76643 20923 142775 106 101 22 1204 128423 64 369 89 8 32 38 145 120 92696 3769 20112 37 36 22 1138 97839 38 417 42 2 25 24 87 66 94785 12252 61023 77 75 15 2158 272458 65 822 99 0 100 43 162 152 93815 28864 132432 132 131 17 1111 172494 52 389 52 0 46 43 165 139 86687 21721 112494 144 131 19 1421 108043 62 466 29 1 45 14 54 38 34553 4821 45109 40 39 20 2833 328107 65 1255 125 3 129 41 159 144 105547 33644 170875 153 144 22 2922 351067 95 1024 95 3 136 45 170 160 213688 42935 214921 220 211 21 1002 158015 29 400 40 0 59 31 119 114 71220 18864 100226 79 78 19 2186 229242 247 719 128 4 63 31 120 119 91721 17939 78876 95 90 21 1035 84207 29 356 73 11 14 30 112 101 111194 325 6940 12 12 18 1417 120445 118 457 72 0 36 16 59 56 51009 13539 49025 63 57 16 3261 324598 110 1402 128 0 113 37 136 133 135777 34538 122037 134 133 20 1587 131069 67 600 61 4 47 30 107 83 51513 12198 53782 69 69 21 1424 204271 42 480 73 0 92 35 130 116 74163 26924 127748 119 119 15 946 116048 64 230 45 0 50 20 75 50 33416 10855 77395 63 61 20 1926 250047 81 651 58 0 41 18 71 61 83305 11932 89324 55 49 23 3352 299775 95 1367 97 9 91 31 120 97 98952 14300 103300 103 101 15 1641 195838 67 564 50 1 111 31 116 98 102372 25515 112283 197 196 18 2035 173260 63 716 37 3 41 21 79 78 37238 2805 10901 16 15 22 2312 254488 83 747 50 10 120 39 150 117 103772 29402 120691 140 136 16 961 92499 32 319 57 0 25 18 71 55 21399 5250 25899 21 21 17 1900 224330 83 612 52 1 131 39 144 132 130115 28608 139296 167 163 24 1254 135781 31 433 98 2 45 14 47 44 24874 8092 52678 32 29 13 1335 74408 67 434 61 4 29 7 28 21 34988 4473 23853 36 35 23 1597 81240 66 503 89 0 58 17 68 50 45549 1572 17306 13 13 5 1645 181633 70 564 48 2 47 30 110 73 64466 14817 89455 96 96 19 2429 271856 103 824 91 1 109 37 147 86 54990 16714 147866 151 151 24 872 95227 34 239 70 0 37 32 111 48 34777 1669 14336 23 23 19 1018 98146 40 459 37 0 15 17 68 48 27114 7768 30059 21 14 20 1314 59194 31 288 247 6 7 24 80 68 37636 6387 22097 20 20 22 1335 139942 42 498 46 0 54 22 88 87 65461 18715 96841 82 72 15 1403 118612 46 454 72 2 54 12 48 43 30080 7936 41907 90 87 19 910 72880 33 376 41 0 14 19 76 67 24094 8643 27080 25 21 25 616 65475 18 225 24 2 16 13 51 46 69008 7294 35885 60 56 21 771 71965 35 252 33 1 32 15 59 56 46090 7185 28313 85 82 19 1376 135131 66 481 87 0 38 15 60 60 34029 8509 36134 41 38 17 1232 108446 60 389 90 1 22 17 68 65 46300 13275 55764 26 25 15 1544 181528 54 609 69 0 32 16 61 60 40662 10737 66956 49 47 21 1230 134019 53 422 51 0 32 18 67 54 28987 8033 47487 46 45 24 1255 121848 39 339 45 0 37 17 64 52 30594 5401 35619 41 41 22 721 81872 45 245 25 0 32 16 64 61 27913 10856 45608 23 23 19 1109 58981 36 384 38 7 0 23 91 61 42744 2154 7721 14 14 20 740 53515 28 212 52 2 5 22 88 81 12934 6117 20634 16 16 21 728 56375 30 229 74 7 10 13 49 40 41385 4820 31931 21 21 19 689 65490 22 224 38 3 27 16 62 40 18653 5615 37754 32 27 22 995 76302 31 333 26 0 29 20 76 68 30976 8702 40557 35 33 14 1613 104011 55 384 67 6 25 22 88 79 63339 15340 94238 42 42 25 2048 98104 54 636 132 2 55 17 66 47 25568 8030 44197 68 68 11 301 30989 14 93 35 0 5 17 68 41 4154 1278 4103 6 6 16 1803 135458 81 581 118 3 43 12 48 29 19474 4236 44144 68 67 19 861 63123 43 304 43 1 34 17 68 60 39067 7196 27640 84 77 17 1451 74914 30 407 64 0 35 23 90 79 65892 6371 28990 30 30 20 628 31774 23 170 48 1 0 17 66 47 4143 1574 4694 0 0 22 1161 81437 38 312 64 0 37 14 54 40 28579 9620 42648 36 36 20 979 65745 53 340 75 0 26 21 77 42 38084 8645 25836 50 48 22 675 56653 45 168 39 0 38 18 68 49 27717 8987 22779 30 29 15 1241 158399 39 443 42 0 23 18 72 57 32928 5544 40820 30 28 23 1049 73624 24 367 93 0 30 17 64 40 19499 6909 32378 33 33 20 1081 91899 35 335 60 0 18 15 59 33 36874 6745 39613 37 33 17 1688 139526 151 364 71 0 28 21 84 77 48259 16724 60865 83 80 20 617 51567 30 206 27 2 21 14 56 45 28207 7025 20107 30 30 25 1656 102538 57 490 79 1 50 15 58 45 45833 9078 48231 51 51 22 705 86678 40 238 44 0 12 15 59 50 29156 4605 39725 19 18 16 1597 150580 77 530 124 0 27 22 83 71 45588 9653 62991 41 39 25 982 99611 35 291 81 0 41 21 81 67 45097 8914 49363 54 54 18 1212 99373 63 397 92 1 12 18 72 62 28394 6700 24552 25 24 19 1143 86230 44 467 42 0 21 17 61 54 18632 5788 31493 25 24 25 435 30837 19 178 10 0 8 4 15 4 2325 593 3439 8 8 23 532 31706 13 175 24 0 26 10 32 25 25139 4506 19555 26 26 24 882 89806 42 299 64 0 27 16 62 40 27975 6382 21228 20 19 21 830 64175 42 260 48 0 37 18 72 59 21792 6928 28893 46 47 21 652 59382 49 227 49 0 29 12 41 24 26263 1514 21425 47 47 22 707 119308 30 239 48 0 32 16 61 58 23686 9238 50276 37 37 21 954 76702 49 333 62 0 35 21 67 42 49303 8204 37643 51 51 18 285 19764 12 75 19 1 10 2 8 4 5752 2416 9927 10 10 13 733 84105 20 261 45 0 17 17 66 63 20055 5432 27184 34 34 22 642 64187 27 238 36 0 10 16 61 54 20154 5576 18475 12 11 23 894 72535 14 329 44 0 17 16 64 39 19540 6095 35873 27 21 15
Names of X columns:
pageviews time_in_rfc logins compendium_views_info compendium_views_pr shared_compendiums blogged_computations compendiums_reviewed feedback_messages_p1 feedback_messages_p120 totsize totrevisions totseconds tothyperlinks totblogs I1
Endogenous Variable (Column Number)
Categorization
none
quantiles
hclust
equal
Number of categories (only if categorization<>none)
Cross-Validation? (only if categorization<>none)
no
yes
Chart options
R Code
library(party) library(Hmisc) par1 <- as.numeric(par1) par3 <- as.numeric(par3) x <- data.frame(t(y)) is.data.frame(x) x <- x[!is.na(x[,par1]),] k <- length(x[1,]) n <- length(x[,1]) colnames(x)[par1] x[,par1] if (par2 == 'kmeans') { cl <- kmeans(x[,par1], par3) print(cl) clm <- matrix(cbind(cl$centers,1:par3),ncol=2) clm <- clm[sort.list(clm[,1]),] for (i in 1:par3) { cl$cluster[cl$cluster==clm[i,2]] <- paste('C',i,sep='') } cl$cluster <- as.factor(cl$cluster) print(cl$cluster) x[,par1] <- cl$cluster } if (par2 == 'quantiles') { x[,par1] <- cut2(x[,par1],g=par3) } if (par2 == 'hclust') { hc <- hclust(dist(x[,par1])^2, 'cen') print(hc) memb <- cutree(hc, k = par3) dum <- c(mean(x[memb==1,par1])) for (i in 2:par3) { dum <- c(dum, mean(x[memb==i,par1])) } hcm <- matrix(cbind(dum,1:par3),ncol=2) hcm <- hcm[sort.list(hcm[,1]),] for (i in 1:par3) { memb[memb==hcm[i,2]] <- paste('C',i,sep='') } memb <- as.factor(memb) print(memb) x[,par1] <- memb } if (par2=='equal') { ed <- cut(as.numeric(x[,par1]),par3,labels=paste('C',1:par3,sep='')) x[,par1] <- as.factor(ed) } table(x[,par1]) colnames(x) colnames(x)[par1] x[,par1] if (par2 == 'none') { m <- ctree(as.formula(paste(colnames(x)[par1],' ~ .',sep='')),data = x) } load(file='createtable') if (par2 != 'none') { m <- ctree(as.formula(paste('as.factor(',colnames(x)[par1],') ~ .',sep='')),data = x) if (par4=='yes') { a<-table.start() a<-table.row.start(a) a<-table.element(a,'10-Fold Cross Validation',3+2*par3,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'',1,TRUE) a<-table.element(a,'Prediction (training)',par3+1,TRUE) a<-table.element(a,'Prediction (testing)',par3+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Actual',1,TRUE) for (jjj in 1:par3) a<-table.element(a,paste('C',jjj,sep=''),1,TRUE) a<-table.element(a,'CV',1,TRUE) for (jjj in 1:par3) a<-table.element(a,paste('C',jjj,sep=''),1,TRUE) a<-table.element(a,'CV',1,TRUE) a<-table.row.end(a) for (i in 1:10) { ind <- sample(2, nrow(x), replace=T, prob=c(0.9,0.1)) m.ct <- ctree(as.formula(paste('as.factor(',colnames(x)[par1],') ~ .',sep='')),data =x[ind==1,]) if (i==1) { m.ct.i.pred <- predict(m.ct, newdata=x[ind==1,]) m.ct.i.actu <- x[ind==1,par1] m.ct.x.pred <- predict(m.ct, newdata=x[ind==2,]) m.ct.x.actu <- x[ind==2,par1] } else { m.ct.i.pred <- c(m.ct.i.pred,predict(m.ct, newdata=x[ind==1,])) m.ct.i.actu <- c(m.ct.i.actu,x[ind==1,par1]) m.ct.x.pred <- c(m.ct.x.pred,predict(m.ct, newdata=x[ind==2,])) m.ct.x.actu <- c(m.ct.x.actu,x[ind==2,par1]) } } print(m.ct.i.tab <- table(m.ct.i.actu,m.ct.i.pred)) numer <- 0 for (i in 1:par3) { print(m.ct.i.tab[i,i] / sum(m.ct.i.tab[i,])) numer <- numer + m.ct.i.tab[i,i] } print(m.ct.i.cp <- numer / sum(m.ct.i.tab)) print(m.ct.x.tab <- table(m.ct.x.actu,m.ct.x.pred)) numer <- 0 for (i in 1:par3) { print(m.ct.x.tab[i,i] / sum(m.ct.x.tab[i,])) numer <- numer + m.ct.x.tab[i,i] } print(m.ct.x.cp <- numer / sum(m.ct.x.tab)) for (i in 1:par3) { a<-table.row.start(a) a<-table.element(a,paste('C',i,sep=''),1,TRUE) for (jjj in 1:par3) a<-table.element(a,m.ct.i.tab[i,jjj]) a<-table.element(a,round(m.ct.i.tab[i,i]/sum(m.ct.i.tab[i,]),4)) for (jjj in 1:par3) a<-table.element(a,m.ct.x.tab[i,jjj]) a<-table.element(a,round(m.ct.x.tab[i,i]/sum(m.ct.x.tab[i,]),4)) a<-table.row.end(a) } a<-table.row.start(a) a<-table.element(a,'Overall',1,TRUE) for (jjj in 1:par3) a<-table.element(a,'-') a<-table.element(a,round(m.ct.i.cp,4)) for (jjj in 1:par3) a<-table.element(a,'-') a<-table.element(a,round(m.ct.x.cp,4)) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable3.tab') } } m bitmap(file='test1.png') plot(m) dev.off() bitmap(file='test1a.png') plot(x[,par1] ~ as.factor(where(m)),main='Response by Terminal Node',xlab='Terminal Node',ylab='Response') dev.off() if (par2 == 'none') { forec <- predict(m) result <- as.data.frame(cbind(x[,par1],forec,x[,par1]-forec)) colnames(result) <- c('Actuals','Forecasts','Residuals') print(result) } if (par2 != 'none') { print(cbind(as.factor(x[,par1]),predict(m))) myt <- table(as.factor(x[,par1]),predict(m)) print(myt) } bitmap(file='test2.png') if(par2=='none') { op <- par(mfrow=c(2,2)) plot(density(result$Actuals),main='Kernel Density Plot of Actuals') plot(density(result$Residuals),main='Kernel Density Plot of Residuals') plot(result$Forecasts,result$Actuals,main='Actuals versus Predictions',xlab='Predictions',ylab='Actuals') plot(density(result$Forecasts),main='Kernel Density Plot of Predictions') par(op) } if(par2!='none') { plot(myt,main='Confusion Matrix',xlab='Actual',ylab='Predicted') } dev.off() if (par2 == 'none') { detcoef <- cor(result$Forecasts,result$Actuals) a<-table.start() a<-table.row.start(a) a<-table.element(a,'Goodness of Fit',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Correlation',1,TRUE) a<-table.element(a,round(detcoef,4)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'R-squared',1,TRUE) a<-table.element(a,round(detcoef*detcoef,4)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'RMSE',1,TRUE) a<-table.element(a,round(sqrt(mean((result$Residuals)^2)),4)) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable1.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Actuals, Predictions, and Residuals',4,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'#',header=TRUE) a<-table.element(a,'Actuals',header=TRUE) a<-table.element(a,'Forecasts',header=TRUE) a<-table.element(a,'Residuals',header=TRUE) a<-table.row.end(a) for (i in 1:length(result$Actuals)) { a<-table.row.start(a) a<-table.element(a,i,header=TRUE) a<-table.element(a,result$Actuals[i]) a<-table.element(a,result$Forecasts[i]) a<-table.element(a,result$Residuals[i]) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable.tab') } if (par2 != 'none') { a<-table.start() a<-table.row.start(a) a<-table.element(a,'Confusion Matrix (predicted in columns / actuals in rows)',par3+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'',1,TRUE) for (i in 1:par3) { a<-table.element(a,paste('C',i,sep=''),1,TRUE) } a<-table.row.end(a) for (i in 1:par3) { a<-table.row.start(a) a<-table.element(a,paste('C',i,sep=''),1,TRUE) for (j in 1:par3) { a<-table.element(a,myt[i,j]) } a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable2.tab') }
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
2 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation