Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1178 476049 2141 474605 2238 470439 2685 461251 4341 454724 5376 455626 4478 516847 6404 525192 4617 522975 3024 518585 1897 509239 2075 512238 1351 519164 2211 517009 2453 509933 3042 509127 4765 500857 4992 506971 4601 569323 6266 579714 4812 577992 3159 565464 1916 547344 2237 554788 1595 562325 2453 560854 2226 555332 3597 543599 4706 536662 4974 542722 5756 593530 5493 610763 5004 612613 3225 611324 2006 594167 2291 595454 1588 590865 2105 589379 2191 584428 3591 573100 4668 567456 4885 569028 5822 620735 5599 628884 5340 628232 3082 612117 2010 595404 2301 597141 1514 593408 1979 590072 2480 579799 3499 574205 4676 572775 5585 572942 5610 619567 5796 625809 6199 619916 3030 587625 1930 565742 2552 557274
Names of X columns:
h w
Method
ward
ward
single
complete
average
mcquitty
median
centroid
# of (top) clusters to display
ALL
ALL
2
3
4
5
6
7
8
9
10
15
20
Horizontal
FALSE
FALSE
TRUE
Triangle
FALSE
FALSE
TRUE
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
par3 <- as.logical(par3) par4 <- as.logical(par4) if (par3 == 'TRUE'){ dum = xlab xlab = ylab ylab = dum } x <- t(y) hc <- hclust(dist(x),method=par1) d <- as.dendrogram(hc) str(d) mysub <- paste('Method: ',par1) bitmap(file='test1.png') if (par4 == 'TRUE'){ plot(d,main=main,ylab=ylab,xlab=xlab,horiz=par3, nodePar=list(pch = c(1,NA), cex=0.8, lab.cex = 0.8),type='t',center=T, sub=mysub) } else { plot(d,main=main,ylab=ylab,xlab=xlab,horiz=par3, nodePar=list(pch = c(1,NA), cex=0.8, lab.cex = 0.8), sub=mysub) } dev.off() if (par2 != 'ALL'){ if (par3 == 'TRUE'){ ylab = 'cluster' } else { xlab = 'cluster' } par2 <- as.numeric(par2) memb <- cutree(hc, k = par2) cent <- NULL for(k in 1:par2){ cent <- rbind(cent, colMeans(x[memb == k, , drop = FALSE])) } hc1 <- hclust(dist(cent),method=par1, members = table(memb)) de <- as.dendrogram(hc1) bitmap(file='test2.png') if (par4 == 'TRUE'){ plot(de,main=main,ylab=ylab,xlab=xlab,horiz=par3, nodePar=list(pch = c(1,NA), cex=0.8, lab.cex = 0.8),type='t',center=T, sub=mysub) } else { plot(de,main=main,ylab=ylab,xlab=xlab,horiz=par3, nodePar=list(pch = c(1,NA), cex=0.8, lab.cex = 0.8), sub=mysub) } dev.off() str(de) } load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Summary of Dendrogram',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Label',header=TRUE) a<-table.element(a,'Height',header=TRUE) a<-table.row.end(a) num <- length(x[,1])-1 for (i in 1:num) { a<-table.row.start(a) a<-table.element(a,hc$labels[i]) a<-table.element(a,hc$height[i]) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable1.tab') if (par2 != 'ALL'){ a<-table.start() a<-table.row.start(a) a<-table.element(a,'Summary of Cut Dendrogram',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Label',header=TRUE) a<-table.element(a,'Height',header=TRUE) a<-table.row.end(a) num <- par2-1 for (i in 1:num) { a<-table.row.start(a) a<-table.element(a,i) a<-table.element(a,hc1$height[i]) 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