R:提取lm()属性
问题描述:
我想从摘要(lm())中提取fstatistic值。到目前为止,我发现的唯一方法是R:提取lm()属性
summary(lm(this_vector ~ that_vector))["fstatistic"][[1]][1]
是否有一种较不详细的方式来获取该单元格值?这个问题有点迂腐,但我认为答案可能会提供一些关于如何使用R列表的有趣信息。
答
尝试以下任一:
summary(lm(this_vector ~ that_vector))$fstatistic[1]
summary(lm(this_vector ~ that_vector))[["fstatistic"]][1]
["fstatistic"]
返回与具有匹配什么是单括号内名称元素的列表,所以你需要[[1]]
拿到的第一个元素。双括号本身返回元素,如使用$
表示法。
答
包broom
包含从各种模型(其中lm
)方便地提取估计和模型拟合统计的功能。为了提取F统计量,使用broom
的glance()
函数。看到这里提到的例子https://github.com/dgrtwo/broom。
在你的情况,这将是
glance(lm(this_vector ~ that_vector))$statistic