<BODY contentEditable=true onload=parent.init()><P>Weerd, </P>
<P>my source code shows that the default is ml:</P>
<P>cd /net/rosie/scratch1/phzwart/CCTBX/Sources/mmtbx/mmtbx/bulk_solvent/</P>
<P>more bulk_solvent_and_scaling.py</P>
<P>master_params = iotbx.phil.parse("""\<BR>&nbsp; bulk_solvent = True<BR>&nbsp;&nbsp;&nbsp; .type = bool<BR>&nbsp; anisotropic_scaling = True<BR>&nbsp;&nbsp;&nbsp; .type = bool<BR>&nbsp; statistical_solvent = False<BR>&nbsp;&nbsp;&nbsp; .type = bool<BR>&nbsp; k_sol_b_sol_grid_search = True<BR>&nbsp;&nbsp;&nbsp; .type = bool<BR>&nbsp; minimization_k_sol_b_sol = True<BR>&nbsp;&nbsp;&nbsp; .type = bool<BR>&nbsp; minimization_b_cart = True<BR>&nbsp;&nbsp;&nbsp; .type = bool<BR>&nbsp; target = ls_wunit_k1 *ml<BR>&nbsp;&nbsp;&nbsp; .type = choice</P>
<P>Pavel developed some heuristics governing&nbsp;optimal settings in&nbsp;certain cases&nbsp;which might come into play here.&nbsp; I suspect it might be due to model size, resolution or other factors. If not, I am puzzled.</P>
<P>Cheers</P>
<P>Peter&nbsp;</P>
<P><BR><BR>----- Original Message ----- <BR>From: Jianghai Zhu &lt;jzhu@cbr.med.harvard.edu&gt; <BR>Date: Sunday, March 25, 2007 5:18 pm <BR>Subject: Re: [phenixbb] ls_wunit_k1 <BR>To: PHENIX user mailing list &lt;phenixbb@phenix-online.org&gt; <BR><BR>&gt; Hi Peter, <BR>&gt; <BR>&gt; You are right.&nbsp; Somehow, the default target function in&nbsp; <BR>&gt; refinement.bulk_solvent_and_scale is LS, not ML. <BR>&gt; <BR>&gt; Jianghai <BR>&gt; <BR>&gt; +++++++++++++++++++++++++++++++ <BR>&gt; Jianghai Zhu, Ph.D <BR>&gt; CBR Institute for Biomedical Research <BR>&gt; Department of Pathology <BR>&gt; Harvard Medical School <BR>&gt; 200 Longwood Ave., Boston, MA 02115 <BR>&gt; Ph: 617-278-3211 <BR>&gt; Fx: 618-278-3232 <BR>&gt; +++++++++++++++++++++++++++++++ <BR>&gt; <BR>&gt; <BR>&gt; <BR>&gt; On Mar 25, 2007, at 8:18 PM, Petrus H Zwart wrote: <BR>&gt; <BR>&gt; &gt; Hi Jianghai, <BR>&gt; &gt; <BR>&gt; &gt; Interesting, I though that by default an ML target is used.&nbsp; <BR>
&gt; &gt; Something must have happened somewhere. <BR>&gt; &gt; <BR>&gt; &gt; Check the parameter file and look for the section on bulk <BR>&gt; solvent&nbsp; <BR>&gt; &gt; scaling ( refinement.bulk_solvent_and_scale ). A target function <BR>&gt; <BR>&gt; &gt; can be defined, you have the choice between ml and ls_wunit_k1, <BR>&gt; ml&nbsp; <BR>&gt; &gt; is (I think) selected by default. <BR>&gt; &gt; <BR>&gt; &gt; There is also a switch present that triggers the refinement of <BR>&gt; bulk&nbsp; <BR>&gt; &gt; solvent parameters using a least squares target if the r value <BR>&gt; is&nbsp; <BR>&gt; &gt; larger than a certain value: <BR>&gt; &gt; <BR>&gt; &gt; ls_bulk_solvent_and_scale_if_r_greater_then = 0.99 <BR>&gt; &gt; <BR>&gt; &gt; Changing the parameter you propose changes the least squares <BR>&gt; target&nbsp; <BR>&gt; &gt; type if the main target selected is ls and has (as far as I&nbsp; <BR>
&gt; &gt; understand) no influence on bulk solvent scaling. If you want, <BR>&gt; you&nbsp; <BR>&gt; &gt; can send a log file for a closer look. <BR>&gt; &gt; <BR>&gt; &gt; HTH <BR>&gt; &gt; <BR>&gt; &gt; Peter <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; ----- Original Message ----- <BR>&gt; &gt; From: Jianghai Zhu &lt;jzhu@cbr.med.harvard.edu&gt; <BR>&gt; &gt; Date: Sunday, March 25, 2007 2:58 pm <BR>&gt; &gt; Subject: [phenixbb] ls_wunit_k1 <BR>&gt; &gt; To: PHENIX user mailing list &lt;phenixbb@phenix-online.org&gt; <BR>&gt; &gt; <BR>&gt; &gt; &gt; Hi, <BR>&gt; &gt; &gt; <BR>&gt; &gt; &gt; I happen to notice that the target function of "Bulk solvent <BR>&gt; &gt; &gt; correction and scaling" in a phenix.refine run is ls_wunit_k1. <BR>&gt; &gt; &gt; Does <BR>&gt; &gt; &gt; that mean LS function for bulk solvent correction and scaling is <BR>
&gt; &gt; &gt; the&nbsp; &gt; default in phenix.refine?&nbsp; If I uncheck ls_wunit_k1 in <BR>&gt; &gt; &gt; refinement.ls_target_names in my parameter file, will <BR>&gt; &gt; &gt; phenix.refine <BR>&gt; &gt; &gt; use the ML function or something bad will happen? <BR>&gt; &gt; <BR>&gt; &gt; &gt; <BR>&gt; &gt; &gt; Jianghai <BR>&gt; &gt; &gt; <BR>&gt; &gt; &gt; +++++++++++++++++++++++++++++++ <BR>&gt; &gt; &gt; Jianghai Zhu, Ph.D <BR>&gt; &gt; &gt; CBR Institute for Biomedical Research <BR>&gt; &gt; &gt; Department of Pathology <BR>&gt; &gt; &gt; Harvard Medical School <BR>&gt; &gt; &gt; 200 Longwood Ave., Boston, MA 02115 <BR>&gt; &gt; &gt; Ph: 617-278-3211 <BR>&gt; &gt; &gt; Fx: 618-278-3232 <BR>&gt; &gt; &gt; +++++++++++++++++++++++++++++++ <BR>&gt; &gt; &gt; <BR>&gt; &gt; &gt; <BR>&gt; &gt; &gt; <BR>&gt; &gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; _______________________________________________ <BR>&gt; &gt; phenixbb mailing list <BR>&gt; &gt; phenixbb@phenix-online.org <BR>
&gt; &gt; http://www.phenix-online.org/mailman/listinfo/phenixbb <BR>&gt; <BR>&gt; </P></BODY>