read.cif issue

Issue #400 resolved
Former user created an issue

Hi! I've been working with bio3d (version 2.2-2) for some months and I'm glad to see the new "read.cif" function in the new version 2.3.0 (btw, congrats for the package). I've started working with this function, but at the moment I've found one PDB entry which I could not read: 4R3I.

read.cif("4R3I.cif") Error in matrix(unlist(value, recursive = FALSE, use.names = FALSE), nrow = nr, : length of 'dimnames' [2] not equal to array extent In addition: Warning message: In read.cif("4R3I.cif") : helix/sheet records will not be parsed in this version of the code

I tried to access it as a local file in two computers using R version 3.3.1 (Ubuntu 16.04.1 LTS and Ubuntu 14.04.02. LTS in each computer). In all cases I received the same output. Many thanks in advance!

Comments (3)

  1. Lars Skjærven

    Thanks for reporting. This is caused by 0's in the 'model' column for this particular PDB. Fixed by counting model id's. See commit 44e1fae.

    You can update using package devtools. See instructions on the bio3d bitbucket front page.

  2. Log in to comment