| 查看: 791 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
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; |
回复此楼» 猜你喜欢
265求调剂
已经有9人回复
-
085600材料与化工调剂
已经有20人回复
-
专硕 351 086100 也是考的材科基 本科也是材料
已经有6人回复
-
085600专硕材料与化工348分求调剂
已经有10人回复
-
085600 295分求调剂
已经有21人回复
-
285求调剂
已经有5人回复
-
一志愿安徽大学0817化学工程与技术，求调剂
已经有9人回复
-
一志愿0817化学工程与技术，求调剂
已经有8人回复
-
271分求调剂学校
已经有3人回复
-
生物学308分求调剂（一志愿华东师大）
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 文章小修 其中一审稿人意见
已经有11人回复
- [求助] sci论文revise, 其中一个拒稿的审稿人意见该如何回复？
已经有14人回复
- 审稿意见如何回复，求助
已经有4人回复
- 求助对JCED审稿意见的理解
已经有7人回复
- 能否在审稿意见中向作者索要程序代码？
已经有39人回复
- SCIENCE CHINA Information Sciences 一个审稿意见为中文一个为英文，如何回复
已经有21人回复
- 第一次发SCI，麻烦高人帮我看看这条审稿意见如何回复
已经有19人回复
- 请大家帮忙看看如何回复审稿意见
已经有4人回复
- 审稿意见回来，一个审稿人让据，编辑说大修，大家帮我看看希望大不
已经有24人回复
paperhunter
荣誉版主 (文学泰斗)
还没想好
-
专家经验: +14 - SEPI: 1
- 应助: 14603 (教授)
- 贵宾: 11.659
- 金币: 345467
- 散金: 9262
- 红花: 1423
- 沙发: 2016
- 帖子: 127554
- 在线: 8880.4小时
- 虫号: 971012
- 注册: 2010-03-14
- 专业: 环境工程
- 管辖: 论文投稿
4楼2012-04-29 23:00:52
myselfsky8159
木虫 (著名写手)
- 应助: 6 (幼儿园)
- 贵宾: 0.01
- 金币: 6879.5
- 散金: 931
- 红花: 6
- 沙发: 2
- 帖子: 2193
- 在线: 349.6小时
- 虫号: 684894
- 注册: 2008-12-30
- 性别: GG
- 专业: 食品科学基础
2楼2012-04-29 20:47:08
他乡客
至尊木虫 (知名作家)
- 应助: 1887 (讲师)
- 金币: 14731.7
- 散金: 70
- 红花: 41
- 帖子: 6254
- 在线: 384.4小时
- 虫号: 1575919
- 注册: 2012-01-13
- 性别: GG
- 专业: 有机合成
3楼2012-04-29 20:59:21
sxlion811
金虫 (正式写手)
- 应助: 12 (小学生)
- 贵宾: 0.128
- 金币: 1079.5
- 散金: 2188
- 红花: 3
- 帖子: 870
- 在线: 445.8小时
- 虫号: 479151
- 注册: 2007-12-14
- 专业: 食品科学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
myselfsky8159: 金币+30, ★★★很有帮助, 谢谢专家，不过还是有些地方不太明白。 2012-05-03 13:08:43
感谢参与，应助指数 +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
