CBCA без опции

Question

мы пытаемся рассчитать CBCA с параметром «Нет». Мы работали с пакетом ChoiceModelR, следуя книге Элеи МакДоннел «R для маркетинговых исследований и аналитики (используйте R!)». Мы нашли там функцию «pred.hb.mnl».

predict.hb.mnl <- function(betadraws, data) {
  data.model <- model.matrix(~ Scan + Eco + price, data = data)#our formula
  data.model <- data.model[,-1]
  nresp <- dim(betadraws)[1]
  ndraws <- dim(bayesLukas$betadraw)[3]
  shares <- array(dim=c(nresp, nrow(data), ndraws))
  for (d in 1:ndraws) {
    for (i in 1:nresp) {
      utility <- data.model %*% betadraws[i,,d]
      shares[i,,d] = exp(utility) / sum(exp(utility))
    }
  }
  shares.agg <- apply(shares, 2:3, mean)
  cbind(share = apply(shares.agg, 1, mean),
        pct = t(apply(shares.agg, 1, quantile, probs = c(0.05, 0.95))), data)
}

Эта функция вычисляет бета-версии для каждого модуля панели, но не учитывает параметр none, который есть в наших данных. Есть ли возможность дополнить указанную выше функцию или есть ли у кого-нибудь другие предложения, как мы можем ее добавить. Большое спасибо!