STATA-FORUM.DE

[ann-kathrin_fisc12]

Hallo,

ich verwende ein mixed logit modell für die Berechnung von Koeffizienten in einer DCE. Meine Variablen sind Effects-coded, außer das Kosten-Attribut zur Berechnung der WTP.

Ich bekomme die Koeffizienten für jeden Level eines Attributs außer dem letzten Level. Diese lasse ich mit dem nlcom Befehl berechnen.

Mit dem Befehl WTP möchte ich die WTP und die CI (95%) für jedes Level meiner DCE berechnen. Dafür gebe ich dem Befehl die Namen der Variablen, die verwendet werden müssen, aus der mixed logit Berechnung. Das funktioniert wunderbar. Meine Frage ist, wie ich den WTP Befehl mit der Berechnung aus nlcom verknüpfe. Wie kann ich in der Syntax wtp angeben, dass die Werte der nlcom Berechnung genutzt werden sollen?

Hat jemand eine Ahnung?

[Staxa]

Keine der genannten Abkürzungen sagt mir was. Bitte doch mal Beispielcode bzw. Output posten, was hier genau gemacht wird und warum Effektkodierung ist in Stata eigl gar nicht explizit notwendig, wie wurde das gemacht?

[ann-kathrin_fisc12]

gmnl choice, rand(v1-v22) group(number) id(rid) nrep(500) mixl dtype(3) seed(1234) -> 1. Syntax

Output:
choice Coefficient Std. err. z P>z [95% conf. interval]

Mean
v1 .12 .1573696 2.33 0.020 .0578835 .6747608
v2 .337 .1554955 2.28 0.023 .0495715 .6591025
v3 -.47 .2205666 -0.51 0.608 -.5455373 .3190679
v4 -.32 .1393158 -2.30 0.022 -.5930877 -.0469796
v5 .0009 .1517926 2.09 0.037 .0191502 .6141663
v6 -.731 .1478734 -1.75 0.080 -.5485673 .0310856
v7 .038 .3032876 0.13 0.899 -.556073 .6327923
v8 -.123 .1181181 -0.85 0.397 -.3316079 .1314066
v9 .111 .0736686 8.82 0.000 .5053547 .7941302
v10 -.909 .1175434 -1.85 0.064 -.4480762 .0126855
v11 .777 .1209047 5.26 0.000 .3988834 .8728212
v12 .653 .1591982 1.80 0.072 -.0252396 .5988057
v13 .907 .4264255 0.37 0.711 -.6777331 .9938241
v14 -.123 .1821828 -1.46 0.144 -.6234036 .0907399
v15 .56 .0957448 6.19 0.000 .4047461 .780059
v16 -.1 .2027354 -0.23 0.820 -.443544 .3511643
v17 .002 .1372006 0.31 0.753 -.2257573 .312059
v18 -.066 .0018181 -18.07 0.000 -.0364145 -.0292878
v19 1.99 .1257158 11.57 0.000 1.208362 1.701159
v20 .10 .1429064 5.56 0.000 .5146153 1.074798
v21 -.07 .1649075 -0.30 0.765 -.3724362 .2739895
v22 -.09 .210808 -0.40 0.687 -.4981315 .3282207

SD
v1 .12 .1573696 2.33 0.020 .0578835 .6747608
v2 ......

(Das ist nur ein Beispieloutput mit veränderten Koeffizienten; SD Tabelle unvollständig dargestellt)

v18 = lin. cod. Kostenattribut
v1-v4; v5-v8; v9-v10; v11-v14; v15-v17; v19-v22 sind jeweils effects codierte Level eines Attributes aus einem Discrete Choice Experiment. Hier fehlt in der berechneten Tabelle jeweils das letzte "Level" eines Attributes/ ein Koeffizient, der mit folgender Syntax berechnet wird:

. nlcom -1*(_b[v1] + _b[v2] + _b[v3] + _b[v4])

_nl_1: -1*(_b[v1] + _b[v2] + _b[v3] + _b[v4])

------------------------------------------------------------------------------
choice | Coefficient Std. err. z P>|z| [95% conf. interval]
-------------+----------------------------------------------------------------
_nl_1 | .133443 .0768008 2.00 0.046 .00228 .2489534
------------------------------------------------------------------------------
Im Anschluss wird folgendes Berechnet für die Willingness to Pay:


. wtp v18 v1

wtp wird für jeden Koeffizienten berechnet.
v18 bleibt konstant.

Outcome ist die berechnete willingess to pay und 95% CI.

-> Mein Problem ist, dass ich das "fehlende Level" per nlcom Befehl nachrechne. Ich müsste diesen Wert in der Syntax wtp wiederverwenden. -> wtp v18 [WERT AUS NLCOM]

[Staxa]

Das Problem hier ist wohl, dass du 2 user written ados nutzt, da weder gmnl noch wtp offizielle Stata Befehle sind. Offenbar erwartet wtp die Ergebnisse in einem ganz bestimmten Format, die nlcom so nicht liefert bzw. im Speichern belässt. Am ehesten würde ich mich an die Autoren von wtp wenden, ob man da etwas machen kann oder wie man den Befehl auch mit nlcom nutzt. Ein ähnliches, aber bisher ungelöstes Problem, zeigt sich etwa hier: https://www.statalist.org/forums/forum/ ... culate-wtp Vielleicht findest du dort einen Hinweis.

[ann-kathrin_fisc12]

Ok danke!