Population EVPI

Data source:

Quenneville SP, Xie X, Brophy JM. The cost-effectiveness of Maze procedures using ablation techniques at the time of mitral valve surgery. Int J Technol Assess Health Care 2009; 25(4):485-96 Link: http://www.ncbi.nlm.nih.gov/pubmed/19818194

SAS:

%let pa=%str(AAAAAAAA\MAZE-model\Inputs for Markov model-new.xls);

proc import

    datafile="&pa"

    out=aa

    DBMS=xls REPLACE;

    getnames=yes;

      RANGE="EVPI graph$a18:e29";

run;

 

 

*/ format data;

data cc;

    set aa ;

      label years_5= 'Five years'

          years_10=  'Ten years'

          years_15= 'Fifteen years'

          years_20= 'Twenty years';

run;

 

proc format;

   value WTP     10000='10K'

                 20000='20K'

                 30000='30K'

                 40000='40K'

                 50000='50K'

                 60000='60K'

                     70000='70K'

                 80000='80K'

                     90000='90K'

                 100000='100K';

     value EVPI  500000='0.5'

                 1000000='1.0'

                 1500000='1.5'

                 2000000='2.0'

                 2500000='2.5'

                 3000000='3.0'

                 3500000='3.5'

                 4000000='4.0';

run;

 

*/ draw graph, href=50000 vref=400000;

goptions reset=all;

proc gplot data=cc;

      plot (years_5 years_10 years_15 years_20)*wtp/

    overlay legend=legend1 haxis=axis1 vaxis=axis2;

    format wtp wtp.;

    format  years_5 years_10 years_15 years_20 EVPI.;

    symbol1  i=spline ci=black    v=square  h=1.5 w=2 l=5;

    symbol2  i=spline ci=red   v=diamond h=1.5 w=2 l=3;

      symbol3  i=spline ci=blue   v=dot     h=1.5 w=2 l=1;

      symbol4  i=spline ci=green  v=star    h=1.5 w=2 l=2;

    axis1 label=(h=2 f=Arial  'Threshold willingness to pay for a QALY (CAD $)') width=1 length=100

          order=(0 to 100000 by 10000) value=(f='Times New Roman' h=16pt) minor=none;

    axis2 label=(h=2  f=Arial  a=90 'Population EVPI (million CAD $)') width=1 length=34.4

          order=(0 to 4000000 by 500000 )  value=(f='Times New Roman' h=16pt) minor=none;

      legend1 across=1 down=4 position=(top left inside)label=(f='Times New Roman' h=16pt 'Time Horizon' justify=center)

          frame value=(f='Times New Roman' h=16pt)  mode=share shape=line(15) shape=SYMBOL(15, 1.2) ;

run;

quit;