Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
0.04374 0.426 0.02182 0.0313 0.02971 0.06545 0.06134 0.626 0.03134 0.04518 0.04368 0.09403 0.05233 0.482 0.02757 0.03858 0.0359 0.0827 0.05492 0.517 0.02924 0.04005 0.03772 0.08771 0.06425 0.584 0.0349 0.04825 0.04465 0.1047 0.04701 0.456 0.02328 0.03526 0.03243 0.06985 0.01608 0.14 0.00779 0.00937 0.01351 0.02337 0.01567 0.134 0.00829 0.00946 0.01256 0.02487 0.02093 0.191 0.01073 0.01277 0.01717 0.03218 0.02838 0.255 0.01441 0.01725 0.02444 0.04324 0.02143 0.197 0.01079 0.01342 0.01892 0.03237 0.02752 0.249 0.01424 0.01641 0.02214 0.04272 0.01259 0.112 0.00656 0.00717 0.0114 0.01968 0.01642 0.154 0.00728 0.00932 0.01797 0.02184 0.01828 0.158 0.01064 0.00972 0.01246 0.03191 0.01503 0.126 0.00772 0.00888 0.01359 0.02316 0.02047 0.192 0.00969 0.012 0.02074 0.02908 0.03327 0.348 0.01441 0.01893 0.0343 0.04322 0.05517 0.542 0.02471 0.03572 0.05767 0.07413 0.03995 0.348 0.01721 0.02374 0.0431 0.05164 0.0381 0.328 0.01667 0.02383 0.04055 0.05 0.04137 0.37 0.02021 0.02591 0.04525 0.06062 0.04351 0.377 0.02228 0.0254 0.04246 0.06685 0.04192 0.364 0.02187 0.0247 0.03772 0.06562 0.01659 0.164 0.00738 0.00948 0.01497 0.02214 0.03767 0.381 0.01732 0.02245 0.0378 0.05197 0.01966 0.186 0.00889 0.01169 0.01872 0.02666 0.01919 0.198 0.00883 0.01144 0.01826 0.0265 0.01718 0.161 0.00769 0.01012 0.01661 0.02307 0.01791 0.168 0.00793 0.01057 0.01799 0.0238 0.01098 0.097 0.00563 0.0068 0.00802 0.01689 0.01015 0.089 0.00504 0.00641 0.00762 0.01513 0.01263 0.111 0.0064 0.00825 0.00951 0.01919 0.00954 0.085 0.00469 0.00606 0.00719 0.01407 0.00958 0.085 0.00468 0.0061 0.00726 0.01403 0.01194 0.107 0.00586 0.0076 0.00957 0.01758 0.02126 0.189 0.01154 0.01347 0.01612 0.03463 0.01851 0.168 0.00938 0.0116 0.01491 0.02814 0.01444 0.131 0.00726 0.00885 0.0119 0.02177 0.01663 0.151 0.00829 0.01003 0.01366 0.02488 0.01495 0.135 0.00774 0.00941 0.01233 0.02321 0.01463 0.132 0.00742 0.00901 0.01234 0.02226 0.01752 0.164 0.01035 0.01024 0.01133 0.03104 0.0176 0.154 0.01006 0.01038 0.01251 0.03017 0.01419 0.126 0.00777 0.00898 0.01033 0.0233 0.01494 0.134 0.00847 0.00879 0.01014 0.02542 0.01608 0.141 0.00906 0.00977 0.01149 0.02719 0.01152 0.103 0.00614 0.0073 0.0086 0.01841 0.01613 0.143 0.00855 0.00776 0.01433 0.02566 0.01681 0.154 0.0093 0.00802 0.014 0.02789 0.02184 0.197 0.01241 0.01024 0.01685 0.03724 0.02033 0.185 0.01143 0.00959 0.01614 0.03429 0.02297 0.21 0.01323 0.01072 0.01677 0.03969 0.02498 0.228 0.01396 0.01219 0.01947 0.04188 0.02719 0.255 0.01483 0.01609 0.02067 0.0445 0.03209 0.307 0.01789 0.01992 0.02454 0.05368 0.03715 0.334 0.02032 0.02302 0.02802 0.06097 0.02293 0.221 0.01189 0.01459 0.01948 0.03568 0.02645 0.265 0.01394 0.01625 0.02137 0.04183 0.03225 0.35 0.01805 0.01974 0.02519 0.05414 0.01861 0.17 0.00975 0.01258 0.01382 0.02925 0.01906 0.165 0.01013 0.01296 0.0134 0.03039 0.01643 0.145 0.00867 0.01108 0.012 0.02602 0.01644 0.145 0.00882 0.01075 0.01179 0.02647 0.01457 0.129 0.00769 0.00957 0.01016 0.02308 0.01745 0.154 0.00942 0.0116 0.01234 0.02827 0.03198 0.313 0.0183 0.0181 0.02428 0.0549 0.03111 0.308 0.01638 0.01759 0.02603 0.04914 0.05384 0.478 0.03152 0.02422 0.03392 0.09455 0.05428 0.497 0.03357 0.02494 0.03635 0.1007 0.03485 0.365 0.01868 0.01906 0.02949 0.05605 0.04978 0.483 0.02749 0.02466 0.03736 0.08247 0.01706 0.152 0.00974 0.00925 0.01345 0.02921 0.02448 0.226 0.01373 0.01375 0.01956 0.0412 0.02442 0.216 0.01432 0.01325 0.01831 0.04295 0.02215 0.206 0.01284 0.01219 0.01715 0.03851 0.03999 0.35 0.02413 0.02231 0.02704 0.07238 0.02199 0.197 0.01284 0.01199 0.01636 0.03852 0.03202 0.263 0.01803 0.01886 0.02455 0.05408 0.03121 0.361 0.01773 0.01783 0.02139 0.0532 0.04024 0.364 0.02266 0.02451 0.02876 0.06799 0.03156 0.296 0.01792 0.01841 0.0219 0.05377 0.02427 0.216 0.01371 0.01421 0.01751 0.04114 0.02223 0.202 0.01277 0.01343 0.01552 0.03831 0.04795 0.435 0.02679 0.03022 0.0351 0.08037 0.03852 0.331 0.02107 0.02493 0.02877 0.06321 0.03759 0.327 0.02073 0.02415 0.02784 0.06219 0.06511 0.58 0.03671 0.04159 0.04683 0.11012 0.06727 0.65 0.03788 0.04254 0.04802 0.11363 0.04313 0.442 0.02297 0.02768 0.03455 0.06892 0.0664 0.634 0.0365 0.04282 0.05114 0.10949 0.07959 0.772 0.04421 0.04962 0.0569 0.13262 0.0419 0.383 0.02383 0.02521 0.03051 0.0715 0.05925 0.637 0.03341 0.03794 0.04398 0.10024 0.03716 0.307 0.02062 0.02321 0.02764 0.06185 0.03272 0.283 0.01813 0.01909 0.02571 0.05439 0.03381 0.307 0.01806 0.02024 0.02809 0.05417 0.03886 0.342 0.02135 0.02174 0.03088 0.06406 0.04689 0.422 0.02542 0.0263 0.03908 0.07625 0.06734 0.659 0.03611 0.03963 0.05783 0.10833 0.09178 0.891 0.05358 0.04791 0.06196 0.16074 0.0617 0.584 0.03223 0.03672 0.05174 0.09669 0.09419 0.93 0.05551 0.05005 0.06023 0.16654 0.01131 0.107 0.00522 0.00659 0.01009 0.01567 0.0103 0.094 0.00469 0.00582 0.00871 0.01406 0.01346 0.126 0.0066 0.00818 0.01059 0.01979 0.01064 0.097 0.00522 0.00632 0.00928 0.01567 0.0145 0.137 0.00633 0.00788 0.01267 0.01898 0.01024 0.093 0.00455 0.00576 0.00993 0.01364 0.03044 0.275 0.01771 0.01815 0.02084 0.05312 0.02286 0.207 0.01192 0.01439 0.01852 0.03576 0.01761 0.155 0.00952 0.01058 0.01307 0.02855 0.02378 0.21 0.01277 0.01483 0.01767 0.03831 0.0168 0.149 0.00861 0.01017 0.01301 0.02583 0.02105 0.209 0.01107 0.01284 0.01604 0.0332 0.01843 0.235 0.00796 0.00832 0.01271 0.02389 0.01458 0.148 0.00606 0.00747 0.01312 0.01818 0.01725 0.175 0.00757 0.00971 0.01652 0.0227 0.01279 0.129 0.00617 0.00744 0.01151 0.01851 0.01299 0.124 0.00679 0.00631 0.01075 0.02038 0.02008 0.221 0.00849 0.01117 0.01734 0.02548 0.01169 0.117 0.00534 0.0063 0.01104 0.01603 0.04479 0.441 0.02587 0.02567 0.0322 0.07761 0.02503 0.231 0.01372 0.0158 0.01931 0.04115 0.02343 0.224 0.01289 0.0142 0.0172 0.03867 0.02362 0.233 0.01235 0.01495 0.01944 0.03706 0.02791 0.246 0.01484 0.01805 0.02259 0.04451 0.02857 0.257 0.01547 0.01859 0.02301 0.04641 0.01033 0.098 0.00538 0.0057 0.00811 0.01614 0.01022 0.09 0.00476 0.00588 0.00903 0.01428 0.01412 0.125 0.00703 0.0082 0.01194 0.0211 0.01516 0.138 0.00721 0.00815 0.0131 0.02164 0.01201 0.106 0.00633 0.00701 0.00915 0.01898 0.01043 0.099 0.0049 0.00621 0.00903 0.01471 0.04932 0.441 0.02683 0.03112 0.03651 0.0805 0.04128 0.379 0.02229 0.02592 0.03316 0.06688 0.04879 0.431 0.02385 0.02973 0.0437 0.07154 0.05279 0.476 0.02896 0.03347 0.04134 0.08689 0.05643 0.517 0.0307 0.0353 0.04451 0.09211 0.03026 0.267 0.01514 0.01812 0.0277 0.04543 0.03273 0.281 0.01713 0.01964 0.02824 0.05139 0.06725 0.571 0.04016 0.04003 0.04464 0.12047 0.03527 0.297 0.02055 0.02076 0.0253 0.06165 0.01997 0.18 0.01117 0.01177 0.01506 0.0335 0.02662 0.228 0.01475 0.01558 0.02006 0.04426 0.02536 0.225 0.01379 0.01478 0.01909 0.04137 0.08143 0.821 0.03804 0.05426 0.08808 0.11411 0.0605 0.618 0.02865 0.04101 0.06359 0.08595 0.07118 0.722 0.03474 0.0458 0.06824 0.10422 0.0717 0.833 0.03515 0.04265 0.0646 0.10546 0.0583 0.784 0.02699 0.03714 0.06259 0.08096 0.11908 1.302 0.05647 0.0794 0.13778 0.16942 0.08684 1.018 0.04284 0.05556 0.08318 0.12851 0.02534 0.241 0.0134 0.01399 0.02056 0.04019 0.02682 0.236 0.01484 0.01405 0.02018 0.04451 0.03087 0.276 0.01659 0.01804 0.02402 0.04977 0.02293 0.223 0.01205 0.01289 0.01771 0.03615 0.04912 0.438 0.0261 0.02161 0.02916 0.0783 0.02852 0.266 0.015 0.01581 0.02157 0.04499 0.03235 0.339 0.0136 0.0165 0.03105 0.04079 0.04009 0.406 0.01579 0.01994 0.04114 0.04736 0.03273 0.325 0.01644 0.01722 0.02931 0.04933 0.03658 0.369 0.01864 0.0194 0.03091 0.05592 0.01756 0.155 0.00967 0.01033 0.01363 0.02902 0.02814 0.272 0.01579 0.01553 0.02073 0.04736 0.02448 0.217 0.0141 0.01426 0.01621 0.04231 0.01242 0.116 0.00696 0.00747 0.00882 0.02089 0.0203 0.197 0.01186 0.0123 0.01367 0.03557 0.02177 0.189 0.01279 0.01272 0.01439 0.03836 0.02018 0.212 0.01176 0.01191 0.01344 0.03529 0.01897 0.181 0.01084 0.01121 0.01255 0.03253 0.01358 0.129 0.00664 0.00786 0.0114 0.01992 0.01484 0.133 0.00754 0.0095 0.01285 0.02261 0.01472 0.133 0.00748 0.00905 0.01148 0.02245 0.01657 0.145 0.00881 0.01062 0.01318 0.02643 0.01503 0.137 0.00812 0.00933 0.01133 0.02436 0.01725 0.155 0.00874 0.01021 0.01331 0.02623 0.01469 0.132 0.00728 0.00886 0.0123 0.02184 0.01574 0.142 0.00839 0.00956 0.01309 0.02518 0.0145 0.131 0.00725 0.00876 0.01263 0.02175 0.02551 0.237 0.01321 0.01574 0.02148 0.03964 0.01831 0.163 0.0095 0.01103 0.01559 0.02849 0.02145 0.198 0.01155 0.01341 0.01666 0.03464 0.01909 0.171 0.00864 0.01223 0.01949 0.02592 0.01795 0.163 0.0081 0.01144 0.01756 0.02429 0.01564 0.136 0.00667 0.0099 0.01691 0.02001 0.0166 0.154 0.0082 0.00972 0.01491 0.0246 0.013 0.117 0.00631 0.00789 0.01144 0.01892 0.01185 0.106 0.00557 0.00721 0.01095 0.01672 0.02574 0.255 0.01454 0.01582 0.01758 0.04363 0.04087 0.405 0.02336 0.02498 0.02745 0.07008 0.02751 0.263 0.01604 0.01657 0.01879 0.04812 0.02308 0.256 0.01268 0.01365 0.01667 0.03804 0.02296 0.241 0.01265 0.01321 0.01588 0.03794 0.01884 0.19 0.01026 0.01161 0.01373 0.03078
Names of X columns:
MDVP:Shimmer MDVP:Shimmer(dB) Shimmer:APQ3 Shimmer:APQ5 MDVP:APQ Shimmer:DDA
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation