<div dir="ltr"><div><div>The "reg_multlin_stats" documentation <br> <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/reg_multlin_stats.shtml">http://www.ncl.ucar.edu/Document/Functions/Contributed/reg_multlin_stats.shtml</a><br>does state that but the code (located in contributed.ncl) has:<br></div><br>undef ("reg_multlin_stats")<br>function reg_multlin_stats(Y[*]:numeric<wbr>, XP:numeric, opt) ; should be float or double<br>;<br>; Nomenclature:<br>; Y - dependent variable (size NY)<br>; missing values (_FillValue) are not allowed. <=============<br>; XP - one [*] or more [*][*] independent variables.<br>; Missing values (_FillValue) are not allowed. <=============<br><br></div><div>[snip]<br><br> if (isatt(Y,"_FillValue") .and. any(ismissing(Y))) then<br> print("reg_multlin_stats: Y has missing values. Not allowed!")<br> exit<br> end if<br><br> if (isatt(XP,"_FillValue") .and. any(ismissing(XP))) then<br> print("reg_multlin_stats: XP has missing values. Not allowed!")<br> exit<br> end if<br></div><div>[snip]<br></div><div>==========<br><br></div><div>Still, there might be a way to use the function. My thought is to manually loop ove ?grid points? skipping any points that have missing values.<br></div><div><br></div><div>How are you using it? Please be clear and concise.<br><br></div><div>I am gone through 22 Jan 2018 so maybe someone else can help.<br></div><div><br><br><br></div><div> <br></div><div><br></div><div> </div><div><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 10, 2018 at 4:00 PM, Anahita Amiri Farahani via ncl-talk <span dir="ltr"><<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@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">Dear all,<div><br></div><div>I wanted to use reg_multlin_stats to calculate linear regression between low cloud and stability, the low cloud and stability are my dependent and independent variables. Stability is calculated based on the potential temperature at the surface and 700 hPa. So we have missing values for some areas. When I read the description of this function, it says <span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">While missing values are allowed, it is recommended that users not input any missing independent values. It just confuses the results. Missing values should be indicated by the </span><a href="https://www.ncl.ucar.edu/Document/Language/fillval.shtml" style="margin:0px;padding:0px;text-decoration-line:none;color:rgb(133,45,133);font-family:verdana,sans-serif;font-size:13.3333px" target="_blank">_FillValue</a><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"> attribute.</span></div><div><br></div><div><br></div><div>But when I use this function I got this error: <span style="font-family:Courier">reg_multlin_stats: XP has missing values. Not allowed!</span></div><div><span style="font-family:Courier"><br></span></div><div>and it popped out from NCL, Any suggestions? <br></div><div><br></div><div>Thanks,</div><div>Ana</div>
<p class="MsoNormal" style="margin:0in 0in 0.0001pt;font-size:12pt;font-family:Cambria"><span style="font-size:19pt"><span></span></span></p>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>