Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_hierarchicalclustering.wasp
Title produced by softwareHierarchical Clustering
Date of computationThu, 24 Apr 2014 10:15:05 -0400
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Apr/24/t1398348967lk5ye6xgn6wiv86.htm/, Retrieved Fri, 17 May 2024 04:19:57 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=234583, Retrieved Fri, 17 May 2024 04:19:57 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact156
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Hierarchical Clustering] [] [2014-04-24 14:15:05] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
5,310	14,343	8,691	10,284	8,200	9,437	3,879	11,554	19,931	12,369	9,576	12,213	8,261	10,401	10,768	7,645	7,371	13,406	10,873	7,290	14,175	12,699
2,389	8,940	4,778	6,233	4,583	7,042	1,459	6,909	14,185	9,301	6,786	9,139	4,400	5,970	9,110	4,782	9,105	12,643	7,129	4,526	12,189	11,632
3,159	9,641	5,794	7,277	5,230	6,104	1,935	7,143	18,415	8,585	6,796	7,076	4,839	9,811	9,770	3,782	8,488	11,750	9,174	3,610	11,163	11,068
4,938	14,371	8,616	10,247	8,138	9,498	3,665	11,089	18,465	12,255	9,226	11,311	8,753	10,378	11,099	7,839	6,729	14,060	11,111	7,853	14,049	12,308
4,003	12,263	7,485	8,691	7,132	8,935	3,597	8,355	17,904	11,736	8,466	11,824	7,662	8,848	10,998	6,594	10,751	13,258	9,048	6,658	14,651	13,040
2,986	9,582	6,483	6,788	5,517	6,727	2,284	8,692	16,283	9,722	7,385	9,228	5,965	7,185	8,884	6,029	9,282	10,785	8,454	5,334	12,066	12,113
8,344	14,455	11,525	10,510	9,036	10,232	6,133	8,650	19,864	12,986	11,184	12,763	10,156	10,821	13,161	6,924	13,702	15,660	11,295	9,881	13,565	16,890
2,705	9,592	5,373	6,486	4,775	6,228	2,135	8,028	14,535	9,196	7,139	9,543	3,994	6,262	8,378	5,304	9,429	8,410	7,388	4,255	12,357	12,084
3,345	10,130	6,866	7,102	5,162	6,317	2,362	8,879	16,720	9,880	7,185	8,563	6,505	8,909	8,943	5,448	9,161	11,236	8,284	5,162	14,029	12,237
4,117	10,880	3,944	7,344	6,482	7,142	3,105	4,650	15,244	10,255	7,824	10,125	3,184	6,593	10,762	5,991	10,654	11,355	8,918	2,186	8,708	10,010
3,184	10,075	6,756	7,024	4,990	6,046	2,270	7,144	15,748	9,610	7,668	9,080	5,683	6,600	8,903	5,310	9,581	9,985	7,361	6,023	12,444	12,470
5,215	12,488	8,394	8,770	7,575	9,215	4,127	10,340	19,583	12,145	9,638	12,228	7,258	9,104	12,520	7,686	11,474	14,793	10,287	6,362	15,059	13,981
4,764	13,859	7,170	9,819	8,002	7,108	3,491	11,102	18,192	11,813	9,433	11,285	6,865	9,706	11,197	8,189	5,829	13,599	11,316	5,952	12,747	11,933
3,314	10,529	4,574	7,018	5,048	5,909	2,399	7,760	15,273	9,408	7,682	9,593	3,463	6,054	8,921	5,253	10,095	8,781	7,360	3,655	12,713	12,863
4,739	11,396	7,596	7,995	6,664	8,475	3,638	9,350	17,295	10,379	8,627	10,605	6,673	8,056	10,551	6,466	10,813	12,842	9,940	6,306	10,255	14,241
2,860	9,775	4,915	7,145	5,278	6,102	2,182	8,220	14,958	9,870	7,525	9,417	5,637	6,176	8,320	5,405	9,134	11,060	9,190	4,249	11,836	11,202
4,091	11,184	7,813	7,240	5,510	6,762	2,804	9,445	16,264	9,889	8,520	9,793	7,988	7,826	9,975	6,767	9,315	12,191	8,952	7,045	12,818	12,449
3,786	12,169	7,619	8,585	8,095	7,948	3,330	8,173	14,649	12,422	7,909	10,877	6,148	8,294	8,252	6,094	10,371	14,186	11,671	6,255	13,031	14,403
4,259	11,062	7,493	7,332	5,553	6,206	3,145	8,763	16,638	9,982	7,925	9,737	6,283	7,329	9,435	6,377	10,276	12,053	8,453	5,770	12,780	12,909
4,654	13,801	7,765	9,531	7,539	8,048	3,457	11,065	18,626	11,608	8,922	10,451	7,818	10,223	10,370	7,271	5,465	13,329	10,457	6,902	13,545	11,098




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net
R Engine error message
Error in array(list(5.31, 14.343, 8.691, 10.284, 8.2, 9.437, 3.879, 11.554,  : 
  length of 'dimnames' [1] not equal to array extent
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 & 1 seconds \tabularnewline
R Server & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
R Engine error message & 
Error in array(list(5.31, 14.343, 8.691, 10.284, 8.2, 9.437, 3.879, 11.554,  : 
  length of 'dimnames' [1] not equal to array extent
Execution halted
\tabularnewline \hline \end{tabular} %Source: https://freestatistics.org/blog/index.php?pk=234583&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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[ROW][C]R Engine error message[/C][C]
Error in array(list(5.31, 14.343, 8.691, 10.284, 8.2, 9.437, 3.879, 11.554,  : 
  length of 'dimnames' [1] not equal to array extent
Execution halted
[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=234583&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=234583&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 time1 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net
R Engine error message
Error in array(list(5.31, 14.343, 8.691, 10.284, 8.2, 9.437, 3.879, 11.554,  : 
  length of 'dimnames' [1] not equal to array extent
Execution halted



Parameters (Session):
par1 = single ; par2 = ALL ; par3 = FALSE ; par4 = FALSE ;
Parameters (R input):
par1 = single ; par2 = ALL ; par3 = FALSE ; par4 = FALSE ;
R code (references can be found in the software module):
par4 <- 'FALSE'
par3 <- 'FALSE'
par2 <- 'ALL'
par1 <- 'single'
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')
}