Scribd
Upload
130 views

Fitting & Interpreting Linear Models in Rinear Models in R

This document discusses linear regression models in R. It begins by explaining how to easily build linear models using the lm function. Key steps include fitting a model to example height data of children and parents using the formula child ~ parent. The summary function is then used to interpret the lm results, which provides important metrics like residuals, coefficient estimates, p-values and R-squared to evaluate the quality and significance of the fitted model. Categorical variables can also be included, with R automatically handling dummy variable coding.

Uploaded by

ReaderRat
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
130 views

Fitting & Interpreting Linear Models in Rinear Models in R

This document discusses linear regression models in R. It begins by explaining how to easily build linear models using the lm function. Key steps include fitting a model to example height data of children and parents using the formula child ~ parent. The summary function is then used to interpret the lm results, which provides important metrics like residuals, coefficient estimates, p-values and R-squared to evaluate the quality and significance of the fitted model. Categorical variables can also be included, with R automatically handling dummy variable coding.

Uploaded by

ReaderRat
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 8

Fitting & Interpreting Linear Models in R

R makes it easy to fit a linear model to your data. The hard part is knowing whether the model
you've built is worth keeping and, if so, figuring out what to do next.
This is a post about linear models in R, how to interpret lm results, and common rules of thumb
to help side-step the most common mistakes.

Building a linear model in R


makes building linear models really easy. Things like dummy variables, categorical features,
interactions, and multiple regression all come very naturally. The centerpiece for linear
regression in R is the lm function.
R

comes with base R, so you don't have to install any packages or import anything special. The
documentation for lm is very extensive, so if you have any questions about using it, just type ?lm
into the R console.
lm

Introduction to lm
For our example linear model, I'm going to use data from the original, or at least one of the
earliest, linear regression models. The dataset consists of heights of children and their parents.
The origin of the term "regression" stems from a 19th century statistician's observation that
children's heights tended to "regress" towards the population mean in relation to their parent's
heights.

galton <- read.csv("http://blog.yhathq.com/static/misc/galton.csv",


header=TRUE, stringsAsFactors=FALSE)
summary(galton)
# child parent
# Min. :61.7 Min. :64.0
# 1st Qu.:66.2 1st Qu.:67.5

#
#
#
#

Median :68.2 Median :68.5


Mean :68.1 Mean :68.3
3rd Qu.:70.2 3rd Qu.:69.5
Max. :73.7 Max. :73.0

head(galton)
# child parent
#1 61.7 70.5
#2 61.7 68.5
#3 61.7 65.5
#4 61.7 64.5
#5 61.7 64.0
#6 62.2 67.5
view raw intro_to_lm_summary.R hosted with by GitHub

Fit the model to the data by creating a formula and passing it to the lm function. In our case we
want to use the parent's height to predict the child's height, so we make the formula (child ~
parent). In other words, we're representing the relationship between parents' heights (X) and
children's heights (y).
We then set the data being used to galton so lm knows what data frame to associate "child" and
"parent" to.
fit <- lm(child ~ parent, data=galton)
fit
#Call:
#lm(formula = child ~ parent, data = galton)
#
#Coefficients:
#(Intercept) parent
# 23.942 0.646

view raw model_lm_summary.R hosted with by GitHub

NOTE: Formulas in R take the form (y ~ x). To add more predictor variables,

just use the + sign. i.e. (y ~ x + z).


Calling summary
We fit a model to our data. That's great! But the important question is, is it
any good?
There are lots of ways to evaluate model fit. lm consolidates some of the
most popular ways into the summary function. You can invoke the summary
function on any model you've fit with lm and get some metrics indicating the
quality of the fit.

summary(fit)
#Call:
#lm(formula = child ~ parent, data = galton)
#
#Residuals:
# Min 1Q Median 3Q Max
#-7.805 -1.366 0.049 1.634 5.926
#
#Coefficients:
# Estimate Std. Error t value Pr(>|t|)
#(Intercept) 23.9415 2.8109 8.52 <2e-16 ***
#parent 0.6463 0.0411 15.71 <2e-16 ***
#--#Signif. codes: 0 *** 0.001 ** 0.01 * 0.05 . 0.1 1
#
#Residual standard error: 2.24 on 926 degrees of freedom
#Multiple R-squared: 0.21, Adjusted R-squared: 0.21
#F-statistic: 247 on 1 and 926 DF, p-value: <2e-16

view raw summary_lm_summary.R hosted with by GitHub

So if you're like I was at first, your reaction was probably something like
"Whoa this is cool...what does it mean?"

Interpreting the output

#
Name
1 Residuals

Description
The residuals are the difference between the actual values
of the variable you're predicting and predicted values from
your regression--y - . For most regressions you want your
residuals to look like a normal distribution when plotted. If
our residuals are normally distributed, this indicates the
mean of the difference between our predictions and the
actual values is close to 0 (good) and that when we miss,
we're missing both short and long of the actual value, and
the likelihood of a miss being far from the actual value gets
smaller as the distance from the actual value gets larger.
Think of it like a dartboard. A good model is going to hit the
bullseye some of the time (but not everytime). When it
doesn't hit the bullseye, it's missing in all of the other
buckets evenly (i.e. not just missing in the 16 bin) and it
also misses closer to the bullseye as opposed to on the

outer edges of the dartboard.


The stars are shorthand for significance levels, with the
number of asterisks displayed according to the p-value
Significance computed. *** for high significance and * for low
2
Stars
significance. In this case, *** indicates that it's unlikely that
no relationship exists b/w heights of parents and heights of
their children.
The estimated coefficient is the value of slope calculated by
the regression. It might seem a little confusing that the
Intercept also has a value, but just think of it as a slope that
Estimated
3
is always multiplied by 1. This number will obviously vary
Coeffecient
based on the magnitude of the variable you're inputting into
the regression, but it's always good to spot check this
number to make sure it seems reasonable.
Measure of the variability in the estimate for the coefficient.
Lower means better but this number is relative to the value
of the coefficient. As a rule of thumb, you'd like this value to
Standard
be at least an order of magnitude less than the coefficient
Error of the
4
estimate.
Coefficient
Estimate
In our example, the std error or the parent variable is 0.04
which is 16x less than the estimate of the coefficient (or 1.6
orders of magnitude greater).
Score that measures whether or not the coefficient for this
t-value of the
variable is meaningful for the model. You probably won't use
5 Coefficient
this value itself, but know that it is used to calculate the pEstimate
value and the significance levels.
Probability the variable is NOT relevant. You want this
number to be as small as possible. If the number is really
Variable p6
small, R will display it in scientific notation. In or example
value
2e-16 means that the odds that parent is meaningless is
about 15000000000000000
The more punctuation there is next to your variables, the
better.
Significance
7
Legend
Blank=bad, Dots=pretty good, Stars=good, More
Stars=very good
8 Residual Std The Residual Std Error is just the standard deviation of your
Error /
residuals. You'd like this number to be proportional to the
Degrees of quantiles of the residuals in #1. For a normal distribution,
Freedom
the 1st and 3rd quantiles should be 1.5 +/- the std error.
The Degrees of Freedom is the difference between the
number of observations included in your training sample

and the number of variables used in your model (intercept


counts as a variable).
Metric for evaluating the goodness of fit of your model.
Higher is better with 1 being the best. Corresponds with the
amount of variability in what you're predicting that is
9 R-squared
explained by the model. In this instance, ~21% of the cause
for a child's height is due to the height their parent.
WARNING: While a high R-squared indicates good
correlation, correlation does not always imply causation.
Performs an F-test on the model. This takes the parameters
of our model (in our case we only have 1) and compares it
to a model that has fewer parameters. In theory the model
with more parameters should fit better. If the model with
more parameters (your model) doesn't perform better than
F-statistic & the model with fewer parameters, the F-test will have a high
1
resulting p- p-value (probability NOT significant boost). If the model with
0
value
more parameters is better than the model with fewer
parameters, you will have a lower p-value.
The DF, or degrees of freedom, pertains to how many
variables are in the model. In our case there is one variable
so there is one degree of freedom.
Categorical Variables
People often wonder how they can include categorical variables in their
regression models. With R this is extremely easy. Just include the categorical
variable in your regression formula and R will take care of the rest. R calls
categorical variables factors. A factor has a set of levels, or possible values.
These levels will show up as variables in the model summary.
Dummy Variable Trap

One very important thing to note is that one of your levels will not appear in
the output. This is because when fitting a regression with a categorical
variable, one option must be left out to avoid overfitting the model. This is
often referred to as the dummy variable trap. In our model, Africa is left out
of the summary but it is still accounted for in the model.

library(reshape2)

phones <- melt(WorldPhones)


names(phones) <- c("year", "continent", "n_phones")
head(phones)
fit <- lm(n_phones ~ year + continent, data=phones)

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy