By downloading this code you agree to the
Source Code Use License (PDF). |
{ version = 3.48; (* of dnaplot.p 2005 May 23}
(* begin module describe.dnaplot *)
(*
name
dnaplot: plot Ri values for one or more large DNA sequences
synopsis
dnaplot(dnain: in, dnaout: out, dnaplotp: in,
positions: in, dnasymbols: in, output: out)
files
dnain: A data input file created by scan. It contains header lines that
begin with asterisks ('*') that are copied to dnaout. Remaining lines
are the data in columns, ending with end of file.
dnaout: output in PostScript format
dnaplotp: Parameter file for dnaplot, which is configured as followed:
frompos topos These two integers are the positions on the
sequence which the graph will represent. If,
instead, the first character on the line is 'r',
then these numbers are read from the positions
file for each sequence.
sCol pCol cCol vCol columns numbers for reading from the dnain file:
sCol = piece number
pCol = column number of piece name
cCol = piece coordinate
vCol = value (Ri evaluation)
numperpg number of graphs per page
numperln number of base pairs per line
bitlower bitupper lower and upper bounds of bits to display
orix oriy x, y origin of plot (in cm)
xaxlength yaxlength length of the x and y axes in cm
showaxis t=true means show coordinate axis to dnaout
xinterval yinterval size of intervals on axes to plot
xsubint ysubint x and y sub intervals to mark
xwidth ywidth width of numbers in characters
xdecimal ydecimal number of decimal places
xticlength xticdx xticdy length of tic mark and shift of number (cm)
yticlength yticdx yticdy length of tic mark and shift of number (cm)
sequencelabel t=true means print sequence number on graphs
xaxislabel the label for the x axis
yaxislabel the label for the y axis
plottype How to draw the plot:
z = lines from zero to value
b = lines from bottom of graph to value
dodash Whether to put vertical dashed lines around
segments of continuous sequence. This is
important for distinguishing between the absence
of sequence and low Ri values, but this often gets
in the way, so it can be turned off:
d = do dashes
n = no dashes
thelinewidth The width of the plot lines in bases.
positions: If the first character on the first line of the parameter file
is an 'r' then this file will be read to determine the positions to graph
for each sequence. The file consists of pairs of integers, one pair per
line, representing the first (frompos) and last (topos) coordinates to be
plotted.
dnasymbols: If the file is not empty, then it contains information on how
and where to plot special symbols to make marks on the graph. Each line
has 5 values:
symbol (first character): the type of symbol to draw:
c = circle
s = square
t = triangle
symboltype (second character): the way to draw the symbol:
s = stroked as a solid line
f = filled
d = dotted line
symbolplacement (third character): where to put the symbol on the graph:
a = use absolute location (given by symbolbits) on graph
r = use relative location (given by symbolbits) from current Ri value
symbolbits (real): the distance in bits
symbolsize (real): diameter of a circle or side of square and triangle
relative to the spacing between graph lines. A value of 1 fits
between the lines.
piece number (integer): The number of the fragment to mark,
as given in the dnain file
piece coordinate (integer): the coordinate on the piece to mark
as given in the dnain file
********************************************************************
* The symbols MUST be in increasing order of position in the plot! *
********************************************************************
Lines that are empty or begin with "*" are ignored.
output: messages to the user
description
dnaplot creates a PostScript graph of information content (or
other values) versus position on a DNA sequence.
examples
dnasymbols
The line:
csa 15 0.5 1 100
means place a circle, stroked, at absolute 15 bits, 0.5 size,
for sequence 1 at coordinate 100.
documentation
see also
scan.p, xyplo.p, dbbk.p
author
Stacy L. Bartram, modified by Tom Schneider
bugs
The program cannot handle negatively numbered base systems because the axis
cannot give decreasing numbers (yet or ever). These simply come out as
blank (or the program will halt, depending on the version).
technical notes
*)
(* end module describe.dnaplot *)
{This manual page was created by makman 1.45}
{created by htmlink 1.62}