如何解决从 lm 的大列表中提取 p 值
我正在尝试从我创建的 lm 大列表中提取 p 值:
<div>
<button (click)="enlargeText()">A+</button>
<button (click)="minimizeText()">A-</button>
</div>
这给我留下了每个组(BayStation)的 24 个 lm 列表。我能够使用以下方法提取相交和斜率:
S_models <- dlply(S,"BayStation",function(df)
lm(Temp ~ Date,data = df))
但是,如果不单独或手动执行,我似乎无法弄清楚如何从这个大列表中获取 p 值的列表或数据框。
解决方法
尝试使用 sapply
/lapply
:
result <- sapply(S_models,function(x) summary(x)$coefficients[,4])
#With `lapply`
#result <- lapply(S_models,4])
,
我们可以使用map
library(purrr)
result <- map(S_models,~ summary(.x)$coefficients[,4])
,
以下仅使用基数 R 并使用函数 pf
获取 p 值。分位数和自由度在 summary
返回的列表中。
sapply(S_models,function(x){
ff <- summary(x)$fstatistic
pf(ff[1],df1 = ff[2],df2 = ff[3],lower.tail = FALSE)
})
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。