应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 学术交流区 » 论文投稿 » 投稿求助 » 审稿意见中的一段SAS程序如何理解
51/1返回列表
查看: 723  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

myselfsky8159

木虫 (著名写手)

[求助] 审稿意见中的一段SAS程序如何理解

审稿意见很中肯，还对数据处理编写了一段SAS程序，请高手帮忙分析下这段程序。
--------------------------------------------------------
data wr18606;
     attrib sitetype length=$30;
     input sitetype n@@;
     do i= 1 to n;
        input value@@;
            if value =0 then do; low=.; high=.3; end;
            else do; low=value; high=value; end;
                output;
         end;
cards;
UpstreamWanzhouCity 12 0 0 0 0 0 0 0 0 0 0 2.3 4.3
OutletWWTP 6 0 0 0 0 2.3 4.3
;
run;
ods trace on;
proc lifereg data=wr18606;
     class sitetype;
         model (low high)=sitetype/distribution=lnormal;
         output out=outlnormal p=p std=stderr q=(.025,.5, .975);
         ods output parameterestimates=pe;
run;

proc glm data=wr18606;
     class sitetype;
         model value=sitetype;
         output out=outglm p=p uclm=u lclm=l;
run;
quit;
proc transpose data=pe out=tpe prefix=e;
     var estimate;
run;
data tpe2;
     set tpe;
         mean1=exp(e1+.5*e4*e4);
         mean2=exp(e1+e2+.5*e4*e4);
run;
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-04-29 15:25:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

他乡客

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
这个不懂，帮你顶顶吧。。。
赞一下 回复此楼
互相帮助，共同进步
3楼2012-04-29 20:59:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

myselfsky8159

木虫 (著名写手)
自己顶个，求高手出招。呵呵。
赞一下 回复此楼
2楼2012-04-29 20:47:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paperhunter

荣誉版主 (文学泰斗)

还没想好

优秀版主优秀版主优秀版主优秀版主优秀版主优秀版主优秀版主

【答案】应助回帖

感谢参与，应助指数 +1
楼主的问题专业性比较强，如果发在专业版块可能应助更有效率。
赞一下 回复此楼
咱也是有组织的人了...
4楼2012-04-29 23:00:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxlion811

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
myselfsky8159: 金币+30, ★★★很有帮助, 谢谢专家，不过还是有些地方不太明白。 2012-05-03 13:08:43
哈哈，我这个程序版的专家终于派上用场了。
/* 这一段读入原始数据， 并且当value等于0时，赋值low和high  */
data wr18606;
     attrib sitetype length=$30;
     input sitetype n@@;
     do i= 1 to n;
        input value@@;
            if value =0 then do; low=.; high=.3; end;
            else do; low=value; high=value; end;
                output;
         end;
cards;
UpstreamWanzhouCity 12 0 0 0 0 0 0 0 0 0 0 2.3 4.3
OutletWWTP 6 0 0 0 0 2.3 4.3
;
run;
/*   做生存分析， 并且输出统计分析结果到数据集pe中 */
ods trace on;
proc lifereg data=wr18606;
     class sitetype;
         model (low high)=sitetype/distribution=lnormal;
         output out=outlnormal p=p std=stderr q=(.025,.5, .975);
         ods output parameterestimates=pe;
run;
/*做方差分析，结果放在数据集outglm中  */
proc glm data=wr18606;
     class sitetype;
         model value=sitetype;
         output out=outglm p=p uclm=u lclm=l;
run;
quit;

/*把生存分析结果pe矩阵进行转置 */
proc transpose data=pe out=tpe prefix=e;
     var estimate;
run;

/*  计算mean1和mean2  */
data tpe2;
     set tpe;
         mean1=exp(e1+.5*e4*e4);
         mean2=exp(e1+e2+.5*e4*e4);
run;


ps：从上面的代码看出，这个SAS水平还可以。
赞一下 回复此楼
开心努力一辈子
5楼2012-05-03 13:03:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭