Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_hierarchicalclustering.wasp
Title produced by softwareHierarchical Clustering
Date of computationWed, 12 Nov 2008 11:30:59 -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/12/t1226514848f4un9qqcpxhxtbz.htm/, Retrieved Tue, 28 May 2024 05:17:27 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=24361, Retrieved Tue, 28 May 2024 05:17:27 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordshierarchical clustering
Estimated Impact146
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F       [Hierarchical Clustering] [workshop 2 Q2 Hie...] [2008-11-12 18:30:59] [74c7506a1ea162af3aa8be25bcd05d28] [Current]
Feedback Forum
2008-11-15 11:57:42 [58d427c57bd46519a715a3a7fea6a80f] [reply
De dendrogram geeft weer: in bepaalde periodes zijn er groepen gemaakt die gelijkaardig zijn . Elke cluster werd opnieuw onderverdeeld tot op1 enkele periode.
2008-11-19 18:43:07 [Bénédicte Soens] [reply
Via deze methode worden er groepen of clusters gevormd, daarbij wordt nagegaan of er met de periodes groepen kunnen worden gevormd. Zo worden er al 2 grote groepen/clusters gevormd (=tijdreeksen opsplitsen met een knooppunt): de eerste telt de eerste 18 maanden en de volgende telt dan de rest. Dus er is een verandering tussen de eerste 18 en laatste maanden. Vervolgens worden er steeds nieuwe opsplitsingen gemaakt tot men op 1 element uitkomt.
2008-11-24 11:53:24 [Julian De Ruyter] [reply
Geen conclusie gegeven, wel een berekening.
Een dendogram geeft gelijkaardige groepen weer per periode. Dit wordt gedaan door middel van clusters (boom-tak verdeling)
38tot46 liggen in dezelfde hoofdcluster (hoofdopsplitsing).
Verder is er geen duidelijk patroon te zien en te verklaren.

Post a new message
Dataseries X:
8,4	7,9	9,1	23,6	6,7
8,4	7,9	9	22,3	6,8
8,6	8,1	9,3	21,8	7,2
8,9	8,2	9,9	20,8	7,6
8,8	8	9,8	19,7	7,6
8,3	7,5	9,4	18,3	7,3
7,5	6,8	8,3	17,4	6,4
7,2	6,5	8	17	6,1
7,5	6,6	8,5	18,1	6,3
8,8	7,6	10,4	23,9	7,1
9,3	8	11,1	25,6	7,5
9,3	8	10,9	25,3	7,4
8,7	7,7	9,9	23,6	7,1
8,2	7,5	9,2	21,9	6,8
8,3	7,6	9,2	21,4	6,9
8,5	7,7	9,5	20,6	7,2
8,6	7,9	9,6	20,5	7,4
8,6	7,8	9,5	20,2	7,3
8,2	7,5	9,1	20,6	6,9
8,1	7,5	8,9	19,7	6,9
8	7,1	9	19,3	6,8
8,6	7,5	10,1	22,8	7,1
8,7	7,5	10,3	23,5	7,2
8,8	7,6	10,2	23,8	7,1
8,5	7,7	9,6	22,6	7
8,4	7,7	9,2	22	6,9
8,5	7,9	9,3	21,7	7
8,7	8,1	9,4	20,7	7,4
8,7	8,2	9,4	20,2	7,5
8,6	8,2	9,2	19,1	7,5
8,5	8,1	9	19,5	7,4
8,3	7,9	9	18,7	7,3
8,1	7,3	9	18,6	7
8,2	6,9	9,8	22,2	6,7
8,1	6,6	10	23,2	6,5
8,1	6,7	9,9	23,5	6,5
7,9	6,9	9,3	21,3	6,5
7,9	7	9	20	6,6
7,9	7,1	9	18,7	6,8
8	7,2	9,1	18,9	6,9
8	7,1	9,1	18,3	6,9
7,9	6,9	9,1	18,4	6,8
8	7	9,2	19,9	6,8
7,7	6,8	8,8	19,2	6,5
7,2	6,4	8,3	18,5	6,1
7,5	6,7	8,4	20,9	6
7,3	6,7	8,1	20,5	5,9
7	6,4	7,8	19,4	5,8
7	6,3	7,9	18,1	5,9
7	6,2	7,9	17	5,9
7,2	6,5	8	17	6,2
7,3	6,8	7,9	17,3	6,3
7,1	6,8	7,5	16,7	6,2
6,8	6,5	7,2	15,5	6
6,6	6,3	6,9	15,3	5,8
6,2	5,9	6,6	13,7	5,5
6,2	5,9	6,7	14,1	5,5
6,8	6,4	7,3	17,3	5,7
6,9	6,4	7,5	18,1	5,8
6,8	6,5	7,2	18,1	5,7
6,7




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

\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
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24361&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]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24361&T=0

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

As an alternative you can also use a 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







Summary of Dendrogram
LabelHeight
10.100000000000001
20.223606797749979
30.264575131106458
40.282842712474619
50.316227766016838
60.316227766016838
70.316227766016839
80.331662479035540
90.331662479035541
100.346410161513775
110.374165738677394
120.386459265903657
130.41097921090005
140.412310562561767
150.463976882780912
160.469041575982343
170.47958315233127
180.5
190.519615242270663
200.547722557505166
210.556274035168819
220.565685424949238
230.569231480695466
240.591607978309961
250.599070478491458
260.60827625302982
270.662891260736982
280.69282032302755
290.695806520791897
300.737587156526468
310.756273476694708
320.778107981368538
330.92107229414546
340.959583847541467
351.01441743544059
361.04383129176540
371.06770782520313
381.08599308836795
391.15012750131221
401.18477405730313
411.21920736167592
421.59341137253781
431.62398355679715
441.78657966302433
452.31201889558532
462.35449770556043
472.49798645884276
482.85703888142444
493.06637949851724
503.33266274194377
513.62951383442793
524.55429912979917
537.08340808870633
547.81315977173791
559.02130808142497
5614.1409332977778
5714.3937719498212
5835.0003797266385
5937.9256528042753
6071.5806833590777

\begin{tabular}{lllllllll}
\hline
Summary of Dendrogram \tabularnewline
Label & Height \tabularnewline
1 & 0.100000000000001 \tabularnewline
2 & 0.223606797749979 \tabularnewline
3 & 0.264575131106458 \tabularnewline
4 & 0.282842712474619 \tabularnewline
5 & 0.316227766016838 \tabularnewline
6 & 0.316227766016838 \tabularnewline
7 & 0.316227766016839 \tabularnewline
8 & 0.331662479035540 \tabularnewline
9 & 0.331662479035541 \tabularnewline
10 & 0.346410161513775 \tabularnewline
11 & 0.374165738677394 \tabularnewline
12 & 0.386459265903657 \tabularnewline
13 & 0.41097921090005 \tabularnewline
14 & 0.412310562561767 \tabularnewline
15 & 0.463976882780912 \tabularnewline
16 & 0.469041575982343 \tabularnewline
17 & 0.47958315233127 \tabularnewline
18 & 0.5 \tabularnewline
19 & 0.519615242270663 \tabularnewline
20 & 0.547722557505166 \tabularnewline
21 & 0.556274035168819 \tabularnewline
22 & 0.565685424949238 \tabularnewline
23 & 0.569231480695466 \tabularnewline
24 & 0.591607978309961 \tabularnewline
25 & 0.599070478491458 \tabularnewline
26 & 0.60827625302982 \tabularnewline
27 & 0.662891260736982 \tabularnewline
28 & 0.69282032302755 \tabularnewline
29 & 0.695806520791897 \tabularnewline
30 & 0.737587156526468 \tabularnewline
31 & 0.756273476694708 \tabularnewline
32 & 0.778107981368538 \tabularnewline
33 & 0.92107229414546 \tabularnewline
34 & 0.959583847541467 \tabularnewline
35 & 1.01441743544059 \tabularnewline
36 & 1.04383129176540 \tabularnewline
37 & 1.06770782520313 \tabularnewline
38 & 1.08599308836795 \tabularnewline
39 & 1.15012750131221 \tabularnewline
40 & 1.18477405730313 \tabularnewline
41 & 1.21920736167592 \tabularnewline
42 & 1.59341137253781 \tabularnewline
43 & 1.62398355679715 \tabularnewline
44 & 1.78657966302433 \tabularnewline
45 & 2.31201889558532 \tabularnewline
46 & 2.35449770556043 \tabularnewline
47 & 2.49798645884276 \tabularnewline
48 & 2.85703888142444 \tabularnewline
49 & 3.06637949851724 \tabularnewline
50 & 3.33266274194377 \tabularnewline
51 & 3.62951383442793 \tabularnewline
52 & 4.55429912979917 \tabularnewline
53 & 7.08340808870633 \tabularnewline
54 & 7.81315977173791 \tabularnewline
55 & 9.02130808142497 \tabularnewline
56 & 14.1409332977778 \tabularnewline
57 & 14.3937719498212 \tabularnewline
58 & 35.0003797266385 \tabularnewline
59 & 37.9256528042753 \tabularnewline
60 & 71.5806833590777 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24361&T=1

[TABLE]
[ROW][C]Summary of Dendrogram[/C][/ROW]
[ROW][C]Label[/C][C]Height[/C][/ROW]
[ROW][C]1[/C][C]0.100000000000001[/C][/ROW]
[ROW][C]2[/C][C]0.223606797749979[/C][/ROW]
[ROW][C]3[/C][C]0.264575131106458[/C][/ROW]
[ROW][C]4[/C][C]0.282842712474619[/C][/ROW]
[ROW][C]5[/C][C]0.316227766016838[/C][/ROW]
[ROW][C]6[/C][C]0.316227766016838[/C][/ROW]
[ROW][C]7[/C][C]0.316227766016839[/C][/ROW]
[ROW][C]8[/C][C]0.331662479035540[/C][/ROW]
[ROW][C]9[/C][C]0.331662479035541[/C][/ROW]
[ROW][C]10[/C][C]0.346410161513775[/C][/ROW]
[ROW][C]11[/C][C]0.374165738677394[/C][/ROW]
[ROW][C]12[/C][C]0.386459265903657[/C][/ROW]
[ROW][C]13[/C][C]0.41097921090005[/C][/ROW]
[ROW][C]14[/C][C]0.412310562561767[/C][/ROW]
[ROW][C]15[/C][C]0.463976882780912[/C][/ROW]
[ROW][C]16[/C][C]0.469041575982343[/C][/ROW]
[ROW][C]17[/C][C]0.47958315233127[/C][/ROW]
[ROW][C]18[/C][C]0.5[/C][/ROW]
[ROW][C]19[/C][C]0.519615242270663[/C][/ROW]
[ROW][C]20[/C][C]0.547722557505166[/C][/ROW]
[ROW][C]21[/C][C]0.556274035168819[/C][/ROW]
[ROW][C]22[/C][C]0.565685424949238[/C][/ROW]
[ROW][C]23[/C][C]0.569231480695466[/C][/ROW]
[ROW][C]24[/C][C]0.591607978309961[/C][/ROW]
[ROW][C]25[/C][C]0.599070478491458[/C][/ROW]
[ROW][C]26[/C][C]0.60827625302982[/C][/ROW]
[ROW][C]27[/C][C]0.662891260736982[/C][/ROW]
[ROW][C]28[/C][C]0.69282032302755[/C][/ROW]
[ROW][C]29[/C][C]0.695806520791897[/C][/ROW]
[ROW][C]30[/C][C]0.737587156526468[/C][/ROW]
[ROW][C]31[/C][C]0.756273476694708[/C][/ROW]
[ROW][C]32[/C][C]0.778107981368538[/C][/ROW]
[ROW][C]33[/C][C]0.92107229414546[/C][/ROW]
[ROW][C]34[/C][C]0.959583847541467[/C][/ROW]
[ROW][C]35[/C][C]1.01441743544059[/C][/ROW]
[ROW][C]36[/C][C]1.04383129176540[/C][/ROW]
[ROW][C]37[/C][C]1.06770782520313[/C][/ROW]
[ROW][C]38[/C][C]1.08599308836795[/C][/ROW]
[ROW][C]39[/C][C]1.15012750131221[/C][/ROW]
[ROW][C]40[/C][C]1.18477405730313[/C][/ROW]
[ROW][C]41[/C][C]1.21920736167592[/C][/ROW]
[ROW][C]42[/C][C]1.59341137253781[/C][/ROW]
[ROW][C]43[/C][C]1.62398355679715[/C][/ROW]
[ROW][C]44[/C][C]1.78657966302433[/C][/ROW]
[ROW][C]45[/C][C]2.31201889558532[/C][/ROW]
[ROW][C]46[/C][C]2.35449770556043[/C][/ROW]
[ROW][C]47[/C][C]2.49798645884276[/C][/ROW]
[ROW][C]48[/C][C]2.85703888142444[/C][/ROW]
[ROW][C]49[/C][C]3.06637949851724[/C][/ROW]
[ROW][C]50[/C][C]3.33266274194377[/C][/ROW]
[ROW][C]51[/C][C]3.62951383442793[/C][/ROW]
[ROW][C]52[/C][C]4.55429912979917[/C][/ROW]
[ROW][C]53[/C][C]7.08340808870633[/C][/ROW]
[ROW][C]54[/C][C]7.81315977173791[/C][/ROW]
[ROW][C]55[/C][C]9.02130808142497[/C][/ROW]
[ROW][C]56[/C][C]14.1409332977778[/C][/ROW]
[ROW][C]57[/C][C]14.3937719498212[/C][/ROW]
[ROW][C]58[/C][C]35.0003797266385[/C][/ROW]
[ROW][C]59[/C][C]37.9256528042753[/C][/ROW]
[ROW][C]60[/C][C]71.5806833590777[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24361&T=1

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

As an alternative you can also use a QR Code:  

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

Summary of Dendrogram
LabelHeight
10.100000000000001
20.223606797749979
30.264575131106458
40.282842712474619
50.316227766016838
60.316227766016838
70.316227766016839
80.331662479035540
90.331662479035541
100.346410161513775
110.374165738677394
120.386459265903657
130.41097921090005
140.412310562561767
150.463976882780912
160.469041575982343
170.47958315233127
180.5
190.519615242270663
200.547722557505166
210.556274035168819
220.565685424949238
230.569231480695466
240.591607978309961
250.599070478491458
260.60827625302982
270.662891260736982
280.69282032302755
290.695806520791897
300.737587156526468
310.756273476694708
320.778107981368538
330.92107229414546
340.959583847541467
351.01441743544059
361.04383129176540
371.06770782520313
381.08599308836795
391.15012750131221
401.18477405730313
411.21920736167592
421.59341137253781
431.62398355679715
441.78657966302433
452.31201889558532
462.35449770556043
472.49798645884276
482.85703888142444
493.06637949851724
503.33266274194377
513.62951383442793
524.55429912979917
537.08340808870633
547.81315977173791
559.02130808142497
5614.1409332977778
5714.3937719498212
5835.0003797266385
5937.9256528042753
6071.5806833590777



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