# NOT RUN {
## Read a PRMTOP file
prmtop <- read.prmtop(system.file("examples/crambin.prmtop", package="bio3d"))
print(prmtop)
## Read a Amber CRD file
crds <- read.crd(system.file("examples/crambin.inpcrd", package="bio3d"))
## Atom selection
ca.inds <- atom.select(prmtop, "calpha")
## Convert to PDB format
pdb <- as.pdb(prmtop, crds, inds=ca.inds)
# }
Run the code above in your browser using DataLab