Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
114971 307836 104263 527070 105531 300074 104241 509846 104919 304198 105141 514258 104782 306122 106018 516922 101281 300414 105866 507561 94545 292133 105944 492622 93248 290616 106379 490243 84031 280244 105082 469357 87486 285179 104915 477580 115867 305486 107026 528379 120327 305957 107306 533590 117008 293886 107051 517945 108811 289441 107922 506174
Names of X columns:
jongerdan25 tussen25en50 ouderdan50 totaal
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