Home » date » 2008 » Nov » 12 »

opdracht8.3 various EDA topics T2

*Unverified author*
R Software Module: rwasp_hierarchicalclustering.wasp (opens new window with default values)
Title produced by software: Hierarchical Clustering
Date of computation: Wed, 12 Nov 2008 05:01:04 -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/12/t1226491364xjv66lw0cibo67y.htm/, Retrieved Wed, 12 Nov 2008 12:02:51 +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/12/t1226491364xjv66lw0cibo67y.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:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
31/08/2001 90.7 78.4 97.8 30/09/2001 94.3 114.6 107.4 31/10/2001 104.6 113.3 117.5 30/11/2001 111.1 117 105.6 31/12/2001 110.8 99.6 97.4 31/01/2002 107.2 99.4 99.5 28/02/2002 99 101.9 98 31/03/2002 99 115.2 104.3 30/04/2002 91 108.5 100.6 31/05/2002 96.2 113.8 101.1 30/06/2002 96.9 121 103.9 31/07/2002 96.2 92.2 96.9 31/08/2002 100.1 90.2 95.5 30/09/2002 99 101.5 108.4 31/10/2002 115.4 126.6 117 30/11/2002 106.9 93.9 103.8 31/12/2002 107.1 89.8 100.8 31/01/2003 99.3 93.4 110.6 28/02/2003 99.2 101.5 104 31/03/2003 108.3 110.4 112.6 30/04/2003 105.6 105.9 107.3 31/05/2003 99.5 108.4 98.9 30/06/2003 107.4 113.9 109.8 31/07/2003 93.1 86.1 104.9 31/08/2003 88.1 69.4 102.2 30/09/2003 110.7 101.2 123.9 31/10/2003 113.1 100.5 124.9 30/11/2003 99.6 98 112.7 31/12/2003 93.6 106.6 121.9 31/01/2004 98.6 90.1 100.6 29/02/2004 99.6 96.9 104.3 31/03/2004 114.3 125.9 120.4 30/04/2004 107.8 112 107.5 31/05/2004 101.2 100 102.9 30/06/2004 112.5 123.9 125.6 31/07/2004 100.5 79.8 107.5 31/08/2004 93.9 83.4 108.8 30/09/2004 116.2 113.6 128.4 31/10/2004 112 112.9 121.1 30/11/2004 106.4 104 119.5 31/12/2004 95.7 109.9 128.7 31/01/2005 96 99 108.7 28/02/2005 95.8 106.3 105.5 31/03/2005 103 128.9 119.8 30/04/2005 102.2 111.1 111.3 31/05/2005 98.4 102.9 110.6 30/06/2005 111.4 130 120.1 31/07/2005 86.6 87 97.5 31/08/2005 91.3 87.5 107.7 30/09/2005 107.9 117.6 127.3 31/10/2005 101.8 103.4 117.2 30/11/2005 104.4 110.8 119.8 31/12/2005 93.4 112.6 116.2 31/01/2006 100.1 102.5 111 28/02/2006 98.5 112.4 112.4 31/03/2006 112.9 135.6 130.6 30/04/2006 101.4 105.1 109.1 31/05/2006 107.1 127.7 118.8 30/06/2006 110.8 137 123.9 31/07/2006 90.3 91 101.6 31/08/2006 95.5 90.5 112.8 30/09/2006 111.4 122.4 128 31/10/2006 113 123.3 129.6 30/11/2006 107.5 124.3 125.8 31/12/2006 95.9 120 119.5 31/01/2007 106.3 118.1 115.7 28/02/2007 105.2 119 113.6 31/03/2007 117.2 142.7 129.7 30/04/2007 106.9 123.6 112 31/05/2007 108.2 129.6 116.8 30/06/2007 113 151.6 127 31/07/2007 97.2 110.4 112.1 31/08/2007 99.9 99.2 114.2 30/09/2007 108.1 130.5 121.1 31/10/2007 118.1 136.2 131.6 30/11/2007 109.1 129.7 125 31/12/2007 93.3 128 120.4 31/01/2008 112.1 121.6 117.7 29/02/2008 111.8 135.8 117.5 31/03/2008 112.5 143.8 120.6 30/04/2008 116.3 147.5 127.5 31/05/2008 110.3 136.2 112.3 30/06/2008 117.1 156.6 124.5 31/07/2008 103.4 123.3 115.2 31/08/2008 96.2 100.4 105.4
 
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'Gwilym Jenkins' @ 72.249.127.135


Summary of Dendrogram
LabelHeight
11.79168992949735
21.94422229282815
32.40416779297055
42.43515913515519
52.535758613688
62.6925824060872
72.73130283500468
82.96984848098390
93.00998364666029
103.18333374323805
113.40293991068268
123.4842503775909
133.59029011408157
143.61109402777765
153.64143007854879
163.92960181467499
174.17255334960276
184.47660633321512
194.60108683643955
204.63896540190065
214.68028947805942
224.75184199499103
235.08428952806283
245.17783738971267
255.1893363725513
265.26214627972185
275.26741047141656
285.28677610146172
295.32183410642805
305.3218731788598
315.32916625821537
325.42678699467211
335.62339142119244
345.96437605691548
356.49628545269212
366.58103450771045
376.70171900091033
386.73058613694655
396.81909084849301
407.20347211236167
417.7150883321705
427.74725803374203
437.8447434629821
448.0466120504973
458.3441425108843
468.45103544032073
478.45990543682378
488.56939387663422
498.91282492113824
509.1933168150248
519.26949089206632
529.3454961405434
539.50012946557253
549.54430973502844
5510.3498792263487
5611.3471868144622
5711.5547179514600
5812.5635840995395
5912.6110661251402
6013.3754496768022
6113.9241416568346
6218.0991828471803
6318.4946525129736
6420.2376117206061
6520.2986371456707
6622.675989995303
6725.0714082712866
6825.2296929088941
6925.7177979187179
7026.2352251129129
7126.6461161159687
7229.3347350238217
7333.2167887544979
7436.5136666427946
7538.5075635155270
7639.8461725290528
7746.8151583875458
7849.0341339285924
7961.0424208159177
8080.6472601342296
81121.547205342734
82152.920783915421
83255.062315074147
84680.616476883125
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Nov/12/t1226491364xjv66lw0cibo67y/1crso1226491262.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Nov/12/t1226491364xjv66lw0cibo67y/1crso1226491262.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Nov/12/t1226491364xjv66lw0cibo67y/28gsx1226491262.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Nov/12/t1226491364xjv66lw0cibo67y/28gsx1226491262.ps (open in new window)


 
Parameters (Session):
par1 = ward ; par2 = ALL ; par3 = TRUE ; par4 = FALSE ;
 
Parameters (R input):
par1 = ward ; par2 = ALL ; par3 = TRUE ; 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