BookmarkSubscribeRSS Feed
ElhamAssadi
Calcite | Level 5

Hi everyone,

Hope I can find a solution for this! I am fitting models to my data mostly broken-line and curvilinear plateau by Proc NLIN procedure. In some models there is a note: Convergence criterion met but a note in the log indicates a possible problem with the model! But it makes the model and gives the estimates but then I have no Std error for my estimates and no 95% confidence limits!

In log I have: (I think my initial parameters are correct)

DER.L not initialized or missing. It will be computed automatically.

NOTE: DER.V not initialized or missing. It will be computed automatically.

NOTE: DER.R not initialized or missing. It will be computed automatically.

NOTE: PROC NLIN grid search time was  0: 0: 0.

NOTE: Convergence criterion met.

NOTE: The (approximate) Hessian is singular.

NOTE: The data set WORK.PPP has 507 observations and 3 variables.

And in some models I have a very big 95% confidence limits which my papers reviewer ask why!?!?

Any idea and comment??

My SAS syntax for model:

proc nlin data=one; *straight broken-line;

parameters L=0.70 U=-1.4 R=0.80;

z1=(x<R)*(R-x);

model y = L + U*(z1);

output out=ppp p=pred;

run;

proc gplot;

title2 '2 linear broken lines';

goptions hpos=35 vpos=35 ftext=swiss;

symbol1 v=dot c=black;

symbol2 i=join v=none c=black;

plot y*x pred*x/overlay;

run;

and results:

ParameterEstimateApprox
Std Error
Approximate 95% Confidence
Limits
L0.690.009490.66370.7163
U-0.90.2324-1.5452-0.2548
R0.8...
3 REPLIES 3
SteveDenham
Jade | Level 19

A singular hessian often means that the model is overspecified, but here I think it is because the broken stick model does not have continuous derivatives at the break point.  You may have to specify the derivatives (in DER. statements) to get around this.

Steve Denham

ElhamAssadi
Calcite | Level 5

Dear Steve,

Thank you very much for your help. Just I have no DER. statement in my SAS syntax! do you mean initial parameters?

Also do you have any idea what is the reason to have quite large confidence limits?

Many thanks again

SteveDenham
Jade | Level 19

No, I meant specifying the derivatives in a DER. statement.  But further thought makes that seem unreasonsable.  I fit the data in the documentation for Example 67.1 Segmented Model to a broken stick model, and found the same problem--Hessian is singular, huge standard error, and an indeterminate correlation matrix.  I think it has to do with a lack of smoothness in this model and that the second derivatives are not continuous.  I also tried using PROC MODEL and the results included a statement that the model was singular and the results were biased.

So I looked at http://www.ats.ucla.edu/stat/sas/faq/nlin_optimal_knots.htm, and found what I think may be useful to you.  It worked on the sample data I tried, and the Hessian error message was no longer showing up.

Steve Denham

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 3 replies
  • 5120 views
  • 3 likes
  • 2 in conversation