Copyright © 1997 by the Regents of the University of California.
plotDihedral: plot dihedral angles
The plotDihedral command reads its input
and generates a PostScript file containing time and
distribution plots of selected dihedral angles.
Command Line Flags
plotDihedral
[ -aeiv ]
[ -b degrees ]
[ -c columns ]
[ -o output-file ]
[ -r rows ]
[ -s dihedral-number ]
[ -t type ]
[ -A angle ]
[ -D color ]
[ -I format ]
[ -O format ]
[ -T color ]
input-files
-
-a
-
Plot all dihedrals.
-
-b degrees
-
Normally, plotDihedral plots a single line for each
dihedral angle sample.
Specifying the -b flag makes it
plot distribution as degrees size bins
rather than individual lines.
The default bin size is zero, which
is the same as using individual lines.
-
-c columns
-
Display columns columns of plots per page.
Note that columns is relative to the portrait page,
regardless of whether the plots are rotated.
The default number of columns is two.
-
-e
-
Generate Encapsulated PostScript. It is the responsibility
of the user to make sure that only one page of output is
generated.
-
-i
-
Interactive mode. Normally, when dihedrals have been selected
using either the -a or -s flag, plotDihedral
will plot the selected dihedrals without displaying the dihedral
selection panel. Specifying the -i flag will make it
always show the panel. Any pre-selected dihedral will be
automatically selected in the panel as well.
-
-o output-file
-
Send PostScript output to output-file. Without this flag,
plotDihedral sends PostScript to standard output.
-
-r rows
-
Display rows rows of plots per page.
Note that rows is relative to the portrait page,
regardless of whether the plots are rotated.
The default number of rows is four.
-
-s dihedral-number
-
If the -a flag is not specified,
select the dihedral-numberth dihedral for plotting;
otherwise,
deselect the dihedral-numberth dihedral
for plotting.
dihedral-number is the zero-based index of dihedrals
found in the CARNAL output file after any type
filtering specified by the -t flag.
The -s flag is only for brave souls; others should
use the selection panel in interactive mode (see the
-i flag).
-
-t type
-
Use only dihedrals of type type from the CARNAL output file.
Currently, the supported values for type are phi,
psi, omega, and any.
Only the type specified by the last -t flag is used.
The default type is any.
-
-v
-
Print informational messages, including
CARNAL output not preceded by `>',
number of atoms, and number of dihedrals.
-
-A angle
-
Specify the angle (in degrees) at the bottom of a Cartesian plot.
The default value is -180, which places zero degrees in the middle
of the plot.
-
-D color
-
Specify the color of the distribution plot.
Only the data portion is shown in color,
i.e., the border will still be black.
color is specified in RGB and must be
in the format (r, g, b);
the parentheses and commas are required elements and
r, g and b must be between 0 and 1 inclusive.
This option will probably be replaced by a more general
color specification scheme.
-
-I format
-
Select input format. Currently, the only supported format
(and hence the default)
is carnal, for CARNAL dihedral analysis output. The two
input-files required for carnal format are the
CARNAL output file and table file.
-
-O format
-
Select output format. Currently, the supported formats
are radial (default) and cartesian. See
below for description of generated plots.
-
-T color
-
Specify the color of the time plot.
Only the data portion is shown in color,
i.e., the border will still be black.
color is specified in RGB and must be
in the format (r, g, b);
the parentheses and commas are required elements and
r, g and b must be between 0 and 1 inclusive.
This option will probably be replaced by a more general
color specification scheme.
-
input-files
-
List of input files to process.
Required and permissible input files depend on the input format.
See -I.
Description
The plotDihedral command proceeds in the following steps:
-
Get Dihedral List
-
Read the output and table files from CARNAL to get
a list of dihedral angle specifications and their values during a trajectory.
-
Filter Dihedral List
-
If the -t flag is specified, the dihedral list is filtered
to select only those of the correct type.
-
Select Dihedrals
-
If no dihedrals are selected by the -a or -s,
or if the -i flag is specified, plotDihedral
displays a panel from which the user can select dihedrals
of interest, as well as change the page layout parameters
such as the number of columns and rows of plots per page,
and distribution plot bin size.
-
Plot Selected Dihedrals
-
Plots of selected dihedral angles are generated in PostScript.
If the -e flag is specified, the output is suitable for
incorporation into other programs such as FrameMaker and Diagram.
The plots generated by plotDihedral have three components:
a title on top,
a time plot on the right,
and a distribution plot on the left.
The title is the dihedral specification as found in CARNAL output
and looks like:
name: atom1 - atom2 = atom3 - atom4
(average / sigma)
where each atom specification looks like:
residue-number:atom-number residue-type atom-name
name is the user-selected name when a dihedral was specified using the
TORSION statement. average and sigma are the average
and standard deviation of the dihedral. These three fields and
residue-type may be absent if no information is available.
Both the time and distribution plots are either radial or Cartesian plots.
The radial time plot proceeds from the center to the perimeter of
the circle, drawing lines between consecutive angle samples.
The radial distribution plot have one of two appearances.
If the distribution plot bin size is zero, a line is drawn
from the center to the perimeter for each angle sample.
If the bin size is non-zero, the circle is divided into
pie slices that subtend the bin size (e.g., if the bin
size is five degrees, then the circle is divided into seventy-two
pie slices). For each pie slice, the number of samples falling
in the slice determines the gray level of the pie slice. A slice
with zero samples is white; a slice with more than half the
total sample is black; all other slices will be a shade of gray
proportional to the number of samples they contain. Any slice
that contains at least one sample is guaranteed to be a visible
shade of gray, i.e., it will not be so light gray as to be
indistinguishable from white.
The Cartesian plots are rectangular instead of circular.
For the Cartesian time plot, the y-axis represents the dihedral
angle, with zero in the middle and 180 and -180 at the top and
bottom respectively (this may be changed using the -A option);
the x-axis represents time and proceeds from left to right.
For the Cartesian distribution plot, the y-axis again represents
the dihedral angle. If the distribution plot bin size is zero,
horizontal lines spanning the plot are drawn for each angle sample.
If the bin size is non-zero, the y-axis is divided into slices that
correspond to the bin size. For each bin with one or more samples,
a rectangle whose width is proportional to the number of samples in
the bin is drawn. A rectangle that horizontally spans the plot
represents a bin that contains 50% or more of the samples.
See Also
CARNAL,
written by Wilson S. Ross, is described in the
AMBER Version 4.1 documentation on pages 234-253.
plotDihedral extension mechanisms.
Conrad Huang, UCSF Computer Graphics Laboratory