sas医学统计案例程序-经管之家官网！

[UseMoney=20]sas医学统计案例程序：差异检验 transpose ods html NPAR1WAY means freq gplot corr macro libname xls
/*本程序都是我跑出来的，拥有 此程序可以让你sas 编程进步飞快
可以登录我的msn查看更多的SAS CODE:
http://cid-7a2c11bd1414b988.spaces.live ... !208.entry
*/

/*先导入每个表单：s1-对照组；s2-间重病例；s3-持重病例;
s11-对照组人员背景；s12-对照组鼻塞,流涕，喷嚏，鼻痒，眼痒，流泪；
s13-28个变量；s14-SF36中8个变量；
SF-36量表。此量表共包括36个条目,分
为8个维度,即生理功能(physical function,PF)、生理职能(re-
sponsibility of physical,RP)、躯体疼痛(body pain,BP)、总体健
康(general health,GH)、活力(vitality,VT)、社会功能(social
function,SF)、情感职能(responsibility of emotion,RE)和精神
健康(mental health,MH)
*/
options validvarname=any;
%macro inputdata;/*导入三个表*/
%do i=1 %to 3;
%let ref="E:\lixin\q&i..xls";/**/
libname lib &ref mixed=yes; *dbmax_text=50;/*虽然有mixed但是要保证exl前八个记录不能全是数值，而且不需要dbmax_text*/
%do j=1 %to 4;/*办法二：LibName xlsLib "H:\ExcelToSas\Demo.xls" mixed=yes ;
Missing M ;
Data Sheet3( drop=_: ) ;
Set xls."Sheet3$"n( reName=( x=_x y=_y ) ) ;
x = input( _x , best. ) ;
y = input( _y , best. ) ;
Run ;
LibName xls clear ;*/
%let sheet="sheet&j$"n;
data biyan.s&i&j;
set lib.&sheet;
run;
%end;
libname lib clear;
%end;
%mend ;
%lixin;
%macro setli;/*按表单合并，共4个*/
%do i=1 %to 4;
%let dui=s1&i;
%let jian=s2&i;
%let chi=s3&i;
%put &dui;
data su&i;
*format class $ ;/*三个组别：对照 间重 持重*/
set biyan.&dui(in=x) biyan.&jian(in=y) biyan.&chi(in=z);
if x=1 then class='A'; /*if不需要加%*/
if y=1 then class='B';
if z=1 then class='C';
run;
%end;
%mend inputdata;
%inputdata;


/*一，两样本或三样本差异检验：持重和间重
在总体分布形式难确定或总体分布为偏态时，不适用T检验，故选用非参数统计方法。使用*/
%macro NPA();/*title显示还有个问题：我想显示把表的名字也显示出来，改如何，因为这个结果看不出来是比较哪个表*/
filename odsout 'e:\lixin\htmls';
ods listing close;
ods html body='CORR&NPA_mapbody.html'
contents='CORR&NPA_contents.html'
frame='相关分析和差异检验_frame.html'
/*nogtitle*/
path=odsout;
/*nogtitle 添加此选项以后，可以自己定义目录名字，具体参考Gmap：Example 5*/
%let title1="间重和持重";
%let title2="对照和持重";
%let title3="对照和间重";
%do i=2 %to 4;
%let s=su&i;
%do j=1 %to 3;
%if &j=1 %then %let zu='A' ; /*注意此处要两个let,而且if 在数据步外，所以要加%if*/
%if &j=2 %then %let zu='B' ;
%if &j=3 %then %let zu='C' ;
/*不同组差异比较*/
proc NPAR1WAY Wilcoxon anova DATA=&s(drop=编号);
class class;
where class^=&zu;
title &&title&j 差异比较;
output out=out;
run;

/*同组相关分析*/
proc corr data=&s(drop=编号) NOSIMPLE;
where class=&zu;
title COLOR=RED height=20pt &&title&j 相关分析;
run;
data plot;
set out;
n+1;
keep _VAR_ PT2_WIL n;
proc print ;/*找出差异显著的项*/
where PT2_WIL<=0.05;
title &&title&j 显著项;/*title可以不用连字符来连接字符串*/
run;
goptions reset=global gunit=pct border cback=white
colors=(blue) htitle=4;
symbol color=RED interpol=none width=2 value=dot height=1;
axis1 label=( angle=90 height=3 "Two-sided P value") ;

proc gplot data=plot;
plot PT2_WIL*_VAR_/vref=(0.05 0.01) vaxis=axis1;
title &&title&j 显著项;
run;
quit;
%end;
%end;
ods html close;
ods listing;
%mend;
%NPA;



/*按性别组间检测*/
%macro sex(zu);
data sex;
merge su1(keep=性别) su4 ;
run;
proc NPAR1WAY Wilcoxon anova DATA=sex ;
class 性别;
where class=&zu;title "&zu";/*宏在双引号才有效*/
output out=wilsex;
run;
proc print data=wilsex;
var PT2_WIL _VAR_;
where PT2_WIL<=0.05;
run;

%mend;
%sex('C');/*参数可以改为B，A*/



%MACRO MEAN;
filename odsout 'e:\lixin\htmls';
ods listing close;
ods html body='mean_mapbody.html'
contents='mean_contents.html'
frame='均值图.html'
/*nogtitle*/
path=odsout;
%let title2="VAS";
%let title3=" RQLQ";
%let title4="SF-36";
%do k=2 %to 4;
%LET DATASET=su&k;
proc means data=&dataset(drop=编号) maxdec=1 mean;/*maxdec=3设定均值小数位数*/
class class;
*title 'mean 求均值';
output OUT=outmeans;
title &&title&k 均值分析;
run;
data meansu2;
set outmeans;
if _type_=1 & _STAT_="MEAN";
drop _type_ _freq_;
run;
PROC transpose data=meansu2 out=trans;
by class;
run;
proc sort data=trans out=bymean;
by _name_;
run;
data chart;
set bymean;format col1 3.1;
label _name_="症状" col1="均值";
run;
goptions reset=global colors=(red) ;
pattern color=blue;
title &&title&k 均值图;
proc gchart data=chart;
vbar3d class/discrete/*分组显示*/
sumvar=col1
inside=mean
group=_NAME_;
run;quit;
%end;
ods html close;
ods listing;
%MEND;
%MEAN;

%macro top();/*RQLQ活动影响最大的*/
filename odsout 'e:\lixin\htmls';
ods listing close;
ods html body='top_mapbody.html'
contents='top_contents.html'
frame='活动前三项.html'
/*nogtitle*/
path=odsout;
data top;
merge su1 su3;
%do i=1 %to 3;
act&i=compress(活动&i,,'ds');
%end;
keep 编号 class act1 act2 act3;/*1 活动2 活动3 活动1*/
by class;
run;
data top1;
set top;
n+1;/*构建一个没有相同值的观察向量n*/
select(class);
when('A') 组别='对照';
when('B') 组别='间重';
when('C') 组别='持重';
end;
run;
proc transpose data=top1
out=toptras(rename=(col1=活动));
var act1-act3 ;/*将此三个变量合并为一列，需要一个没有相同值的向量n，*/
by class n 组别 ;run;

proc freq data=toptras order=freq ;
tables 活动*组别/NOCOL NOROW;
*where class='B';
title 'RQLQ活动影响排名' ;
run;
ods html close;
ods listing;
%mend;
%top;

[/UseMoney]
爱萌 金币+5 魅力+15 经验+10 奖励公布了程序,2009-5-21 0:42:15