Model run by stephane.hess using Apollo 0.2.9 on R 4.0.5 for Darwin. www.ApolloChoiceModelling.com Model name : MDCEV_with_outside_good Model description : MDCEV model on time use data, alpha-gamma profile with outside good and socio-demographics Model run at : 2023-05-10 22:27:57 Estimation method : bfgs Model diagnosis : successful convergence Optimisation diagnosis : Maximum found hessian properties : Negative definitive maximum eigenvalue : -0.000194 Number of individuals : 447 Number of rows in database : 2826 Number of modelled outcomes : 2826 Number of cores used : 1 Model without mixing LL(start) : -21594.92 LL at equal shares, LL(0) : NA LL at observed shares, LL(C) : NA LL(final) : -15007.36 Rho-squared vs equal shares : Not applicable Adj.Rho-squared vs equal shares : Not applicable Rho-squared vs observed shares : Not applicable Adj.Rho-squared vs observed shares : Not applicable AIC : 30044.72 BIC : 30133.92 Estimated parameters : 15 Time taken (hh:mm:ss) : 00:00:10.04 pre-estimation : 00:00:0.25 estimation : 00:00:6.68 initial estimation : 00:00:6.38 estimation after rescaling : 00:00:0.3 post-estimation : 00:00:3.11 Iterations : 67 initial estimation : 64 estimation after rescaling : 3 Unconstrained optimisation. Estimates: Estimate s.e. t.rat.(0) Rob.s.e. Rob.t.rat.(0) alpha_base -14.7609 71.86163 -0.2054 0.12283 -120.172 gamma_work 4.8990 0.29747 16.4689 0.28123 17.420 gamma_school 3.0985 0.58589 5.2885 0.54811 5.653 gamma_shopping 0.4264 0.02697 15.8126 0.03088 13.808 gamma_private 0.6200 0.05038 12.3075 0.06763 9.168 gamma_leisure 2.0957 0.11401 18.3815 0.12633 16.590 delta_work -3.7172 0.07070 -52.5777 0.11544 -32.201 delta_school -7.4182 0.22459 -33.0308 0.34571 -21.458 delta_shopping -3.8044 0.04154 -91.5828 0.05212 -72.999 delta_private -4.2786 0.04785 -89.4163 0.06188 -69.148 delta_leisure -3.4001 0.04423 -76.8731 0.05645 -60.227 delta_work_FT 1.3249 0.08157 16.2421 0.13484 9.825 delta_work_wknd -2.8609 0.14288 -20.0229 0.19540 -14.641 delta_school_young 2.3441 0.25750 9.1032 0.40177 5.834 delta_leisure_wknd 0.2949 0.07139 4.1303 0.07397 3.987 sig 1.0000 NA NA NA NA Overview of choices for model component : outside work school shopping private leisure Times available 2826.00 2826.00 2826.00 2826.0 2826.00 2826.00 Observations in which chosen 2826.00 1139.00 85.00 783.0 535.00 1174.00 Average consumption when available 18.21 2.89 0.11 0.5 0.44 1.85 Average consumption when chosen 18.21 7.17 3.52 1.8 2.34 4.45 Classical covariance matrix: alpha_base gamma_work gamma_school gamma_shopping gamma_private alpha_base 5164.093488 -0.014792 -0.008291 -0.001349 -0.002184 gamma_work -0.014792 0.088490 -4.7924e-04 -3.032e-05 -5.857e-05 gamma_school -0.008291 -4.7924e-04 0.343272 -4.525e-06 -1.887e-05 gamma_shopping -0.001349 -3.032e-05 -4.525e-06 7.2727e-04 3.078e-07 gamma_private -0.002184 -5.857e-05 -1.887e-05 3.078e-07 0.002538 gamma_leisure -0.006809 -5.8250e-04 -1.7293e-04 -2.710e-06 -1.059e-05 delta_work 0.005651 -0.002152 1.2331e-04 4.549e-06 6.641e-06 delta_school 0.005585 -1.0685e-04 -2.2872e-04 2.067e-06 7.219e-07 delta_shopping 0.005596 -6.181e-05 3.816e-06 -1.6228e-04 -6.371e-07 delta_private 0.005613 -5.705e-05 1.852e-06 3.014e-07 -2.1534e-04 delta_leisure 0.005591 -4.349e-05 3.477e-05 -2.693e-06 -8.976e-07 delta_work_FT -1.7952e-04 -0.004915 -1.0582e-04 -6.964e-06 -7.118e-06 delta_work_wknd 2.1221e-04 0.006179 -7.759e-05 1.055e-05 6.500e-06 delta_school_young -2.701e-05 7.568e-05 -0.003040 -4.638e-06 -4.827e-07 delta_leisure_wknd 2.4181e-04 8.532e-05 -6.774e-05 1.083e-05 7.623e-06 gamma_leisure delta_work delta_school delta_shopping delta_private alpha_base -0.006809 0.005651 0.005585 0.005596 0.005613 gamma_work -5.8250e-04 -0.002152 -1.0685e-04 -6.181e-05 -5.705e-05 gamma_school -1.7293e-04 1.2331e-04 -2.2872e-04 3.816e-06 1.852e-06 gamma_shopping -2.710e-06 4.549e-06 2.067e-06 -1.6228e-04 3.014e-07 gamma_private -1.059e-05 6.641e-06 7.219e-07 -6.371e-07 -2.1534e-04 gamma_leisure 0.012999 4.638e-05 -1.806e-05 -2.059e-05 -1.487e-05 delta_work 4.638e-05 0.004998 3.8084e-04 3.9198e-04 3.8867e-04 delta_school -1.806e-05 3.8084e-04 0.050438 3.8367e-04 3.8733e-04 delta_shopping -2.059e-05 3.9198e-04 3.8367e-04 0.001726 3.8697e-04 delta_private -1.487e-05 3.8867e-04 3.8733e-04 3.8697e-04 0.002290 delta_leisure -8.8791e-04 5.1000e-04 3.7852e-04 3.8834e-04 3.8100e-04 delta_work_FT -2.779e-05 -0.004565 2.160e-06 -7.660e-06 -1.094e-05 delta_work_wknd -2.842e-05 -8.7597e-04 3.903e-05 3.423e-06 2.608e-05 delta_school_young 1.450e-07 2.841e-05 -0.050218 4.725e-06 -6.547e-06 delta_leisure_wknd -5.3943e-04 -3.6577e-04 5.464e-05 2.392e-05 4.254e-05 delta_leisure delta_work_FT delta_work_wknd delta_school_young delta_leisure_wknd alpha_base 0.005591 -1.7952e-04 2.1221e-04 -2.701e-05 2.4181e-04 gamma_work -4.349e-05 -0.004915 0.006179 7.568e-05 8.532e-05 gamma_school 3.477e-05 -1.0582e-04 -7.759e-05 -0.003040 -6.774e-05 gamma_shopping -2.693e-06 -6.964e-06 1.055e-05 -4.638e-06 1.083e-05 gamma_private -8.976e-07 -7.118e-06 6.500e-06 -4.827e-07 7.623e-06 gamma_leisure -8.8791e-04 -2.779e-05 -2.842e-05 1.450e-07 -5.3943e-04 delta_work 5.1000e-04 -0.004565 -8.7597e-04 2.841e-05 -3.6577e-04 delta_school 3.7852e-04 2.160e-06 3.903e-05 -0.050218 5.464e-05 delta_shopping 3.8834e-04 -7.660e-06 3.423e-06 4.725e-06 2.392e-05 delta_private 3.8100e-04 -1.094e-05 2.608e-05 -6.547e-06 4.254e-05 delta_leisure 0.001956 -3.144e-06 -3.6992e-04 -7.159e-06 -0.001715 delta_work_FT -3.144e-06 0.006653 -0.001224 -5.774e-05 -1.489e-06 delta_work_wknd -3.6992e-04 -0.001224 0.020415 -1.384e-05 0.001190 delta_school_young -7.159e-06 -5.774e-05 -1.384e-05 0.066308 6.147e-06 delta_leisure_wknd -0.001715 -1.489e-06 0.001190 6.147e-06 0.005097 Robust covariance matrix: alpha_base gamma_work gamma_school gamma_shopping gamma_private alpha_base 0.015087 0.006720 0.006149 7.4962e-04 5.6947e-04 gamma_work 0.006720 0.079092 0.011724 -3.0702e-04 0.001014 gamma_school 0.006149 0.011724 0.300429 -6.2123e-04 7.5428e-04 gamma_shopping 7.4962e-04 -3.0702e-04 -6.2123e-04 9.5379e-04 1.3060e-04 gamma_private 5.6947e-04 0.001014 7.5428e-04 1.3060e-04 0.004573 gamma_leisure 0.003727 0.003811 0.001502 4.8374e-04 -4.5304e-04 delta_work -0.001252 -0.002909 -0.002683 2.3990e-04 -7.3746e-04 delta_school 0.010629 0.009100 0.045071 4.1163e-04 1.0307e-04 delta_shopping 8.1896e-04 -9.9808e-04 0.001006 -2.4733e-04 -2.4934e-04 delta_private 7.5242e-04 0.001046 0.001523 -7.425e-05 -1.2222e-04 delta_leisure 8.900e-05 -0.001758 0.002269 -5.746e-05 -2.1638e-04 delta_work_FT -0.001022 -0.012015 0.003376 -1.3694e-04 4.7884e-04 delta_work_wknd 0.002976 0.020238 0.001382 8.685e-05 6.9013e-04 delta_school_young -0.010313 -0.004020 -0.046995 -1.4011e-04 2.4370e-04 delta_leisure_wknd -6.404e-05 9.0878e-04 0.001588 -6.766e-05 3.241e-05 gamma_leisure delta_work delta_school delta_shopping delta_private alpha_base 0.003727 -0.001252 0.010629 8.1896e-04 7.5242e-04 gamma_work 0.003811 -0.002909 0.009100 -9.9808e-04 0.001046 gamma_school 0.001502 -0.002683 0.045071 0.001006 0.001523 gamma_shopping 4.8374e-04 2.3990e-04 4.1163e-04 -2.4733e-04 -7.425e-05 gamma_private -4.5304e-04 -7.3746e-04 1.0307e-04 -2.4934e-04 -1.2222e-04 gamma_leisure 0.015958 8.2695e-04 0.005251 -6.8343e-04 -6.2042e-04 delta_work 8.2695e-04 0.013326 0.004196 5.7882e-04 -4.0312e-04 delta_school 0.005251 0.004196 0.119517 -3.5682e-04 8.7800e-04 delta_shopping -6.8343e-04 5.7882e-04 -3.5682e-04 0.002716 7.1246e-04 delta_private -6.2042e-04 -4.0312e-04 8.7800e-04 7.1246e-04 0.003829 delta_leisure -0.001664 8.539e-05 3.5050e-04 8.8867e-04 0.001031 delta_work_FT -0.001236 -0.013172 -0.005040 -2.0933e-04 4.0090e-04 delta_work_wknd 0.003070 9.7600e-04 0.011725 -0.001205 -2.5944e-04 delta_school_young -0.004576 -0.008306 -0.119801 4.4853e-04 -0.001116 delta_leisure_wknd 4.1977e-04 1.5148e-04 -1.7539e-04 2.4589e-04 4.567e-06 delta_leisure delta_work_FT delta_work_wknd delta_school_young delta_leisure_wknd alpha_base 8.900e-05 -0.001022 0.002976 -0.010313 -6.404e-05 gamma_work -0.001758 -0.012015 0.020238 -0.004020 9.0878e-04 gamma_school 0.002269 0.003376 0.001382 -0.046995 0.001588 gamma_shopping -5.746e-05 -1.3694e-04 8.685e-05 -1.4011e-04 -6.766e-05 gamma_private -2.1638e-04 4.7884e-04 6.9013e-04 2.4370e-04 3.241e-05 gamma_leisure -0.001664 -0.001236 0.003070 -0.004576 4.1977e-04 delta_work 8.539e-05 -0.013172 9.7600e-04 -0.008306 1.5148e-04 delta_school 3.5050e-04 -0.005040 0.011725 -0.119801 -1.7539e-04 delta_shopping 8.8867e-04 -2.0933e-04 -0.001205 4.4853e-04 2.4589e-04 delta_private 0.001031 4.0090e-04 -2.5944e-04 -0.001116 4.567e-06 delta_leisure 0.003187 5.8629e-04 -8.7438e-04 -2.9211e-04 -0.001780 delta_work_FT 5.8629e-04 0.018182 -0.006396 0.008273 -3.0017e-04 delta_work_wknd -8.7438e-04 -0.006396 0.038180 -0.005667 2.4376e-04 delta_school_young -2.9211e-04 0.008273 -0.005667 0.161418 0.001891 delta_leisure_wknd -0.001780 -3.0017e-04 2.4376e-04 0.001891 0.005471 Classical correlation matrix: alpha_base gamma_work gamma_school gamma_shopping gamma_private alpha_base 1.000000 -6.9195e-04 -1.9691e-04 -6.9631e-04 -6.0319e-04 gamma_work -6.9195e-04 1.000000 -0.002750 -0.003780 -0.003909 gamma_school -1.9691e-04 -0.002750 1.000000 -2.8640e-04 -6.3932e-04 gamma_shopping -6.9631e-04 -0.003780 -2.8640e-04 1.000000 2.2660e-04 gamma_private -6.0319e-04 -0.003909 -6.3932e-04 2.2660e-04 1.000000 gamma_leisure -8.3107e-04 -0.017175 -0.002589 -8.8145e-04 -0.001845 delta_work 0.001112 -0.102327 0.002977 0.002386 0.001865 delta_school 3.4604e-04 -0.001599 -0.001738 3.4120e-04 6.381e-05 delta_shopping 0.001875 -0.005002 1.5678e-04 -0.144864 -3.0447e-04 delta_private 0.001632 -0.004008 6.607e-05 2.3358e-04 -0.089334 delta_leisure 0.001759 -0.003306 0.001342 -0.002258 -4.0283e-04 delta_work_FT -3.063e-05 -0.202559 -0.002214 -0.003166 -0.001732 delta_work_wknd 2.067e-05 0.145374 -9.2691e-04 0.002737 9.0311e-04 delta_school_young -1.460e-06 9.8796e-04 -0.020152 -6.6781e-04 -3.721e-05 delta_leisure_wknd 4.713e-05 0.004018 -0.001619 0.005627 0.002120 gamma_leisure delta_work delta_school delta_shopping delta_private alpha_base -8.3107e-04 0.001112 3.4604e-04 0.001875 0.001632 gamma_work -0.017175 -0.102327 -0.001599 -0.005002 -0.004008 gamma_school -0.002589 0.002977 -0.001738 1.5678e-04 6.607e-05 gamma_shopping -8.8145e-04 0.002386 3.4120e-04 -0.144864 2.3358e-04 gamma_private -0.001845 0.001865 6.381e-05 -3.0447e-04 -0.089334 gamma_leisure 1.000000 0.005754 -7.0523e-04 -0.004348 -0.002726 delta_work 0.005754 1.000000 0.023985 0.133469 0.114889 delta_school -7.0523e-04 0.023985 1.000000 0.041125 0.036042 delta_shopping -0.004348 0.133469 0.041125 1.000000 0.194681 delta_private -0.002726 0.114889 0.036042 0.194681 1.000000 delta_leisure -0.176079 0.163094 0.038106 0.211361 0.180021 delta_work_FT -0.002989 -0.791507 1.1792e-04 -0.002261 -0.002804 delta_work_wknd -0.001744 -0.086715 0.001216 5.7670e-04 0.003815 delta_school_young 4.938e-06 0.001560 -0.868347 4.4176e-04 -5.3133e-04 delta_leisure_wknd -0.066271 -0.072463 0.003408 0.008064 0.012453 delta_leisure delta_work_FT delta_work_wknd delta_school_young delta_leisure_wknd alpha_base 0.001759 -3.063e-05 2.067e-05 -1.460e-06 4.713e-05 gamma_work -0.003306 -0.202559 0.145374 9.8796e-04 0.004018 gamma_school 0.001342 -0.002214 -9.2691e-04 -0.020152 -0.001619 gamma_shopping -0.002258 -0.003166 0.002737 -6.6781e-04 0.005627 gamma_private -4.0283e-04 -0.001732 9.0311e-04 -3.721e-05 0.002120 gamma_leisure -0.176079 -0.002989 -0.001744 4.938e-06 -0.066271 delta_work 0.163094 -0.791507 -0.086715 0.001560 -0.072463 delta_school 0.038106 1.1792e-04 0.001216 -0.868347 0.003408 delta_shopping 0.211361 -0.002261 5.7670e-04 4.4176e-04 0.008064 delta_private 0.180021 -0.002804 0.003815 -5.3133e-04 0.012453 delta_leisure 1.000000 -8.7147e-04 -0.058536 -6.2854e-04 -0.543124 delta_work_FT -8.7147e-04 1.000000 -0.104999 -0.002749 -2.5568e-04 delta_work_wknd -0.058536 -0.104999 1.000000 -3.7628e-04 0.116657 delta_school_young -6.2854e-04 -0.002749 -3.7628e-04 1.000000 3.3436e-04 delta_leisure_wknd -0.543124 -2.5568e-04 0.116657 3.3436e-04 1.000000 Robust correlation matrix: alpha_base gamma_work gamma_school gamma_shopping gamma_private alpha_base 1.000000 0.19452 0.09133 0.19761 0.068556 gamma_work 0.194519 1.00000 0.07606 -0.03535 0.053337 gamma_school 0.091328 0.07606 1.00000 -0.03670 0.020349 gamma_shopping 0.197608 -0.03535 -0.03670 1.00000 0.062531 gamma_private 0.068556 0.05334 0.02035 0.06253 1.000000 gamma_leisure 0.240222 0.10727 0.02169 0.12399 -0.053030 delta_work -0.088281 -0.08961 -0.04240 0.06729 -0.094463 delta_school 0.250314 0.09359 0.23785 0.03855 0.004409 delta_shopping 0.127935 -0.06810 0.03522 -0.15367 -0.070746 delta_private 0.098998 0.06009 0.04491 -0.03886 -0.029207 delta_leisure 0.012835 -0.11072 0.07331 -0.03295 -0.056676 delta_work_FT -0.061709 -0.31683 0.04567 -0.03288 0.052511 delta_work_wknd 0.124006 0.36828 0.01290 0.01439 0.052227 delta_school_young -0.208983 -0.03557 -0.21341 -0.01129 0.008969 delta_leisure_wknd -0.007048 0.04369 0.03918 -0.02962 0.006479 gamma_leisure delta_work delta_school delta_shopping delta_private alpha_base 0.24022 -0.08828 0.250314 0.12793 0.09900 gamma_work 0.10727 -0.08961 0.093594 -0.06810 0.06009 gamma_school 0.02169 -0.04240 0.237854 0.03522 0.04491 gamma_shopping 0.12399 0.06729 0.038554 -0.15367 -0.03886 gamma_private -0.05303 -0.09446 0.004409 -0.07075 -0.02921 gamma_leisure 1.00000 0.05671 0.120244 -0.10381 -0.07937 delta_work 0.05671 1.00000 0.105149 0.09621 -0.05644 delta_school 0.12024 0.10515 1.000000 -0.01980 0.04104 delta_shopping -0.10381 0.09621 -0.019805 1.00000 0.22094 delta_private -0.07937 -0.05644 0.041045 0.22094 1.00000 delta_leisure -0.23332 0.01310 0.017959 0.30205 0.29516 delta_work_FT -0.07255 -0.84620 -0.108126 -0.02979 0.04805 delta_work_wknd 0.12439 0.04327 0.173569 -0.11837 -0.02146 delta_school_young -0.09016 -0.17908 -0.862523 0.02142 -0.04490 delta_leisure_wknd 0.04492 0.01774 -0.006859 0.06379 9.9787e-04 delta_leisure delta_work_FT delta_work_wknd delta_school_young delta_leisure_wknd alpha_base 0.01284 -0.06171 0.12401 -0.208983 -0.007048 gamma_work -0.11072 -0.31683 0.36828 -0.035574 0.043686 gamma_school 0.07331 0.04567 0.01290 -0.213406 0.039180 gamma_shopping -0.03295 -0.03288 0.01439 -0.011292 -0.029619 gamma_private -0.05668 0.05251 0.05223 0.008969 0.006479 gamma_leisure -0.23332 -0.07255 0.12439 -0.090155 0.044922 delta_work 0.01310 -0.84620 0.04327 -0.179081 0.017739 delta_school 0.01796 -0.10813 0.17357 -0.862523 -0.006859 delta_shopping 0.30205 -0.02979 -0.11837 0.021421 0.063786 delta_private 0.29516 0.04805 -0.02146 -0.044899 9.9787e-04 delta_leisure 1.00000 0.07702 -0.07926 -0.012879 -0.426303 delta_work_FT 0.07702 1.00000 -0.24276 0.152707 -0.030096 delta_work_wknd -0.07926 -0.24276 1.00000 -0.072186 0.016865 delta_school_young -0.01288 0.15271 -0.07219 1.000000 0.063647 delta_leisure_wknd -0.42630 -0.03010 0.01687 0.063647 1.000000 20 worst outliers in terms of lowest average per choice prediction: ID Avg prob per choice 2146576 1.310627e-06 9902059 1.643701e-05 9880000 2.639109e-05 6292582 3.777632e-05 3051038 6.564045e-05 8415029 6.887219e-05 8530000 7.160639e-05 9230000 7.424581e-05 7560000 8.724992e-05 8494828 9.341737e-05 1226238 9.510117e-05 5436896 9.791518e-05 7177352 9.979044e-05 6689995 1.057281e-04 1566942 1.153969e-04 3730000 1.300519e-04 2942411 1.303263e-04 9216479 1.404902e-04 6490000 1.466017e-04 2119561 1.511202e-04 Changes in parameter estimates from starting values: Initial Estimate Difference alpha_base 0.000 -14.7609 -14.7609 gamma_work 1.000 4.8990 3.8990 gamma_school 1.000 3.0985 2.0985 gamma_shopping 1.000 0.4264 -0.5736 gamma_private 1.000 0.6200 -0.3800 gamma_leisure 1.000 2.0957 1.0957 delta_work 0.000 -3.7172 -3.7172 delta_school 0.000 -7.4182 -7.4182 delta_shopping 0.000 -3.8044 -3.8044 delta_private 0.000 -4.2786 -4.2786 delta_leisure 0.000 -3.4001 -3.4001 delta_work_FT 0.000 1.3249 1.3249 delta_work_wknd 0.000 -2.8609 -2.8609 delta_school_young 0.000 2.3441 2.3441 delta_leisure_wknd 0.000 0.2949 0.2949 sig 1.000 1.0000 0.0000 Settings and functions used in model definition: apollo_control -------------- Value modelName "MDCEV_with_outside_good" modelDescr "MDCEV model on time use data, alpha-gamma profile with outside good and socio-demographics" indivID "indivID" 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 alpha_base 14.7608888 gamma_work 4.8990258 gamma_school 3.0984826 gamma_shopping 0.4264321 gamma_private 0.6200096 gamma_leisure 2.0956962 delta_work 3.7172460 delta_school 7.4182089 delta_shopping 3.8043597 delta_private 4.2786317 delta_leisure 3.4000778 delta_work_FT 1.3248514 delta_work_wknd 2.8608873 delta_school_young 2.3441172 delta_leisure_wknd 0.2948808 Scaling used in computing Hessian --------------------------------- Value alpha_base 14.7608915 gamma_work 4.8990249 gamma_school 3.0984826 gamma_shopping 0.4264320 gamma_private 0.6200096 gamma_leisure 2.0956978 delta_work 3.7172370 delta_school 7.4182123 delta_shopping 3.8043560 delta_private 4.2786312 delta_leisure 3.4000873 delta_work_FT 1.3248517 delta_work_wknd 2.8608874 delta_school_young 2.3441173 delta_leisure_wknd 0.2948808 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() ### Define individual alternatives alternatives = c("outside", "work", "school", "shopping", "private", "leisure") ### Define availabilities avail = list(outside = 1, work = 1, school = 1, shopping = 1, private = 1, leisure = 1) ### Define continuous consumption for individual alternatives continuousChoice = list(outside = t_outside/60, work = t_a02/60, school = t_a03/60, shopping = t_a04/60, private = t_a05/60, leisure = t_leisure/60) ### Define utilities for individual alternatives V = list() V[["outside"]] = 0 V[["work"]] = delta_work + delta_work_FT * occ_full_time + delta_work_wknd * weekend V[["school"]] = delta_school + delta_school_young * (age<=30) V[["shopping"]] = delta_shopping V[["private"]] = delta_private V[["leisure"]] = delta_leisure + delta_leisure_wknd*weekend ### Define alpha parameters alpha = list(outside = 1 /(1 + exp(-alpha_base)), work = 1 /(1 + exp(-alpha_base)), school = 1 /(1 + exp(-alpha_base)), shopping = 1 /(1 + exp(-alpha_base)), private = 1 /(1 + exp(-alpha_base)), leisure = 1 /(1 + exp(-alpha_base))) ### Define gamma parameters gamma = list(work = gamma_work, school = gamma_school, shopping = gamma_shopping, private = gamma_private, leisure = gamma_leisure) ### Define costs for individual alternatives cost = list(outside = 1, work = 1, school = 1, shopping = 1, private = 1, leisure = 1) ### Define settings for MDCEV model mdcev_settings <- list(alternatives = alternatives, avail = avail, continuousChoice = continuousChoice, utilities = V, alpha = alpha, gamma = gamma, sigma = sig, cost = cost, budget = 24) ### Compute probabilities using MDCEV model P[["model"]] = apollo_mdcev(mdcev_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) }