включить столбец с 0/1 до пересечения в линейной фиксированной модели

Question

Давайте предположим, что у меня есть кадр данных, который выглядит следующим образом:

|--------------|-----------|-----------|-------|-----|-----------|
|coffeeshop_id | americano | capuccino | latte | tea | mean_price|
|--------------|-----------|-----------|-------|-----|-----------|
|15476         |1          | 0         | 0     | 1   | 2.35      |
|--------------|-----------|-----------|-------|-----|-----------|
|13466         |0          | 1         | 0     | 1   | 2.45      |
|--------------|-----------|-----------|-------|-----|-----------|

Я хочу построить линейную регрессию и сделать американо эталоном в перехвате. Как я могу это сделать?

Answer 1

Вы можете использовать relevel, чтобы установить новый опорный уровень для подгонки модели.

# create some data
set.seed(1)
x <- factor(sample(letters[1:3], size=20, replace=T))
y <- rnorm(20)
my.df <- data.frame(x=x, y=y)

# model fitted with "a" as the reference level
summary(lm(y~x, data=my.df))

# set "b" as the reference level
my.df$x <- relevel(my.df$x, ref='b')
summary(lm(y~x, data=my.df))