Model run by stephane.hess using Apollo 0.2.9 on R 4.0.5 for Darwin. www.ApolloChoiceModelling.com Model name : CNL Model description : CNL model with socio-demographics on mode choice SP data Model run at : 2023-05-10 20:20:02 Estimation method : bfgs Model diagnosis : successful convergence Optimisation diagnosis : Maximum found hessian properties : Negative definitive maximum eigenvalue : -4.440449 Number of individuals : 500 Number of rows in database : 7000 Number of modelled outcomes : 7000 Number of cores used : 1 Model without mixing LL(start) : -4830.94 LL at equal shares, LL(0) : -8196.02 LL at observed shares, LL(C) : -6706.94 LL(final) : -4742.25 Rho-squared vs equal shares : 0.4214 Adj.Rho-squared vs equal shares : 0.419 Rho-squared vs observed shares : 0.2929 Adj.Rho-squared vs observed shares : 0.2904 AIC : 9524.5 BIC : 9661.57 Estimated parameters : 20 Time taken (hh:mm:ss) : 00:00:39.6 pre-estimation : 00:00:0.75 estimation : 00:00:18.99 initial estimation : 00:00:17.52 estimation after rescaling : 00:00:1.47 post-estimation : 00:00:19.85 Iterations : 33 initial estimation : 30 estimation after rescaling : 3 Unconstrained optimisation. Estimates: Estimate s.e. t.rat.(0) t.rat(1) Rob.s.e. Rob.t.rat.(0) Rob.t.rat.(1) asc_car 0.000000 NA NA NA NA NA NA asc_bus 0.299289 0.445857 0.6713 -1.572 0.420733 0.7114 -1.665 asc_air -0.646373 0.292126 -2.2126 -5.636 0.285269 -2.2658 -5.771 asc_rail -1.563495 0.287744 -5.4336 -8.909 0.278678 -5.6104 -9.199 asc_bus_shift_female 0.273647 0.105276 2.5993 -6.899 0.111518 2.4538 -6.513 asc_air_shift_female 0.240335 0.076638 3.1360 -9.912 0.075900 3.1665 -10.009 asc_rail_shift_female 0.178849 0.067296 2.6577 -12.202 0.067973 2.6312 -12.081 b_tt_car -0.010391 7.1322e-04 -14.5688 -1416.657 7.0558e-04 -14.7266 -1432.004 b_tt_bus -0.016555 0.001241 -13.3375 -818.976 0.001199 -13.8093 -847.944 b_tt_air -0.013073 0.001943 -6.7299 -521.507 0.001956 -6.6836 -517.914 b_tt_rail -0.005460 0.001284 -4.2522 -783.036 0.001253 -4.3592 -802.744 b_tt_shift_business -0.006010 5.1448e-04 -11.6817 -1955.402 4.7225e-04 -12.7264 -2130.264 b_access -0.017829 0.001983 -8.9918 -513.325 0.001894 -9.4112 -537.265 b_cost -0.057587 0.002219 -25.9532 -476.631 0.002234 -25.7755 -473.367 b_cost_shift_business 0.023261 0.002199 10.5778 -444.162 0.002050 11.3476 -476.489 cost_income_elast -0.621005 0.030935 -20.0747 -52.401 0.031203 -19.9024 -51.951 b_no_frills 0.000000 NA NA NA NA NA NA b_wifi 0.774396 0.043847 17.6614 -5.145 0.043516 17.7958 -5.184 b_food 0.337786 0.038453 8.7844 -17.221 0.040823 8.2744 -16.222 lambda_fastPT 0.401188 0.037274 10.7632 -16.065 0.041319 9.7095 -14.492 lambda_groundPT 0.528452 0.043296 12.2055 -10.891 0.045230 11.6837 -10.426 alpha_rail_fastPT 0.492901 0.053623 9.1919 -9.457 0.051644 9.5442 -9.819 Structure for CNL model component : car (alpha) bus (alpha) air (alpha) rail (alpha) lambda car nest 1.000 0.000 0.000 0.0000 1.0000 groundPT nest 0.000 1.000 0.000 0.5071 0.5285 fastPT nest 0.000 0.000 1.000 0.4929 0.4012 Overview of choices for CNL model component : car bus air rail Times available 5446.00 6314.00 5264.00 6118.00 Times chosen 1946.00 358.00 1522.00 3174.00 Percentage chosen overall 27.80 5.11 21.74 45.34 Percentage chosen when available 35.73 5.67 28.91 51.88 Classical covariance matrix: asc_bus asc_air asc_rail asc_bus_shift_female asc_bus 0.198789 0.036424 0.032361 -0.005655 asc_air 0.036424 0.085338 0.050177 -0.001679 asc_rail 0.032361 0.050177 0.082797 -0.001772 asc_bus_shift_female -0.005655 -0.001679 -0.001772 0.011083 asc_air_shift_female -0.001462 -0.002844 -0.001887 0.003092 asc_rail_shift_female -0.001744 -0.001738 -0.002342 0.003366 b_tt_car 1.0331e-04 1.4072e-04 1.4838e-04 -1.386e-06 b_tt_bus -4.5297e-04 3.323e-05 4.904e-05 -1.852e-06 b_tt_air -2.366e-05 -3.3537e-04 -7.031e-06 -9.493e-07 b_tt_rail 1.284e-05 3.191e-06 -2.3368e-04 -1.222e-06 b_tt_shift_business 4.004e-06 9.549e-06 1.484e-05 4.581e-07 b_access -2.015e-05 -2.7134e-04 -1.879e-05 4.745e-07 b_cost -3.231e-05 -2.501e-05 3.227e-05 -1.066e-05 b_cost_shift_business 4.414e-05 -3.684e-05 -9.832e-06 1.084e-05 cost_income_elast -6.0902e-04 9.1816e-04 3.4153e-04 -1.7148e-04 b_wifi 3.6890e-04 -0.001199 -0.001609 8.078e-05 b_food -3.414e-05 -8.6770e-04 -0.001162 9.811e-06 lambda_fastPT 1.679e-05 -6.5425e-04 -7.0757e-04 1.0202e-04 lambda_groundPT 4.7313e-04 -4.663e-05 -9.3948e-04 1.5926e-04 alpha_rail_fastPT -3.9997e-04 -1.9817e-04 5.3192e-04 -2.541e-05 asc_air_shift_female asc_rail_shift_female b_tt_car b_tt_bus asc_bus -0.001462 -0.001744 1.0331e-04 -4.5297e-04 asc_air -0.002844 -0.001738 1.4072e-04 3.323e-05 asc_rail -0.001887 -0.002342 1.4838e-04 4.904e-05 asc_bus_shift_female 0.003092 0.003366 -1.386e-06 -1.852e-06 asc_air_shift_female 0.005873 0.003765 -1.047e-06 -1.494e-06 asc_rail_shift_female 0.003765 0.004529 -5.129e-07 -3.161e-07 b_tt_car -1.047e-06 -5.129e-07 5.087e-07 2.034e-07 b_tt_bus -1.494e-06 -3.161e-07 2.034e-07 1.541e-06 b_tt_air -1.807e-06 -9.815e-07 5.111e-08 1.485e-07 b_tt_rail -1.448e-06 2.158e-07 2.486e-08 3.704e-08 b_tt_shift_business -7.693e-07 -3.664e-07 -1.965e-08 -3.777e-08 b_access -1.886e-06 -2.359e-06 6.436e-08 5.947e-08 b_cost -6.843e-06 -3.207e-06 6.136e-07 9.924e-07 b_cost_shift_business 1.492e-06 1.018e-06 -4.208e-07 -7.698e-07 cost_income_elast 1.868e-05 7.312e-06 -5.671e-07 2.473e-06 b_wifi 1.0093e-04 3.347e-05 -8.127e-06 -1.180e-05 b_food 6.104e-05 4.890e-06 -3.659e-06 -4.488e-06 lambda_fastPT 2.317e-05 2.955e-05 -3.955e-06 -6.944e-06 lambda_groundPT 6.984e-05 1.318e-05 -6.392e-06 -1.338e-05 alpha_rail_fastPT -6.848e-05 -3.261e-06 3.739e-06 4.231e-06 b_tt_air b_tt_rail b_tt_shift_business b_access asc_bus -2.366e-05 1.284e-05 4.004e-06 -2.015e-05 asc_air -3.3537e-04 3.191e-06 9.549e-06 -2.7134e-04 asc_rail -7.031e-06 -2.3368e-04 1.484e-05 -1.879e-05 asc_bus_shift_female -9.493e-07 -1.222e-06 4.581e-07 4.745e-07 asc_air_shift_female -1.807e-06 -1.448e-06 -7.693e-07 -1.886e-06 asc_rail_shift_female -9.815e-07 2.158e-07 -3.664e-07 -2.359e-06 b_tt_car 5.111e-08 2.486e-08 -1.965e-08 6.436e-08 b_tt_bus 1.485e-07 3.704e-08 -3.777e-08 5.947e-08 b_tt_air 3.774e-06 -4.613e-08 -3.011e-08 1.288e-06 b_tt_rail -4.613e-08 1.649e-06 -7.618e-08 -1.708e-07 b_tt_shift_business -3.011e-08 -7.618e-08 2.647e-07 6.041e-08 b_access 1.288e-06 -1.708e-07 6.041e-08 3.932e-06 b_cost 1.131e-06 7.200e-07 -1.958e-07 1.002e-06 b_cost_shift_business -2.912e-07 -5.125e-07 7.255e-07 3.654e-08 cost_income_elast -7.037e-06 -2.594e-06 -5.059e-07 -6.611e-06 b_wifi -1.115e-05 -8.348e-06 -2.181e-06 -1.320e-05 b_food -5.920e-06 -3.363e-06 -1.226e-06 -5.811e-06 lambda_fastPT 4.384e-07 -2.617e-06 4.638e-07 1.228e-06 lambda_groundPT -1.165e-05 -3.865e-06 -1.121e-07 -1.409e-05 alpha_rail_fastPT 1.461e-05 2.603e-06 4.079e-06 1.918e-05 b_cost b_cost_shift_business cost_income_elast b_wifi asc_bus -3.231e-05 4.414e-05 -6.0902e-04 3.6890e-04 asc_air -2.501e-05 -3.684e-05 9.1816e-04 -0.001199 asc_rail 3.227e-05 -9.832e-06 3.4153e-04 -0.001609 asc_bus_shift_female -1.066e-05 1.084e-05 -1.7148e-04 8.078e-05 asc_air_shift_female -6.843e-06 1.492e-06 1.868e-05 1.0093e-04 asc_rail_shift_female -3.207e-06 1.018e-06 7.312e-06 3.347e-05 b_tt_car 6.136e-07 -4.208e-07 -5.671e-07 -8.127e-06 b_tt_bus 9.924e-07 -7.698e-07 2.473e-06 -1.180e-05 b_tt_air 1.131e-06 -2.912e-07 -7.037e-06 -1.115e-05 b_tt_rail 7.200e-07 -5.125e-07 -2.594e-06 -8.348e-06 b_tt_shift_business -1.958e-07 7.255e-07 -5.059e-07 -2.181e-06 b_access 1.002e-06 3.654e-08 -6.611e-06 -1.320e-05 b_cost 4.923e-06 -3.197e-06 -1.825e-05 -4.792e-05 b_cost_shift_business -3.197e-06 4.836e-06 9.943e-08 1.558e-05 cost_income_elast -1.825e-05 9.943e-08 9.5695e-04 1.6662e-04 b_wifi -4.792e-05 1.558e-05 1.6662e-04 0.001923 b_food -2.095e-05 4.659e-06 1.1556e-04 0.001002 lambda_fastPT -2.786e-05 2.411e-05 2.527e-05 2.2136e-04 lambda_groundPT -4.289e-05 2.276e-05 -1.5147e-04 5.9444e-04 alpha_rail_fastPT 2.962e-05 7.910e-06 -2.1151e-04 -6.7298e-04 b_food lambda_fastPT lambda_groundPT alpha_rail_fastPT asc_bus -3.414e-05 1.679e-05 4.7313e-04 -3.9997e-04 asc_air -8.6770e-04 -6.5425e-04 -4.663e-05 -1.9817e-04 asc_rail -0.001162 -7.0757e-04 -9.3948e-04 5.3192e-04 asc_bus_shift_female 9.811e-06 1.0202e-04 1.5926e-04 -2.541e-05 asc_air_shift_female 6.104e-05 2.317e-05 6.984e-05 -6.848e-05 asc_rail_shift_female 4.890e-06 2.955e-05 1.318e-05 -3.261e-06 b_tt_car -3.659e-06 -3.955e-06 -6.392e-06 3.739e-06 b_tt_bus -4.488e-06 -6.944e-06 -1.338e-05 4.231e-06 b_tt_air -5.920e-06 4.384e-07 -1.165e-05 1.461e-05 b_tt_rail -3.363e-06 -2.617e-06 -3.865e-06 2.603e-06 b_tt_shift_business -1.226e-06 4.638e-07 -1.121e-07 4.079e-06 b_access -5.811e-06 1.228e-06 -1.409e-05 1.918e-05 b_cost -2.095e-05 -2.786e-05 -4.289e-05 2.962e-05 b_cost_shift_business 4.659e-06 2.411e-05 2.276e-05 7.910e-06 cost_income_elast 1.1556e-04 2.527e-05 -1.5147e-04 -2.1151e-04 b_wifi 0.001002 2.2136e-04 5.9444e-04 -6.7298e-04 b_food 0.001479 1.607e-05 2.6305e-04 -3.5153e-04 lambda_fastPT 1.607e-05 0.001389 -2.246e-05 9.7579e-04 lambda_groundPT 2.6305e-04 -2.246e-05 0.001875 -9.2818e-04 alpha_rail_fastPT -3.5153e-04 9.7579e-04 -9.2818e-04 0.002875 Robust covariance matrix: asc_bus asc_air asc_rail asc_bus_shift_female asc_bus 0.177016 0.031231 0.030635 -0.004713 asc_air 0.031231 0.081378 0.048657 -9.3249e-04 asc_rail 0.030635 0.048657 0.077662 -9.3277e-04 asc_bus_shift_female -0.004713 -9.3249e-04 -9.3277e-04 0.012436 asc_air_shift_female 8.2547e-04 -0.002310 -0.002398 0.003264 asc_rail_shift_female -7.9086e-04 -0.001148 -0.001598 0.003735 b_tt_car 9.518e-05 1.3598e-04 1.4001e-04 7.577e-07 b_tt_bus -4.0540e-04 4.006e-05 4.284e-05 -3.523e-06 b_tt_air -4.727e-05 -3.2725e-04 -9.992e-06 1.019e-05 b_tt_rail 7.235e-08 -1.441e-06 -2.1600e-04 -6.766e-07 b_tt_shift_business 1.095e-05 1.855e-05 2.487e-05 1.043e-06 b_access 5.016e-05 -2.4901e-04 -3.300e-05 -1.181e-05 b_cost -4.226e-05 -2.178e-05 1.097e-05 -1.528e-05 b_cost_shift_business 8.367e-05 -2.633e-05 2.902e-05 1.258e-05 cost_income_elast -0.001381 9.4450e-04 8.5459e-04 1.1572e-04 b_wifi -4.7817e-04 -0.001038 -0.001152 8.693e-07 b_food 5.8153e-04 -6.0373e-04 -9.0033e-04 2.3761e-04 lambda_fastPT -4.4838e-04 -0.001338 -7.3633e-04 1.4584e-04 lambda_groundPT 0.001666 0.001166 -4.7380e-04 3.8619e-04 alpha_rail_fastPT -8.1181e-04 -0.001104 -1.9726e-04 -7.679e-05 asc_air_shift_female asc_rail_shift_female b_tt_car b_tt_bus asc_bus 8.2547e-04 -7.9086e-04 9.518e-05 -4.0540e-04 asc_air -0.002310 -0.001148 1.3598e-04 4.006e-05 asc_rail -0.002398 -0.001598 1.4001e-04 4.284e-05 asc_bus_shift_female 0.003264 0.003735 7.577e-07 -3.523e-06 asc_air_shift_female 0.005761 0.003890 -1.811e-06 -8.256e-06 asc_rail_shift_female 0.003890 0.004620 1.426e-06 -1.061e-06 b_tt_car -1.811e-06 1.426e-06 4.978e-07 2.138e-07 b_tt_bus -8.256e-06 -1.061e-06 2.138e-07 1.437e-06 b_tt_air -1.124e-07 4.975e-06 6.848e-08 2.169e-07 b_tt_rail 4.140e-06 1.797e-06 5.554e-08 1.104e-07 b_tt_shift_business -7.036e-08 5.619e-07 -2.327e-10 -4.113e-08 b_access -4.309e-06 -2.443e-06 6.204e-08 -8.929e-08 b_cost -1.786e-05 -8.249e-06 5.979e-07 1.040e-06 b_cost_shift_business 1.205e-05 1.193e-05 -3.897e-07 -8.527e-07 cost_income_elast 2.5715e-04 1.9379e-04 5.438e-07 5.166e-06 b_wifi -4.205e-05 -8.592e-05 -7.087e-06 -9.068e-06 b_food -1.4095e-04 -9.296e-05 -3.656e-06 -6.568e-06 lambda_fastPT 6.031e-05 2.0247e-04 -4.232e-06 -6.136e-06 lambda_groundPT 5.495e-05 -1.5076e-04 -4.477e-06 -1.616e-05 alpha_rail_fastPT -2.930e-05 1.9762e-04 1.163e-06 3.824e-06 b_tt_air b_tt_rail b_tt_shift_business b_access asc_bus -4.727e-05 7.235e-08 1.095e-05 5.016e-05 asc_air -3.2725e-04 -1.441e-06 1.855e-05 -2.4901e-04 asc_rail -9.992e-06 -2.1600e-04 2.487e-05 -3.300e-05 asc_bus_shift_female 1.019e-05 -6.766e-07 1.043e-06 -1.181e-05 asc_air_shift_female -1.124e-07 4.140e-06 -7.036e-08 -4.309e-06 asc_rail_shift_female 4.975e-06 1.797e-06 5.619e-07 -2.443e-06 b_tt_car 6.848e-08 5.554e-08 -2.327e-10 6.204e-08 b_tt_bus 2.169e-07 1.104e-07 -4.113e-08 -8.929e-08 b_tt_air 3.826e-06 2.089e-08 -5.610e-08 1.208e-06 b_tt_rail 2.089e-08 1.569e-06 -1.029e-07 -6.281e-08 b_tt_shift_business -5.610e-08 -1.029e-07 2.230e-07 5.354e-08 b_access 1.208e-06 -6.281e-08 5.354e-08 3.589e-06 b_cost 9.779e-07 8.229e-07 -2.039e-07 1.058e-06 b_cost_shift_business -1.559e-07 -7.019e-07 5.480e-07 -1.400e-07 cost_income_elast -5.489e-06 -4.915e-06 -5.178e-07 -3.763e-06 b_wifi -1.130e-05 -9.431e-06 -2.652e-06 -8.950e-06 b_food -9.769e-06 -4.923e-06 -1.299e-06 -4.878e-06 lambda_fastPT 8.459e-06 -4.042e-06 -5.797e-07 4.277e-06 lambda_groundPT -1.534e-05 -1.421e-06 -4.992e-07 -2.055e-05 alpha_rail_fastPT 1.775e-05 2.105e-06 4.401e-06 2.017e-05 b_cost b_cost_shift_business cost_income_elast b_wifi asc_bus -4.226e-05 8.367e-05 -0.001381 -4.7817e-04 asc_air -2.178e-05 -2.633e-05 9.4450e-04 -0.001038 asc_rail 1.097e-05 2.902e-05 8.5459e-04 -0.001152 asc_bus_shift_female -1.528e-05 1.258e-05 1.1572e-04 8.693e-07 asc_air_shift_female -1.786e-05 1.205e-05 2.5715e-04 -4.205e-05 asc_rail_shift_female -8.249e-06 1.193e-05 1.9379e-04 -8.592e-05 b_tt_car 5.979e-07 -3.897e-07 5.438e-07 -7.087e-06 b_tt_bus 1.040e-06 -8.527e-07 5.166e-06 -9.068e-06 b_tt_air 9.779e-07 -1.559e-07 -5.489e-06 -1.130e-05 b_tt_rail 8.229e-07 -7.019e-07 -4.915e-06 -9.431e-06 b_tt_shift_business -2.039e-07 5.480e-07 -5.178e-07 -2.652e-06 b_access 1.058e-06 -1.400e-07 -3.763e-06 -8.950e-06 b_cost 4.992e-06 -3.298e-06 -1.718e-05 -4.503e-05 b_cost_shift_business -3.298e-06 4.202e-06 3.981e-08 1.492e-05 cost_income_elast -1.718e-05 3.981e-08 9.7360e-04 1.3039e-04 b_wifi -4.503e-05 1.492e-05 1.3039e-04 0.001894 b_food -2.039e-05 5.316e-06 3.531e-05 0.001055 lambda_fastPT -3.244e-05 2.330e-05 2.1631e-04 2.4141e-04 lambda_groundPT -4.605e-05 2.312e-05 -2.1548e-04 5.4293e-04 alpha_rail_fastPT 2.157e-05 1.198e-05 -3.523e-05 -5.1586e-04 b_food lambda_fastPT lambda_groundPT alpha_rail_fastPT asc_bus 5.8153e-04 -4.4838e-04 0.001666 -8.1181e-04 asc_air -6.0373e-04 -0.001338 0.001166 -0.001104 asc_rail -9.0033e-04 -7.3633e-04 -4.7380e-04 -1.9726e-04 asc_bus_shift_female 2.3761e-04 1.4584e-04 3.8619e-04 -7.679e-05 asc_air_shift_female -1.4095e-04 6.031e-05 5.495e-05 -2.930e-05 asc_rail_shift_female -9.296e-05 2.0247e-04 -1.5076e-04 1.9762e-04 b_tt_car -3.656e-06 -4.232e-06 -4.477e-06 1.163e-06 b_tt_bus -6.568e-06 -6.136e-06 -1.616e-05 3.824e-06 b_tt_air -9.769e-06 8.459e-06 -1.534e-05 1.775e-05 b_tt_rail -4.923e-06 -4.042e-06 -1.421e-06 2.105e-06 b_tt_shift_business -1.299e-06 -5.797e-07 -4.992e-07 4.401e-06 b_access -4.878e-06 4.277e-06 -2.055e-05 2.017e-05 b_cost -2.039e-05 -3.244e-05 -4.605e-05 2.157e-05 b_cost_shift_business 5.316e-06 2.330e-05 2.312e-05 1.198e-05 cost_income_elast 3.531e-05 2.1631e-04 -2.1548e-04 -3.523e-05 b_wifi 0.001055 2.4141e-04 5.4293e-04 -5.1586e-04 b_food 0.001667 -5.780e-05 2.8688e-04 -3.6835e-04 lambda_fastPT -5.780e-05 0.001707 -8.113e-05 0.001085 lambda_groundPT 2.8688e-04 -8.113e-05 0.002046 -0.001031 alpha_rail_fastPT -3.6835e-04 0.001085 -0.001031 0.002667 Classical correlation matrix: asc_bus asc_air asc_rail asc_bus_shift_female asc_bus 1.000000 0.279650 0.25224 -0.120488 asc_air 0.279650 1.000000 0.59693 -0.054588 asc_rail 0.252244 0.596934 1.00000 -0.058497 asc_bus_shift_female -0.120488 -0.054588 -0.05850 1.000000 asc_air_shift_female -0.042781 -0.127018 -0.08559 0.383269 asc_rail_shift_female -0.058118 -0.088412 -0.12095 0.475173 b_tt_car 0.324882 0.675419 0.72299 -0.018462 b_tt_bus -0.818499 0.091637 0.13730 -0.014176 b_tt_air -0.027319 -0.590974 -0.01258 -0.004642 b_tt_rail 0.022428 0.008506 -0.63247 -0.009043 b_tt_shift_business 0.017457 0.063533 0.10024 0.008457 b_access -0.022787 -0.468443 -0.03293 0.002273 b_cost -0.032661 -0.038581 0.05055 -0.045640 b_cost_shift_business 0.045020 -0.057346 -0.01554 0.046835 cost_income_elast -0.044156 0.101603 0.03837 -0.052654 b_wifi 0.018870 -0.093635 -0.12756 0.017499 b_food -0.001991 -0.077244 -0.10499 0.002423 lambda_fastPT 0.001010 -0.060085 -0.06597 0.025999 lambda_groundPT 0.024509 -0.003687 -0.07541 0.034941 alpha_rail_fastPT -0.016729 -0.012650 0.03447 -0.004501 asc_air_shift_female asc_rail_shift_female b_tt_car b_tt_bus asc_bus -0.042781 -0.058118 0.32488 -0.818499 asc_air -0.127018 -0.088412 0.67542 0.091637 asc_rail -0.085588 -0.120945 0.72299 0.137304 asc_bus_shift_female 0.383269 0.475173 -0.01846 -0.014176 asc_air_shift_female 1.000000 0.730097 -0.01915 -0.015703 asc_rail_shift_female 0.730097 1.000000 -0.01069 -0.003784 b_tt_car -0.019153 -0.010687 1.00000 0.229761 b_tt_bus -0.015703 -0.003784 0.22976 1.000000 b_tt_air -0.012136 -0.007508 0.03689 0.061572 b_tt_rail -0.014716 0.002498 0.02714 0.023240 b_tt_shift_business -0.019510 -0.010582 -0.05356 -0.059148 b_access -0.012412 -0.017681 0.04551 0.024162 b_cost -0.040242 -0.021479 0.38773 0.360330 b_cost_shift_business 0.008851 0.006878 -0.26828 -0.282038 cost_income_elast 0.007879 0.003512 -0.02570 0.064395 b_wifi 0.030036 0.011343 -0.25989 -0.216778 b_food 0.020713 0.001890 -0.13342 -0.094033 lambda_fastPT 0.008110 0.011779 -0.14878 -0.150089 lambda_groundPT 0.021049 0.004523 -0.20700 -0.248930 alpha_rail_fastPT -0.016664 -9.0380e-04 0.09777 0.063564 b_tt_air b_tt_rail b_tt_shift_business b_access asc_bus -0.027319 0.022428 0.017457 -0.022787 asc_air -0.590974 0.008506 0.063533 -0.468443 asc_rail -0.012578 -0.632468 0.100238 -0.032933 asc_bus_shift_female -0.004642 -0.009043 0.008457 0.002273 asc_air_shift_female -0.012136 -0.014716 -0.019510 -0.012412 asc_rail_shift_female -0.007508 0.002498 -0.010582 -0.017681 b_tt_car 0.036889 0.027142 -0.053564 0.045509 b_tt_bus 0.061572 0.023240 -0.059148 0.024162 b_tt_air 1.000000 -0.018492 -0.030128 0.334494 b_tt_rail -0.018492 1.000000 -0.115312 -0.067101 b_tt_shift_business -0.030128 -0.115312 1.000000 0.059216 b_access 0.334494 -0.067101 0.059216 1.000000 b_cost 0.262302 0.252715 -0.171513 0.227827 b_cost_shift_business -0.068161 -0.181493 0.641227 0.008380 cost_income_elast -0.117093 -0.065302 -0.031788 -0.107783 b_wifi -0.130859 -0.148265 -0.096682 -0.151856 b_food -0.079246 -0.068114 -0.061960 -0.076214 lambda_fastPT 0.006054 -0.054683 0.024188 0.016613 lambda_groundPT -0.138543 -0.069518 -0.005034 -0.164122 alpha_rail_fastPT 0.140287 0.037803 0.147866 0.180392 b_cost b_cost_shift_business cost_income_elast b_wifi asc_bus -0.03266 0.045020 -0.044156 0.01887 asc_air -0.03858 -0.057346 0.101603 -0.09364 asc_rail 0.05055 -0.015539 0.038369 -0.12756 asc_bus_shift_female -0.04564 0.046835 -0.052654 0.01750 asc_air_shift_female -0.04024 0.008851 0.007879 0.03004 asc_rail_shift_female -0.02148 0.006878 0.003512 0.01134 b_tt_car 0.38773 -0.268285 -0.025704 -0.25989 b_tt_bus 0.36033 -0.282038 0.064395 -0.21678 b_tt_air 0.26230 -0.068161 -0.117093 -0.13086 b_tt_rail 0.25271 -0.181493 -0.065302 -0.14827 b_tt_shift_business -0.17151 0.641227 -0.031788 -0.09668 b_access 0.22783 0.008380 -0.107783 -0.15186 b_cost 1.00000 -0.655148 -0.265918 -0.49254 b_cost_shift_business -0.65515 1.000000 0.001462 0.16162 cost_income_elast -0.26592 0.001462 1.000000 0.12284 b_wifi -0.49254 0.161622 0.122843 1.00000 b_food -0.24549 0.055092 0.097147 0.59412 lambda_fastPT -0.33682 0.294184 0.021918 0.13544 lambda_groundPT -0.44646 0.239054 -0.113093 0.31313 alpha_rail_fastPT 0.24892 0.067078 -0.127507 -0.28623 b_food lambda_fastPT lambda_groundPT alpha_rail_fastPT asc_bus -0.001991 0.001010 0.024509 -0.016729 asc_air -0.077244 -0.060085 -0.003687 -0.012650 asc_rail -0.104990 -0.065972 -0.075410 0.034473 asc_bus_shift_female 0.002423 0.025999 0.034941 -0.004501 asc_air_shift_female 0.020713 0.008110 0.021049 -0.016664 asc_rail_shift_female 0.001890 0.011779 0.004523 -9.0380e-04 b_tt_car -0.133424 -0.148782 -0.206997 0.097769 b_tt_bus -0.094033 -0.150089 -0.248930 0.063564 b_tt_air -0.079246 0.006054 -0.138543 0.140287 b_tt_rail -0.068114 -0.054683 -0.069518 0.037803 b_tt_shift_business -0.061960 0.024188 -0.005034 0.147866 b_access -0.076214 0.016613 -0.164122 0.180392 b_cost -0.245492 -0.336816 -0.446465 0.248925 b_cost_shift_business 0.055092 0.294184 0.239054 0.067078 cost_income_elast 0.097147 0.021918 -0.113093 -0.127507 b_wifi 0.594115 0.135445 0.313129 -0.286229 b_food 1.000000 0.011214 0.157999 -0.170482 lambda_fastPT 0.011214 1.000000 -0.013920 0.488196 lambda_groundPT 0.157999 -0.013920 1.000000 -0.399787 alpha_rail_fastPT -0.170482 0.488196 -0.399787 1.000000 Robust correlation matrix: asc_bus asc_air asc_rail asc_bus_shift_female asc_bus 1.00000 0.260207 0.26129 -0.100459 asc_air 0.26021 1.000000 0.61206 -0.029312 asc_rail 0.26129 0.612057 1.00000 -0.030014 asc_bus_shift_female -0.10046 -0.029312 -0.03001 1.000000 asc_air_shift_female 0.02585 -0.106680 -0.11339 0.385669 asc_rail_shift_female -0.02765 -0.059197 -0.08434 0.492674 b_tt_car 0.32062 0.675559 0.71206 0.009630 b_tt_bus -0.80374 0.117148 0.12824 -0.026348 b_tt_air -0.05744 -0.586459 -0.01833 0.046707 b_tt_rail 1.3730e-04 -0.004032 -0.61883 -0.004844 b_tt_shift_business 0.05509 0.137678 0.18894 0.019804 b_access 0.06294 -0.460758 -0.06251 -0.055923 b_cost -0.04496 -0.034170 0.01762 -0.061317 b_cost_shift_business 0.09702 -0.045029 0.05080 0.055013 cost_income_elast -0.10519 0.106110 0.09828 0.033255 b_wifi -0.02612 -0.083607 -0.09500 1.7913e-04 b_food 0.03386 -0.051842 -0.07914 0.052193 lambda_fastPT -0.02579 -0.113521 -0.06395 0.031651 lambda_groundPT 0.08756 0.090340 -0.03759 0.076565 alpha_rail_fastPT -0.03736 -0.074927 -0.01371 -0.013334 asc_air_shift_female asc_rail_shift_female b_tt_car b_tt_bus asc_bus 0.025850 -0.02765 0.320616 -0.80374 asc_air -0.106680 -0.05920 0.675559 0.11715 asc_rail -0.113388 -0.08434 0.712060 0.12824 asc_bus_shift_female 0.385669 0.49267 0.009630 -0.02635 asc_air_shift_female 1.000000 0.75394 -0.033819 -0.09074 asc_rail_shift_female 0.753945 1.00000 0.029738 -0.01302 b_tt_car -0.033819 0.02974 1.000000 0.25274 b_tt_bus -0.090736 -0.01302 0.252743 1.00000 b_tt_air -7.5713e-04 0.03742 0.049616 0.09249 b_tt_rail 0.043552 0.02110 0.062845 0.07351 b_tt_shift_business -0.001963 0.01750 -6.9843e-04 -0.07264 b_access -0.029968 -0.01897 0.046413 -0.03931 b_cost -0.105341 -0.05432 0.379293 0.38847 b_cost_shift_business 0.077477 0.08562 -0.269447 -0.34697 cost_income_elast 0.108580 0.09137 0.024698 0.13810 b_wifi -0.012732 -0.02905 -0.230820 -0.17383 b_food -0.045490 -0.03350 -0.126919 -0.13420 lambda_fastPT 0.019232 0.07209 -0.145154 -0.12387 lambda_groundPT 0.016008 -0.04904 -0.140291 -0.29796 alpha_rail_fastPT -0.007475 0.05630 0.031927 0.06177 b_tt_air b_tt_rail b_tt_shift_business b_access asc_bus -0.057438 1.3730e-04 0.055087 0.06294 asc_air -0.586459 -0.004032 0.137678 -0.46076 asc_rail -0.018330 -0.618827 0.188942 -0.06251 asc_bus_shift_female 0.046707 -0.004844 0.019804 -0.05592 asc_air_shift_female -7.5713e-04 0.043552 -0.001963 -0.02997 asc_rail_shift_female 0.037416 0.021101 0.017505 -0.01897 b_tt_car 0.049616 0.062845 -6.9843e-04 0.04641 b_tt_bus 0.092492 0.073513 -0.072645 -0.03931 b_tt_air 1.000000 0.008526 -0.060733 0.32599 b_tt_rail 0.008526 1.000000 -0.173991 -0.02647 b_tt_shift_business -0.060733 -0.173991 1.000000 0.05985 b_access 0.325992 -0.026471 0.059847 1.00000 b_cost 0.223772 0.294065 -0.193276 0.24995 b_cost_shift_business -0.038879 -0.273393 0.566065 -0.03605 cost_income_elast -0.089929 -0.125754 -0.035139 -0.06365 b_wifi -0.132745 -0.173036 -0.129064 -0.10856 b_food -0.122341 -0.096283 -0.067367 -0.06307 lambda_fastPT 0.104665 -0.078092 -0.029709 0.05464 lambda_groundPT -0.173392 -0.025086 -0.023370 -0.23985 alpha_rail_fastPT 0.175745 0.032548 0.180462 0.20617 b_cost b_cost_shift_business cost_income_elast b_wifi asc_bus -0.04496 0.09702 -0.10519 -0.02612 asc_air -0.03417 -0.04503 0.10611 -0.08361 asc_rail 0.01762 0.05080 0.09828 -0.09500 asc_bus_shift_female -0.06132 0.05501 0.03326 1.7913e-04 asc_air_shift_female -0.10534 0.07748 0.10858 -0.01273 asc_rail_shift_female -0.05432 0.08562 0.09137 -0.02905 b_tt_car 0.37929 -0.26945 0.02470 -0.23082 b_tt_bus 0.38847 -0.34697 0.13810 -0.17383 b_tt_air 0.22377 -0.03888 -0.08993 -0.13275 b_tt_rail 0.29406 -0.27339 -0.12575 -0.17304 b_tt_shift_business -0.19328 0.56606 -0.03514 -0.12906 b_access 0.24995 -0.03605 -0.06365 -0.10856 b_cost 1.00000 -0.72019 -0.24648 -0.46313 b_cost_shift_business -0.72019 1.00000 6.2242e-04 0.16730 cost_income_elast -0.24648 6.2242e-04 1.00000 0.09603 b_wifi -0.46313 0.16730 0.09603 1.00000 b_food -0.22352 0.06353 0.02772 0.59364 lambda_fastPT -0.35140 0.27511 0.16778 0.13426 lambda_groundPT -0.45569 0.24934 -0.15269 0.27585 alpha_rail_fastPT 0.18695 0.11315 -0.02186 -0.22954 b_food lambda_fastPT lambda_groundPT alpha_rail_fastPT asc_bus 0.03386 -0.02579 0.08756 -0.037362 asc_air -0.05184 -0.11352 0.09034 -0.074927 asc_rail -0.07914 -0.06395 -0.03759 -0.013706 asc_bus_shift_female 0.05219 0.03165 0.07657 -0.013334 asc_air_shift_female -0.04549 0.01923 0.01601 -0.007475 asc_rail_shift_female -0.03350 0.07209 -0.04904 0.056296 b_tt_car -0.12692 -0.14515 -0.14029 0.031927 b_tt_bus -0.13420 -0.12387 -0.29796 0.061766 b_tt_air -0.12234 0.10466 -0.17339 0.175745 b_tt_rail -0.09628 -0.07809 -0.02509 0.032548 b_tt_shift_business -0.06737 -0.02971 -0.02337 0.180462 b_access -0.06307 0.05464 -0.23985 0.206172 b_cost -0.22352 -0.35140 -0.45569 0.186954 b_cost_shift_business 0.06353 0.27511 0.24934 0.113153 cost_income_elast 0.02772 0.16778 -0.15269 -0.021861 b_wifi 0.59364 0.13426 0.27585 -0.229545 b_food 1.00000 -0.03427 0.15537 -0.174719 lambda_fastPT -0.03427 1.00000 -0.04341 0.508485 lambda_groundPT 0.15537 -0.04341 1.00000 -0.441502 alpha_rail_fastPT -0.17472 0.50849 -0.44150 1.000000 20 worst outliers in terms of lowest average per choice prediction: ID Avg prob per choice 400 0.2057383 464 0.2468192 181 0.2623595 146 0.2670334 317 0.2733098 276 0.2772089 370 0.2807020 447 0.2859821 367 0.2879349 161 0.2892147 59 0.2933226 434 0.2937566 379 0.2970106 307 0.2971931 293 0.2974712 73 0.2975922 186 0.2994407 259 0.2999630 441 0.3017757 467 0.3053905 Changes in parameter estimates from starting values: Initial Estimate Difference asc_car 0.000000 0.000000 0.000000 asc_bus 0.286548 0.299289 0.012741 asc_air -0.903342 -0.646373 0.256969 asc_rail -2.092649 -1.563495 0.529155 asc_bus_shift_female 0.340176 0.273647 -0.066529 asc_air_shift_female 0.268177 0.240335 -0.027843 asc_rail_shift_female 0.189615 0.178849 -0.010767 b_tt_car -0.013107 -0.010391 0.002717 b_tt_bus -0.021266 -0.016555 0.004711 b_tt_air -0.016578 -0.013073 0.003504 b_tt_rail -0.007051 -0.005460 0.001591 b_tt_shift_business -0.006234 -0.006010 2.2369e-04 b_access -0.021153 -0.017829 0.003323 b_cost -0.076190 -0.057587 0.018603 b_cost_shift_business 0.033381 0.023261 -0.010120 cost_income_elast -0.613795 -0.621005 -0.007211 b_no_frills 0.000000 0.000000 0.000000 b_wifi 1.026713 0.774396 -0.252316 b_food 0.422069 0.337786 -0.084283 lambda_fastPT 1.000000 0.401188 -0.598812 lambda_groundPT 1.000000 0.528452 -0.471548 alpha_rail_fastPT 0.500000 0.492901 -0.007099 Settings and functions used in model definition: apollo_control -------------- Value modelName "CNL" modelDescr "CNL model with socio-demographics on mode choice SP data" indivID "ID" outputDirectory "output/" debug "FALSE" nCores "1" workInLogs "FALSE" seed "13" mixing "FALSE" HB "FALSE" noValidation "FALSE" noDiagnostics "FALSE" calculateLLC "TRUE" panelData "TRUE" analyticGrad "TRUE" analyticGrad_manualSet "FALSE" overridePanel "FALSE" preventOverridePanel "FALSE" noModification "FALSE" Hessian routines attempted -------------------------- numerical second derivative of LL (using numDeriv) Scaling in estimation --------------------- Value asc_bus 0.299289420 asc_air 0.646373098 asc_rail 1.563497075 asc_bus_shift_female 0.273647349 asc_air_shift_female 0.240334530 asc_rail_shift_female 0.178848711 b_tt_car 0.010390738 b_tt_bus 0.016555182 b_tt_air 0.013073470 b_tt_rail 0.005460034 b_tt_shift_business 0.006009987 b_access 0.017829193 b_cost 0.057586544 b_cost_shift_business 0.023261216 cost_income_elast 0.621005828 b_wifi 0.774395461 b_food 0.337785910 lambda_fastPT 0.401188489 lambda_groundPT 0.528451251 alpha_rail_fastPT 0.492899453 Scaling used in computing Hessian --------------------------------- Value asc_bus 0.299289375 asc_air 0.646372936 asc_rail 1.563494625 asc_bus_shift_female 0.273647337 asc_air_shift_female 0.240334520 asc_rail_shift_female 0.178848694 b_tt_car 0.010390793 b_tt_bus 0.016555200 b_tt_air 0.013073460 b_tt_rail 0.005460028 b_tt_shift_business 0.006009982 b_access 0.017829183 b_cost 0.057587032 b_cost_shift_business 0.023261126 cost_income_elast 0.621005465 b_wifi 0.774396451 b_food 0.337785884 lambda_fastPT 0.401187796 lambda_groundPT 0.528451531 alpha_rail_fastPT 0.492900543 apollo_probabilities ---------------------- function(apollo_beta, apollo_inputs, functionality="estimate"){ ### Attach inputs and detach after function exit apollo_attach(apollo_beta, apollo_inputs) on.exit(apollo_detach(apollo_beta, apollo_inputs)) ### Create list of probabilities P P = list() ### Create alternative specific constants and coefficients using interactions with socio-demographics asc_bus_value = asc_bus + asc_bus_shift_female * female asc_air_value = asc_air + asc_air_shift_female * female asc_rail_value = asc_rail + asc_rail_shift_female * female b_tt_car_value = b_tt_car + b_tt_shift_business * business b_tt_bus_value = b_tt_bus + b_tt_shift_business * business b_tt_air_value = b_tt_air + b_tt_shift_business * business b_tt_rail_value = b_tt_rail + b_tt_shift_business * business b_cost_value = ( b_cost + b_cost_shift_business * business ) * ( income / mean_income ) ^ cost_income_elast ### List of utilities: these must use the same names as in cnl_settings, order is irrelevant V = list() V[["car"]] = asc_car + b_tt_car_value * time_car + b_cost_value * cost_car V[["bus"]] = asc_bus_value + b_tt_bus_value * time_bus + b_access * access_bus + b_cost_value * cost_bus V[["air"]] = asc_air_value + b_tt_air_value * time_air + b_access * access_air + b_cost_value * cost_air + b_no_frills * ( service_air == 1 ) + b_wifi * ( service_air == 2 ) + b_food * ( service_air == 3 ) V[["rail"]] = asc_rail_value + b_tt_rail_value * time_rail + b_access * access_rail + b_cost_value * cost_rail + b_no_frills * ( service_rail == 1 ) + b_wifi * ( service_rail == 2 ) + b_food * ( service_rail == 3 ) ### Specify nests for CNL model cnlNests = list(car=1,groundPT=lambda_groundPT,fastPT=lambda_fastPT) ### Specify nest allocation parameters for alternatives included in multiple nests alpha_rail_groundPT = 1 - alpha_rail_fastPT ### Specify tree structure, showing membership in nests (one row per nest, one column per alternative) cnlStructure = matrix(0, nrow=length(cnlNests), ncol=length(V)) cnlStructure[1,] = c( 1, 0, 0, 0 ) # car cnlStructure[2,] = c( 0, 1, 0, alpha_rail_groundPT) # groundPT cnlStructure[3,] = c( 0, 0, 1, alpha_rail_fastPT ) # fastPT ### Define settings for CNL model cnl_settings <- list( alternatives = c(car=1, bus=2, air=3, rail=4), avail = list(car=av_car, bus=av_bus, air=av_air, rail=av_rail), choiceVar = choice, utilities = V, cnlNests = cnlNests, cnlStructure = cnlStructure ) ### Compute probabilities using CNL model P[["model"]] = apollo_cnl(cnl_settings, functionality) ### Take product across observation for same individual P = apollo_panelProd(P, apollo_inputs, functionality) ### Prepare and return outputs of function P = apollo_prepareProb(P, apollo_inputs, functionality) return(P) }