如何提取p值

Question

我已经尝试了一段时间，在我的数据集的第一个观察值和199个观察值中的其余观察值之间进行回归。我使用了lapply函数，回归结果以列表的形式存储在环境中。我的目标是只获取p_values列表作为数据帧，并确定有多少观测值小于0.05。任何帮助都将不胜感激！

## Here are the code I am using right now.
myre1 <- apply(2:ncol(muscle), function(x) lm(muscle[,1] ~ muscle[,x], data = muscle))
myre2 <- lapply(muscle[,-1], function(x) lm(muscle$GIR ~ x))

## To extract the coefficient
myre3 <- lapply(2:ncol(muscle), function(x) coefficients(lm(muscle[,1] ~ muscle[,x], data = muscle)))
myre4 <- lapply(muscle[,-1], function(x) coefficients(lm(muscle$GIR ~ x)))

Answer 1

尝尝这个。您将获得系数和p.values

sum1 <- lapply(mre1, function(x) summary(x)$coefficients)
sum2 <- lapply(mre2, function(x) summary(x)$coefficients)

您可能需要将sum1和sum2转换为数据帧。

编辑：

仅提取p值

sum1 <- lapply(mre1, function(x) summary(x)$coefficients[10:12])
sum2 <- lapply(mre2, function(x) summary(x)$coefficients[7:8])