par7 <- '0' par6 <- '0.5' par5 <- '0.5' par4 <- '0.5' par3 <- '1' par2 <- '1' par1 <- '1' par3 <- as.numeric(par3) par4 <- as.numeric(par4) par5 <- as.numeric(par5) par6 <- as.numeric(par6) par7 <- as.numeric(par7) library(cluster) if (par2 == 'flexible') { if (par3 == 1) pm <- c(par4) if (par3 == 3) pm <- c(par4,par5,par6) if (par3 == 4) pm <- c(par4,par5,par6,par7) ag <- agnes(t(y),metric=par1,method=par2,par.method=pm) } else { ag <- agnes(t(y),metric=par1,method=par2) } mysub <- paste('Method: ',par2) summary(ag) bitmap(file='test1.png') plot(ag,which.plots=2,main=main,sub=mysub,xlab=xlab,ylab=ylab) dev.off() bitmap(file='test2.png') plot(ag,which.plots=1,main='Banner',sub=mysub,xlab=ylab,ylab=xlab) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Agglomerative Nesting (Hierarchical Clustering)',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Agglomerative Coefficient',header=TRUE) a<-table.element(a,ag$ac) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|