Model pro čokoládu se všemi náhodnými koeficienty
Teď, když máme effects coding uložený v datech chocolate, můžeme přistoupit k sestavení modelu, kde jsou všechny koeficienty normálně rozdělené. K tomu potřebujeme vytvořit vektor rpar, který předáme funkci mlogit(). Je to trochu složitější, takže jsem kód připravil za tebe – spusť ho a podívej se, jak funguje. Pak napíšeš samotné volání funkce mlogit().
Toto cvičení je součástí kurzu
Modely výběru pro marketing v R
Pokyny k cvičení
- První dva vstupy jsou vzorec modelu
Selection ~ 0 + Brand + Type + Pricea datachocolate. - Dalším vstupem by mělo být
rpar = my_rpar, které funkcimlogit()říká, které koeficienty mají být normálně rozdělené. - Posledním vstupem by mělo být
panel = TRUE.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# create my_rpar vector
choc_m2 <- mlogit(Selection ~ 0 + Brand + Type + Price, data=chocolate)
my_rpar <- rep("n", length(choc_m2$coef))
names(my_rpar) <- names(choc_m2$coef)
my_rpar
# fit model with random coefficients
choc_m7 <- mlogit(___)