计算测定分项系数的R包?

计算测定分项系数的R包?,r,asbio,R,Asbio,有人知道在多元回归中计算部分R^2的R包吗?我已经尝试了来自包asbio的命令partial.R2,但它给出了错误消息,即使是来自所提供文档的示例 非常感谢。我发现包lmSupport中的命令lm.sumSquares通过偏相关和半偏相关提供。来自John Neter、Michael H Kutner、William Wasserman、Christopher J.Nachtsheim的“应用线性统计模型”的数据 第274页第7.4节: # body fat example from Neter

有人知道在多元回归中计算部分R^2的R包吗?我已经尝试了来自包asbio的命令partial.R2,但它给出了错误消息,即使是来自所提供文档的示例


非常感谢。

我发现包lmSupport中的命令lm.sumSquares通过偏相关和半偏相关提供。

来自John Neter、Michael H Kutner、William Wasserman、Christopher J.Nachtsheim的“应用线性统计模型”的数据

第274页第7.4节:

# body fat example from Neter et al.  via rhelp archives:

bf.dat <- read.table(text="x1 x2 x3 y
1 19.5 43.1 29.1 11.9
2 24.7 49.8 28.2 22.8
3 30.7 51.9 37.0 18.7
4 29.8 54.3 31.1 20.1
5 19.1 42.2 30.9 12.9
6 25.6 53.9 23.7 21.7
7 31.4 58.5 27.6 27.1
8 27.9 52.1 30.6 25.4
9 22.1 49.9 23.2 21.3
10 25.5 53.5 24.8 19.3
11 31.1 56.6 30.0 25.4
12 30.4 56.7 28.3 27.2
13 18.7 46.5 23.0 11.7
14 19.7 44.2 28.6 17.8
15 14.6 42.7 21.3 12.8
16 29.5 54.4 30.1 23.9
17 27.7 55.3 25.7 22.6
18 30.2 58.6 24.6 25.4
19 22.7 48.2 27.1 14.8
20 25.2 51.0 27.5 21.1 ", header=TRUE)

library(rms)  # will also load Hmisc
fit <- ols(y ~ x1 + x2, data=bf.dat) 

plt <- plot(anova(fit), what='partial R2')
plt
#         x2          x1 
#0.066955220 0.007010427 
#Neter等人通过rhelp档案提供的体脂示例:

bf.dat可能对您有用?我发现包lmSupport中的命令lm.sumSquares通过部分和半部分相关性提供。