if (FALSE) {
## 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