<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p style="margin-top:0;margin-bottom:0">Hi there,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">I am using NCL 6.3.0 and following example 1 (using the same data) from here:
<a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/reg_multlin_stats.shtml" class="OWAAutoLink" id="LPlnk405831" previewremoved="true">
https://www.ncl.ucar.edu/Document/Functions/Contributed/reg_multlin_stats.shtml</a><br>
</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">For some reason my p-value vector is <span>( 0.1488179, 0.006746806, 0.007730724 )</span> which differs somewhat from the example 1 result. I have checked against R's multiple linear regression function which matches
the NCL example 1 exactly. Is it possible version 6.3.0 from 6.2.0 introduced a bug in the calculation of reg_multin_stats p-values? It seems odd that my stderr and tvals match exactly to the example (and R) just not the p-values.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Thanks,</p>
<p style="margin-top:0;margin-bottom:0">Peter<br>
</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">--------------------------------<br>
</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">ANOVA outout is as follows:<br>
</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0"></p>
<div>Number Of Attributes: 31<br>
_FillValue : 9.96921e+36<br>
long_name : multiple regression coefficients<br>
model : Yest = b(0) + b(1)*X1 + b(2)*X2 + ...+b(M)*XM<br>
N : 17<br>
NP : 2<br>
M : 3<br>
bstd : ( 0, 0.5029772, 0.4922954 )<br>
SST : 13239.72<br>
SSE : 423.3723<br>
SSR : 12816.36<br>
MST : 827.4825<br>
MSE : 30.24088<br>
MSE_dof : 14<br>
MSR : 6408.178<br>
RSE : 5.499171<br>
RSE_dof : 13<br>
F : 211.9045<br>
F_dof : ( 2, 14 )<br>
F_pval : 3.419033e-11<br>
r2 : 0.9680232<br>
r : 0.9838817<br>
r2a : 0.9634551<br>
fuv : 0.03197682<br>
Yest : <ARRAY of 17 elements><br>
Yavg : 288.4<br>
Ystd : 28.766<br>
Xavg : ( 55.35294, 30.9 )<br>
Xstd : ( 11.68028, 1.172071 )<br>
stderr : ( 100.8641, 0.3945341, 3.931675 )<br>
tval : ( -1.521964, 3.139713, 3.073079 )<br>
pval : ( 0.1488179, 0.006746806, 0.007730724 )<br>
(0) -153.5115<br>
(1) 1.238724<br>
(2) 12.08235</div>
<br>
<p></p>
</div>
</body>
</html>