source('/home/pw/wessanet/cretab')



myrfcuid = ''

x <- array(list(70.8,66.67,69.6,66.33,69.87,64.33,67.47,64,67.6,63.33,67.13,61.33,66.27,64.67,66.73,63,68.07,60.67,67.8,63.67,64.8,60.67,64.6,61.67,64.2,62.33,64.2,60.33,63.67,59.67,61,60.33,59.67,59.33,59.67,58.67,59.8,58.67,60.73,59.33,59.4,57.33,58.07,59.33,57.47,56,70.73,53.67,72.87,58.67,66,49.33,66.07,71.33,66,70.33,66.27,69,64,66,63.67,66,63.73,63.33,63.33,65.33,63.53,64.33,63.53,64,62.87,61.67,59.53,63.67,62.8,64.67,60.8,61.67,59.8,62,56.67,61.33,57.67,63.67,58.4,61.33,55.47,62.33,56.2,59.67,71.53,59.33,68.67,61.67,65.67,58.67,66.73,58,67.33,56.67,66.73,59.67,66.87,58,65.8,57,64.73,57.67,65.47,58.67,63.6,55.33,64.07,56,64.67,55.67,63.73,53.33,62.53,53.67,61.93,51,62.67,47,62.8,4.33,61.33,70,62.6,68.67,59.13,67.67,61.27,66,59.47,65.67,57.87,65.67,59.73,63.67,61.4,63.67,58.8,64,58.33,62,57.47,62,57.13,61.67,55,61.67,51.53,63.33,72.73,61,73,62.33,70.8,60.33,70.07,60.33,71.67,60.67,71.07,57.67,70.67,58.33,70.73,58,70.73,57.33,68.6,56.67,69.6,58,66.47,55.33,67.07,55.67,68.67,54.67,66.93,56.33,65.93,55,68.87,55,66.53,54.67,65.8,54.33,66.6,49,66,48.33,65,49.67,66.8,43.67,65.6,6.33,66,3,65.67,66.67,64.67,67.33,65.07,65.33,64.67,66,65.07,65.67,65.2,66.67,64.87,65.67,63.47,65,62.6,64.67,64.07,66.67,63.73,63.67,64.67,63.33,61.6,63.67,61.6,63.33,60.47,63.67,61.27,63,63,61.67,61.47,61.33,60.87,60.67,61.67,60,62.87,61.67,62.4,61.33,59.73,58.67,60.13,60.33,58.8,59.67,59.6,59.33,58.93,59.67,60.13,61,58.2,61,58.27,60,58.27,60,55.07,58.67,53.87,58.33,52.33,58,47.2,56.33,37.93,54.67,72.73,55.33,70.07,54,70.67,52.67,72.07,44,68.8,65.67,68.8,65,67.47,66.33,66.73,64,66.53,62.33,66,61.33,67.6,63,66,63.67,66,62,66.53,61.33,65.8,64.67,64.27,62.67,64.67,64,64.6,61,64.13,60.67,65.47,59.67,62.93,60.33,63.53,56.67,62.13,56.67,63.87,54.33,64.67,51,63.33,51,63.13,47,62.8,68,62.4,65,62.4,64,62.6,64,61.47,64,62.2,62,63,61,61.8,60,59.73,60,60.33,62,60.13,60,59.53,59,59,61,55.93,60,41.87,60,36.33,58,71.67,58,71.47,60,70.47,58,69.53,59,70.73,56,69.93,54,68.73,51,67.53,47),dim=c(2,189),dimnames=list(c('X1','X2'),1:189))
 y <- array(NA,dim=c(2,189),dimnames=list(c('X1','X2'),1:189))
 for (i in 1:dim(x)[1])
 {
 	for (j in 1:dim(x)[2])
 	{
 		y[i,j] <- as.numeric(x[i,j])
 	}
 }
par7 = '0'
par6 = '0.5'
par5 = '0.5'
par4 = '0.5'
par3 = '1'
par2 = '1'
par1 = '1'
ylab = 'height'
xlab = 'cases'
main = 'Dendrogram'
par7 <- '0'
par6 <- '0.5'
par5 <- '0.5'
par4 <- '0.5'
par3 <- '1'
par2 <- '1'
par1 <- '1'
#'GNU S' R Code compiled by R2WASP v. 1.2.327 (Sun, 16 Jul 2017 10:30:11 +0200)
#Author: root
#To cite this work: Wessa, P., (2017), Agglomerative Nesting (v1.0.5) in Free Statistics Software (v$_version), Office for Research Development and Education, URL https://www.wessa.net/rwasp_agglomerativehierarchicalclustering.wasp/
#Source of accompanying publication: Office for Research, Development, and Education
#
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)
postscript(file="/home/pw/wessanet/rcomp/tmp/16obw1638205031.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
plot(ag,which.plots=2,main=main,sub=mysub,xlab=xlab,ylab=ylab)
dev.off()
postscript(file="/home/pw/wessanet/rcomp/tmp/2ly301638205031.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
plot(ag,which.plots=1,main='Banner',sub=mysub,xlab=ylab,ylab=xlab)
dev.off()

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="/home/pw/wessanet/rcomp/tmp/3pyon1638205031.tab") 

try(system("convert /home/pw/wessanet/rcomp/tmp/16obw1638205031.ps /home/pw/wessanet/rcomp/tmp/16obw1638205031.png",intern=TRUE))
try(system("convert /home/pw/wessanet/rcomp/tmp/2ly301638205031.ps /home/pw/wessanet/rcomp/tmp/2ly301638205031.png",intern=TRUE))