• Exemplo: Extrair coeficientes de regressão de lm() em R

Pode utilizar os seguintes métodos para extrair os coeficientes de regressão da função lm() no R:

Método 1: Extrair apenas os coeficientes de regressão

 model$coefficients 

Método 2: Extrair coeficientes de regressão com erro padrão, estatística T, & valores P

 summary(model)$coefficients 

O exemplo seguinte mostra como utilizar estes métodos na prática.

Exemplo: Extrair coeficientes de regressão de lm() em R

Suponha que ajustamos o seguinte modelo de regressão linear múltipla em R:

 #criar quadro de dados df <- data. frame (rating=c(67, 75, 79, 85, 90, 96, 97), points=c(8, 12, 16, 15, 22, 28, 24), assists=c(4, 6, 6, 5, 3, 8, 7), rebounds=c(1, 4, 3, 3, 2, 6, 7)) #ajustar modelo de regressão linear múltipla model <- lm(rating ~ points + assists + rebounds, data=df) 

Podemos utilizar o resumo() para visualizar o resumo completo do modelo de regressão:

 #ver resumo do modelo resumo(modelo) Chamada: lm(fórmula = classificação ~ pontos + assistências + ressaltos, dados = df) Resíduos: 1 2 3 4 5 6 7 -1.5902 -1.7181 0.2413 4.8597 -1.0201 -0.6082 -0.1644 Coeficientes: Estimativa Erro Std. valor t Pr(>Códigos dos sinais: 0 '***' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 ' ' 1 Erro-padrão residual: 3,193 em 3 graus de liberdade Quadrado R múltiplo: 0,9589, Quadrado R ajustado: 0,9179 Estatística F: 23,35 em 3 e 3 DF, valor p: 0,01396 

Para visualizar apenas os coeficientes de regressão, podemos utilizar model$coefficients da seguinte forma:

 #Ver apenas os coeficientes de regressão do modelo model$coefficients (Interceção) pontos assistências ressaltos 66.435519 1.215203 -2.596789 2.820224 

Podemos utilizar estes coeficientes para escrever a seguinte equação de regressão ajustada:

Classificação = 66,43551 + 1,21520(pontos) - 2,59678(assistências) + 2,82022(ressaltos)

Para visualizar os coeficientes de regressão juntamente com os seus erros padrão, estatísticas t e valores p, podemos utilizar summary(model)$coefficients da seguinte forma:

 #ver coeficientes de regressão com erros padrão, estatística t e valores p summary(model)$coefficients Estimate Std. Error t value Pr(> 

Também podemos aceder a valores específicos nesta saída.

Por exemplo, podemos utilizar o seguinte código para aceder ao ficheiro valor de p para o pontos variável:

 #ver o valor p para a variável pontos summary(model)$coefficients[" points ", " Pr(> 

Ou podemos utilizar o seguinte código para aceder ao valor p de cada um dos coeficientes de regressão:

 #visualizar o valor p para todas as variáveis summary(model)$coefficients[, " Pr(> 

Os valores de p são apresentados para cada coeficiente de regressão no modelo.

É possível utilizar uma sintaxe semelhante para aceder a qualquer um dos valores na saída da regressão.