| 查看: 789 | 回复: 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; |
回复此楼» 猜你喜欢
一志愿0817化学工程与技术，求调剂
已经有8人回复
-
271分求调剂学校
已经有3人回复
-
285求调剂
已经有4人回复
-
生物学308分求调剂（一志愿华东师大）
已经有7人回复
-
296材料专硕求调剂
已经有18人回复
-
材料与化工306分找调剂
已经有18人回复
-
302求调剂
已经有9人回复
-
一志愿北京科技大学材料工程085601，求调剂
已经有11人回复
-
283分材料与化工求调剂
已经有19人回复
-
326求调剂
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 文章小修 其中一审稿人意见
已经有11人回复
- [求助] sci论文revise, 其中一个拒稿的审稿人意见该如何回复？
已经有14人回复
- 审稿意见如何回复，求助
已经有4人回复
- 求助对JCED审稿意见的理解
已经有7人回复
- 能否在审稿意见中向作者索要程序代码？
已经有39人回复
- SCIENCE CHINA Information Sciences 一个审稿意见为中文一个为英文，如何回复
已经有21人回复
- 第一次发SCI，麻烦高人帮我看看这条审稿意见如何回复
已经有19人回复
- 请大家帮忙看看如何回复审稿意见
已经有4人回复
- 审稿意见回来，一个审稿人让据，编辑说大修，大家帮我看看希望大不
已经有24人回复
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
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
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
chenjung0825
木虫 (小有名气)
- 应助: 32 (小学生)
- 金币: 2213
- 散金: 8
- 红花: 2
- 帖子: 221
- 在线: 360.6小时
- 虫号: 1176821
- 注册: 2010-12-24
- 专业: 生物大分子结构与功能
6楼2012-05-03 13:18:06
sesame_oil
荣誉版主 (文坛精英)
- SEPI: 1
- 应助: 4518 (副教授)
- 贵宾: 2.947
- 金币: 35328.4
- 散金: 385
- 红花: 188
- 沙发: 2
- 帖子: 18585
- 在线: 536.5小时
- 虫号: 1188734
- 注册: 2011-01-12
- 专业: 计算机应用技术
- 管辖: 论文投稿
7楼2012-05-03 13:18:27
