On Tue, Aug 3, 2010 at 11:21 AM, Joseph Noel <span dir="ltr"><<a href="mailto:noel@salk.edu">noel@salk.edu</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Thanks so much for the suggestions. They are really helpful! A few questions relative to Phenix. Is there a way to check with Phenix the thin resolution shells? I created Free Rs in Phenix using the thin resolution shells options but removed the 2000 limit and instead used 5%. Maybe I over-did it relative to your reply. Just wondering how to check each of these thin shells. Alternatively, does Phenix ensure one has an adequate number of test reflections in each "thin shell" if this option is set? If so, what should I used for the maximum number of reflections and/or %?<br>
</blockquote><div><br></div><div>There isn't a good way to check the shells in the GUI - but all of the pieces are there, and I've been meaning to add some utilities for examining reflection files anyway. �For now, on the command line, run this:</div>
<div><br></div><div>iotbx.r_free_flags_accumulation data.mtz</div><div><br></div><div>and it will print out something like this:</div><div><br></div><div><div><div>Number of work/free reflections by resolution:</div><div>
�� � � � � � � � � � � � � � � � � � � � � work �free �%free</div><div>��bin �1: 135.8135 - � 4.3073 [7861/7861] �7663 � 198 � 2.5%</div><div>��bin �2: � 4.3073 - � 3.4187 [7682/7682] �7484 � 198 � 2.6%</div><div>��bin �3: � 3.4187 - � 2.9866 [7645/7645] �7444 � 201 � 2.6%</div>
<div>��bin �4: � 2.9866 - � 2.7135 [7605/7605] �7406 � 199 � 2.6%</div><div>��bin �5: � 2.7135 - � 2.5190 [7593/7593] �7395 � 198 � 2.6%</div><div>��bin �6: � 2.5190 - � 2.3704 [7618/7618] �7416 � 202 � 2.7%</div><div>��bin �7: � 2.3704 - � 2.2517 [7535/7535] �7337 � 198 � 2.6%</div>
<div>��bin �8: � 2.2517 - � 2.1537 [7573/7573] �7374 � 199 � 2.6%</div><div>��bin �9: � 2.1537 - � 2.0708 [7556/7556] �7354 � 202 � 2.7%</div><div>��bin 10: � 2.0708 - � 1.9993 [7547/7547] �7348 � 199 � 2.6%</div><div>�� � � � � � � � � � � � � � � � �overall 74221 �1994 � 2.6%</div>
</div></div><div><br></div><div>This is a real example, using the default settings in the reflection file editor (20 thin shells for the test set). �It looks much better than what I had originally done with this dataset five years ago using a different program, where the last shell had only 0.2% of reflections flagged. �However, the code that does this isn't particularly sophisticated, so I would recommend double-checking the output for your data. �I don't think you're going to do any harm by sticking with 5% and removing the absolute limit.</div>
<div><br></div><div>-Nat</div></div>