<div>Hi you all!</div>
<div>&nbsp;</div>
<div>I have found myself the answer to my inquiry.</div>
<div>&nbsp;</div>
<div>For this may help someone else, the way to estimate a linear trend (or regression)&nbsp; on one-monthly series from a multi-year multi-monthly data file, wihtout having to seperate the chosen month interannual series into a single file, is just to perform the regression calculus for that month!</div>

<div>&nbsp;</div>
<div>I help myself from other previous prompts from the GRADSUSRS Archive.</div>
<div>&nbsp;</div>
<div>Here it is the script that can do that:</div>
<div>&nbsp;</div>
<div>&#39;sdfopen <a href="http://file.nc">file.nc</a>&#39;<br>&#39;set lev &#39;k&#39;&#39;<br>&#39;set lat -90 0&#39;<br>&#39;set lon&nbsp; -180 180&#39;<br>&#39;set t 1 276&#39;<br>&#39;define yvar = z&#39; *this is the variable on which to estimate the linear trend</div>

<div><br>* the lineal variable 1, 2, ...,&nbsp;276 to&nbsp;estimate &nbsp;m&nbsp;and c&nbsp;of y = mx + c<br>* y = mx + c<br>* y = amp*x + ave</div>
<div>&#39;set lat -60&#39;<br>&#39;set lon&nbsp; -180&#39;<br>&#39;set t 1 276&#39;<br>&#39;xvar = yvar*0 + 1&#39;<br>&#39;xvar = const(xvar,1,-u)&#39;<br>&#39;xvar = 277 - sum(xvar,t+0,t=276)&#39;</div>
<div>&#39;set lat -90 0&#39;<br>&#39;set lon&nbsp; -180 180&#39;<br>&#39;set t 1&#39;</div>
<div>&#39;define p = sum(xvar,t=10,t=276,12)&#39; * starttime=10, supposing you want to estimate October linear trend<br>&#39;define q = sum(yvar,t=10,t=276,12)&#39;<br>&#39;define r = sum(xvar*yvar,t=10,t=276,12)&#39;<br>
&#39;define s = sum(xvar*xvar,t=10,t=276,12)&#39;</div>
<div>&#39;define rn = ave(xvar*yvar,t=10,t=276,12)&#39;<br>&#39;define n = r/rn&#39;</div>
<div>&#39;define amp =&nbsp; (n*r - p*q)/(n*s - p*p)&#39;&nbsp; *this is m<br>&#39;define ave = (q - amp*p)/n&#39;</div>
<div>&#39;set t 1 276&#39;<br>&#39;monmask 10&#39;&nbsp;&nbsp;&nbsp; * this is for month October, you need the <a href="http://monmas.gs">monmas.gs</a> script from the Archive<br>&#39;define trend=amp*xvar*mmask + ave&#39; * this is the final linear trend for Oct alone<br>
</div>
<div>&nbsp;</div>
<div>That`s all!!</div>
<div>&nbsp;</div>
<div>Bye.</div>
<div>&nbsp;</div>
<div>Eduardo.</div>
<div>&nbsp;</div>
<div><br><br>&nbsp;</div>
<div class="gmail_quote">2008/7/2 Eduardo Agosta Scarel &lt;<a href="mailto:eduardo.agosta@gmail.com">eduardo.agosta@gmail.com</a>&gt;:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div class="gmail_quote">
<div>Dear Grads&#39; Users,</div>
<div>&nbsp;</div>
<div>I am trying to extimate linear trend&nbsp;in an interannual series of one-monthly data from a multy-year monthly data file.</div>
<div>I can estimate properly the linear trend using the script below, if the&nbsp;data file has the one-monthly data series previously separately. The problem is when using the&nbsp;full monthly multy-year data file.</div>
<div>&nbsp;</div>
<div>How can I just select one given month from the data , for instance June, and estimate the linear trend for the June-series alone, having a multy-year data file with all the months of the year?<br clear="all"></div>
<div>I have tried using the function monmask (with argument 6), but I cannot solve it properly.</div>
<div>&nbsp;</div>
<div>Any suggestion I will much appreciate.</div>
<div>&nbsp;</div>
<div>Eduardo.<br></div>
<div>&nbsp;</div>
<div>reinit&#39;<br>&#39;sdfopen d:\anterior\era\<a href="http://mes.nc/" target="_blank">mes.nc</a>&#39;</div>
<div>&#39;set clopts -1 -1 0.16&#39;<br>&#39;set mproj sps&#39;</div>
<div>k=70</div>
<div><br>&#39;set grads off&#39;<br>&#39;set lev &#39;k&#39;&#39;<br>&#39;set mproj sps&#39;<br>&#39;set lat -90 0&#39;<br>&#39;set lon&nbsp; -180 180&#39;</div>
<div><br>&#39;set lat -60&#39;<br>&#39;set lon&nbsp; -180&#39;<br>&#39;set t 1 132&#39;</div>
<div>* contruyo variable lineal 1, 2, ..., 276 to para calcular m y c de y = mx + c<br>&#39;xvar = 100*t*0 + 1&#39;<br>&#39;xvar = const(xvar,1,-u)&#39;<br>&#39;xvar = 133 - sum(xvar,t+0,t=132)&#39;</div>
<div>&nbsp;</div>
<div>&#39;set lat -90 0&#39;<br>&#39;set lon&nbsp; -180 180&#39;<br>&#39;set t 1&#39;</div>
<div><br>&#39;define coeff = tregr(xvar,100*t,t=1, t=132)&#39;<br>&#39;define tave = ave(100*t, t=1, t=132)&#39;<br>&#39;define xvarave = ave(xvar, t=1, t=132)&#39;<br>&#39;define trend = coeff*(xvar - xvarave) + tave&#39;</div>
</div></blockquote></div>