[metapost] Attn: Taco

Mia Champion champion at broadinstitute.org
Thu Jul 30 19:02:03 CEST 2009


Hi,

I am not at all familiar with the metapost language and have an output 
file from a program that I am using for genome analysis (see below).
In order for the data to be linked into other output files from the 
program, the metapost script below must be converted into individual pdf 
files.
I have tried a couple of different metapost converters online just to 
see what the output would look like and keep getting errors that quit 
the conversion and seem to
have general trouble with the code.  I would really appreciate it if a 
metapost expert could give me advice regarding this code and what might 
be the best conversion
program to use in order to generate pdf image files (the image outputs 
are phylogenetic trees).

Thanks so much,

Mia


beginfig(1)
u:=1in;
pair p[];
pair mid[];
p0 = (0.000000u,0.000000u);
p1 = (0.875000u,2.333333u);
p2 = (1.750000u,0.000000u);
p3 = (2.625000u,4.666667u);
p4 = (3.500000u,0.000000u);
p5 = (4.375000u,2.333333u);
p6 = (5.250000u,0.000000u);
p7 = (3.500000u,7.000000u);
p8 = (7.000000u,0.000000u);

% draw line
draw p[0] -- ( xpart p[0], ypart p[1] ) -- p[1] ;
label.urt( btex 0 etex, p[0]);
xpart mid[0] = xpart(p[0]);
ypart mid[0] = (ypart(p[0])+ypart(p[1]))/2;
label.rt( btex $l = 6$ etex, mid[0]);
label.rt( btex $\lambda=0.001700$ etex, mid[0] + (0,-0.150000u));
draw p[1] -- ( xpart p[1], ypart p[3] ) -- p[3] ;
label.urt( btex 0 \small{(0.500, 0.500)} etex, p[1]);
xpart mid[1] = xpart(p[1]);
ypart mid[1] = (ypart(p[1])+ypart(p[3]))/2;
label.rt( btex $l = 81$ etex, mid[1]);
label.rt( btex $\lambda=0.001700$ etex, mid[1] + (0,-0.150000u));
draw p[2] -- ( xpart p[2], ypart p[1] ) -- p[1] ;
label.urt( btex 0 etex, p[2]);
xpart mid[2] = xpart(p[2]);
ypart mid[2] = (ypart(p[2])+ypart(p[1]))/2;
label.rt( btex $l = 6$ etex, mid[2]);
label.rt( btex $\lambda=0.001700$ etex, mid[2] + (0,-0.150000u));
draw p[3] -- ( xpart p[3], ypart p[7] ) -- p[7] ;
label.urt( btex 0 \small{(0.500, 0.500)} etex, p[3]);
xpart mid[3] = xpart(p[3]);
ypart mid[3] = (ypart(p[3])+ypart(p[7]))/2;
label.rt( btex $l = 6$ etex, mid[3]);
label.rt( btex $\lambda=0.001700$ etex, mid[3] + (0,-0.150000u));
draw p[4] -- ( xpart p[4], ypart p[5] ) -- p[5] ;
label.urt( btex 0 etex, p[4]);
xpart mid[4] = xpart(p[4]);
ypart mid[4] = (ypart(p[4])+ypart(p[5]))/2;
label.rt( btex $l = 17$ etex, mid[4]);
label.rt( btex $\lambda=0.001700$ etex, mid[4] + (0,-0.150000u));
draw p[5] -- ( xpart p[5], ypart p[3] ) -- p[3] ;
label.urt( btex 0 \small{(0.500, 0.500)} etex, p[5]);
xpart mid[5] = xpart(p[5]);
ypart mid[5] = (ypart(p[5])+ypart(p[3]))/2;
label.rt( btex $l = 70$ etex, mid[5]);
label.rt( btex $\lambda=0.001700$ etex, mid[5] + (0,-0.150000u));
draw p[6] -- ( xpart p[6], ypart p[5] ) -- p[5] ;
label.urt( btex 0 etex, p[6]);
xpart mid[6] = xpart(p[6]);
ypart mid[6] = (ypart(p[6])+ypart(p[5]))/2;
label.rt( btex $l = 17$ etex, mid[6]);
label.rt( btex $\lambda=0.001700$ etex, mid[6] + (0,-0.150000u));
draw p[7] -- ( p[7] + (0,0.333333u) );
label.urt( btex 1 \small{(0.015, 0.186)} etex, p[7]);
draw p[8] -- ( xpart p[8], ypart p[7] ) -- p[7] ;
label.urt( btex 0 etex, p[8]);
xpart mid[8] = xpart(p[8]);
ypart mid[8] = (ypart(p[8])+ypart(p[7]))/2;
label.rt( btex $l = 93$ etex, mid[8]);
label.rt( btex $\lambda=0.001700$ etex, mid[8] + (0,-0.150000u));

% add name
labeloffset := 10; dotlabel.bot(btex chimp etex, p[0]);
labeloffset := 10; dotlabel.bot(btex human etex, p[2]);
labeloffset := 10; dotlabel.bot(btex mouse etex, p[4]);
labeloffset := 10; dotlabel.bot(btex rat etex, p[6]);
labeloffset := 10; dotlabel.bot(btex dog etex, p[8]);

% Mark duplication

% annotation
label( btex  etex, (0.1u, 7.700000u));
label( btex ENSF00000002057 etex, (0.1u, 7.500000u));
label(btex Max p-value = 0.000 etex, (0.1u, 7.300000u));
endfig;


beginfig(2)
u:=1in;
pair p[];
pair mid[];
p0 = (0.000000u,0.000000u);
p1 = (0.875000u,2.333333u);
p2 = (1.750000u,0.000000u);
p3 = (2.625000u,4.666667u);



More information about the metapost mailing list