Existem várias formas de rodar várias regressões por categoria no R. Vou te mostrar como fazer com as funções base do R e com o dplyr
. Como exemplo, vamos usar a base de dados mtcars
.
Suponha que você queira rodar a regressão mpg ~ disp + hp
para cada nível da variável cyl
de mtcars
(são 3 categorias).
Primeiramente, você pode usar a função split()
para construir uma lista com três data.frames
diferentes, um para cada categoria:
data.frame_por_categoria <- split(mtcars, mtcars$cyl)
Agora, basta usar lapply()
para aplicar a regressão em cada data.frame
:
modelos <- lapply(data.frame_por_categoria, function(x) lm(mpg ~ disp + hp, data = x))
O resultado, modelos
é uma lista com as três regressões. Para acessar o primeiro modelo:
modelos[[1]]
Call:
lm(formula = mpg ~ disp + hp, data = x)
Coefficients:
(Intercept) disp hp
43.04006 -0.11954 -0.04609
Também é possível fazer a mesma coisa com o pacote dplyr
.
Você manda agrupar pela categoria e depois utiliza função do()
para rodar a regressão, colocando um ponto .
onde o data.frame
precisaria entrar:
library(dplyr)
resultado <- mtcars %>% group_by(cyl) %>% do(modelo = lm(mpg ~ disp + hp, data = .))
O resultado
da operação é um data.frame
com uma coluna chamada modelo, e cada elemento dessa coluna é a regressão. Para acessar o primeiro modelo:
resultado$modelo[[1]]
Call:
lm(formula = mpg ~ disp + hp, data = .)
Coefficients:
(Intercept) disp hp
43.04006 -0.11954 -0.04609
X
sãoX1
,X2
,X3
eX4
ou1
,2
,3
e4
?