Zróżnicowanie preferencji dla innych cech

Modelowanie wyborów w marketingu w R

Elea McDonnell Feit

Assistant Professor of Marketing, Drexel University

Inny sposób kodowania czynników

Kodowanie zero-jedynkowe (to, czego używaliśmy do tej pory)

  seat4 seat5
2     0     0
4     1     0
5     0     1

Kodowanie efektów (lepsze dla modeli hierarchicznych)

   seat2  seat4
2      1      0
4      0      1
5     -1     -1

Zmiana kodowania czynnika

contrasts(sportscar$seat) <- contr.sum(levels(sportscar$seat))

dimnames(contrasts(sportscar$seat))[[2]] <- levels(sportscar$seat)[1:2]

contrasts(sportscar$seat)

Zróżnicowanie wszystkich współczynników

my_rpar <- c("n", "n", "n", "n", "n")

m3 <- mlogit(choice ~ 0 + seat + trans + convert + price,
             data = sportscar)
names(my_rpar) <- names(m3$coefficients)

my_rpar

      seat2       seat4 transmanual  convertyes       price 
        "n"         "n"         "n"         "n"         "n"

m8 <- mlogit(choice ~ 0 + seat + trans + convert + price, 
             data = sportscar, 
             panel = TRUE, rpar = my_rpar)

Parametry modelu hierarchicznego

m8 <- mlogit(choice ~ 0 + seat + trans + convert + price, 
              data = sportscar, 
             panel = TRUE, rpar = my_rpar)
plot(m8, par = c("seat2", "seat4"))

rozkłady parametrów kodowanych efektami

Współczynnik dla poziomu bazowego

m8$coef[1:2]

     seat2      seat4 
-0.1686077 -0.1912170

-sum(m8$coef[1:2])

[1] 0.3598248

Przetestujmy to na danych `chocolate`!

Modelowanie wyborów w marketingu w R