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 computationSun, 09 Nov 2008 01:56:48 -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/09/t12262210920d6lpoik4k83rcy.htm/, Retrieved Sun, 19 May 2024 03:02:49 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=22675, Retrieved Sun, 19 May 2024 03:02:49 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsreeks 1; reeks 2;reeks 3;reeks 4
Estimated Impact206
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F     [Kendall tau Correlation Matrix] [Hypothesis Testin...] [2008-11-03 18:32:17] [6743688719638b0cb1c0a6e0bf433315]
- RMPD  [Bivariate Kernel Density Estimation] [Various EDA topic...] [2008-11-08 11:46:03] [6743688719638b0cb1c0a6e0bf433315]
F RMPD      [Hierarchical Clustering] [Various EDA metho...] [2008-11-09 08:56:48] [9b05d7ef5dbcfba4217d280d9092f628] [Current]
F RMPD        [Box-Cox Normality Plot] [Various EDA Q4] [2008-11-09 20:31:09] [6743688719638b0cb1c0a6e0bf433315]
- RMPD        [Maximum-likelihood Fitting - Normal Distribution] [Various eda q5] [2008-11-09 20:36:57] [6743688719638b0cb1c0a6e0bf433315]
Feedback Forum
2008-11-23 12:50:05 [Pieter Broos] [reply
Om tijdreeksen in clusters op te delen wordt gebruik gemaakt van een dendogram, of boomstructuur.
De bedoeling hiervan is om inzicht te krijgen in welke observaties van de tijdreeks gelijkaardig zijn. De puntjes in het dendrogram zijn observaties. Hier worden clusters van gemaakt en zo verder onderverdeeld. Deze methode wordt meestal gebruikt voor niet- tijdreeksen maar eerder in de marketing sector om bijvoorbeeld te zien welke producten samen horen.

Post a new message
Dataseries X:
569	117,8	77,9	85,7
580	113,5	60	61,9
578	121,2	99,5	104,9
565	130,4	95	107,9
547	115,2	105,6	95,6
555	117,9	102,5	79,8
562	110,7	93,3	94,8
561	107,6	97,3	93,7
555	124,3	127	108,1
544	115,1	111,7	96,9
537	112,5	96,4	88,8
543	127,9	133	106,7
594	117,4	72,2	86,8
611	119,3	95,8	69,8
613	130,4	124,1	110,9
611	126	127,6	105,4
594	125,4	110,7	99,2
595	130,5	104,6	84,4
591	115,9	112,7	87,2
589	108,7	115,3	91,9
584	124	139,4	97,9
573	119,4	119	94,5
567	118,6	97,4	85
569	131,3	154	100,3
621	111,1	81,5	78,7
629	124,8	88,8	65,8
628	132,3	127,7	104,8
612	126,7	105,1	96
595	131,7	114,9	103,3
597	130,9	106,4	82,9
593	122,1	104,5	91,4
590	113,2	121,6	94,5
580	133,6	141,4	109,3
574	119,2	99	92,1
573	129,4	126,7	99,3
573	131,4	134,1	109,6
620	117,1	81,3	87,5
626	130,5	88,6	73,1
620	132,3	132,7	110,7
588	140,8	132,9	111,6
566	137,5	134,4	110,7
557	128,6	103,7	84
561	126,7	119,7	101,6
549	120,8	115	102,1
532	139,3	132,9	113,9
526	128,6	108,5	99
511	131,3	113,9	100,4
499	136,3	142	109,5
555	128,8	97,7	93,1
565	133,2	92,2	77
542	136,3	128,8	108
527	151,1	134,9	119,9
510	145	128,2	105,9
514	134,4	114,8	78,2
517	135,7	117,9	100,3
508	128,7	119,1	102,2
493	129,2	120,7	97
490	138,6	129,1	101,3
469	132,7	117,6	89,2
478	132,5	129,2	93,3
528	135,2	99,1	86,4




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time5 seconds
R Server'George Udny Yule' @ 72.249.76.132

\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 & 5 seconds \tabularnewline
R Server & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=22675&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]5 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=22675&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=22675&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 time5 seconds
R Server'George Udny Yule' @ 72.249.76.132







Summary of Dendrogram
LabelHeight
13.10644491340182
25.27446679769624
36.78527818147495
46.9217049922689
58.1154174261094
68.22800097228968
78.66833317310772
89.35467797414748
99.53362470417207
109.73755616158387
1110.1158291800524
1210.6995327000762
1311.0837719211467
1411.1176927402039
1511.1269942032878
1611.2836164415492
1711.7175082675456
1812.7935987209316
1913.5029626378806
2013.6532047519987
2113.9061137633776
2214.2912560679599
2315.244999180059
2416.0269525882465
2516.8193341128595
2617.1662459495371
2719.4935455419710
2820.0022498734517
2920.0336167449125
3021.0595438480890
3122.0762446720005
3223.7975771603291
3323.9556004851981
3424.2371643624292
3524.9908236827359
3625.0168668958763
3725.6682683482934
3825.7766400098992
3926.3351708478073
4028.1081516181605
4128.9168438491848
4231.1928358628244
4332.9363646128943
4433.6722116314038
4535.0108618388809
4636.5513987414996
4738.9268424596249
4842.7107786497974
4949.2672534002783
5050.4039449846778
5160.0205150909532
5275.5401160570481
53106.383508431499
54113.264347553177
55127.052901228084
56141.116772317358
57221.984795688692
58338.522129849007
59358.319593075369
60976.283725306065

\begin{tabular}{lllllllll}
\hline
Summary of Dendrogram \tabularnewline
Label & Height \tabularnewline
1 & 3.10644491340182 \tabularnewline
2 & 5.27446679769624 \tabularnewline
3 & 6.78527818147495 \tabularnewline
4 & 6.9217049922689 \tabularnewline
5 & 8.1154174261094 \tabularnewline
6 & 8.22800097228968 \tabularnewline
7 & 8.66833317310772 \tabularnewline
8 & 9.35467797414748 \tabularnewline
9 & 9.53362470417207 \tabularnewline
10 & 9.73755616158387 \tabularnewline
11 & 10.1158291800524 \tabularnewline
12 & 10.6995327000762 \tabularnewline
13 & 11.0837719211467 \tabularnewline
14 & 11.1176927402039 \tabularnewline
15 & 11.1269942032878 \tabularnewline
16 & 11.2836164415492 \tabularnewline
17 & 11.7175082675456 \tabularnewline
18 & 12.7935987209316 \tabularnewline
19 & 13.5029626378806 \tabularnewline
20 & 13.6532047519987 \tabularnewline
21 & 13.9061137633776 \tabularnewline
22 & 14.2912560679599 \tabularnewline
23 & 15.244999180059 \tabularnewline
24 & 16.0269525882465 \tabularnewline
25 & 16.8193341128595 \tabularnewline
26 & 17.1662459495371 \tabularnewline
27 & 19.4935455419710 \tabularnewline
28 & 20.0022498734517 \tabularnewline
29 & 20.0336167449125 \tabularnewline
30 & 21.0595438480890 \tabularnewline
31 & 22.0762446720005 \tabularnewline
32 & 23.7975771603291 \tabularnewline
33 & 23.9556004851981 \tabularnewline
34 & 24.2371643624292 \tabularnewline
35 & 24.9908236827359 \tabularnewline
36 & 25.0168668958763 \tabularnewline
37 & 25.6682683482934 \tabularnewline
38 & 25.7766400098992 \tabularnewline
39 & 26.3351708478073 \tabularnewline
40 & 28.1081516181605 \tabularnewline
41 & 28.9168438491848 \tabularnewline
42 & 31.1928358628244 \tabularnewline
43 & 32.9363646128943 \tabularnewline
44 & 33.6722116314038 \tabularnewline
45 & 35.0108618388809 \tabularnewline
46 & 36.5513987414996 \tabularnewline
47 & 38.9268424596249 \tabularnewline
48 & 42.7107786497974 \tabularnewline
49 & 49.2672534002783 \tabularnewline
50 & 50.4039449846778 \tabularnewline
51 & 60.0205150909532 \tabularnewline
52 & 75.5401160570481 \tabularnewline
53 & 106.383508431499 \tabularnewline
54 & 113.264347553177 \tabularnewline
55 & 127.052901228084 \tabularnewline
56 & 141.116772317358 \tabularnewline
57 & 221.984795688692 \tabularnewline
58 & 338.522129849007 \tabularnewline
59 & 358.319593075369 \tabularnewline
60 & 976.283725306065 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=22675&T=1

[TABLE]
[ROW][C]Summary of Dendrogram[/C][/ROW]
[ROW][C]Label[/C][C]Height[/C][/ROW]
[ROW][C]1[/C][C]3.10644491340182[/C][/ROW]
[ROW][C]2[/C][C]5.27446679769624[/C][/ROW]
[ROW][C]3[/C][C]6.78527818147495[/C][/ROW]
[ROW][C]4[/C][C]6.9217049922689[/C][/ROW]
[ROW][C]5[/C][C]8.1154174261094[/C][/ROW]
[ROW][C]6[/C][C]8.22800097228968[/C][/ROW]
[ROW][C]7[/C][C]8.66833317310772[/C][/ROW]
[ROW][C]8[/C][C]9.35467797414748[/C][/ROW]
[ROW][C]9[/C][C]9.53362470417207[/C][/ROW]
[ROW][C]10[/C][C]9.73755616158387[/C][/ROW]
[ROW][C]11[/C][C]10.1158291800524[/C][/ROW]
[ROW][C]12[/C][C]10.6995327000762[/C][/ROW]
[ROW][C]13[/C][C]11.0837719211467[/C][/ROW]
[ROW][C]14[/C][C]11.1176927402039[/C][/ROW]
[ROW][C]15[/C][C]11.1269942032878[/C][/ROW]
[ROW][C]16[/C][C]11.2836164415492[/C][/ROW]
[ROW][C]17[/C][C]11.7175082675456[/C][/ROW]
[ROW][C]18[/C][C]12.7935987209316[/C][/ROW]
[ROW][C]19[/C][C]13.5029626378806[/C][/ROW]
[ROW][C]20[/C][C]13.6532047519987[/C][/ROW]
[ROW][C]21[/C][C]13.9061137633776[/C][/ROW]
[ROW][C]22[/C][C]14.2912560679599[/C][/ROW]
[ROW][C]23[/C][C]15.244999180059[/C][/ROW]
[ROW][C]24[/C][C]16.0269525882465[/C][/ROW]
[ROW][C]25[/C][C]16.8193341128595[/C][/ROW]
[ROW][C]26[/C][C]17.1662459495371[/C][/ROW]
[ROW][C]27[/C][C]19.4935455419710[/C][/ROW]
[ROW][C]28[/C][C]20.0022498734517[/C][/ROW]
[ROW][C]29[/C][C]20.0336167449125[/C][/ROW]
[ROW][C]30[/C][C]21.0595438480890[/C][/ROW]
[ROW][C]31[/C][C]22.0762446720005[/C][/ROW]
[ROW][C]32[/C][C]23.7975771603291[/C][/ROW]
[ROW][C]33[/C][C]23.9556004851981[/C][/ROW]
[ROW][C]34[/C][C]24.2371643624292[/C][/ROW]
[ROW][C]35[/C][C]24.9908236827359[/C][/ROW]
[ROW][C]36[/C][C]25.0168668958763[/C][/ROW]
[ROW][C]37[/C][C]25.6682683482934[/C][/ROW]
[ROW][C]38[/C][C]25.7766400098992[/C][/ROW]
[ROW][C]39[/C][C]26.3351708478073[/C][/ROW]
[ROW][C]40[/C][C]28.1081516181605[/C][/ROW]
[ROW][C]41[/C][C]28.9168438491848[/C][/ROW]
[ROW][C]42[/C][C]31.1928358628244[/C][/ROW]
[ROW][C]43[/C][C]32.9363646128943[/C][/ROW]
[ROW][C]44[/C][C]33.6722116314038[/C][/ROW]
[ROW][C]45[/C][C]35.0108618388809[/C][/ROW]
[ROW][C]46[/C][C]36.5513987414996[/C][/ROW]
[ROW][C]47[/C][C]38.9268424596249[/C][/ROW]
[ROW][C]48[/C][C]42.7107786497974[/C][/ROW]
[ROW][C]49[/C][C]49.2672534002783[/C][/ROW]
[ROW][C]50[/C][C]50.4039449846778[/C][/ROW]
[ROW][C]51[/C][C]60.0205150909532[/C][/ROW]
[ROW][C]52[/C][C]75.5401160570481[/C][/ROW]
[ROW][C]53[/C][C]106.383508431499[/C][/ROW]
[ROW][C]54[/C][C]113.264347553177[/C][/ROW]
[ROW][C]55[/C][C]127.052901228084[/C][/ROW]
[ROW][C]56[/C][C]141.116772317358[/C][/ROW]
[ROW][C]57[/C][C]221.984795688692[/C][/ROW]
[ROW][C]58[/C][C]338.522129849007[/C][/ROW]
[ROW][C]59[/C][C]358.319593075369[/C][/ROW]
[ROW][C]60[/C][C]976.283725306065[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=22675&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=22675&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
13.10644491340182
25.27446679769624
36.78527818147495
46.9217049922689
58.1154174261094
68.22800097228968
78.66833317310772
89.35467797414748
99.53362470417207
109.73755616158387
1110.1158291800524
1210.6995327000762
1311.0837719211467
1411.1176927402039
1511.1269942032878
1611.2836164415492
1711.7175082675456
1812.7935987209316
1913.5029626378806
2013.6532047519987
2113.9061137633776
2214.2912560679599
2315.244999180059
2416.0269525882465
2516.8193341128595
2617.1662459495371
2719.4935455419710
2820.0022498734517
2920.0336167449125
3021.0595438480890
3122.0762446720005
3223.7975771603291
3323.9556004851981
3424.2371643624292
3524.9908236827359
3625.0168668958763
3725.6682683482934
3825.7766400098992
3926.3351708478073
4028.1081516181605
4128.9168438491848
4231.1928358628244
4332.9363646128943
4433.6722116314038
4535.0108618388809
4636.5513987414996
4738.9268424596249
4842.7107786497974
4949.2672534002783
5050.4039449846778
5160.0205150909532
5275.5401160570481
53106.383508431499
54113.264347553177
55127.052901228084
56141.116772317358
57221.984795688692
58338.522129849007
59358.319593075369
60976.283725306065



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