Plotting partials of multiple regression in R (avPlots)

You can just use function avPlots() from ‘car’ package:

My multiple regression model follows:

mreg=lm(activity~rain+temperature)

That provided these results:

Call:
lm(formula = activity ~ rain + temperature)

Residuals:
Min        1Q    Median        3Q       Max
-0.016638 -0.001452  0.001483  0.003645  0.010533

Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.190e+00  1.916e-02 114.322 1.52e-15 ***
rain        9.569e-04  2.905e-05  32.945 1.08e-10 ***
temperature 2.273e-03  8.618e-04   2.637    0.027 *

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.00819 on 9 degrees of freedom
Multiple R-squared: 0.9923,    Adjusted R-squared: 0.9906
F-statistic: 578.5 on 2 and 9 DF,  p-value: 3.12e-10

than I just use:

library(car)  #a native library of R (you don’t need to install it 😀 )

avPlots(mreg)

That will kindly plot it as:

Partial plots from the multiple regression

Partial plots from the multiple regression

You can find a another example here

Thats it! 😉

Anúncios

Sobre Conrado

Meu interesse é no estudo do comportamento social e comunicação de lagartos. "Infófilo", entusiasta do software livre, da licença creative commons e do livre acesso da produção artística e intelectual.

Publicado em abril 25, 2013, em Dicas, Graphics e marcado como , . Adicione o link aos favoritos. Deixe um comentário.

Deixe um comentário

Preencha os seus dados abaixo ou clique em um ícone para log in:

Logotipo do WordPress.com

Você está comentando utilizando sua conta WordPress.com. Sair / Alterar )

Imagem do Twitter

Você está comentando utilizando sua conta Twitter. Sair / Alterar )

Foto do Facebook

Você está comentando utilizando sua conta Facebook. Sair / Alterar )

Foto do Google+

Você está comentando utilizando sua conta Google+. Sair / Alterar )

Conectando a %s

%d blogueiros gostam disto: