[有偿编程] 群设计的样本量估计公式，哪里错了？ [推广有奖]

3关注
0粉丝

已卖：15份资源

硕士生

70%

还不是VIP/贵宾

威望: 0 级
论坛币: 766 个
通用积分: 5.2403
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 10427 点
帖子: 120
精华: 0
在线时间: 242 小时
注册时间: 2010-5-9
最后登录: 2025-11-9

楼主

hwb5258289 发表于 2023-9-28 09:42:09 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

/*SAS 9.4软件实现*/
*ｔ为每个个体的观察时间;
*INDEX为指示变量;
*index = 1时表示求群数，此时需要赋值ｘ为群样本量;
*index = 2时表示求群样本量，此时需赋值ｘ为群数;
*ｔ为每个观察对象被观察的时间；ｋ为研究所需群数;

运行提示：ERROR: (execution) Matrix has not been set to a value.

proc IML;
start CRCR(a,s,lamda1,lamda2,CV,t,index,x,power);
error=0;
if (ａ<0 | a>1)|(ｓ^=1 & s^=2) |(lamda1<0|lamda1>1)|(lamda2<0|lamda2>1)|(index^=1 & index^=2)|(power>100|power<1) then error=1;
if error = 1 then stop ;
cha=abs(lamda1-lamda2);
y=1;
do until (pw>=power);
y=y+1;
if(index=1) then do;
n=x;
k=y;
end;
if (index=2) then do;
k=x;
n=y;
end;
if =1+cv##2#(lamda1##2+lamda2##2)#N#t/(lamda1+lamda2);
pw= 100#CDF('normal',sqrt(n#t#k#(lamda1-lamda2)##2)/sqrt((lamda1+lamda2)#IF)-probit(1-a/s));
end;
pw= round(pw,0.001);
print
a[label="alpha"]
s[label="1or2sided"]
lamda1[label="lamda1"]
lamda2[label="lamda2"]
cha[label="lamda1-lamda2"]
CV[label="CV"]
t[label="t"]
k[label="k"]
N[label="N"]
PW[label="POWER(%)"];
FINISH CRCR;
RUN CRCR(0.05,2,0.005,0.01,0.25,2,1,1000,80);
QUIT;