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的glance()函数。看到这里提到的例子https://github.com/dgrtwo/broom。

在你的情况，这将是

glance(lm(this_vector ~ that_vector))$statistic