Создание фрейма данных с использованием извлеченных коэффициентов из списка

Я запускаю почти 40 логистических регрессий, используя один и тот же набор независимых переменных, но 40 разных зависимых переменных. Я хотел бы извлечь коэффициенты из 40 линий регрессии и создать фрейм данных для построения точечного графика. Я думаю, что лучше визуализировать это, чем помещать все линии регрессии в одну таблицу.

Я не мог понять, как извлечь влияние независимой переменной (один и тот же тип коэффициентов из каждого уравнения)

Признателен, если вы можете мне помочь!

Вот воспроизводимый пример.

set.seed(10) 
y <- matrix(rnorm(10000 * 14), ncol = 14) 
x <- matrix(rnorm(10000 * 2), ncol = 2) 
res <- lapply(1:14, function(i) lm(y[, i] ~ x))

r list extract coefficients
person redoksss    schedule 03.06.2014    source источник
comment
res <- apply(y, 2, function(z) lm(z~x)) будет немного эффективнее, так как apply был специально разработан для операций со столбцами.   -  person Señor O    schedule 04.06.2014

Ответы (1)


arrow_upward
1
arrow_downward

Попробуйте использовать do.call и rbind в своем списке res

> do.call(rbind, lapply(res, function(x) coef(x)))
        (Intercept)            x1            x2
 [1,]  1.440115e-03 -0.0198232209 -0.0005720764
 [2,] -2.227644e-02 -0.0134155339 -0.0092420757
 [3,]  3.535811e-03 -0.0284229117  0.0140198529
 [4,] -2.031279e-02  0.0032004789 -0.0036719760
 [5,] -1.127532e-02 -0.0004463859 -0.0116754425
 [6,] -1.369851e-02  0.0174797415  0.0112791379

Это матрица, используйте data.frame, чтобы получить data.frame

person Jilber Urbina    schedule 03.06.2014