實在是想不出來，為啥每個數(shù)據(jù)集數(shù)值都一樣��？求大神指點迷津

是 否 +2 論壇幣

k人 參與回答

經(jīng)管之家送您一份

應(yīng)屆畢業(yè)生專屬福利!

求職就業(yè)群

趙安豆老師微信：zhaoandou666

經(jīng)管之家聯(lián)合CDA

送您一個全額獎學(xué)金名額~ !

立即領(lǐng)取

感謝您參與論壇問題回答

經(jīng)管之家送您兩個論壇幣！

+2 論壇幣

/*group=3*/
%macro sim ;
%let j=1 ;
%do j=1 %to 100;**隨機抽取樣本100次***；
%let i=0 ;
%do i= 0 %to 2 ;
DATA SIMU&i.&j ;
d=0.1+&i*0.4;  ***因變量***；
DO n = 1 to 67  ;/*67sample size，樣本含量67例*/
R=1.002+9.023/(1+exp(10.47*(0.5-d)));**效應(yīng)值***；
v&i=R+3*rannor(1); /*標(biāo)準(zhǔn)差不變,蒙特卡羅模擬效應(yīng)值*/
output ;
end;
run ;
%end ;
proc sql ; create table new&j as
select mean(v0) as m0,mean(v1)as m1 ,mean(v2) as m2
from simu&i.&j;
output;
quit;
proc append data=new&j base=a3 ; run;
我想要d=0.1、0.5、0.9時的效應(yīng)值被模擬67例的均數(shù)，放入新的數(shù)據(jù)集&newj,共抽取100次樣本，所以應(yīng)該每種劑量100個均值，可我為什么每次生成的simu&i&j數(shù)據(jù)集都是相同的��？求大神指點？

掃碼加我 拉你入群

請注明：姓名-公司-職位

以便審核進(jìn)群資格，未注明則拒絕

分享0 收藏0 回帖

關(guān)鍵詞：指點迷津 數(shù)據(jù)集 proc sql output Select

求解，跪求，高人啊，快快光顧指導(dǎo)吧

差距太大嗎？為啥沒人回貼

能不能告訴我怎樣才有人回帖？急啊，畢業(yè)論文急需��！

rannor(1) 用的是同一個 random seed 1, 所以每個循環(huán)跑出來的結(jié)果都一樣。只需要把 1 改成 -1， 讓 SAS 根據(jù)系統(tǒng)時鐘來選 random seed 就可以了。

那d=0.5和0.1的時候數(shù)據(jù)都一樣啊，這是怎么回事

高貴林 發(fā)表于 2016-3-26 09:06
rannor(1) 用的是同一個 random seed 1, 所以每個循環(huán)跑出來的結(jié)果都一樣。只需要把 1 改成 -1， 讓 SAS 根 ...

那d=0.5和d=0.1時數(shù)據(jù)都一樣啊，這是怎么回事？先跪謝一下，上一次你對我的回答，這次再拜

沒時間仔細(xì)看，不過是不是每次循環(huán)應(yīng)該用不同的seed？

lipeng0809 發(fā)表于 2016-3-26 23:32
那d=0.5和d=0.1時數(shù)據(jù)都一樣啊，這是怎么回事？先跪謝一下，上一次你對我的回答，這次再拜

你的程序有點問題，幫你改了一下：

%macro sim ;
%do j=1 %to 100;
      DATA SIMU&j ;
          do n = 1 to 67  ;
              %do i= 0 %to 2 ;
                  d=0.1+&i*0.4;
                  R=1.002+9.023/(1+exp(10.47*(0.5-d)));
                  v&i=R+3*rannor(-1);   
                     /* 如果想讓每次運行結(jié)果完全一樣，可以把 rannor(-1) 改成 rannor(&j) */
                  %end;
                  output;
          end;
      run;

      proc sql ; create table new&j as
      select mean(v0) as m0,mean(v1)as m1 ,mean(v2) as m2
      from SIMU&j;
      quit;

      proc append data=new&j base=a3;
      run;
%end;
%mend;

%sim

高貴林 發(fā)表于 2016-3-27 09:33
你的程序有點問題，幫你改了一下：

%macro sim ;

高手，為啥要把n=1 to 67 寫道%do i=0 % to 2的上邊，要是反過來會怎么樣？生成的simu&j數(shù)據(jù)庫,v0應(yīng)該對應(yīng)d=0.1,v1對應(yīng)d=0.5,v2對應(yīng)d=0.9,時的隨機R值。為什么只有d=0.9?是只保存了最后的d的緣故嗎？再次感謝您對我的幫助。