FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_hierarchicalclustering.wasp
Title produced by softwareHierarchical Clustering
Date of computationThu, 13 Nov 2008 08:39:16 -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/2008/Nov/13/t1226590790c12pnjgur4uziw0.htm/, Retrieved Sat, 18 May 2024 23:42:17 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=24653, Retrieved Sat, 18 May 2024 23:42:17 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F       [Hierarchical Clustering] [] [2008-11-13 15:39:16] [1768685c15539a739b6b33586be71b78] [Current]
Feedback Forum
2008-11-18 13:16:56 [Julie Govaerts] [reply] 
maakt in periodes groepen die gelijkaardig zijn = de tijdsreeks wordt gesplitst

dit is eerder een exploratief instrument (bv voor in de marketing = welke producten horen samen in 1 groep?)
hier zie je dat er 2 grote vertakkingen zijn die nog verder vertakken in kleinere groepjes van gegevens
2008-11-24 14:04:01 [Sören Van Donink] [reply] 
Deze berekening toont verschillende groepen in gelijke peridoen

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
103.0	101.5	103.2
103.5	101.3	102.9
103.9	101.1	102.6
104.4	101.1	102.3
103.5	101.4	102.1
104.7	101.4	101.5
101.6	101.3	101.1
102.3	101.2	101.9
101.2	101.2	101.5
101.4	101.4	101.5
102.7	101.1	102.0
101.5	101.1	101.5
102.1	101.0	101.6
100.5	101.0	101.5
101.7	101.0	101.9
103.6	100.6	102.0
103.2	100.6	102.0
101.7	100.7	102.2
105.6	100.8	102.2
103.7	100.9	102.3
104.3	100.8	102.0
105.9	100.8	101.8
104.8	100.8	101.8
106.0	100.8	101.6
105.6	100.9	101.9
106.7	101.0	102.3
103.3	101.0	102.4
103.0	101.5	102.1
101.6	101.1	102.3
101.7	101.0	102.2

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001
R Framework error message
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.

\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 & 3 seconds \tabularnewline
R Server & 'Herman Ole Andreas Wold' @ 193.190.124.10:1001 \tabularnewline
R Framework error message & 
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.
 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24653&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]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Herman Ole Andreas Wold' @ 193.190.124.10:1001[/C][/ROW]
[ROW][C]R Framework error message[/C][C]
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.
[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24653&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=24653&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 time3 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001
R Framework error message
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.






Summary of Dendrogram
LabelHeight
10.173205080756886
20.223606797749979
30.282842712474623
40.316227766016838
50.327356117197582
60.399999999999991
70.41231056256176
80.412310562561777
90.424264068711924
100.425529493977513
110.435889894354065
120.509901951359277
130.572280807449562
140.592935994404436
150.616441400296897
160.630746517572738
170.678232998312532
180.761420060860425
190.957267660893179
201.04291197594972
211.28793589582476
221.32984860661636
231.44836779444046
241.85695257445782
251.8860898377495
263.27733684673851
274.14266950887764
2811.0568565864898
2924.3683561071989

\begin{tabular}{lllllllll}
\hline
Summary of Dendrogram \tabularnewline
Label & Height \tabularnewline
1 & 0.173205080756886 \tabularnewline
2 & 0.223606797749979 \tabularnewline
3 & 0.282842712474623 \tabularnewline
4 & 0.316227766016838 \tabularnewline
5 & 0.327356117197582 \tabularnewline
6 & 0.399999999999991 \tabularnewline
7 & 0.41231056256176 \tabularnewline
8 & 0.412310562561777 \tabularnewline
9 & 0.424264068711924 \tabularnewline
10 & 0.425529493977513 \tabularnewline
11 & 0.435889894354065 \tabularnewline
12 & 0.509901951359277 \tabularnewline
13 & 0.572280807449562 \tabularnewline
14 & 0.592935994404436 \tabularnewline
15 & 0.616441400296897 \tabularnewline
16 & 0.630746517572738 \tabularnewline
17 & 0.678232998312532 \tabularnewline
18 & 0.761420060860425 \tabularnewline
19 & 0.957267660893179 \tabularnewline
20 & 1.04291197594972 \tabularnewline
21 & 1.28793589582476 \tabularnewline
22 & 1.32984860661636 \tabularnewline
23 & 1.44836779444046 \tabularnewline
24 & 1.85695257445782 \tabularnewline
25 & 1.8860898377495 \tabularnewline
26 & 3.27733684673851 \tabularnewline
27 & 4.14266950887764 \tabularnewline
28 & 11.0568565864898 \tabularnewline
29 & 24.3683561071989 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24653&T=1

[TABLE]
[ROW][C]Summary of Dendrogram[/C][/ROW]
[ROW][C]Label[/C][C]Height[/C][/ROW]
[ROW][C]1[/C][C]0.173205080756886[/C][/ROW]
[ROW][C]2[/C][C]0.223606797749979[/C][/ROW]
[ROW][C]3[/C][C]0.282842712474623[/C][/ROW]
[ROW][C]4[/C][C]0.316227766016838[/C][/ROW]
[ROW][C]5[/C][C]0.327356117197582[/C][/ROW]
[ROW][C]6[/C][C]0.399999999999991[/C][/ROW]
[ROW][C]7[/C][C]0.41231056256176[/C][/ROW]
[ROW][C]8[/C][C]0.412310562561777[/C][/ROW]
[ROW][C]9[/C][C]0.424264068711924[/C][/ROW]
[ROW][C]10[/C][C]0.425529493977513[/C][/ROW]
[ROW][C]11[/C][C]0.435889894354065[/C][/ROW]
[ROW][C]12[/C][C]0.509901951359277[/C][/ROW]
[ROW][C]13[/C][C]0.572280807449562[/C][/ROW]
[ROW][C]14[/C][C]0.592935994404436[/C][/ROW]
[ROW][C]15[/C][C]0.616441400296897[/C][/ROW]
[ROW][C]16[/C][C]0.630746517572738[/C][/ROW]
[ROW][C]17[/C][C]0.678232998312532[/C][/ROW]
[ROW][C]18[/C][C]0.761420060860425[/C][/ROW]
[ROW][C]19[/C][C]0.957267660893179[/C][/ROW]
[ROW][C]20[/C][C]1.04291197594972[/C][/ROW]
[ROW][C]21[/C][C]1.28793589582476[/C][/ROW]
[ROW][C]22[/C][C]1.32984860661636[/C][/ROW]
[ROW][C]23[/C][C]1.44836779444046[/C][/ROW]
[ROW][C]24[/C][C]1.85695257445782[/C][/ROW]
[ROW][C]25[/C][C]1.8860898377495[/C][/ROW]
[ROW][C]26[/C][C]3.27733684673851[/C][/ROW]
[ROW][C]27[/C][C]4.14266950887764[/C][/ROW]
[ROW][C]28[/C][C]11.0568565864898[/C][/ROW]
[ROW][C]29[/C][C]24.3683561071989[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24653&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=24653&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.173205080756886
20.223606797749979
30.282842712474623
40.316227766016838
50.327356117197582
60.399999999999991
70.41231056256176
80.412310562561777
90.424264068711924
100.425529493977513
110.435889894354065
120.509901951359277
130.572280807449562
140.592935994404436
150.616441400296897
160.630746517572738
170.678232998312532
180.761420060860425
190.957267660893179
201.04291197594972
211.28793589582476
221.32984860661636
231.44836779444046
241.85695257445782
251.8860898377495
263.27733684673851
274.14266950887764
2811.0568565864898
2924.3683561071989


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')
}