[求助]怎样在SAS中生成二分类两水平logistic随机效应数据？谢谢

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我是这样编的，但是P可能大于1，恳请高手指点，谢谢

data aa;
 do r=1 to 10;
 rj=rannor(12345);
   do f=1 to 20;
   treat=(ranuni(12345)>0.5);
   fij=1+1*treat+rj;
 p1=exp(fij)/(1+(exp(fij)));
   eij=rannor(0);
   ez=sqrt((1-p1)*p1)*eij;
   p=p1+ez;
   y=ranbin(0,1,p);
   output;
   end;
   end;
   run;

[此贴子已经被作者于2008-4-29 10:39:09编辑过]

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：logistic ogistic logisti logist logis 数据 SAS 效应 随机 logistic

  请问产生p1后计算p的目的是什么？ 

  fij=1+1*treat+rj;p1=exp(fij)/(1+(exp(fij)));

  eij=rannor(0);ez=sqrt((1-p1)*p1)*eij;p=p1+ez;

  ez是在反映过度离散的情况吗，而且在模拟产生y的时候，p是没有被用到的

  y=ranbin(0,1,p1);

我是直接按照书上多水平logistic模型的公式编的， ez是水平1的误差

yij=（pai）ij+ez=exp(fij)/(1+(exp(fij)))+ez;

fij=1+1*treat+rj;

ez=sqrt((1-p1)*p1)*eij

其中，ez服从二项分布

 

由于自变量和因变量为非线性关系，所以第一水平误差项不用加，只添加第二水平误差项即可。