Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
106.22 103.7 107.97 106.31 103.75 108.13 107.38 103.85 108.54 109.31 104.02 109.86 110.82 104.13 109.75 111.22 104.17 109.99 110.66 104.18 112.01 110.76 104.2 111.96 110.69 104.5 111.41 111.08 104.78 112.11 110.97 104.88 111.67 110.24 104.89 111.95 112.51 104.9 112.31 111.52 104.95 113.26 112.13 105.24 113.5 112.23 105.35 114.43 112.92 105.44 115.02 111.89 105.46 115.1 111.99 105.47 117.11 111.51 105.48 117.52 112.33 105.75 116.1 112.04 106.1 116.39 112.09 106.19 116.01 111.41 106.23 116.74 112.61 106.24 116.68 113.14 106.25 117.45 113.65 106.35 117.8 114.26 106.48 119.37 114.4 106.52 118.9 114.93 106.55 119.05 114.86 106.55 120.46 114.95 106.56 120.99 116.17 106.89 119.86 114.6 107.09 120.18 114.62 107.24 119.81 113.82 107.28 120.15 115.02 107.3 119.8 115.18 107.31 120.27 115.59 107.47 120.71 116.6 107.35 121.87 117.07 107.31 121.87 116.96 107.32 121.92 116.66 107.32 123.72 116.07 107.34 124.38 116.04 107.53 123.21 115.81 107.72 123.17 116.22 107.75 122.95 115.85 107.79 123.46 116.43 107.81 123.24 117.39 107.9 123.86 119.17 107.8 124.28 119.24 107.86 124.78 120.03 107.8 125.19 119.34 107.74 125.46 118.49 107.75 127.6 118.59 107.83 127.8 117.5 107.8 126.63 117.56 107.81 127.06 118.25 107.86 126.77 118.01 107.83 127.05
Names of X columns:
V&D K&S HCR
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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation