最尤法の数値計算結果がどうみても正しくない

考えられる原因は：

1. 尤度関数のコーディングに誤りがある
2. 最適化手法の指定に不適切な点が隠れている（誤りではないと思うのだが、不等式制約が加わると最適化の計算はかなり困難になるらしい。）
3. 推定したいパラメータの数に対し、データ量が不十分（ここに原因を求めても仕方がないのだが、生存時間分析ではレコード数ではなくイベント数が重要で、そもそもサンプルサイズが少ない中でイベント率が20％弱となると、適用できる解析手法も限られるという、ごく当たり前の原理である。）

といったところであろうか。

3.については対処のしようがないので、1.と2.の可能性を再度模索するよりない。