Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
-0.35 -0.76 -0.57 -0.69 -0.16 0.48 0.00 -0.02 -0.19 -0.53 -0.35 -0.56 -0.19 -0.84 0.28 0.55 -0.12 0.15 0.26 0.29 -0.04 0.70 -0.13 -0.47 0.53 0.10 -0.39 -1.19 -0.19 1.72 0.09 0.22 -0.30 0.13 1.33 3.46 -0.55 -1.30 -0.73 -0.52 -0.38 -1.50 -0.67 -0.86 -0.62 -1.00 -0.74 -0.06 -0.26 0.28 0.26 0.86 0.25 1.04 1.30 0.64 0.18 1.90 0.23 -0.35 2.22 1.74 0.43 -1.16 -0.62 0.34 -0.40 2.36 -0.56 -0.85 -0.58 0.24 -0.09 -0.60 -0.19 -0.49 -0.39 0.18 -0.65 -0.81 -0.34 -0.65 -0.30 -0.06 -0.57 -1.07 -0.71 -0.33 -0.35 -0.64 -0.51 -0.53 -0.53 -0.62 -0.70 -0.53 -0.59 1.18 -0.59 0.60 -0.61 -0.35 -0.46 1.74 -0.54 -0.94 -0.73 -0.60 -0.67 -0.70 -0.83 -0.24 4.39 1.24 3.58 -0.74 0.04 0.65 -0.24 -0.73 -0.06 -0.59 0.74 0.81 0.33 -0.89 -0.52 -1.23 1.87 2.30 2.98 -0.11 -0.55 -0.29 -0.52 0.50
Names of X columns:
Population Density TotalKWH KWHpercapita
Method
complete
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
par4 <- 'FALSE' par3 <- 'FALSE' par2 <- 'ALL' par1 <- 'complete' 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