Home » date » 2007 » Nov » 05 » attachments

Various EDA Q2

R Software Module: rwasp_hierarchicalclustering.wasp (opens new window with default values)
Title produced by software: Hierarchical Clustering
Date of computation: Mon, 05 Nov 2007 09:10:19 -0700
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2007/Nov/05/eblhjavttfxck1u1194278904.htm/, Retrieved Mon, 05 Nov 2007 17:08:27 +0100
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
100,70 101,30 97,90 97,60 96,50 96,40 96,60 97,00 96,60 96,40 95,50 94,70 91,80 89,30 89,30 85,90 87,00 83,30 85,90 81,50 88,00 85,00 87,90 84,80 89,20 87,50 90,90 89,00 91,60 90,00 90,20 89,60 89,10 87,40 87,50 84,80 86,30 81,90 86,00 81,10 84,40 79,10 86,10 80,50 91,00 88,50 92,70 90,90 88,00 84,90 84,30 80,00 82,20 76,50 80,80 75,40 79,40 73,50 80,20 74,30 82,20 77,70 82,20 77,90 81,20 76,70 82,10 77,20 88,10 86,00 88,50 86,90 92,10 92,00 98,60 101,70 100,90 104,50 100,60 101,70 101,10 100,60 102,10 100,30 103,60 102,50 102,80 101,00 108,30 108,60 104,00 103,40 106,10 106,40 106,30 106,60 109,00 108,90 111,00 110,50 113,70 114,00 112,70 112,80 110,30 109,60 114,50 116,00 119,30 124,60 121,80 129,00 125,40 131,50 129,70 138,60 129,40 138,10 134,50 146,30 141,20 157,60 141,40 158,40 152,20 176,30 167,70 199,90 173,30 210,40 168,70 202,60 172,60 207,10 169,80 202,00 172,00 203,40 179,40 216,30 174,60 207,30 172,50 203,50 172,60 204,40 176,30 203,70 178,90 205,70 179,60 208,00 179,90 209,30 180,30 208,70 180,90 206,50 177,70 204,50
 
Text written by user:
 
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 compuational 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
10.0999999999999943
20.0999999999999943
30.141421356237312
40.200000000000003
50.210018769324862
60.282842712474613
70.41231056256177
80.412310562561774
90.509901951359277
100.509901951359277
110.583095189484521
120.646704523705155
130.707106781186552
140.72111025509282
150.728010988928055
160.761577310586398
170.772184168686542
180.809287012133409
190.824621125123547
200.905538513813748
210.984885780179618
220.984885780179618
230.98994949366117
241.13054235081598
251.13137084989847
261.14017542509914
271.20558339393895
281.21118527806878
291.21574746538694
301.25299640861417
311.25299640861418
321.30904051501827
331.32822761018399
341.35143368061546
351.56204993518133
361.61245154965970
371.86224923923909
382.00997512422418
392.01872442950477
402.03159715193742
412.15406592285381
422.30438975736918
432.34988124600521
442.68446690830628
453.03005445557119
463.36303527026687
473.48172419096934
483.48907312167021
493.81999017291256
503.83678102332728
514.10466231591478
524.13081681110551
534.38292140016223
544.57132749413426
555.19129076390141
565.32334293741009
577.66974696830948
588.052632723089
598.38949465318321
609.22083661977003
619.5772051114983
629.92924084560596
6310.5939848167218
6412.292429605067
6513.0404975979163
6618.8867727923453
6722.5616760628438
6828.1258809942865
6930.7546475326866
7030.7597876840117
7133.9898174724017
7236.0787240676087
7340.4085200367006
7499.437909506743
75110.073975553826
76117.599074034966
77447.285243173008
78754.270793994751
793016.3342992471
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2007/Nov/05/eblhjavttfxck1u1194278904/1mwhc1194279015.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2007/Nov/05/eblhjavttfxck1u1194278904/1mwhc1194279015.ps (open in new window)


 
Parameters:
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:

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