elety "NA" produces NA in pdb$atom$elety

[Lars Skjærven]

pdb = read.pdb("1c75")
any(is.na(pdb$atom$elety))
[1] TRUE

[Lars Skjærven]

occurs with atom name "NA" in the PDB file.

[Xinqiu Yao]

Fixed by using empty string "" instead of "NA" as NA values. Temporarily checked with 1c75 example and it was correct:

 > pdb = read.pdb("1C75")
  Note: Accessing on-line PDB file
  HEADER    ELECTRON TRANSPORT                      09-FEB-00   1C75               
   PDB has ALT records, taking A only, rm.alt=TRUE
> any(is.na(pdb$atom$elety))
[1] FALSE

[Lars Skjærven]

looks good. thanks !