FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_hierarchicalclustering.wasp
Title produced by softwareHierarchical Clustering
Date of computationWed, 12 Nov 2008 04:08:29 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Nov/12/t1226488137s0r7w258znjwyri.htm/, Retrieved Mon, 20 May 2024 08:12:41 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=24126, Retrieved Mon, 20 May 2024 08:12:41 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact164

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F       [Hierarchical Clustering] [Q1 part 3] [2008-11-12 11:08:29] [b5110a3ab194da7214bdf478e0a05dbd] [Current]
Feedback Forum
2008-11-23 14:40:58 [Chi-Kwong Man] [reply] 
Geeft een error aan, dus waarschijnlijk iets fout gelopen met gegevens in te geven?

Hierarchical Clustering: welke observaties kan men plaatsen in één cluster? Dit is puur exploratief.

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
4.43	109.1	2187	88.3
4.61	113.8	1852	88.6
4.54	97.4	1570	91
4.2	72.5	1851	91.5
4.08	82.7	1954	95.4
3.95	88.9	1828	98.7
4.19	105.9	2251	99.9
4.23	100.8	2277	98.6
3.89	94	2085	100.3
3.92	105	2282	100.2
4.14	58.5	2266	100.4
4.24	87.6	1878	101.4
4.08	113.1	2267	103
4.37	112.5	2069	109.1
4.43	89.6	1746	111.4
4.3	74.5	2299	114.1
4.27	82.7	2360	121.8
4.06	90.1	2214	127.6
3.96	109.4	2825	129.9
4.21	96	2355	128
4.31	89.2	2333	123.5
4.35	109.1	3016	124
4.25	49.1	2155	127.4
4.06	92.9	2172	127.6
4	107.7	2150	128.4
3.87	103.5	2533	131.4
3.71	91.1	2058	135.1
3.63	79.8	2160	134
3.48	71.9	2259	144.5
3.6	82.9	2498	147.3
3.66	90.1	2695	150.9
3.45	100.7	2799	148.7
3.3	90.7	2945	141.4
3.14	108.8	2930	138.9
3.21	44.1	2318	139.8
3.12	93.6	2540	145.6
3.14	107.4	2570	147.9
3.4	96.5	2669	148.5
3.42	93.6	2450	151.1
3.29	76.5	2842	157.5
3.49	76.7	3439	167.5
3.52	84	2677	172.3
3.81	103.3	2979	173.5
4.03	88.5	2257	187.5
3.98	99	2842	205.5
4.1	105.9	2546	195.1
3.96	44.7	2455	204.5
3.83	94	2293	204.5
3.72	107.1	2379	201.7
3.82	104.8	2478	207
3.76	102.5	2054	206.6
3.98	77.7	2272	210.6
4.14	85.2	2351	211.1
4	91.3	2271	215
4.13	106.5	2542	223.9
4.28	92.4	2304	238.2
4.46	97.5	2194	238.9
4.63	107	2722	229.6
4.49	51.1	2395	232.2
4.41	98.6	2146	222.1
4.5	102.2	1894	221.6
4.39	114.3	2548	227.3
4.33	99.4	2087	221
4.45	72.5	2063	213.6
4.17	92.3	2481	243.4
4.13	99.4	2476	253.8
4.33	85.9	2212	265.3
4.47	109.4	2834	268.2
4.63	97.6	2148	268.5
4.9	104.7	2598	266.9

Tables (Output of Computation)





\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 & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
R Framework error message & 
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.
 \tabularnewline
R Engine error message & 
Error in dimnames(data) <- dimnames : 
  length of 'dimnames' [1] not equal to array extent
Calls: array
Execution halted

 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24126&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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[ROW][C]R Framework error message[/C][C]
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.
[/C][/ROW]
[ROW][C]R Engine error message[/C][C]
Error in dimnames(data) <- dimnames : 
  length of 'dimnames' [1] not equal to array extent
Calls: array
Execution halted

[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24126&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=24126&T=0

As an alternative you can also use a QR Code:  
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'Gwilym Jenkins' @ 72.249.127.135
R Framework error message
The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.
R Engine error message
Error in dimnames(data) <- dimnames : 
  length of 'dimnames' [1] not equal to array extent
Calls: array
Execution halted


Figures (Output of Computation)

Input Parameters & R Code

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