r - cr.setup 以 rms 为单位

Question

有人可以解释一下包中的cr.setup功能是如何rms工作的吗？我似乎无法弄清楚它如何重新映射初始数据以及这种重新映射对延续比率模型有何用处，并且帮助和示例没有那么有用。网上也找不到其他解释。

Answer 1

在他的优秀著作“回归建模策略”中，Harrell 在他关于序数逻辑回归的章节末尾有三页专门讨论了延续比率模型。cr.setup通过复制某些行并为各种比较创建层标记来支持“欺骗普通逻辑回归”的过程：Y >= 0; Y>=1, ... Y>=K-1并创建适当的响应变量来表示特定层的“结果”。看看他的 cr.setup 的第一个例子：

y <- c(NA, 10, 21, 32, 32)
> cr.setup(y)
$y
[1] NA  1  0  1  0  0  0  0

$cohort
[1] <NA>  all   all   y>=21 all   y>=21 all   y>=21
Levels: all y>=21

$subs
[1] 1 2 3 3 4 4 5 5

$reps
[1] 1 1 2 2 2

如果有 3 个级别的非 N Y，则只有 2 个级别的新结果。y 向量是新结果。向量元素是原始数据的subs索引。reps向量告诉软件需要多少次复制。您可以通过以下示例了解如何在实践中使用它：

combinations <- expand.grid(cohort=levels(cohort), sex=levels(sex))
combinations