Statistics 板


LINE

※ 引述《Brunhilde (不见。)》之铭言: : 刚学macro的语法 : 有点问题想请教一下 : 我有设定了 : %do i=0 %to 1; : 可是只跑出了0的部份 : 以下为全部语法 谢谢 : data y n;set HY.taidd96_06480487; : if ICD9_2 in ('493','490','491','492','494','496') then output y; : else output n; : run; : %macro aa; : data L0;Set HY.taipollcli;day=date; ^^^^ : %do i=0 %to 1; : L0&i.PM10=lag&i(PM10); : L0&i.SO2=lag&i(SO2); : L0&i.NO2=lag&i(NO2); : L0&i.CO=lag&i(CO); : L0&i.O3=lag&i(O3); : if 1996 le year(date) le 2006 then output; : if temp<23 then output; : drop PM10 CO SO2 NO2 O3 stno; : format day date9.; : run; ^^^^ data 和run 是一组的 也就是说,虽然你在最後面加了对应%do的%end 但是sas在这边就会把你的%do卡掉 你下面的程式我没有细看,所以不确定你想做什麽 如果要把整个内容用i=0 和1各执行一便的话 可以把%do那一行放到data上面那一行试试看 : proc sort data= L0; : by day;run; : proc sort data=N;by IN_DATE;run; : data b;set N; day=IN_date; : BY IN_DATE; : IF FIRST.IN_DATE THEN DO; : NO=0; : END; : NO+1; : IF LAST.IN_DATE; : if 1996 le year(IN_date) le 2006 then output; : format day date9.; : run; : data ev; merge b(in=e) L0;by day; : if e=1; : run; : data event; set ev; : eventday=In_date; : format eventday date9.; : Run; : data ref1; set event; : day=eventday+7; EVEN=2; : drop L0&i.PM10 L0&i.CO L0&i.SO2 L0&i.NO2 L0&i.O3 RH temp; : format day date9.;run; : proc sort data=ref1; : by day;run; : data ref11; merge ref1(in=e) L0; : by day;if e=1;case=0; : ;run; : data ref2; set event; : day=eventday-7; EVEN=2; : drop L0&i.PM10 L0&i.CO L0&i.SO2 L0&i.NO2 L0&i.O3 RH temp; : format day date9.;run; : proc sort data=ref2; by day;run; : data ref22; merge ref2(in=e) L0; by day;if e=1;case=0; : run; : data event0;set event; : drop L0&i.PM10 L0&i.CO L0&i.SO2 L0&i.NO2 L0&i.O3 RH temp; : day=eventday; EVEN=1;format day date9.;run; : data event00;merge event0(in=e) L0; by day; : ;if e=1;case=1;run; : data test; set event00 ref11 ref22; run; : proc sort data=test;by day;run; : data test1;merge test(in=e) b;by day; : if e=1; : run; : data test2;set test1; : a&i.PM10=L0&i.PM10/24.99 ; : a&i.SO2=L0&i.SO2/2.42; : a&i.NO2=L0&i.NO2/8.67; : a&i.CO=L0&i.CO/0.37; : a&i.O3=L0&i.O3/12.2; : Run; : /*单变项ties case(0) TIES=BRESLOW weight*/ : proc phreg data=test2; : model even*case(0)=a&i.PM10 /rl TIES=BRESLOW;strata eventday; weight no;run; : proc phreg data=test2; : model even*case(0)=a&i.SO2 /rl TIES=BRESLOW; strata eventday; weight no;run; : proc phreg data=test2; : model even*case(0)=a&i.NO2 /rl TIES=BRESLOW ;strata eventday; weight no;run; : proc phreg data=test2; : model even*case(0)=a&i.CO /rl TIES=BRESLOW; strata eventday; weight no;run; : proc phreg data=test2; : model even*case(0)=a&i.O3 / rl TIES=BRESLOW; strata eventday; weight no; run; : %end; : %mend aa; : %aa; -- (C)lass 【 分组讨论区 】 11 卡漫梦工厂 卡通, 漫画, 动画 [Ckmike] 4 C_J_Artist 作家 Σ日本漫画家专区 [Ckmike] 10 TakahasiShin 作家 ◎****高桥真作品讨论版**** kenshin528 http://kenshinsasstory.blogspot.com/ SAS Story --



※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.119.47.118
1F:推 Brunhilde:thanks 11/04 13:50







like.gif 您可能会有兴趣的文章
icon.png[问题/行为] 猫晚上进房间会不会有憋尿问题
icon.pngRe: [闲聊] 选了错误的女孩成为魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一张
icon.png[心得] EMS高领长版毛衣.墨小楼MC1002
icon.png[分享] 丹龙隔热纸GE55+33+22
icon.png[问题] 清洗洗衣机
icon.png[寻物] 窗台下的空间
icon.png[闲聊] 双极の女神1 木魔爵
icon.png[售车] 新竹 1997 march 1297cc 白色 四门
icon.png[讨论] 能从照片感受到摄影者心情吗
icon.png[狂贺] 贺贺贺贺 贺!岛村卯月!总选举NO.1
icon.png[难过] 羡慕白皮肤的女生
icon.png阅读文章
icon.png[黑特]
icon.png[问题] SBK S1安装於安全帽位置
icon.png[分享] 旧woo100绝版开箱!!
icon.pngRe: [无言] 关於小包卫生纸
icon.png[开箱] E5-2683V3 RX480Strix 快睿C1 简单测试
icon.png[心得] 苍の海贼龙 地狱 执行者16PT
icon.png[售车] 1999年Virage iO 1.8EXi
icon.png[心得] 挑战33 LV10 狮子座pt solo
icon.png[闲聊] 手把手教你不被桶之新手主购教学
icon.png[分享] Civic Type R 量产版官方照无预警流出
icon.png[售车] Golf 4 2.0 银色 自排
icon.png[出售] Graco提篮汽座(有底座)2000元诚可议
icon.png[问题] 请问补牙材质掉了还能再补吗?(台中半年内
icon.png[问题] 44th 单曲 生写竟然都给重复的啊啊!
icon.png[心得] 华南红卡/icash 核卡
icon.png[问题] 拔牙矫正这样正常吗
icon.png[赠送] 老莫高业 初业 102年版
icon.png[情报] 三大行动支付 本季掀战火
icon.png[宝宝] 博客来Amos水蜡笔5/1特价五折
icon.pngRe: [心得] 新鲜人一些面试分享
icon.png[心得] 苍の海贼龙 地狱 麒麟25PT
icon.pngRe: [闲聊] (君の名は。雷慎入) 君名二创漫画翻译
icon.pngRe: [闲聊] OGN中场影片:失踪人口局 (英文字幕)
icon.png[问题] 台湾大哥大4G讯号差
icon.png[出售] [全国]全新千寻侘草LED灯, 水草

请输入看板名称,例如:Tech_Job站内搜寻

TOP