read.crd(file, verbose = TRUE)
nrow(atom)
with TRUE values indicating a C-alpha read.pdb
for more details. For a description of CHARMM CARD (CRD) format see:
write.crd
, read.pdb
,
atom.select
, write.pdb
,
read.dcd
, read.fasta.pdb
,
read.fasta
pdb <- read.pdb("1bg2")
crdfile <- tempfile()
write.crd(pdb, file=crdfile)
crd <- read.crd(crdfile)
ca.inds <- which(crd$calpha)
crd$atom[ca.inds[1:20],c("x","y","z")]
# write.pdb(crd, file=tempfile())
Run the code above in your browser using DataLab