R：提取lm（）属性

Question

我想从摘要（lm（））中提取fstatistic值。到目前为止，我找到的唯一方法是

summary(lm(this_vector ~ that_vector))["fstatistic"][[1]][1]

是否有更简洁的方法来获得该单元格值？问题有点迂腐，但我认为答案可能会提供一些有关如何使用R列表的有趣信息。

Answer 1

尝试其中任何一种：

summary(lm(this_vector ~ that_vector))$fstatistic[1]
summary(lm(this_vector ~ that_vector))[["fstatistic"]][1]

["fstatistic"]返回一个列表，其中的元素名称与单括号内的名称相匹配，因此您需要[[1]]来获取第一个元素。双括号返回元素本身，使用$表示法。

Answer 2

包broom包含从各种模型（其中lm）中方便地提取估计和模型拟合统计的功能。要提取F统计量，请使用broom的{​​{1}}函数。请参阅此处提到的示例https://github.com/dgrtwo/broom。

在你的情况下，那将是

glance()