<div class="gmail_quote">
<div>Hi Mark,</div>
<div> </div>
<div>Thank you very  much for your help. I got my answer. </div>
<div> </div>
<div>I chose protein as the group for calculating SASA and residue1 for output and i got the SASA of residue 1 vs time. This gave me the correct mean SASA of residue1 and also its std deviation. </div>
<div>I repeated the same way for all the residues and got SASA per residue with std deviation of each.</div>
<div> </div>
<div>Thanks again!!</div>
<div> </div>
<div>Regards,</div>
<div>Harpreet</div>
<div> </div>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">On 14/12/2011 4:01 PM, Harpreet Basra wrote:<br>&gt; Hi<br>&gt; I want to calculate the standard deviation of SASA of each residue of<br>
&gt; a peptide. When using -or option for complete peptide in g_sas we get<br>&gt; residue vs sasa only.<br><br>I don&#39;t follow what result you are describing as unsuitable for you needs.<br><br>&gt; So i tried by creating an index file for each residue separately<br>
&gt; (residue1.ndx, residue2.ndx and so forth) and then by calculating by<br>&gt; using -o option i got time vs sasa for each residue, thinking that<br>&gt; this will give me average sasa per residue as well as its standard<br>
&gt; deviation. But I see that the average value obtained by this method<br>&gt; for each residue is different from the one i got by using -or option<br>&gt; (and calculating for complete peptide). Why is it so? Doesnt -or<br>
&gt; option gives average sasa over entire trajectory per residue? Can<br>&gt; anyone suggest any way to compute the standard deviation?<br>&gt;<br><br>It&#39;s hard to comment when we don&#39;t know which groups you chose for which<br>
calculation. Using an index file to change *only* the output group and<br>using g_sas -or should definitely work, but I don&#39;t see why custom<br>groups should be needed at all.<br><br>Mark<br><br></blockquote></div>