Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*Unverified author*

R Software Module

rwasp_hierarchicalclustering.wasp

Title produced by software

Hierarchical Clustering

Date of computation

Thu, 28 Feb 2013 06:27:11 -0500

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2013/Feb/28/t1362050915qcx5dni9904yxi6.htm/, Retrieved Wed, 01 May 2024 23:06:55 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=207184, Retrieved Wed, 01 May 2024 23:06:55 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

139

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-       [Hierarchical Clustering] [] [2013-02-28 11:27:11] [d41d8cd98f00b204e9800998ecf8427e] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

Transcript Cluster ID	FC(Y vs. C)	FC(R vs. C)
10608222	1.06733	10.4944
10417371	1.23429	9.51651
10608377	1.13728	7.20195
10608608	1.1611	6.95525
10608293	1.06539	6.89676
10608368	1.07486	6.52938
10602221	1.2787	6.42874
10608373	1.12006	6.22858
10608606	1.10482	6.13692
10427035	2.31741	6.10738
10608327	1.08794	5.98788
10608567	1.03463	5.37065
10608484	1.15616	5.34158
10436773	1.152	5.28816
10550208	1.04269	5.08896
10608480	1.21904	5.07814
10563937	1.24442	5.07291
10560744	-1.01175	4.99595
10560780	1.04843	4.989
10542221	1.06874	4.98554
10608365	1.12771	4.81329
10560752	1.00131	4.69144
10560730	1.05448	4.65383
10550193	1.12201	4.59927
10608506	1.16904	4.52743
10560797	1.05355	4.30597
10560728	-1.00511	4.2809
10608209	1.05347	4.21098
10560754	1.01775	4.18404
10550760	1.0558	4.18311
10550786	1.00374	4.18221
10550770	1.06302	4.13593
10447617	1.47076	4.13098
10560787	-1.01573	3.95536
10560785	1.01978	3.93917
10410475	1.23699	3.89159
10550189	1.01413	3.85637
10560732	1.02649	3.85156
10560789	1.05526	3.83015
10598083	1.50704	3.82926
10608394	1.04898	3.81541
10550768	1.04519	3.79488
10560742	1.05534	3.67149
10550782	1.04072	3.64983
10608342	1.01953	3.63535
10608282	1.15024	3.63399
10608308	1.0638	3.53606
10550765	-1.00985	3.43764
10608424	1.03751	3.2463
10608350	1.05683	3.19688
10564161	1.2697	3.17248
10608260	-1.04651	3.14887
10608371	1.06556	3.12941
10608613	-1.02234	3.02691
10608361	1.00989	2.996
10564183	1.39435	2.95423
10373651	1.16809	2.88128
10419122	1.06206	2.78474
10608557	1.04444	2.77212
10559928	-1.04919	2.76363
10608385	1.0821	2.76124
10403011	1.28639	2.75861
10519811	1.23328	2.75457
10608630	1.1636	2.73136
10476106	3.02386	2.70675
10608247	1.04045	2.67949
10564177	1.13581	2.6682
10608625	1.17744	2.6602
10419119	-1.02395	2.65758
10608628	1.05529	2.65529
10608407	-1.07634	2.64345
10608420	1.00118	2.64342
10414241	1.03838	2.6001
10397346	1.68135	2.57466
10482772	1.65121	2.49602
10356271	1.71695	2.48235
10545247	1.13327	2.4617
10608410	1.01284	2.44346
10403988	-1.03938	2.44335
10608521	1.0868	2.43989
10369221	1.2831	2.4373
10608460	1.03499	2.41123
10608488	1.0399	2.40984
10361091	1.54953	2.40405
10608273	1.14343	2.40081
10419116	-1.00649	2.38695
10608482	-1.02915	2.38046
10608464	1.07306	2.3774
10454782	1.57835	2.3467
10608523	1.04175	2.34023
10399428	-1.20318	2.28624
10414706	1.1119	2.28572
10477475	1.69519	2.26072
10417452	1.0091	2.25425
10608302	1.0552	2.22408
10356291	1.13846	2.20589
10369225	1.21737	2.20259
10375574	1.28263	2.19915
10369223	1.23348	2.1982
10351037	2.79559	2.1981
10375608	1.43577	2.18805
10570634	1.19935	2.18613
10544523	1.11729	2.17428
10522581	1.09762	2.16598
10536155	1.09762	2.16598
10572527	1.18764	2.16516
10347915	1.00373	2.15741
10408130	-1.09195	2.15606
10536048	1.06208	2.13553
10582860	1.13045	2.10787
10582877	1.13045	2.10787
10412517	1.03838	2.10741
10514329	-1.02814	2.10408
10608313	1.02634	2.10311
10419125	1.03993	2.10084
10404008	-1.06166	2.08805
10536147	1.07207	2.07893
10577388	1.20778	2.07715
10418169	1.35383	2.05858
10522569	1.07473	2.05842
10522577	1.12771	2.055
10536159	1.12771	2.055
10449741	1.72868	2.05217
10608317	1.00181	2.05015
10560801	-1.015	2.04454
10608517	1.00287	2.04347
10563114	1.89506	2.03573
10356274	1.08035	2.02431
10579060	1.42277	2.0197
10414234	1.02299	2.01817
10435501	1.37385	2.01217
10414228	1.10496	2.00709
10551226	1.68876	2.00487
10551025	1.46934	-2.0122
10545210	5.6231	-2.03272
10419096	-4.55133	-2.07281
10538903	1.21031	-2.08504
10403028	2.34398	-2.10831
10403054	1.12304	-2.12912
10351873	-1.09944	-2.2125
10466172	1.51757	-2.22388
10403079	1.26492	-2.27145
10545175	1.80601	-2.30621
10403063	1.88736	-2.31881
10403060	1.7118	-2.34578
10545233	2.8099	-2.4483
10403018	1.15321	-2.4952
10403015	1.05557	-2.53764
10538882	1.77366	-2.87171
10403048	1.38896	-2.88242
10531126	1.34735	-3.44677
10403031	1.52322	-3.5702
10403043	-1.00889	-3.70188
10545196	2.15942	-4.58852
10545184	2.15529	-4.68987
10545194	1.37134	-4.92176
10538880	-1.10245	-5.16504
10545198	1.61942	-5.3434
10545187	1.76138	-5.59846
10403034	1.09055	-5.68889
10403069	1.42707	-9.07851
10545215	-2.8595	-10.4461

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	0 seconds
R Server	'Gwilym Jenkins' @ jenkins.wessa.net
R Engine error message	Error in array(list(Transcript, Cluster, ID, FC(Y, vs., C), FC(R, vs., : object 'Transcript' not found Execution halted

\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 & 0 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ jenkins.wessa.net \tabularnewline
R Engine error message & Error in array(list(Transcript, Cluster, ID, FC(Y, vs., C), FC(R, vs.,  : 
  object 'Transcript' not found
Execution halted
 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=207184&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]0 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ jenkins.wessa.net[/C][/ROW]
[ROW][C]R Engine error message[/C][C]Error in array(list(Transcript, Cluster, ID, FC(Y, vs., C), FC(R, vs.,  : 
  object 'Transcript' not found
Execution halted
[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=207184&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=207184&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 Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	0 seconds
R Server	'Gwilym Jenkins' @ jenkins.wessa.net
R Engine error message	Error in array(list(Transcript, Cluster, ID, FC(Y, vs., C), FC(R, vs., : object 'Transcript' not found Execution halted

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

par4 <- 'FALSE'
par3 <- 'FALSE'
par2 <- 'ALL'
par1 <- 'ward'
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')
}

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code