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
Dealing with date and timeHi, I am using Stata 16 and seem to have problem extracting year, month and day from a date-time co…
Interpolate to end of month instead of mid month data?Hello, all my data are for the end of month. But I have one variable where the data is listed mid mo…
Meta-analysis: how to use fixed-effect and random-effect model for different subgroupsHi all, I am new to STATA. I was previously told that subgroup analysis by sex can be performed usin…
strtrim not producing desired resultsHi Stata Users, I am using strtrim command in Stata 15. The goal is to remove blank spaces at the b…
LinktestHey guys, I have a question about the output of the linktest. I know that a significant _hatsq in c…
Subscribe to:
Post Comments (Atom)
0 Response to Estimating parameters for each respondent in the given model
Post a Comment