read.cif 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)
-
-
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.
-
- changed status to resolved
- Log in to comment
Thanks for the very helpful report - we will dig into this and post back soon.