FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_hierarchicalclustering.wasp
Title produced by softwareHierarchical Clustering
Date of computationSat, 23 Jan 2010 11:54:54 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Jan/23/t12642729387zyvb4ulcak3uw2.htm/, Retrieved Sun, 05 May 2024 06:11:00 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=72372, Retrieved Sun, 05 May 2024 06:11:00 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact168

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Hierarchical Clustering] [VAC hierarchical ...] [2008-12-17 08:40:02] [379d6c32f73e3218fd773d79e4063d07]
-  M      [Hierarchical Clustering] [Dendogram] [2010-01-23 18:54:54] [f32a893c5a60da9308cd5d37e6977c4f] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
124.1	188.5
124.4	188.6
115.7	191.9
108.3	193.5
102.3	194.9
104.6	194.9
104	196.2
103.5	196.2
96	198
96.6	198.6
95.4	201.3
92.1	203.5
93	204.1
90.4	204.8
93.3	206.5
97.1	207.8
111	208.6
114.1	209.7
113.3	210
111	211.7
107.2	212.4
118.3	213.7
134.1	214.8
139	216.4
116.7	217.5
112.5	218.6
122.8	220.4
130	221.8
125.6	222.5
123.8	223.4
135.8	225.5
136.4	226.5
135.3	227.8
149.5	228.5
159.6	229.1
161.4	229.9
175.2	230.8
199.5	231.9
245	236
257.8	237.5

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 1 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=72372&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=72372&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=72372&T=0

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135






Summary of Dendrogram
LabelHeight
10.316227766016847
20.5
30.848528137423849
40.854400374531754
51.08166538263920
61.16619037896906
71.92314716194683
81.96977156035924
92.01246117974981
102.31570927086406
112.49885023811262
122.86122877876728
133.09999999999999
143.51993819543067
153.77069804686231
163.92299682538441
174.12310562561767
184.28840107403477
195.15460958754395
206.58191037700063
217.29633877774428
227.63246027910915
237.6508634172364
247.89726286167676
2512.1391547839476
2612.8875909308140
2714.0766601627723
2818.0034363557983
2919.5109805940994
3022.6642080010320
3123.6505843432484
3232.2133099825052
3353.096970120444
3454.1660599548719
3599.264220758254
36120.353472977802
37200.380024343742
38284.858535685259
39588.966483108861

\begin{tabular}{lllllllll}
\hline
Summary of Dendrogram \tabularnewline
Label & Height \tabularnewline
1 & 0.316227766016847 \tabularnewline
2 & 0.5 \tabularnewline
3 & 0.848528137423849 \tabularnewline
4 & 0.854400374531754 \tabularnewline
5 & 1.08166538263920 \tabularnewline
6 & 1.16619037896906 \tabularnewline
7 & 1.92314716194683 \tabularnewline
8 & 1.96977156035924 \tabularnewline
9 & 2.01246117974981 \tabularnewline
10 & 2.31570927086406 \tabularnewline
11 & 2.49885023811262 \tabularnewline
12 & 2.86122877876728 \tabularnewline
13 & 3.09999999999999 \tabularnewline
14 & 3.51993819543067 \tabularnewline
15 & 3.77069804686231 \tabularnewline
16 & 3.92299682538441 \tabularnewline
17 & 4.12310562561767 \tabularnewline
18 & 4.28840107403477 \tabularnewline
19 & 5.15460958754395 \tabularnewline
20 & 6.58191037700063 \tabularnewline
21 & 7.29633877774428 \tabularnewline
22 & 7.63246027910915 \tabularnewline
23 & 7.6508634172364 \tabularnewline
24 & 7.89726286167676 \tabularnewline
25 & 12.1391547839476 \tabularnewline
26 & 12.8875909308140 \tabularnewline
27 & 14.0766601627723 \tabularnewline
28 & 18.0034363557983 \tabularnewline
29 & 19.5109805940994 \tabularnewline
30 & 22.6642080010320 \tabularnewline
31 & 23.6505843432484 \tabularnewline
32 & 32.2133099825052 \tabularnewline
33 & 53.096970120444 \tabularnewline
34 & 54.1660599548719 \tabularnewline
35 & 99.264220758254 \tabularnewline
36 & 120.353472977802 \tabularnewline
37 & 200.380024343742 \tabularnewline
38 & 284.858535685259 \tabularnewline
39 & 588.966483108861 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=72372&T=1

[TABLE]
[ROW][C]Summary of Dendrogram[/C][/ROW]
[ROW][C]Label[/C][C]Height[/C][/ROW]
[ROW][C]1[/C][C]0.316227766016847[/C][/ROW]
[ROW][C]2[/C][C]0.5[/C][/ROW]
[ROW][C]3[/C][C]0.848528137423849[/C][/ROW]
[ROW][C]4[/C][C]0.854400374531754[/C][/ROW]
[ROW][C]5[/C][C]1.08166538263920[/C][/ROW]
[ROW][C]6[/C][C]1.16619037896906[/C][/ROW]
[ROW][C]7[/C][C]1.92314716194683[/C][/ROW]
[ROW][C]8[/C][C]1.96977156035924[/C][/ROW]
[ROW][C]9[/C][C]2.01246117974981[/C][/ROW]
[ROW][C]10[/C][C]2.31570927086406[/C][/ROW]
[ROW][C]11[/C][C]2.49885023811262[/C][/ROW]
[ROW][C]12[/C][C]2.86122877876728[/C][/ROW]
[ROW][C]13[/C][C]3.09999999999999[/C][/ROW]
[ROW][C]14[/C][C]3.51993819543067[/C][/ROW]
[ROW][C]15[/C][C]3.77069804686231[/C][/ROW]
[ROW][C]16[/C][C]3.92299682538441[/C][/ROW]
[ROW][C]17[/C][C]4.12310562561767[/C][/ROW]
[ROW][C]18[/C][C]4.28840107403477[/C][/ROW]
[ROW][C]19[/C][C]5.15460958754395[/C][/ROW]
[ROW][C]20[/C][C]6.58191037700063[/C][/ROW]
[ROW][C]21[/C][C]7.29633877774428[/C][/ROW]
[ROW][C]22[/C][C]7.63246027910915[/C][/ROW]
[ROW][C]23[/C][C]7.6508634172364[/C][/ROW]
[ROW][C]24[/C][C]7.89726286167676[/C][/ROW]
[ROW][C]25[/C][C]12.1391547839476[/C][/ROW]
[ROW][C]26[/C][C]12.8875909308140[/C][/ROW]
[ROW][C]27[/C][C]14.0766601627723[/C][/ROW]
[ROW][C]28[/C][C]18.0034363557983[/C][/ROW]
[ROW][C]29[/C][C]19.5109805940994[/C][/ROW]
[ROW][C]30[/C][C]22.6642080010320[/C][/ROW]
[ROW][C]31[/C][C]23.6505843432484[/C][/ROW]
[ROW][C]32[/C][C]32.2133099825052[/C][/ROW]
[ROW][C]33[/C][C]53.096970120444[/C][/ROW]
[ROW][C]34[/C][C]54.1660599548719[/C][/ROW]
[ROW][C]35[/C][C]99.264220758254[/C][/ROW]
[ROW][C]36[/C][C]120.353472977802[/C][/ROW]
[ROW][C]37[/C][C]200.380024343742[/C][/ROW]
[ROW][C]38[/C][C]284.858535685259[/C][/ROW]
[ROW][C]39[/C][C]588.966483108861[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=72372&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=72372&T=1

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Summary of Dendrogram
LabelHeight
10.316227766016847
20.5
30.848528137423849
40.854400374531754
51.08166538263920
61.16619037896906
71.92314716194683
81.96977156035924
92.01246117974981
102.31570927086406
112.49885023811262
122.86122877876728
133.09999999999999
143.51993819543067
153.77069804686231
163.92299682538441
174.12310562561767
184.28840107403477
195.15460958754395
206.58191037700063
217.29633877774428
227.63246027910915
237.6508634172364
247.89726286167676
2512.1391547839476
2612.8875909308140
2714.0766601627723
2818.0034363557983
2919.5109805940994
3022.6642080010320
3123.6505843432484
3232.2133099825052
3353.096970120444
3454.1660599548719
3599.264220758254
36120.353472977802
37200.380024343742
38284.858535685259
39588.966483108861


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = ward ; par2 = ALL ; par3 = FALSE ; par4 = FALSE ;
Parameters (R input):
par1 = ward ; par2 = ALL ; par3 = FALSE ; par4 = FALSE ;
R code (references can be found in the software module):
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')
}