Help! SAS NLMIXED. The final Hessian matrix is not positive definite?
- From: "sphuang via MathKB.com" <u39394@uwe>
- Date: Wed, 28 Nov 2007 20:03:45 GMT
Hi, there,
I am currently running sas PROC NLMIXED to analyze my data. My sample are 27
animals. Their running speed were reperatedly measured at 9 temperature
treatments. The running speed is tested whether it can be predicted by
temperatues using a logistic-expotential function. Therefore, the function I
wrote in SAS code is an logistic-expotential.
My sas code is listed as following text(A). The problem I faced is that there
is a warning in log page, which is as following text(B). Could anyone give me
any suggestion to patch it on? any suggestions would be highly appreciated.
many many thanks.
Shu-Ping
text(A), My sas code:
proc nlmixed data= sprint_30s method=firo ;
parms b1=0.04 b6=4.5 b2=0.2 b4=0.2 U1=0 U2=0 VAR=0.1 t11=0.1 t22=0 t12=0;
bounds VAR>=0;
LOW=1/(0.015+U1+b6*exp(-b2*(Temp-2))); HIGH=1-exp(b4*(Temp-43.1)+U2);
pred=b1*LOW*HIGH;
VAR1=t11*t11;
cov12=t11*t12;
VAR2=t12*t12+t22*t22;
model speed_m~normal (pred,VAR);
random U1 U2 ~normal ([0,0],[VAR1, cov12,VAR2]) subject=ID;
run;
text B:
WARNING: The final Hessian matrix is not positive definite, and therefore the
estimated
covariance matrix is not full rank and may be unreliable. The
variance of some
parameter estimates is zero or some parameters are linearly related
to other
parameters.
--
Message posted via http://www.mathkb.com
.
- Follow-Ups:
- Re: Help! SAS NLMIXED. The final Hessian matrix is not positive definite?
- From: Herman Rubin
- Re: Help! SAS NLMIXED. The final Hessian matrix is not positive definite?
- From: vontressms
- Re: Help! SAS NLMIXED. The final Hessian matrix is not positive definite?
- Prev by Date: Re: What number in my formula is missing for finding 95% intervals?
- Next by Date: Re: Any test about whether data are independent?
- Previous by thread: Tests for distributional assumptions in Tobit models
- Next by thread: Re: Help! SAS NLMIXED. The final Hessian matrix is not positive definite?
- Index(es):
Relevant Pages
|
Loading