An example of 3-D plot

Author: Xuanqian Xie

This plot was prepared using data from a HTA report of Subglottic Secretion Drainage Endotracheal Tubes (Xie et al 2011).
 

 

 

 

options nodate nonumber pageno=1 ps=200 ls=100 formdlim='-';
options mprint mlogic symbolgen;

%let pa=%str(C:\Documents and Settings…..\rr_c_vap_winbugs.xls);
proc import
    datafile="&pa"
    out=cc
    DBMS=EXCEL2000 REPLACE;
    getnames=yes;
    sheet="Sheet1";
run;

proc contents data=cc;
run;

proc univariate data=cc normal;
   var  rr cost_vap ;
run;

%macro cea(risk_1000d=, vet_days=);
data cea;
set cc;
price_tube=18.5;
pre_vap= &risk_1000d* &vet_days;

cost_ctl = cost_vap* pre_vap;
cost_ssd = cost_vap* pre_vap*rr + price_tube;
IC = cost_ssd - cost_ctl;
run;

proc univariate data=cea normal;
   var  IC;
   OUTPUT out=IC  pctlpre=P_ pctlpts=2.5,50,97.5;
run;

data icf;
set IC;
risk=&risk_1000d;
days=&vet_days;
run;

%mend;

%cea (risk_1000d=0.0083, vet_days=10);



data CEA_results;
run;

%MACRO Update;
%do dd= 5 %to 15;
%do nn= 0 %to 100 %by 5;
%let bb=&nn/10000;
%put &bb;
%put &nn;
%cea (risk_1000d=&bb, vet_days=&dd);

/*updat dataset*/
data CEA_results;
    set CEA_results icf;
run;

%end;
%end;
%Mend;

%Update;

 /* Define a title for the plot */
title1 '';

 /* Create the plot */
proc g3d data=CEA_results;
   plot risk*days=P_50 / grid rotate=315;
run;
quit;

proc g3d data=CEA_results;
scatter risk*days=P_50;
run;
 

Prev productPlot the data precisely?
Next product:No next product