That works well, thanks Eric. Is there any way to reassign the ramaProbs to the bfactor column on the command line, without using the attribute calculator? That way I can set the whole sequence of commands as an alias.

Cheers,
Oliver.
On Aug 26, 2014, at 1:45 PM, Eric Pettersen <pett@cgl.ucsf.edu> wrote:

On Aug 26, 2014, at 10:27 AM, Elaine Meng <meng@cgl.ucsf.edu> wrote:

Hi Oliver,
The more direct way I tried used Attribute Calculator (in menu under Tools.. Structure Analysis) to calculate attribute "bfactor" for "atoms" using the formula:

residue.ramaProb

(after assignment using Ramachandran).  It currently generates an error, which seems buglike.  I could create a new atom attribute with that formula, but not overwrite "bfactor". We'll look into that .

The problem here is that bfactors must be numeric, but some residues have a ramaProb of "None" (e.g. they are terminal residues or non-peptides and lack either a phi or a psi).  You can assign such residues a ramaProb of -1.0 with this command:

setattr r ramaProb -1.0 :/^ramaProb

You could then use the Attribute Calculator to assign the ramaProbs to the bfactors.

--Eric

                        Eric Pettersen
                        UCSF Computer Graphics Lab
                        http://www.cgl.ucsf.edu