关于本站
人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!
获取电子版《CDA一级教材》
完整电子版已上线CDA网校,累计已有10万+在读~ 教材严格按考试大纲编写,适合CDA考生备考,也适合业务及数据分析岗位的从业者提升自我。
TOP热门关键词
免费学术公开课,扫码加入![]() |
我需要修改下面的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 | 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
您可能感兴趣的文章
- SAS软件培训 ... | SAS help 中 catmod 例28.8程序中 ...
- SAS软件培训 ... | 求助,在sas9.13的mdc里运行mixe ...
- SAS软件培训 ... | SAS 9.2 9.3 32/64位迅雷快传下载 ...
- SAS软件培训 ... | 送一个9.1.3的sid和永久破解sash ...
- SAS软件培训 ... | Moody's- Impact of Japan's Dis ...
- SAS软件培训 ... | 北美SAS base 11月5号
- SAS软件培训 ... | 好书! Using SAS in Financial R ...
- SAS软件培训 ... | What’s New in SAS 9.2.pdf
人气文章
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。



