Home Â» The Difference Between glm and lm in R

# The Difference Between glm and lm in R

The programming language R offers the following functions for fitting linear models:

1. lm â€“ Used to fit linear models

This function uses the following syntax:

lm(formula, data, â€¦)

where:

• formula: The formula for the linear model (e.g. y ~ x1 + x2)
• data: The name of the data frame that contains the data

2. glm â€“ Used to fit generalized linear models

This function uses the following syntax:

glm(formula, family=gaussian, data, â€¦)

where:

• formula: The formula for the linear model (e.g. y ~ x1 + x2)
• family: The statistical family to use to fit the model. Default is gaussian but other options include binomial, Gamma, and poisson among others.
• data: The name of the data frame that contains the data

Note that the only difference between these two functions is the family argument included in the glm() function.

If you use lm() or glm() to fit a linear regression model,Â they will produce the exact same results.

However, the glm() function can also be used to fit more complex models like:

The following examples show how to use the lm() function and glm() function in practice.

### Example of Using the lm() Function

The following code shows how to fit a linear regression model using the lm() function:

```#fit multiple linear regression model
model #view model summary
summary(model)

Call:
lm(formula = mpg ~ disp + hp, data = mtcars)

Residuals:
Min      1Q  Median      3Q     Max
-4.7945 -2.3036 -0.8246  1.8582  6.9363

Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 30.735904   1.331566  23.083  ```

### Examples of Using the glm() Function

The following code shows how to fit the exact same linear regression model using the glm() function:

```#fit multiple linear regression model
model #view model summary
summary(model)

Call:
glm(formula = mpg ~ disp + hp, data = mtcars)

Deviance Residuals:
Min       1Q   Median       3Q      Max
-4.7945  -2.3036  -0.8246   1.8582   6.9363

Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 30.735904   1.331566  23.083  ```

Notice that the coefficient estimates and standard errors of the coefficient estimates are the exact same as those produced by the lm() function.

Note that we can also use the glm() function to fit a logistic regression model by specifying family=binomial as follows:

```#fit logistic regression model
model #view model summary
summary(model)

Call:
glm(formula = am ~ disp + hp, family = binomial, data = mtcars)

Deviance Residuals:
Min       1Q   Median       3Q      Max
-1.9665  -0.3090  -0.0017   0.3934   1.3682

Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept)  1.40342    1.36757   1.026   0.3048
disp        -0.09518    0.04800  -1.983   0.0474 *
hp           0.12170    0.06777   1.796   0.0725 .
---
Signif. codes:  0 â€˜***â€™ 0.001 â€˜**â€™ 0.01 â€˜*â€™ 0.05 â€˜.â€™ 0.1 â€˜ â€™ 1

(Dispersion parameter for binomial family taken to be 1)

Null deviance: 43.230  on 31  degrees of freedom
Residual deviance: 16.713  on 29  degrees of freedom
AIC: 22.713

Number of Fisher Scoring iterations: 8
```

We can also use the glm() function to fit a Poisson regression model by specifying family=poisson as follows:

```#fit Poisson regression model
model #view model summary
summary(model)

Call:
glm(formula = am ~ disp + hp, family = poisson, data = mtcars)

Deviance Residuals:
Min       1Q   Median       3Q      Max
-1.1266  -0.4629  -0.2453   0.1797   1.5428

Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept)  0.214255   0.593463   0.361  0.71808
disp        -0.018915   0.007072  -2.674  0.00749 **
hp           0.016522   0.007163   2.307  0.02107 *
---
Signif. codes:  0 â€˜***â€™ 0.001 â€˜**â€™ 0.01 â€˜*â€™ 0.05 â€˜.â€™ 0.1 â€˜ â€™ 1

(Dispersion parameter for poisson family taken to be 1)

Null deviance: 23.420  on 31  degrees of freedom
Residual deviance: 10.526  on 29  degrees of freedom
AIC: 42.526

Number of Fisher Scoring iterations: 6
```