How to force R to use a specified factor level as reference in a regression?

by

See the relevel() function. Here is an example:

set.seed(123)
x <- rnorm(100)
DF <- data.frame(x = x,
                 y = 4 + (1.5*x) + rnorm(100, sd = 2),
                 b = gl(5, 20))
head(DF)
str(DF)

m1 <- lm(y ~ x + b, data = DF)
summary(m1)

Now alter the factor b in DF by use of the relevel() function:

DF <- within(DF, b <- relevel(b, ref = 3))
m2 <- lm(y ~ x + b, data = DF)
summary(m2)

The models have estimated different reference levels.

> coef(m1)
(Intercept)           x          b2          b3          b4          b5 
  3.2903239   1.4358520   0.6296896   0.3698343   1.0357633   0.4666219 
> coef(m2)
(Intercept)           x          b1          b2          b4          b5 
 3.66015826  1.43585196 -0.36983433  0.25985529  0.66592898  0.09678759

More Related Contents:

  1. Linear Regression and group by in R
  2. Fitting a linear model with multiple LHS
  3. `lm` summary not display all factor levels
  4. How does predict.lm() compute confidence interval and prediction interval?
  5. Adding a regression line on a ggplot
  6. Fast pairwise simple linear regression between variables in a data frame
  7. Linear Regression with a known fixed intercept in R
  8. lme4::lmer reports “fixed-effect model matrix is rank deficient”, do I need a fix and how to?
  9. R: lm() result differs when using `weights` argument and when using manually reweighted data
  10. How `poly()` generates orthogonal polynomials? How to understand the “coefs” returned?
  11. predict.lm() with an unknown factor level in test data
  12. Regression (logistic) in R: Finding x value (predictor) for a particular y value (outcome)
  13. Messy plot when plotting predictions of a polynomial regression using lm() in R
  14. how to debug “factor has new levels” error for linear model and prediction [duplicate]
  15. Add regression line equation and R^2 on graph
  16. Dummify character column and find unique values [duplicate]
  17. Extract regression coefficient values
  18. Is there a better alternative than string manipulation to programmatically build formulas?
  19. What does the capital letter “I” in R linear regression formula mean?
  20. Creating dummy variables in R data.table
  21. read.csv doesn’t seem to detect factors in R 4.0.0
  22. R error “sum not meaningful for factors”
  23. Make Frequency Histogram for Factor Variables
  24. Add extra level to factors in dataframe
  25. predict.lm() in a loop. warning: prediction from a rank-deficient fit may be misleading
  26. One-Hot Encoding in [R] | Categorical to Dummy Variables [duplicate]
  27. lm function in R does not give coefficients for all factor levels in categorical data
  28. Ordering of points in R lines plot
  29. Add Regression Plane to 3d Scatter Plot in Plotly
  30. Double clustered standard errors for panel data

Leave a Comment