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