1
R实现部分确定系数
有没有人有建议或一揽子计划来计算部分确定系数? 可以将部分确定系数定义为无法在简化模型中解释的变化百分比,但可以由完整模型中指定的预测变量解释。该系数用于深入了解一个或多个其他预测变量在更完全指定的回归模型中是否有用。 在估计您的两个模型并为它们生成ANOVA表后,部分r ^ 2的计算相对简单。局部r ^ 2的计算为: (减少SSE-减少SSEfull)/减少SSE 我已经编写了这个相对简单的函数,可以为多元线性回归模型计算该函数。我不熟悉R中的其他模型结构,其中该功能可能无法很好地发挥作用: partialR2 <- function(model.full, model.reduced){ anova.full <- anova(model.full) anova.reduced <- anova(model.reduced) sse.full <- tail(anova.full$"Sum Sq", 1) sse.reduced <- tail(anova.reduced$"Sum Sq", 1) pR2 <- (sse.reduced - sse.full) / sse.reduced return(pR2) } 对于完成该任务的更健壮的功能和/或以上代码的更有效实现的任何建议或技巧将不胜感激。
9
r
regression
anova