Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1 1595 17 60720 319369 0 5565 47 94355 493408 0 601 6 60720 319210 1 188 11 77655 381180 1 7146 105 134028 297978 0 1135 17 62285 290476 1 450 8 59325 292136 1 34 8 60630 314353 1 133 14 65990 339445 1 119 6 59118 303677 0 2053 53 100423 397144 0 4036 63 100269 424898 1 0 0 60720 315380 1 4655 393 60720 341570 1 131 11 61808 308989 1 1766 73 60438 305959 1 312 14 58598 318690 1 448 40 59781 323361 1 115 13 60945 318903 1 60 10 61124 314049 1 0 0 60720 315380 1 364 35 47705 298466 1 0 0 60720 315380 0 1442 118 121173 438493 0 1389 9 57530 378049 1 149 5 62041 313332 1 2212 14 60720 319864 0 7489 322 136996 430866 1 0 0 60720 315380 1 402 26 84990 332743 0 7419 29 63255 286963 1 0 0 60720 315380 1 307 15 58856 331955 0 1134 162 146216 527021 0 7561 87 82425 364304 1 5131 71 86111 292154 1 9 6 60835 314987 1 2264 24 55174 272713 1 40949 481 129352 1398893 0 4980 91 113521 409642 1 272 8 112995 429112 1 757 19 45689 382712 1 1424 94 131116 374943 1 0 0 60720 315380 0 378 52 71873 297413 1 1061 5 70415 304252 1 0 0 60720 315380 1 203 9 61938 309596 0 3689 25 101792 377305 1 567 11 103345 386212 0 27664 282 74996 657954 1 1686 40 68696 368186 0 1164 17 57320 269753 1 15824 64 63838 585715 1 6575 159 151352 1071292 0 37597 248 37884 992426 1 314 5 60720 315371 0 1932 36 72954 249898 0 5692 72 64239 357312 1 0 0 60720 315380 1 0 0 60720 315380 0 6533 75 126910 364839 0 2126 138 21001 230621 1 14 1 61016 315877 0 3192 60 59831 125390 1 47 3 60720 314882 0 4891 217 210568 370837 0 531 10 71641 330068 1 2332 4 72680 324385 1 0 0 60720 315380 0 16275 109 47818 1443586 1 9252 48 79510 253537 0 29790 309 81767 4321023 1 0 0 60720 315380 0 2300 29 64873 352108 1 267 32 57640 330059 0 382 21 59365 340968 1 397 28 60646 317736 1 11920 182 31080 209458 0 54660 302 126942 1491348 1 5 2 60720 314887 1 0 0 60720 315380 1 0 0 60720 315380 0 2647 37 62920 353058 1 0 0 60720 315380 0 0 0 60720 315380 1 94 7 60793 314533 1 3 2 60698 315354 1 454 51 63261 302187 1 227 9 76644 336639 0 0 0 60720 315380 1 206 9 53110 296702 1 25830 115 245546 1073089 1 1528 38 60326 146494 1 271 2 60720 325249 1 138 6 69817 331420 1 0 0 60720 315380 1 278 26 61404 314922 1 282 30 62452 320016 1 571 80 41477 280398 1 2253 99 63593 452469 1 290 2 58790 301164 1 78 16 62700 317330 1 20 2 60805 315576 1 0 0 60720 315380 1 18073 180 27284 -7170 1 1866 29 56225 322331 0 42634 156 157214 1629616 1 249 62 54323 292754 1 422 28 57935 318056 1 2675 79 59017 355178 1 965 3 73490 204325 1 0 0 60720 315380 1 621 16 68005 317046 0 0 0 60720 315380 1 365 8 54820 309560 0 3122 34 94670 414462 0 12988 119 82340 857217 0 5336 81 112477 697458 0 3160 108 108094 530670 1 2489 40 67804 238125 1 0 0 60720 315380 0 6984 45 80570 741409 1 2001 52 95551 393343 1 15236 90 77440 372631 0 0 0 60720 315380 1 530 60 73433 317291 1 0 0 60720 315380 1 35624 223 157278 306275 0 1383 11 73221 317892 1 875 14 67000 334280 0 0 0 60720 315380 0 0 0 60720 315380 0 0 0 60720 315380 1 72 3 60398 314210 0 265 31 60720 306948 1 0 0 60720 315380 1 335 33 64175 320398 0 0 0 60720 315380 1 0 0 60720 315380 1 4525 226 93811 501749 0 3045 58 27330 202055 1 0 0 60720 315380 1 0 0 60720 315380 0 638 14 60370 333210 0 0 0 60720 315380 1 607 11 60436 322340 1 1558 63 55637 369448 0 1324 117 67440 291841 1 0 0 60720 315380 0 0 0 60720 315380 1 611 8 59190 296919 1 0 0 60720 315380 1 923 8 58620 309038 1 661 3 65920 246541 0 2397 14 74020 289513 1 366 53 61808 344425 1 0 0 60720 315380 1 135 2 62065 314210 0 1659 29 107577 480382 1 316 9 60505 315009 1 0 0 60720 315380 1 309 11 56535 312878 0 49 8 64107 322031 1 0 0 60720 315380 0 4519 49 102129 597793 1 0 0 60720 315380 1 837 69 61262 315688 1 5119 49 39039 378525 1 1280 117 69465 312378 1 2564 22 130140 403560 0 2045 20 97890 510834 0 2234 16 77200 214215 1 975 23 90534 235133 0 1136 32 48522 343613 1 453 21 81125 365959 1 0 0 60720 315380 1 61 14 60720 314551 1 368 17 61686 303230 0 0 0 60720 315380 1 4901 84 121920 469107 1 540 14 103960 354228 1 0 0 60720 315380 1 0 0 60720 315380 0 0 0 60720 315380 1 2 9 60735 315394 1 36 4 61564 312412 1 776 9 64230 333505 1 84738 588 -26007 223193 1 0 0 60720 315380 1 3 4 60761 315656 1 529 9 63870 296261 1 0 0 60720 315380 1 405 7 60845 336425 1 972 34 71642 359335 1 0 0 60720 315380 1 2099 64 106611 308636 0 3437 47 48022 158492 0 0 0 60720 315380 1 0 0 60720 315380 0 22330 84 79801 711969 0 0 0 60720 315380 0 0 0 60720 315380 1 483 27 60830 306268 1 0 0 60720 315380 1 2239 21 88590 442882 0 2949 41 82903 378509 1 0 0 60720 315380 1 365 83 87192 346611 1 2461 57 55792 314289 0 21950 519 114337 856956 0 3294 23 75832 217193 1 141 6 61630 315366 1 572 16 58580 307930 1 13326 102 165548 702380 1 2284 33 76403 194493 0 10 2 61656 316155 1 0 0 60720 315380 0 1414 198 70184 330546 1 1975 35 118881 394510 1 43 3 60887 312846 1 0 0 60720 315380 1 844 78 60925 296139 1 304 15 62969 295580 1 458 11 58625 297765 1 18562 155 102313 377934 1 0 0 60720 315380 1 7123 109 288170 638830 1 622 73 73007 304376 1 174 1 64820 307424 1 2220 22 301670 644190 1 121 24 56178 295370 0 11819 85 106113 574339 1 0 0 60720 315380 1 125 12 60798 310201 1 1182 13 70694 327007 1 1503 39 56364 343466 1 0 0 60720 315380 1 0 0 60720 315380 1 30 4 62045 318098 0 3310 33 75230 448243 1 554 11 79285 325738 0 0 0 60720 315380 1 468 21 60720 312161 1 4917 42 52811 243650 1 3256 66 35250 407159 1 0 0 60720 315380 1 125 13 59734 317698 1 0 0 60720 315380 1 22 2 60722 312502 0 0 0 60720 315380 1 514 52 78780 322378 0 0 0 60720 315380 1 0 0 60720 315380 0 0 0 60720 315380 1 0 0 60720 315380 1 0 0 60720 315380 1 0 0 60720 315380 1 10327 91 88577 640273 1 13718 61 96448 345783 1 3748 9 50350 652925 1 14416 361 49857 439798 0 1526 25 69351 278990 0 666 171 117869 339836 1 2844 67 72683 240897 1 0 0 60720 315380 1 368 17 61167 297141 1 0 0 60720 315380 1 333 13 70811 331323 1 26 1 60896 313880 1 0 0 60720 315380 1 0 0 60720 315380 1 1303 61 69863 309422 1 20 3 60938 315245 1 2384 97 61348 405972 0 203 27 50804 300962 0 71 1 60745 316176 1 53 14 59506 302409 1 562 13 58456 283587 1 622 7 60950 263276 1 645 5 60720 312075 0 1763 24 61600 308336 1 0 0 60720 315380 1 317 9 63915 298700 1 1 6 60719 315372 1 275 4 59500 318745 1 0 0 60720 315380 1 936 30 67939 408881 1 8568 37 32168 786690 0 11528 70 -14545 -83265 0 0 0 60720 315380 1 738 10 60720 321376 1 0 0 60720 315380 1 592 7 64270 276898 1 126 10 60951 315547 1 2 1 60743 315487 1 0 0 60720 315380 0 18014 239 -1710 1405225 0 0 0 60720 315380 0 0 0 60720 315380 0 37704 452 60448 983660 0 63 4 65688 318574 1 1431 34 106885 310768 1 94 3 61360 312887 1 192 4 65276 312339 1 32 1 59988 314964 1 6869 105 117520 379983 1 0 0 60720 315380 1 0 0 60720 315380 1 1 4 60722 315398 0 0 0 60720 315380 1 2328 94 82732 253588 0 209 13 64016 316647 1 28 2 60890 315688 1 176 28 68136 310670 1 1920 15 79420 165404 0 87550 458 153198 4111912 1 520 10 58650 291650 0 0 0 60720 315380 1 1013 2 69770 253468 1 15 8 60831 315688 1 587 12 59595 325699 1 5371 66 87720 446211 1 0 0 60720 315380 1 1012 47 114768 368078 1 0 0 60720 315380 1 876 59 138971 352850 0 162556 1081 213118 6282154 1 43556 247 32648 227132 0 3425 12 83620 283910 1 810 43 74015 236761 1 0 0 60720 315380 0 0 0 60720 315380 0 2365 79 191778 550608 0 1261 33 76114 307528 1 0 0 60720 315380 1 1585 107 93099 355864 1 16189 295 116384 358589 1 0 0 60720 315380 0 0 0 60720 315380 1 10579 92 110309 375195 1 474 20 61977 288985 1 0 0 60720 315380 1 3642 86 90262 458343 1 0 0 60720 315380 0 472 20 80045 269587 1 98 2 61490 315236 1 3999 30 51252 42754 1 0 0 60720 315380 1 621 48 60720 308256 1 0 0 60720 315380 1 0 0 60720 315380 0 30 4 60798 313164 1 0 0 60720 315380 1 746 23 71561 269661 0 0 0 60720 315380 1 0 0 60720 315380 0 0 0 60720 315380 1 0 0 60720 315380 0 4150 150 134759 518365 0 4658 158 156608 233773 1 814 18 39625 301881 1 1002 15 56750 298568 0 496 17 87390 325479 1 389 14 58990 325506 1 12679 110 48020 984885 1 400 152 60720 313267 1 53 6 60349 315793 1 0 0 60720 315380 0 5109 97 67038 215362 1 576 18 113761 314073 1 438 11 58320 298096 1 165 24 62841 325176 1 4069 145 79804 207393 0 23 5 62555 314806 1 1285 107 99489 341340 0 4677 75 90131 426280 0 24811 505 95350 929118 1 167 2 64245 307322 1 0 0 60720 315380 1 4628 62 69159 387475 1 226 5 65745 291787 0 1765 20 77623 247060 1 460 15 63346 329784 1 36 3 60894 315834 1 989 11 58930 304555 0 3055 43 60247 376641 0 0 0 60720 315380 1 0 0 60720 315380 0 0 0 60720 315380 1 0 0 60720 315380 1 0 0 60720 315380 1 13253 947 90829 357760 1 24 8 59818 315637 1 0 0 60720 315380 0 0 0 60720 315380 1 0 0 60720 315380 1 0 0 60720 315380 1 131 6 59661 327071 1 514 15 102725 377516 1 3366 80 108479 299243 0 9327 101 87419 387699 1 384 22 54683 309836 1 17821 373 122844 444477 1 397 10 62710 322327 1 897 5 60720 314913 1 218 21 60720 315553 0 3369 28 87161 688779 1 5702 61 101481 321896 1 8636 855 128294 301607 1 534 7 62620 304485 1 0 0 60720 315380 0 0 0 60720 315380 1 726 3 69980 231861 1 1380 26 60982 347385 1 180 11 59635 316386 0 7285 115 188873 491303 1 880 74 80791 261216 1 1 2 60727 315388 1 0 0 60720 315380 1 96 5 60379 313729 1 1889 45 37527 358649 0 45187 353 234817 1926517 1 288 4 60510 296656 1 1270 26 69206 275311 1 6526 26 55830 -42143 0 0 0 60720 315380 1 226 12 72835 343929 1 694 8 68060 367655 1 0 0 60720 315380 1 249 72 56726 313491
Names of X columns:
group costs trades dividends wealth
Endogenous Variable (Column Number)
Categorization
1
none
quantiles
hclust
equal
Number of categories (only if categorization<>none)
Cross-Validation? (only if categorization<>none)
0
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