求救:SAS宏程序
发布:jiao_taishan | 分类:SAS软件培训
关于本站
人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!
获取电子版《CDA一级教材》
完整电子版已上线CDA网校,累计已有10万+在读~ 教材严格按考试大纲编写,适合CDA考生备考,也适合业务及数据分析岗位的从业者提升自我。
TOP热门关键词
本人第一次用SAS编程,做了3周了,下面的程序是统计分析最后的一段了,但是出问题了,实在想不出原因,麻烦各位帮忙看看。赋值%letsamplesize=3000;%letgroups=298;%letgama00_hat=5091.25;%letgama01_hat=3002.76;% ...
免费学术公开课,扫码加入![]() |
赋值
%let samplesize=3000; %let groups=298;
%let gama00_hat=5091.25; %let gama01_hat=3002.76; %let gama02_hat=1200.02;%let gama03_hat=3456.93;
%let gama40_hat=773.48;%let gama41_hat=-837.86; %let gama42_hat=-748.40; %let gama43_hat=1339.29;
%let gama50_hat=7668.16; %let gama51_hat=6169.1;%let gama52_hat=8481.79; %let gama53_hat=-2298.41;
%let gama60_hat=2361.84; %let gama61_hat=1150.01; %let gama62_hat=1853.28; %let gama63_hat=-2032.24;
%let gama80_hat=406.75;%let gama81_hat=293.52;%let gama82_hat=-193.78; %let gama83_hat=-255.55;
%let beta1_hat=-67.9889; %let beta2_hat=-369.49;%let beta3_hat=-648.17;
程序
proc sql;
create table data100 as select w1,w2,w3,country from sasuser.stand order by country;
create table sasuser.data_level_2 as select distinct country,w1=1 as w1,w2=1 as w2, w3=1 as w3 from sasuser.stand;
quit;
宏程序
%macro ml_bootstrap;
%do sample=1 %to 100;
%put samp=&sample;
proc surveyselect out=data1 data=sasuser.new_level_1_resd method=URS n=&samplesize outhits noprint;
data newdata1;
merge data1 sasuser.stand1;
proc surveyselect out=data2 data=sasuser.new_level_2_resd method=URS n=&groups outhits noprint;
data newdata2;
merge data2 sasuser.data_level_2;
new_alpha=&gama00_hat+&gama01_hat*w1+&gama02_hat*w2+&gama03_hat*w3+new_u0;
new_beta4=&gama40_hat+&gama41_hat*w1+&gama42_hat*w2+&gama43_hat*w3+new_u1;
new_beta5=&gama50_hat+&gama51_hat*w1+&gama52_hat*w2+&gama53_hat*w3+new_u2;
new_beta6=&gama60_hat+&gama61_hat*w1+&gama62_hat*w2+&gama63_hat*w3+new_u3;
new_beta8=&gama80_hat+&gama81_hat*w1+&gama82_hat*w2+&gama83_hat*w3+new_u5;
data bootstrap_data (keep=country new_y1 new_alpha new_beta4 new_beta5 new_beta6 new_beta8 x1 x2 x3 x4 x5 x6 x8 w1 w2 w3 new_resd_1);
merge newdata1 newdata2; by country;
new_y1=new_alpha+new_beta4*x4+new_beta5*x5+new_beta6*x6+new_beta8*x8+&beta1_hat*x1+&beta2_hat*x2+&beta3_hat*x3+new_resd_1;
proc mixed method=reml noinfo noitprint covtest;
class country;
model new_y1=w1|x4 w1|x5 w1|x6 w1|x8 w2|x4 w2|x5 w2|x6 w2|x8 w3|x4 w3|x5 w3|x6 w3|x8 x1 x2 x3
/solution ddfm=bw notest outp=predicted;
random int x4 x5 x6 x8/solution type=un subject=country;
ods output solutionf=reg_coef covparms=resd_cov;
run;
%put &syslast;
%if &syslast=work.reg_coef %then %do; 请问这句程序是什么意思?&syslast是个什么变量?
proc transpose data=reg_coef out=fixed_p; var estimate; id effect; 这句程序没有运行,什么原因呢?
「经管之家」APP:经管人学习、答疑、交友,就上经管之家!
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
您可能感兴趣的文章
人气文章
1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责;
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。



