2 weeks ago I asked if anyone had a mathematica script of filter.math
edited to plot the results of a direct rotation function. Several people
wanted any results I came up with so here you go. Turns out to be pretty
simple if you understand Mathematica scripting language (which I don't). I
still haven't figured out how to get a simple calculation of the mean
plotted as a straight line through the graph. No one volunteered it
presumably because it's real simple. Okay I'll go RTFM if I can just
figure out which of the 100 of so books on Mathematica is going to cover
this.
direct_rotation.math
(* file xtalmr/filter.math *)
(* These are Mathematica instructions to make a plot of the *)
(* PCs of the refined structures as a function of rotation *)
(* function peak index *)
(* change line 20 to PlotRange->{{0, number +2 }, {0., max +0.02 }}, *)
(* if you want a 0-> max on the y axis *)
Off[General::spell1]
(* change the name of the file "filter.list" *)
(*===>*) data = ReadList["search.dat",
{Number,Number,Number,Number}];
column4= Transpose[data][[4]];
number=Dimensions[column4][[1]];
max=Max[column4];
ListPlot[column4,
AxesLabel->{"Solution Number","Linear Correlation Coefficient"},
Ticks->{Range[0, number , 50],Range[-0.04, max ,0.01]},
PlotLabel->"Direct Rotation Function",AxesOrigin->{0,-0.04},PlotJoined->True,
PlotRange->Automatic,
AxesStyle->{PostScript
["/Times-Roman findfont 13 scalefont setfont"]}];
Fool with this 'Range[-0.04, max ,0.01]' to print (or not) negative values
of the Linear Correlation Coefficient and be sure the x axis doesn't
overlay values making the peak numbers invisible. Range[0.0, max ,0.01]
and AxisOrigin->{0,0} will plot only the positive values.
Now, when I can figure out how to plot 3d graphs like DeLano & Brünger did
in their recent Acta Cryst paper (D51, 744) on this function I'll be
really happy.
--
Michael Skidmore (619) 534-1894
UCSD, Dept of Pathology mskidmore at ucsd.edu
9500 Gilman Drive, Dept 0612 <www-acs.ucsd.edu/~mskidmor/>
La Jolla CA 92093-0612
"Standard orbit, Mr. Sulu"