# NOT RUN {
# Redundant testing excluded
# Color a distance matrix
pdb <- read.pdb( system.file("examples/1hel.pdb", package="bio3d") )
d <- dm(pdb,"calpha")
plot(d, color.palette=bwr.colors)
plot(d,
resnum.1 = pdb$atom[pdb$calpha,"resno"],
color.palette = mono.colors,
xlab="Residue Number", ylab="Residue Number")
# }
Run the code above in your browser using DataLab