见SAS代码转成R代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

如题，我想通过模拟数据来看看Power与所要检验的均值差的关系图（power vs mean difference）， 已经通过SAS做好模拟数据和计算出了两样本均值差与模拟的Power数据集。 现50币求R高手来将其使用R实现：

%let n1 =8;           /* group sizes*/

%let n2 = 8;

%let NumSamples= 3000; /* number of simulated samples */

dataPowerSim(drop=i);

callstreaminit(321) ;

do Delta = 0 to2 by 0.1;

   do SampleID = 1 to &NumSamples;

      c = 1;                            /* Group 1 */

      do i = 1 to &n1;

         x = rand("Normal", 0,1);      /* x ~ N(0,1) */

         output;

      end;

      c = 2;                            /* Group 2 */

      do i = 1 to &n2;

         x = rand("Normal", Delta,1); /* x ~ N(Delta, 1) */

         output;

      end;

   end;

end;

run;

/* 2. Compute(pooled) t test for each sample */

%ODSOff

proc ttest data=PowerSim;

   by Delta SampleID;

   class c;

   var x;

   ods outputttests=TTests(where=(method="Pooled"));

run;

%ODSOn

/* Constructindicator var for obs that reject H0 at 0.05 significance */

data Results;

   set TTests;

   RejectH0 = (Probt <= 0.05);

run;

/* 3. Computeproportion: (# that reject H0)/NumSamples and CI */

%ODSoff

proc freqdata=Results;

   by Delta;

   tables RejectH0 / nocum binomial(level='1');

   output out=Est binomial;

run;

%ODSOn

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：sas代码 R代码 significance difference Simulated

Done by myself

thanks for your sharing, xie xie