<div dir="ltr">On Fri, Jul 11, 2014 at 3:28 PM, Jouko Virtanen <span dir="ltr">&lt;<a href="mailto:jouko@umich.edu" target="_blank">jouko@umich.edu</a>&gt;</span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">

<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">Thank Nathaniel.  I am still confused.  Do you intend me to copy and paste the script that you included to a file called score.py.  When I did that I got the following error message<div>

<br></div><div>Syntax error: expected &quot;=&quot;, found &quot;__future__&quot; (file &quot;score.py&quot;, line 1)</div>
<div>Sorry: score.py is not a valid parameter file.</div></div></blockquote><div><br></div><div>That&#39;s the kind of error message I&#39;d expect if you ran &quot;phenix.refine score.py&quot; instead of &quot;phenix.python score.py&quot;.</div>

<div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div>What are the inputs?  I would think that at a minimum I would have to input a pdb file, mtz file, fasta file, and an estimate for the rmsd.</div>

</div></blockquote><div><br></div><div>You would definitely need a PDB file and MTZ file - those will just be additional command-line arguments.  The script I sent you does not require a sequence or estimated RMSD; it also assumes that the model is already placed correctly in the unit cell.  Phaser does require the sequence and RMSD, however; is that what you really wanted to use?  If you actually need to place the model you might as well just run the full Phaser MR_AUTO mode; otherwise you can use the MR_RNP mode.  I had expected to be able to do it like this:</div>

<div><br></div><div>phenix.phaser model.pdb data.mtz seq.fa model_rmsd=1.0 phaser.mode=MR_RNP<br></div><div><br></div><div>although that crashes, which I guess may be a bug.  There are other ways to get the same information but it is not clear what exactly you&#39;re trying to do.</div>

<div><br></div><div>-Nat</div></div></div></div>