Python中蒙特卡罗生存分析的最佳实践

Question

我正在使用python创建一个蒙特卡罗模拟，它正在寻找不同类别的夫妻之间的离婚率模型。

总结详情：

• 婚姻可以是＆＃34;安排婚姻＆＃34;或者基于爱情 婚姻＆＃34;概率分配给两者。

• 在每种婚姻类型中，有多个子组的婚姻，为简单起见，我们说有两个：＆＃34;十几岁婚姻＆＃34;，＆＃34;非青少年婚姻＆＃34 ;，概率分配给两者。

• 他们可以离婚有两个原因：＆＃34;由于不忠而离婚＆＃34;，＆＃34;由于其他原因而离婚＆＃34;及其各自的概率。

数据的结构如下并且已经存在多年：

我需要建立一张生存表，详细说明他们离婚的月份和原因。

目前，我只是为每对新人（n = 10000）运行一个循环进行数千次试验。

通过（1）根据各自的概率随机分配给每个组，（2）掷骰子是否离婚以及出于什么原因，然后（3）将结果附加到生存表，记录所有相关信息。

这让我感到困扰，因为我觉得这不是解决这个问题的最佳方法，但我似乎无法弄清楚是否有更优雅，计算效率更高的方法来解决这个问题。我真诚地感谢任何人都可以提供的见解。