我正在尝试为具有伽马脆弱性的 cox 模型的似然函数编写 R 代码。我知道 R 中有一些包可以很容易地做到这一点。但是我想为我自己的数据而不是使用包编写具有伽马脆弱性的 cox 模型的自定义似然函数。
ID Start_Time End_Time Gap_Time Status X1 X3 X3
1 0 10 10 1 12 A 150
1 11 90 79 0 16 B 601
2 0 45 45 1 25 A 278
2 46 69 23 1 18 A 895
2 70 90 20 0 23 C 246
3 0 14 14 1 9 A 516
3 15 34 19 1 17 D 385
3 35 61 26 1 10 B 520
3 62 85 23 1 5 C 669
3 86 90 4 0 32 C 553
由于我认为自己是 R 和统计学的高级初学者,因此我需要这方面的帮助来学习如何在 R 中逐步完成它。我想知道是否有任何可用的代码可以参考?
我感谢任何帮助和建议。