- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 4511 个
- 通用积分
- 76.8946
- 学术水平
- 171 点
- 热心指数
- 204 点
- 信用等级
- 133 点
- 经验
- 24578 点
- 帖子
- 583
- 精华
- 1
- 在线时间
- 1865 小时
- 注册时间
- 2008-5-25
- 最后登录
- 2025-4-28
|
- %macro comment(); /*这是一种比较特殊的注释方式*/
- Picture: LFT Safety Panel, Baseline versus Study
- (using template and sgrender procedure)
- Software: SAS 9.4 TS Level 1M2
- Name: justsoso
- Date: 20180305
- From: Getting Started with the Graph Template Language in SAS:
- Examples, Tips, and Techniques for Creating Custom Graphs
- History: 20180305;
- %mend;
- /*--Fig_5_5_2 LFT Lattice--*/
- data GTL_GS_LFT_Lattice (keep=drug alat biltot alkph asat
- palat pbiltot palkph pasat visitnum);
- label alat="ALAT (/ULN)";
- label biltot="BILTOT (/ULN)";
- label alkph="ALKPH (/ULN)";
- label asat="ASAT (/ULN)";
- visitnum=1;
- do i= 1 to 100;
- palat = min (4, 2.5 * (abs(rannor(123))) / 3.0);
- pbiltot = min (4, 2.5 * (abs(rannor(123))) / 3.0);
- palkph = min (4, 2.5 * (abs(rannor(123))) / 3.0);
- pasat = min (4, 2.5 * (abs(rannor(123))) / 3.0);
- alat = min (4, 2.5 * (abs(rannor(345))) / 3.0);
- biltot = min (4, 2.5 * (abs(rannor(345))) / 3.0);
- alkph = min (4, 2.5 * (abs(rannor(345))) / 3.0);
- asat = min (4, 2.5 * (abs(rannor(345))) / 3.0);
- j = rannor(345);
- if j > 0 then drug = "A";
- else drug="B";
- output;
- end;
- visitnum=2;
- do i= 1 to 100;
- palat = min (4, 2.5 * (abs(rannor(789))) / 3.0);
- pbiltot = min (4, 2.5 * (abs(rannor(789))) / 3.0);
- palkph = min (4, 2.5 * (abs(rannor(789))) / 3.0);
- pasat = min (4, 2.5 * (abs(rannor(789))) / 3.0);
- alat = min (4, 2.5 * (abs(rannor(567))) / 3.5);
- biltot = min (4, 2.5 * (abs(rannor(567))) / 3.5);
- alkph = min (4, 2.5 * (abs(rannor(567))) / 3.5);
- asat = min (4, 2.5 * (abs(rannor(567))) / 3.5);
- j = rannor(567);
- if j > 0 then drug = "A";
- else drug="B";
- output;
- end;
- visitnum=3;
- do i= 1 to 100;
- palat = min (4, 2.5 * (abs(rannor(321))) / 3.0);
- pbiltot = min (4, 2.5 * (abs(rannor(321))) / 3.0);
- palkph = min (4, 2.5 * (abs(rannor(321))) / 3.0);
- pasat = min (4, 2.5 * (abs(rannor(321))) / 3.0);
- alat = min (4, 2.5 * (abs(rannor(975))) / 2.5);
- biltot = min (4, 2.5 * (abs(rannor(975))) / 2.5);
- alkph = min (4, 2.5 * (abs(rannor(975))) / 2.5);
- asat = min (4, 2.5 * (abs(rannor(975))) / 2.5);
- j = rannor(975);
- if j > 0 then drug = "A";
- else drug="B";
- output;
- end;
- run;
- proc format;
- value wk
- 1='1 Week'
- 2='3 Months'
- 3='6 Months';
- value lab
- 1='ALAT'
- 2='Bilirubin Total'
- 3='Alk Phosphatase'
- 4='ASAT';
- value $trt
- "A"="Drug A (n=240)"
- "B"="Drug B (n=195)";
- run;
- data GTL_GS_LFT_Lattice2 (keep=visitnum drug labtest result pre);
- format visitnum wk. labtest lab. drug $trt.;
- label pre='Baseline (/ULN) *' result='Study (/ULN)';
- set GTL_GS_LFT_Lattice;
- pre=palat;
- labtest=1;
- result=alat;
- output;
- pre=pbiltot;
- labtest=2;
- result=biltot;
- output;
- pre=palkph;
- labtest=3;
- result=alkph;
- output;
- pre=pasat;
- labtest=4;
- result=asat;
- output;
- run;
- ods _all_ close; /*关闭所有ods destinations,避免不必要的输出,节省计算机资源*/
- /*ods listing gpath='C:\Users\Administrator\Desktop\test20180305' image_dpi=300;*/
- ods listing gpath='C:\Users\Administrator\Desktop\test20180305';
- /*打开ods listing这一destination,指定图片保存路径,★★★★★需提前建立*/
- /*ods graphics on / imagename="picture" outputfmt=pdf;*/
- /*图片名称picture,输出类型pdf*/
- ods graphics off;
- /*LFT Safety Panel, Baseline versus Study*/
- proc template;
- define statgraph Fig_5_5_2_LFT_Lattice;
- begingraph;
- entrytitle 'LFT Safety Panel, Baseline vs. Study';
- layout datalattice rowvar=visitnum columnvar=labtest /
- headerlabeldisplay=value
- rowaxisopts=(linearopts=(viewmin=0 viewmax=4))
- columnaxisopts=(linearopts=(viewmin=0 viewmax=4));
- layout prototype;
- scatterplot x=pre y=result / group=drug;
- referenceline x=1 / lineattrs=(pattern=dash);
- referenceline x=1.5 / lineattrs=(pattern=dash);
- referenceline x=2 / lineattrs=(pattern=dash);
- referenceline y=1 / lineattrs=(pattern=dash);
- referenceline y=1.5 / lineattrs=(pattern=dash);
- referenceline y=2 / lineattrs=(pattern=dash);
- endlayout;
- endlayout;
- endgraph;
- end;
- run;
- proc sgrender data=GTL_GS_LFT_Lattice2 template=Fig_5_5_2_LFT_Lattice;
- run;
- ods _all_ close;
- /*ods html; 最好将ods destination恢复至默认的html*/
复制代码
|
|