<div dir="ltr"><div class="gmail_extra"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear Sir Dennis and Sir Will,</div><div><br></div><div>Many thanks for the help!</div><div>I got it working already.</div><div><br></div><div><br></div><div>Sincerely,</div><div><br></div><div>Lyndz</div></div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Thu, Jun 22, 2017 at 10:56 PM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>[SNIP]<br> dimsst = dimsizes(sst)<br> ntim = dimsst(0)<br> nlat = dimsst(1)<br> mlon = dimsst(2)<br><br>;*****************************<wbr>***********************<br>;Calculate Regression coefficient via regCoef_n<br>;*****************************<wbr>***********************<br><br> rc = regCoef_n(data,sst,0,0)<br> printVarSummary(rc) ;(nlat,mlon) <br> ;[1D arrays]: nptxy,rstd,yintercept,tval<br><br>;;Variable: rc<br>;;Type: float<br>;;Total Size: 259200 bytes<br>;; 64800 values<br>;;Number of Dimensions: 2<br>;;Dimensions and sizes: [180] x [360]<br>;;Coordinates: <br>;;Number Of Attributes: 5<br>;; _FillValue : -1e+30<br>;; nptxy : <ARRAY of 64800 elements><br>;; rstd : <ARRAY of 64800 elements><br>;; yintercept : <ARRAY of 64800 elements><br>;; tval : <ARRAY of 64800 elements><br><br>;*****************************<wbr>***********************<br>;Assign coordinates<br>;*****************************<wbr>***********************<br> rc!0 = "lat"<br> rc!1 = "lon"<br> rc&lat = sst&lat<br> rc&lon = sst&lon<span class=""><br><br>;*****************************<wbr>***********************<br>;Test for significance<br>;*****************************<wbr>***********************<br></span> df = onedtond(rc@nptxy-2, (/nlat,mlon/)) ; degrees of freedom <br> tval = onedtond(rc@tval , (/nlat,mlon/)) ; t-statistic <br></div>[SNIP]<br><br><br></div>The latter is what Will was suggesting.<br></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 22, 2017 at 4:28 AM, Lyndon Mark Olaguera <span dir="ltr"><<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.<wbr>com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Sir Will and NCL experts,<div><br></div><div>I still get the same error <b>(fatal:where: dimension sizes of parameter 0 and parameter 2 do not match)</b>. I am not sure if I am implementing this correctly.</div><div><br></div><div>I changed this part in my code based on the suggestion suggestion:</div><div><br></div><div class="gmail_extra"><div class="gmail_extra">;*****************************<wbr>***********************</div><div class="gmail_extra">;Test for significance</div><div class="gmail_extra">;*****************************<wbr>***********************</div><span><div class="gmail_extra"> df = rc@nptxy-2 ; degrees of freedom </div><div class="gmail_extra"> tval = rc@tval ; t-statistic </div><div class="gmail_extra"> b = tval ; b must be same size as tval (and df) </div><div class="gmail_extra"> b = 0.5</div><div class="gmail_extra"> prob = betainc(df/(df+tval^2),df/2.0,<wbr>b)</div></span><div class="gmail_extra"> printVarSummary(prob);dim=1 size=64800</div><div class="gmail_extra"><br></div><div class="gmail_extra"><b> prob_test =onedtond(prob,(/2,64800/))</b></div><div class="gmail_extra"><b> printVarSummary(prob_test)</b></div><div class="gmail_extra"><br></div><div class="gmail_extra"> rc95 = where(prob_test.lt.0.95,rc@_Fi<wbr>llValue,rc)</div><div class="gmail_extra"> copy_VarCoords(rc, rc95)</div><div class="gmail_extra"> printVarSummary(rc95)</div></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><div class="gmail_extra">I'll appreciate any suggestion on how I can do this correctly.</div><div class="gmail_extra"><br></div><div class="gmail_extra">Sincerely,</div><div class="gmail_extra"><br clear="all"><div><div class="m_-6182335580252815879m_-4990892993453598119gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><font face="comic sans ms, sans-serif"><b><u>Lyndz</u></b></font></div></div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote"><div><div class="m_-6182335580252815879h5">On Thu, Jun 22, 2017 at 3:34 PM, Will Hobbs <span dir="ltr"><<a href="mailto:will.hobbs@utas.edu.au" target="_blank">will.hobbs@utas.edu.au</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div class="m_-6182335580252815879h5">
<div bgcolor="white" lang="EN-US">
<div class="m_-6182335580252815879m_-4990892993453598119gmail-m_2220355447108225520WordSection1">
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri">Change the dimensions of ‘prob’ so that it matches the dimensions of rc.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri">You can do this using the onedtond() function; you’ll find a complete description of this function on the NCL web page.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri">Will<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri"><u></u> <u></u></span></p>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(181,196,223);padding:3pt 0cm 0cm">
<p class="MsoNormal"><b><span style="font-family:Calibri;color:black">From: </span>
</b><span style="font-family:Calibri;color:black"><<a href="mailto:ncl-talk-bounces@ucar.edu" target="_blank">ncl-talk-bounces@ucar.edu</a>> on behalf of Lyndon Mark Olaguera <<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.c<wbr>om</a>><br>
<b>Date: </b>Thursday, 22 June 2017 at 4:28 PM<br>
<b>To: </b>"<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
<b>Subject: </b>[ncl-talk] Question about regressing a time series and a gridded data using regCoef_n<u></u><u></u></span></p>
</div><div><div class="m_-6182335580252815879m_-4990892993453598119gmail-h5">
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">Dear NCL experts,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I am trying to calculate the linear regression coefficient between a single time series and a gridded sst using the regCoef_n in ncl. I am using ncl v.6.4.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Attached is my script and the text file that I am regressing. Here is the link to the sst data:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">[Link to the data]<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><a href="https://drive.google.com/file/d/0B9faET7Bc2o8TjNSWnd1ODREcVU/view?usp=sharing" target="_blank">https://drive.google.com/file/<wbr>d/0B9faET7Bc2o8TjNSWnd1ODREcVU<wbr>/view?usp=sharing</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">What I wanted to do:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">1. I want to extract and plot significant coefficients (above 95% confidence level), so I have this line in the script:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"> df = rc@nptxy-2 ; degrees of freedom <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> tval = rc@tval ; t-statistic <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> b = tval ; b must be same size as tval (and df) <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> b = 0.5 <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><b> prob = betainc(df/(df+tval^2),df/2.0,<wbr>b) </b><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><b> rc95 = where(prob.lt.0.95,rc@_FillVal<wbr>ue,rc)</b><u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I encountered an error after this:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><b>fatal:where: condition variable (parameter 0) dimension mismatch with parameter 2</b><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I understand that the <b>"prob"</b> variable has only one dimension and the "<b>rc"</b> has 2 dimensions (lat & lon). <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Is it possible to assign coordinates to "prob" similar to "rc" other than the copyVarCoords?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I'll appreciate any suggestions on how to do this correctly.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Sincerely,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Lyndz<u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div></div></div>
<p style="font-size:10pt;line-height:10pt;font-family:Calibri,sans-serif"><br>
<br>
University of Tasmania Electronic Communications Policy (December, 2014). <br>
This email is confidential, and is for the intended recipient only. Access, disclosure, copying, distribution, or reliance on any of it by anyone outside the intended recipient organisation is prohibited and may be a criminal offence. Please delete if obtained
in error and email confirmation to the sender. The views expressed in this email are not necessarily the views of the University of Tasmania, unless clearly intended otherwise.
</p>
</div>
<br></div></div>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div>