这种数据包含有 右删失 数据,需要先把似然函数写出来,几何分布是离散分布,需要查一下书写出似然函数,然后优化这个函数,可以用 R中的 optimize,几何分布只有一个参数需要估计.
似然函数可在 陈家鼎 生存分析与可靠性 北京大学出版社 2005(第五章第三节) 查到或参考
(Kalbfleisch, J. D. and Prentice, R. L. (2002). The Statistical Analysis of Failure Time Data(Second Edition). New York: John Wiley,6.2节)