打包以获得具有随机效应的生存函数的估计？

Question

我想估算一个包含随机（和固定）效应的函数中的累积危害。 估计应该是半参数的（用样条估计，压电常数指数等）

我尝试了R软件包frailtypack，但是additivePenal函数不能很好地工作（它不能与分层项收敛）。

有人建议使用R函数来获取生存函数（带有分层项）和随机效应的估计吗？

非常感谢！