inst/doc/mxl_models.R

In logitr: Logit Models w/Preference & WTP Space Utility Parameterizations

## ----setup, include=FALSE, message=FALSE, warning=FALSE-----------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  warning = FALSE,
  message = FALSE,
  fig.retina = 3,
  comment = "#>"
)

library(logitr)

# Read in results from already estimated models so that the
# examples aren't actually run when building this page,
# otherwise it'll take much longer to build
mxl_pref <- readRDS(here::here('inst', 'extdata', 'mxl_pref.Rds'))
mxl_wtp  <- readRDS(here::here('inst', 'extdata', 'mxl_wtp.Rds'))
mxl_pref_cor <- readRDS(here::here('inst', 'extdata', 'mxl_pref_cor.Rds'))

## ----eval=FALSE---------------------------------------------------------------
# library("logitr")
# 
# set.seed(456)
# 
# mxl_pref <- logitr(
#   data     = yogurt,
#   outcome  = 'choice',
#   obsID    = 'obsID',
#   panelID  = 'id',
#   pars     = c('price', 'feat', 'brand'),
#   randPars = c(feat = 'n', brand = 'n'),
#   numMultiStarts = 10
# )

## -----------------------------------------------------------------------------
summary(mxl_pref)

## -----------------------------------------------------------------------------
wtp_mxl_pref <- wtp(mxl_pref, scalePar =  "price")
wtp_mxl_pref

## ----eval=FALSE---------------------------------------------------------------
# set.seed(6789)
# 
# mxl_wtp <- logitr(
#   data       = yogurt,
#   outcome    = 'choice',
#   obsID      = 'obsID',
#   panelID    = 'id',
#   pars       = c('feat', 'brand'),
#   scalePar   = 'price',
#   randPars   = c(feat = 'n', brand = 'n'),
#   numMultiStarts = 10,
#   startVals = wtp_mxl_pref$Estimate
# )

## -----------------------------------------------------------------------------
summary(mxl_wtp)

## -----------------------------------------------------------------------------
wtpCompare(mxl_pref, mxl_wtp, scalePar = 'price')

## ----eval=FALSE---------------------------------------------------------------
# library("logitr")
# 
# set.seed(456)
# 
# mxl_pref_cor <- logitr(
#   data     = yogurt,
#   outcome  = 'choice',
#   obsID    = 'obsID',
#   panelID  = 'id',
#   pars     = c('price', 'feat', 'brand'),
#   randPars = c(feat = 'n', brand = 'n'),
#   numMultiStarts = 10,
#   correlation = TRUE
# )

## -----------------------------------------------------------------------------
summary(mxl_pref_cor)