Re: [ccp4bb] [phenixbb] Calculate average B-factor?

2014-10-06 Thread Tim Gruene
Dear Jose, the question came up again because I did not receive an answer to my question. The thread discussed benefits and malefits of PDB vs. mmCIF, which was not my question. This time, Nat Echolls gave a very reasonable answer (at least for phenix) on the phenixbb, i.e., that there are no

Re: [ccp4bb] [phenixbb] Calculate average B-factor?

2014-10-06 Thread Dale Tronrud
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 On 10/6/2014 4:20 AM, Tim Gruene wrote: Dear Jose, the question came up again because I did not receive an answer to my question. The thread discussed benefits and malefits of PDB vs. mmCIF, which was not my question. This time, Nat Echolls

Re: [ccp4bb] [phenixbb] Calculate average B-factor?

2014-10-05 Thread Tim Gruene
Hi all, reading this beauty I would like to ask a question to the respective developers: Will the PDB format remain the working format for the users and only upon deposition will it be converted to PDBml for archiving purposes, or are the refinement programs (et al.) going to abandon PDB, too?

Re: [ccp4bb] [phenixbb] Calculate average B-factor?

2014-10-05 Thread Jose Manuel Duarte
Well, if you simply replace that beauty by this one: grep ^ATOM filename.cif | awk '{print $15}' | awk '{s+=$1;} END {print s/NR;}' You will achieve exactly the same result (the b-factors are in the 15th field of the _atom_site section in deposited mmCIF files). I'm not an expert in awk,

Re: [ccp4bb] [phenixbb] Calculate average B-factor?

2014-10-05 Thread Frances C. Bernstein
mmCIF is a very general format with tag-value pairs, and loops so that tags do not need to be repeated endlessly. It was designed so that there is the flexibility of defining new terms easily and presenting the data in any order and with any kind of spacing. I understand that there are 10+

Re: [ccp4bb] [phenixbb] Calculate average B-factor?

2014-10-05 Thread Kay Diederichs
On Sun, 5 Oct 2014 13:13:14 +0200, Jose Manuel Duarte jose.dua...@psi.ch wrote: Well, if you simply replace that beauty by this one: grep ^ATOM filename.cif | awk '{print $15}' | awk '{s+=$1;} END {print s/NR;}' the problem with this is that it will break if there is any B-factor of 100 or

Re: [ccp4bb] [phenixbb] Calculate average B-factor?

2014-10-05 Thread Kay Diederichs
On Sun, 5 Oct 2014 18:54:16 +0100, Kay Diederichs kay.diederi...@uni-konstanz.de wrote: grep ^ATOM filename.cif | awk '{print $15}' | awk '{s+=$1;} END {print s/NR;}' the problem with this is that it will break if there is any B-factor of 100 or higher, because then the blank will be missing

Re: [ccp4bb] [phenixbb] Calculate average B-factor?

2014-10-05 Thread Jose Manuel Duarte
Thanks Frances for the explanation. Indeed mmCIF format is a lot more complicated and grep can be a dangerous tool to use with them. But for most cases it can do the job and thus it maintains some sort of backwards compatibility. I can't agree more that using specialised tools (for either PDB