|
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]
Re: st: survival analysis - iterations "not concave"
Maximizing likelihood is like climbing a hill to the top. "Not concave"
simply means that when Stata calculates the numerical derivatives around
the current point, it might appear to be in a valley or a saddle, with
no clear direction to the top.
P
Simon Oertel wrote:
Dear Stata list members,
I am using parametric models for a survival analysis. To figure out
which model fits my data best, I used the common procedures (theoretical
implications, graphical plot of the cumulative hazard-rate against the
cox-snell residuals, AIC values). For calculating the AIC values of the
different (standard) models (exponential, weibull, gompertz,
log-logistic, log-normal, gamma) I first used the "streg" command and
fit each of them to my data. Finally, I found the model that will fit my
data best (weibull) - but while running the iteration-processes for the
log-logistic, log-normal and gamma model, STATA tells me for the first
(and for the log-logistic also for the second) iteration that it is "not
concave".
Also if I will not use a log-logistic, log-normal or gamma model for my
data, I am really interested in figuring out what this message
implicates. I suggested that it is not a big problem, as long as the
following iterations (especially the last one) will be concave. I think
that the message "not concave" only tells me that during the iteration
process STATA found a local maximum (what depending on the local maximum
could still be a problem)- but again I am not sure about that.
I appreciate any answer and help.
Thanks
Simon
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
--
E. Paul Wileyto, Ph.D.
Assistant Professor of Biostatistics
Tobacco Use Research Center
School of Medicine, U. of Pennsylvania
3535 Market Street, Suite 4100
Philadelphia, PA 19104-3309
215-746-7147
Fax: 215-746-7140
[email protected]
http://mail.med.upenn.edu/~epw/
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/