Как извлечь p-значения

Question

Я некоторое время пытался регрессировать между первым наблюдением и остальными из 199 наблюдений моего набора данных. Я использовал функцию lapply, и результат регрессии сохраняется в виде списка в среде. Моя цель - получить только список значений p_value в качестве фрейма данных и определить, сколько наблюдений меньше 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])