[Pengpeng Ye]

I totally missed this follow-up question.

Dear Carlo,

Thanks for your reply. Are the codes above appropriate for the questions proposed by Elena?

Best wishes,
Pengpeng

[Pengpeng Ye]

Dear Carlo,

I am also conducting a Cost-Effectiveness Analysis of data from a randomised controlled trial, I have gotten to estimating the ICER (via a mixed multilevel model). The next step is to have my results on a Cost-Effectiveness Plane.
As I can see from your post above, I need to bootstrap the cost and effect and replicate 5,000.
1) is it possible to do it in one line of code (can I estimate 5000 ICERs, or do I need to estimate costs and effects separately) ?
2) is it possible to then plot them across the 4 quadrants?
3) is it possible to also estimate the probability of CE at different willingness to pay thresholds?
And to do all of this in Stata, or do I need to use other software?

Thank you very much.

Elena

Dear Elena,

I found the answer of your 3rd question from Carlo Lazzaro in other post named "Code for Cost-Effectiveness Acceptability Curve" started by Marta Wilson-Barthes.

I think that all your questions could be done in Stata.

// from Carlo Lazzaro
clear
range Threshold_values 0 200000 1000

g Delta_C_b=runiform()*1000
g Delta_E_b=runiform()*10
g double acc=.

forval i = 1/1000 {
qui g t`i' =Threshold_values[`i']* Delta_E_b> Delta_C_b
qui sum t`i', meanonly
replace acc=r(mean) in `i'
drop t`i'
 }

twoway (scatter acc Threshold_values)

[Carlo Lazzaro]

Pengpeng:
yes, I think they are (at forst glance, at least).
Those interested in advanced economic evaluation of health care programmes, can find the following link useful: https://www.herc.ox.ac.uk/downloads/...linical-trials

[Phyu Zin]

Hi Pengpeng

I got errors in running the below " local ratio"
Is it because of the ' ' or " "?
local ratio=`cdiff' / `ediff' di `cdiff' " " `ediff' " " `ratio' return scalar cdf=`cdiff' return scalar edf=`ediff' return scalar r=`ratio' restore

[Renira Corinne Angeles]

Dear Carlo and all Statalist members

This thread has been very helpful for me. I am doing an economic evaluation with trial data, and I need to plot a cost-effectiveness plane. I have followed Carlo's code (post #2). So now I have: one .dta file with boostrapped effect-difference and another with bootstrapped cost-difference. How do I make a cost-effectiveness plane figure (scatterplot)? Shouldn't the boostrapped effect- and costdiff be on the same file?

Thanks in advanced.

Renira

[Carlo Lazzaro]

Renira:
yes, you should have both incremental cos and effectiveness in the same dataset.