# NOT RUN {
eq_d <- c("ASC1 * 1 + B11_dur * dur_1" , "ASC2 * 1 + B12_dur * dur_2",
"ASC3 * 1 + B13_dur * dur_3 + B20_cost * cost_3 + B53_parkman * PbAvl_3",
"ASC4 * 1 + B14_dur * dur_4 + B20_cost * cost_4 + B34_serv * servIdx_4 + B44_stop * stopUs1R1_4")
parl <- c(paste0("ASC", 1:4), paste0("B1", 1:4, "_dur"), "B20_cost", "B53_parkman", "B34_serv",
"B44_stop")
obj <- get_par(parl, eq_d)
ffor <- obj$cheqs0
res <- MNlogitf(ffor, separatenmm=FALSE, transform=FALSE)
ff <- f_create(res$formula, data="data", fixed=1)
# }
Run the code above in your browser using DataLab