SAS 循环问题。。。-经管之家官网!

人大经济论坛-经管之家 收藏本站
您当前的位置> 软件培训>>

SAS软件培训

>>

SAS 循环问题。。。

SAS 循环问题。。。

发布:jzhang04 | 分类:SAS软件培训

关于本站

人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!

获取电子版《CDA一级教材》

完整电子版已上线CDA网校,累计已有10万+在读~ 教材严格按考试大纲编写,适合CDA考生备考,也适合业务及数据分析岗位的从业者提升自我。

完整电子版已上线CDA网校,累计已有10万+在读~ 教材严格按考试大纲编写,适合CDA考生备考,也适合业务及数据分析岗位的从业者提升自我。

有个关于SAS循环的问题求大家帮忙。我很少用SAS,主要用R。但是下面的分析用R我没有搞定。我需要修改下面的SAS程序。目前的程序,我每次只能运行一个plots,通过设置“ifplots=6;”。如果我运行一个plots,结果如下: ...
免费学术公开课,扫码加入


有个关于SAS循环的问题求大家帮忙。我很少用SAS,主要用R。但是下面的分析用R我没有搞定。
我需要修改下面的SAS程序。目前的程序,我每次只能运行一个plots,通过设置“if plots=6;”。如果我运行一个plots,结果如下:
Parameter Estimate Standard Error DF t Value Pr > |t| Alpha Lower Upper Gradient
a -35.7527 52.2230 6 -0.68 0.5191 0.05 -163.54 92.0324 -0.00002
b 0.01587 0.02638 6 0.60 0.5694 0.05 -0.04868 0.08042 -0.03539
alpha 0.3721 0.2290 6 1.63 0.1553 0.05 -0.1882 0.9325 -8.6E-6
我想修改一下这个程序,可以同时运行很多个plots,如10个plots (1,2,3,4,5,6,7,8,9,10) <我的实际数据有上千个plots>. 最后得到一个总结性的结果表,如下(添加了一列“plots”):

Parameter

Estimate

Standard
Error

DF

t Value

Pr>|t|

Alpha

Lower

Upper

Gradient

plots

a

-35.7527

52.223

6

-0.68

0.5191

0.05

-163.54

92.0324

-0.00002

1

b

0.01587

0.02638

6

0.6

0.5694

0.05

-0.04868

0.08042

-0.03539

1

alpha

0.3721

0.229

6

1.63

0.1553

0.05

-0.1882

0.9325

-8.60E-06

1

a

-35.7527

52.223

6

-0.68

0.5191

0.05

-163.54

92.0324

-0.00002

2

b

0.01587

0.02638

6

0.6

0.5694

0.05

-0.04868

0.08042

-0.03539

2

alpha

0.3721

0.229

6

1.63

0.1553

0.05

-0.1882

0.9325

-8.60E-06

2

a

-35.7527

52.223

6

-0.68

0.5191

0.05

-163.54

92.0324

-0.00002

3

b

0.01587

0.02638

6

0.6

0.5694

0.05

-0.04868

0.08042

-0.03539

3

alpha

0.3721

0.229

6

1.63

0.1553

0.05

-0.1882

0.9325

-8.60E-06

3

a

-35.7527

52.223

6

-0.68

0.5191

0.05

-163.54

92.0324

-0.00002

4

b

0.01587

0.02638

6

0.6

0.5694

0.05

-0.04868

0.08042

-0.03539

4

alpha

0.3721

0.229

6

1.63

0.1553

0.05

-0.1882

0.9325

-8.60E-06

4


下面是我的SAS程序:
********
proc import out=trees datafile='C:\GNMM_AW\p.gnmmData.dbh91n.ASRD.STD.xls'
dbms=xls replace;
run;
data trees;
set trees;
if plots=6;
*proc print data=trees;
run;
******** FIT GENERALIZED NONLINEAR MIXED MODEL************* ;
proc nlmixed data=trees; parms a=0.1 b=-0.002;
eta = a + b*year;
p_ann_mort = exp(eta)/(1+exp(eta));
p_ann_surv = 1-p_ann_mort;
p_surv = p_ann_surv**yr_interval;
p_mort = 1-p_surv;
mean_mort = live_off*p_mort;
ll=lgamma(1/alpha + n_dead) - lgamma(1 + n_dead) - lgamma(1/alpha)
+ n_dead*log(alpha*mean_mort)
- (1/alpha + n_dead)*log(1 + alpha*mean_mort); * log likelihood function for negative binomial regression ;
model n_dead ~ general(ll);
predict live_off*p_mort out=predobs1;
run;
不知道哪位高手是否有时间帮我修改一下这段程序。
多谢了~
「经管之家」APP:经管人学习、答疑、交友,就上经管之家!
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
本文关键词:

本文论坛网址:https://bbs.pinggu.org/thread-1270572-1-1.html

人气文章

1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责;
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
数据分析师 人大经济论坛 大学 专业 手机版
联系客服
值班时间:工作日(9:00--18:00)