[gmx-users] Caution: Bug in g_cluster/gmx_cluster

[Roman Affentranger]

Dear all

I found a bug in g_cluster/gmx_cluster (gmx-3.2.1, but also present in the 
3.1.4 version):

(All line numbers refer to gmx_cluster.c of gmx-3.2.1).
The cluster centers, defined as the structure with the lowest RMSD to all 
cluster members, are not correctly calculated by g_cluster (of gmx-3.2.1 
and 3.1.4). The problem is that the function 'analyze_clusters', which 
determines the cluster centers, is called (line 1376) after the function 
'plot_clusters' (line 1364), where one half of the RMSD matrix is replaced 
by the cluster indices. Thus, not only the cluster centers are incorrectly 
determined, but also the average RMSD values reported in the log-file 
(average RMSD of all pairs of structures in cluster as well as average RMSD 
of all cluster members to cluster center).

As a solution to the problem one can simply shift the call of 
'plot_clusters' right after the call of 'analyze clusters' in gmx_cluster.c.

Another solution would be to replace the following line of the function 
'analyze_clusters' (line 846)

	r += rmsd[structure[i1]][structure[i]];

with

	if (i<i1)
	   r += rmsd[structure[i]][structure[i1]];
	else
	   r += rmsd[structure[i1]][structure[i]];


Apart from that, I think that the clustering method denoted "linkage" does 
not represent a full linkage, but rather a single linkage clustering.

Greetings

Roman

-----------------------------------------------------------------
Roman Affentranger
Swiss Federal Institute of Technology Zurich
Institute of Biochemistry
Schafmattstrasse 18
ETH Hoenggerberg HPM G 9.3
8093 Zurich
Switzerland
Phone: +41 1 632 31 39