R 罚款及罚款;灰色回归误差

R 罚款及罚款;灰色回归误差,r,regression,cox-regression,R,Regression,Cox Regression,如果仅使用1个cov执行此FG回归,则在没有任何参数tf定义的情况下,不会获得任何错误: > str(fu) num [1:87] 15 27 3 10 16 13 65 64 5 29 ... > str(tdm) Factor w/ 2 levels "0","1": 2 2 1 2 2 2 2 2 2 2 ... > str(gfr) Factor w/ 4 levels "1","2","3","4": 2 4 2 1 4 3 3 2 3 2 ... > s

如果仅使用1个cov执行此FG回归,则在没有任何参数tf定义的情况下,不会获得任何错误:

> str(fu)
 num [1:87] 15 27 3 10 16 13 65 64 5 29 ...
> str(tdm)
 Factor w/ 2 levels "0","1": 2 2 1 2 2 2 2 2 2 2 ...
> str(gfr)
 Factor w/ 4 levels "1","2","3","4": 2 4 2 1 4 3 3 2 3 2 ...
> str(exitus1)
 Factor w/ 3 levels "0","1","2": 3 1 1 2 1 3 1 1 1 2 ...

> fu
[1]  15  27   3  10  16  13  65  64   5  29 103  72 121 162  22         139 149
[18] 144  26  38  24  19 196 212 215  12 293  48  22 284   2  31   2   5
[35] 319   4  11  15 305 366   5 388  32 386  34  15   4 436 431 436  12
[52] 488 493 489 486  10 516 530 496   7 529 533 586   5 611  13   5 616
[69] 625 596 674 653 646 204 669  20   7 683  20   6 691  34  55 571  41
[86]  50 117
> exitus1
 [1] 2 0 0 1 0 2 0 0 0 1 0 0 0 0 0 0 0 0 0 0 2 2 0 0 0 1 0 0 1 0 2 0 1 1
[35] 0 1 1 2 0 0 1 0 1 0 2 1 1 0 0 0 1 0 0 0 0 1 0 0 0 1 0 0 0 1 0 1 1 0
[69] 0 0 0 0 0 0 0 2 2 0 1 1 0 0 0 0 1 1 2
Levels: 0 1 2
> gfr
 [1] 2 4 2 1 4 3 3 2 3 2 2 4 3 2 4 4 3 3 2 1 2 3 3 1 3 2 3 2 1 2 2 2 1 1
[35] 2 3 1 3 4 3 3 1 2 2 1 3 2 3 2 1 2 3 3 1 1 1 2 1 2 4 3 3 1 1 1 4 3 1
[69] 4 3 3 3 3 1 3 3 1 2 2 1 4 3 2 1 3 3 3
Levels: 1 2 3 4
> tdm
 [1] 1 1 0 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 0 1 0 0 1 1
[35] 0 1 1 1 1 1 0 1 1 1 0 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
[69] 1 1 1 1 0 1 1 1 1 0 1 1 1 1 1 1 1 1 1
Levels: 0 1

fg1Hi,您是否可以包括来自的功能,例如最小数据集?(您在数据集中包含了一个片段,但仍然不足以重现代码。)嗨!我刚刚编辑了文本,希望这能有所帮助。非常感谢。看起来您遗漏了该页面中可能值得重读的几个部分,特别是“生成最小数据集”部分和“复制数据”部分。
> str(fu)
 num [1:87] 15 27 3 10 16 13 65 64 5 29 ...
> str(tdm)
 Factor w/ 2 levels "0","1": 2 2 1 2 2 2 2 2 2 2 ...
> str(gfr)
 Factor w/ 4 levels "1","2","3","4": 2 4 2 1 4 3 3 2 3 2 ...
> str(exitus1)
 Factor w/ 3 levels "0","1","2": 3 1 1 2 1 3 1 1 1 2 ...

> fu
[1]  15  27   3  10  16  13  65  64   5  29 103  72 121 162  22         139 149
[18] 144  26  38  24  19 196 212 215  12 293  48  22 284   2  31   2   5
[35] 319   4  11  15 305 366   5 388  32 386  34  15   4 436 431 436  12
[52] 488 493 489 486  10 516 530 496   7 529 533 586   5 611  13   5 616
[69] 625 596 674 653 646 204 669  20   7 683  20   6 691  34  55 571  41
[86]  50 117
> exitus1
 [1] 2 0 0 1 0 2 0 0 0 1 0 0 0 0 0 0 0 0 0 0 2 2 0 0 0 1 0 0 1 0 2 0 1 1
[35] 0 1 1 2 0 0 1 0 1 0 2 1 1 0 0 0 1 0 0 0 0 1 0 0 0 1 0 0 0 1 0 1 1 0
[69] 0 0 0 0 0 0 0 2 2 0 1 1 0 0 0 0 1 1 2
Levels: 0 1 2
> gfr
 [1] 2 4 2 1 4 3 3 2 3 2 2 4 3 2 4 4 3 3 2 1 2 3 3 1 3 2 3 2 1 2 2 2 1 1
[35] 2 3 1 3 4 3 3 1 2 2 1 3 2 3 2 1 2 3 3 1 1 1 2 1 2 4 3 3 1 1 1 4 3 1
[69] 4 3 3 3 3 1 3 3 1 2 2 1 4 3 2 1 3 3 3
Levels: 1 2 3 4
> tdm
 [1] 1 1 0 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 0 1 0 0 1 1
[35] 0 1 1 1 1 1 0 1 1 1 0 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
[69] 1 1 1 1 0 1 1 1 1 0 1 1 1 1 1 1 1 1 1
Levels: 0 1
fg1<- crr(fu, exitus1, gfr, failcode=1, cencode=0)
summary(fg1)