Home » date » 2008 » Nov » 14 »

Hierarchical clustering own time series

*The author of this computation has been verified*
R Software Module: rwasp_hierarchicalclustering.wasp (opens new window with default values)
Title produced by software: Hierarchical Clustering
Date of computation: Thu, 13 Nov 2008 18:11:20 -0700
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2008/Nov/14/t1226625233ax7v76n8w6ptzgp.htm/, Retrieved Fri, 14 Nov 2008 01:13:53 +0000
 
BibTeX entries for LaTeX users:
@Manual{KEY,
    author = {{YOUR NAME}},
    publisher = {Office for Research Development and Education},
    title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2008/Nov/14/t1226625233ax7v76n8w6ptzgp.htm/},
    year = {2008},
}
@Manual{R,
    title = {R: A Language and Environment for Statistical Computing},
    author = {{R Development Core Team}},
    organization = {R Foundation for Statistical Computing},
    address = {Vienna, Austria},
    year = {2008},
    note = {{ISBN} 3-900051-07-0},
    url = {http://www.R-project.org},
}
 
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
 
Feedback Forum:

Post a new message
 
Original text written by user:
 
IsPrivate?
No (this computation is public)
 
User-defined keywords:
Hierarchical clustering own time series
 
Dataseries X:
» Textbox « » Textfile « » CSV «
1,1372 9924 789,6 1,1139 10371 773,3 1,1222 10846 804,3 1,1692 10413 817,8 1,1702 10709 836,7 1,2286 10662 721,8 1,2613 10570 760,8 1,2646 10297 841,4 1,2262 10635 1045,6 1,1985 10872 949,2 1,2007 10296 850,1 1,2138 10383 957,4 1,2266 10431 851,8 1,2176 10574 913,9 1,2218 10653 888 1,249 10805 973,8 1,2991 10872 927,6 1,3408 10625 833 1,3119 10407 879,5 1,3014 10463 797,3 1,3201 10556 834,5 1,2938 10646 735,1 1,2694 10702 835 1,2165 11353 892,8 1,2037 11346 697,2 1,2292 11451 821,1 1,2256 11964 732,7 1,2015 12574 797,6 1,1786 13031 866,3 1,1856 13812 826,3 1,2103 14544 778,6 1,1938 14931 779,2 1,202 14886 951 1,2271 16005 692,3 1,277 17064 841,4 1,265 15168 857,3 1,2684 16050 760,7 1,2811 15839 841,2 1,2727 15137 810,3 1,2611 14954 1007,4 1,2881 15648 931,3 1,3213 15305 931,2 1,2999 15579 855,8 1,3074 16348 858,4 1,3242 15928 925,9 1,3516 16171 930,7 1,3511 15937 1035,6 1,3419 15713 979,2 1,3716 15594 942,6 1,3622 15683 843,9 1,3896 16438 854,3 1,4227 17032 1029,8 1,4684 17696 944 1,457 17745 856,4 1,4718 19394 1059,4 1,4748 20148 959,3 1,5527 20108 941,5 1,575 18584 1026,4 1,5557 18441 921,3 1,5553 18391 968 1,577 19178 1129
 
Output produced by software:

Enter (or paste) a matrix (table) containing all data (time) series. Every column represents a different variable and must be delimited by a space or Tab. Every row represents a period in time (or category) and must be delimited by hard returns. The easiest way to enter data is to copy and paste a block of spreadsheet cells. Please, do not use commas or spaces to seperate groups of digits!


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24


Summary of Dendrogram
LabelHeight
17.20415440145477
28.75751581271771
320.8061109061737
421.6002342663222
536.6510201234564
643.7818006528968
754.0394992282497
855.1697106413474
956.3027535675654
1061.7172112218302
1168.4170300448653
1275.0179542339033
1381.8752814083104
1486.1664345672216
1588.3457044389256
1690.0933779855101
1794.638546640384
1898.608332196074
19100.373104614533
20104.678621892390
21110.068572824444
22111.553035372776
23118.439516387373
24121.428539319593
25153.661395965953
26163.015157902912
27185.441419946744
28189.608862107693
29190.35504249452
30191.098354855530
31198.278994155723
32222.992943829071
33226.936491263613
34229.934653655664
35241.852680526349
36265.413187480134
37270.218921631927
38279.873134015257
39282.577280868231
40401.58162712244
41449.644667330783
42462.134926752361
43537.208741112394
44743.11661299062
45755.905159584658
46805.676549061245
47827.060977594431
48937.063339539848
491192.85921242141
501216.63528345156
511573.50571062893
521773.79679191880
532308.31009356899
542929.56777876145
554103.92258652658
565717.74468143114
576586.50171807535
5811217.4750954007
5928917.0981165906
60134628.812085719
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Nov/14/t1226625233ax7v76n8w6ptzgp/1e7a11226625075.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Nov/14/t1226625233ax7v76n8w6ptzgp/1e7a11226625075.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Nov/14/t1226625233ax7v76n8w6ptzgp/2wgkb1226625075.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Nov/14/t1226625233ax7v76n8w6ptzgp/2wgkb1226625075.ps (open in new window)


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




Copyright

Creative Commons License

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License.

Software written by Ed van Stee & Patrick Wessa


Disclaimer

Information provided on this web site is provided "AS IS" without warranty of any kind, either express or implied, including, without limitation, warranties of merchantability, fitness for a particular purpose, and noninfringement. We use reasonable efforts to include accurate and timely information and periodically update the information, and software without notice. However, we make no warranties or representations as to the accuracy or completeness of such information (or software), and we assume no liability or responsibility for errors or omissions in the content of this web site, or any software bugs in online applications. Your use of this web site is AT YOUR OWN RISK. Under no circumstances and under no legal theory shall we be liable to you or any other person for any direct, indirect, special, incidental, exemplary, or consequential damages arising from your access to, or use of, this web site.


Privacy Policy

We may request personal information to be submitted to our servers in order to be able to:

  • personalize online software applications according to your needs
  • enforce strict security rules with respect to the data that you upload (e.g. statistical data)
  • manage user sessions of online applications
  • alert you about important changes or upgrades in resources or applications

We NEVER allow other companies to directly offer registered users information about their products and services. Banner references and hyperlinks of third parties NEVER contain any personal data of the visitor.

We do NOT sell, nor transmit by any means, personal information, nor statistical data series uploaded by you to third parties.

We carefully protect your data from loss, misuse, alteration, and destruction. However, at any time, and under any circumstance you are solely responsible for managing your passwords, and keeping them secret.

We store a unique ANONYMOUS USER ID in the form of a small 'Cookie' on your computer. This allows us to track your progress when using this website which is necessary to create state-dependent features. The cookie is used for NO OTHER PURPOSE. At any time you may opt to disallow cookies from this website - this will not affect other features of this website.

We examine cookies that are used by third-parties (banner and online ads) very closely: abuse from third-parties automatically results in termination of the advertising contract without refund. We have very good reason to believe that the cookies that are produced by third parties (banner ads) do NOT cause any privacy or security risk.

FreeStatistics.org is safe. There is no need to download any software to use the applications and services contained in this website. Hence, your system's security is not compromised by their use, and your personal data - other than data you submit in the account application form, and the user-agent information that is transmitted by your browser - is never transmitted to our servers.

As a general rule, we do not log on-line behavior of individuals (other than normal logging of webserver 'hits'). However, in cases of abuse, hacking, unauthorized access, Denial of Service attacks, illegal copying, hotlinking, non-compliance with international webstandards (such as robots.txt), or any other harmful behavior, our system engineers are empowered to log, track, identify, publish, and ban misbehaving individuals - even if this leads to ban entire blocks of IP addresses, or disclosing user's identity.

FreeStatistics.org is powered by