<div dir="ltr"><div style="font-size:12.8px">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" target="_blank" style="color:rgb(133,45,133);margin:0px;padding:0px;text-decoration-line:none;font-family:verdana,sans-serif;font-size:13.3333px">_FillValue</a><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"> attribute.</span></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">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 style="font-size:12.8px"><span style="font-family:Courier"><br></span></div><div style="font-size:12.8px"><span style="font-family:Courier"><br></span></div><div style="font-size:12.8px"><span style="font-family:Courier">Here is my code:</span></div><div style="font-size:12.8px"><span style="font-family:Courier"><br></span></div><div style="font-size:12.8px"><span style="font-family:Courier"><br></span></div><div style="font-size:12.8px"><div><div>ana=addfile("clouds_E_total.<wbr>nc","r")</div><div>lat=ana->lat</div><div>lon=ana->lon</div><div>lev=ana->lev</div><div>CLDLOW=ana->CLDLOW</div><div>CLDMED=ana->CLDMED</div><div>CLDHGH=ana->CLDHGH</div><div>CLDTOT=ana->CLDTOT</div><div>CLOUD=ana->CLOUD</div><div>delete(ana)</div></div><div><br></div><div><br></div><div><div>ana=addfile("stability_E_<wbr><a href="http://total.nc">total.nc</a>","r")</div><div>theta_ANN=ana->theta_ANN</div><div>theta_summer=ana->theta_<wbr>summer1</div><div>theta_winter=ana->theta_<wbr>summer1</div><div>theta_fall=ana->theta_fall1</div><div>theta_spring=ana->theta_<wbr>spring1</div><div>delete(ana)</div></div><div><br></div><div><br></div><div><div><br></div><div><br></div><div>CLDLOW_1d=ndtooned(CLDLOW)</div><div>theta_ANN_1d=ndtooned(theta_<wbr>ANN)</div><div><br></div><div><br></div><div>opt = True</div><div>opt@print_anova = True     ; optional</div><div>opt@print_data  = True</div><div><br></div><div>b   =  reg_multlin_stats(CLDLOW_1d,<wbr>theta_ANN_1d,opt) ; partial regression coef</div></div><div><br></div><div><br></div><div>what do you suggest?</div><div><br></div><div>Thanks</div><div>Ana</div></div></div>