Hello everyone; this model has estimated the alpha and gamma parameters from the experimental data using the MLE estimation procedure. my study consisted of 121 respondents. I wanted these parameters for each individual, but I am unable to get this. Can anyone help to understand how I can estimate parameters (alpha and gamma) for each respondent? please help me
thank you
* define prospect theory model following Harrison and Rutstrom (2008)
program define MLcpt
args lnf alpha gamma
tempvar prob0l prob1l prob0r prob1r y0 y1 y2 y3
tempvar euL euR euDiff euRatio tmp
quietly {
generate double `tmp' = (($ML_y2^`gamma')+((1-$ML_y2)^`gamma'))
replace `tmp' = `tmp'^(1/`gamma')
generate double `prob0l' = ($ML_y2^`gamma')/`tmp'
generate double `prob1l' = ((1-$ML_y2)^(`gamma'))/`tmp'
replace `tmp' = (($ML_y3^`gamma')+((1-$ML_y3)^`gamma'))
replace `tmp' = `tmp'^(1/`gamma')
generate double `prob0r' = ($ML_y3^`gamma')/`tmp'
generate double `prob1r' = ((1-$ML_y3)^(`gamma'))/`tmp'
generate double `y0' = .
replace `y0' = ( $ML_y4)^(`alpha') if $ML_y4>=0
replace `y0' = -`lambda'*(-$ML_y4)^(`alpha') if $ML_y4<0
generate double `y1' = .
replace `y1' = ( $ML_y5)^(`alpha') if $ML_y5>=0
replace `y1' = -`lambda'*(-$ML_y5)^(`alpha') if $ML_y5<0
generate double `y2' = .
replace `y2' = ( $ML_y6)^(`alpha') if $ML_y6>=0
replace `y2' = -`lambda'*(-$ML_y6)^(`alpha') if $ML_y6<0
generate double `y3' = .
replace `y3' = ( $ML_y7)^(`alpha') if $ML_y7>=0
replace `y3' = -`lambda'*(-$ML_y7)^(`alpha') if $ML_y7<0
gen double `euL'=(`prob0l'*`y0')+(`prob1l'*`y1')
gen double `euR'=(`prob0r'*`y2')+(`prob1r'*`y3')
generate double `euDiff' = `euR' - `euL'
replace `lnf' = ln(normal( `euDiff')) if $ML_y1==1
replace `lnf' = ln(normal(-`euDiff')) if $ML_y1==0
}
end
ml model lf HL_cpt0 (alpha: choices prob1l prob2l prize0 prize1 prize2 prize3= )(gamma: ), cluster(ID) technique(nr) maximize
Related Posts with Estimating parameters for each respondent in the given model
please help. making bar graph more simple with "over" optionsHi, guys.
I'm a newcomers with stata, and there is some trouble in making bar graph.
My data has 4…
Please help cant work out how to t-test. Im new to stataHi guys I'm struggling to compare 2 different regions with a t-test.
As seen from the screenshots I…
Adding rows under a variableStata command for adding rows lets say 10 rows under a variable which has fixed observations like 40…
Using cii proportions with loopHello,
I am trying to compute confidence intervals for proportions (number of cases / total) on eac…
Pseudo Panel Data and Mediation AnalysisDear Stata Experts
I'm a PhD student, and I'd like to know how to use the codes to perform the Medi…
Subscribe to:
Post Comments (Atom)
0 Response to Estimating parameters for each respondent in the given model
Post a Comment