[Chimera-users] representative fits in chimera fitmap search

Elaine Meng meng at cgl.ucsf.edu
Fri Mar 18 12:48:26 PDT 2016

Hi Jan!
I’ll have to leave the specific answer to Tom G, but if it is indeed just keeping the first solution in a cluster:

I believe that the reasoning is that when the fits are “non-unique" within small shifts or rotations, they would effectively be equivalent to one another and likely to converge on the same local minimum.  However, you could try reducing the uniqueness criteria (see clusterAngle and clusterShift options)...


...and/or change other fitmap parameters for more thorough local optimization (number of steps, initial step size, convergence criterion) to see if they have any substantive effect on the results.
Elaine C. Meng, Ph.D.                       
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab
Department of Pharmaceutical Chemistry
University of California, San Francisco

On Mar 18, 2016, at 10:19 AM, Jan Kosinski <kosinski at embl.de> wrote:

> Hi,
> In the fitmap command global search, how the unique fit is selected from the cluster of similar fits?
> Ideally, from each cluster, I would like to obtain the fit the gave the best cross-correlation. But, I think, based on the fit_search code in FitMap/search.py, what is happening is that the first unique fit ever encountered is added to the fit list:
>        close = b.close_transforms(ptf)
>        if len(close) == 0:
>            transforms = [M.multiply_matrices(ptf, mtv) for mtv in mtv_list]
>            stats['hits'] = 1
>            f = Fit(models, transforms, volume, stats)
>            f.ptf = ptf
>            flist.append(f)
>            b.add_transform(ptf)
>            fo[id(ptf)] = f
>        else:
>            s = fo[id(close[0])].stats
>            s['hits'] += 1
> and any subsequent fit that is close to the first, would be discarded, even if it gives better cross-correlation.
> Do I understand correctly that the unique fit is not necessarily the best scoring fit from the cluster?
> Thanks in advance for clarification,
> Jan

More information about the Chimera-users mailing list