Question

我想知道是否只能显示R Summary输出的特定组件。从下面的输出中，我只想显示F统计量（或仅下面的行-F-统计量：2和10 DF上的0.834，p值：0.4624）。我实际上正在处理估算的数据并尝试从一系列摘要输出中提取此值。

谢谢！

 > hanes=with(nhanes, lm(bmi~hyp+chl))
 > summary(hanes)
 Call:
 lm(formula = bmi ~ hyp + chl)

 Residuals:
 Min     1Q Median     3Q    Max 
-6.425 -3.296  0.035  3.123  7.632 

 Coefficients:
 Estimate Std. Error t value Pr(>|t|)   
 (Intercept) 20.10564    5.70716   3.523  0.00551 **
 hyp         -0.68459    3.39596  -0.202  0.84428   
 chl          0.03783    0.03054   1.239  0.24370   
 ---
 Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

 Residual standard error: 4.66 on 10 degrees of freedom
 (12 observations deleted due to missingness)
 Multiple R-squared:  0.143,    Adjusted R-squared:  -0.02846 
 F-statistic: 0.834 on 2 and 10 DF,  p-value: 0.4624

使用RAN中的NHANES2数据更新可重复示例：

  > imp10=mice(nhanes2, m=10, seed=11111)

  > hanes=with(imp10, lm(bmi~hyp+chl))

  > summary(hanes)

上面的

产生了10个插补中每个插补的摘要：

  ## summary of imputation 1 :

  Call:
  lm(formula = bmi ~ hyp + chl)

  Residuals:
  Min      1Q  Median      3Q     Max 
  -7.4597 -3.1938  0.1403  2.8153  7.0753 

  Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
  Intercept) 19.30988    4.47439   4.316 0.000279 ***
  hyp2        -0.38202    2.41568  -0.158 0.875789    
  chl          0.04265    0.02398   1.779 0.089101 .  
  ---
  Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

  Residual standard error: 4.332 on 22 degrees of freedom
  Multiple R-squared:  0.1425,  Adjusted R-squared:  0.06456 
  F-statistic: 1.828 on 2 and 22 DF,  p-value: 0.1843

fs = summary（hanes）$ fstatistic 代码仍会生成完整的摘要。再次感谢你！

Answer 1

如下我的评论：

要检索您可能使用的 F-statistic ：

fs <- summary(mode)$fstatistic

与P值相关联：

pva <- anova(mode)$"Pr(>F)"[1]

然后与sprintf等组合在一起或者只使用保存字符向量中所有输出的capture.output，其中每个元素对应一行;你只需要感兴趣的行，在这种情况下是18：

capture.output(summary(mode))[18]
[1] "F-statistic: 30.19 on 1 and 30 DF,  p-value: 5.766e-06"

我在R。

中使用了mtcars数据集的拟合模型

在R摘要中显示对象的子集

1 个答案: