请教这段SAS代码哪里出错了？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

运行之后出现问题，确实不是很懂，求教！

PROC IML;

start MTE2co ( a, Lelr, Uelr, rm, ca, power );

error = 0;

if ( a>1 | a<0 ) then do; error =1; print “error”

“Test significance level must be in 0-1”; end;

If ( Lelr>Uelr ) then do; error=1; print “error”

“Lower equivalence limit ratio must be <=Upper equivalence

limit ratio”; end;

if ( Uelr < Lelr ) then do; error=1; print “error”

“Upper equivalence limit ratio must be >=Lowerequivalence

limit ratio”; end;

if ( rm <= Lelr | rm >=Uelr) then do; error=1;

print “error”“Expected ratio must be inLeld-Ueld”;

end;

if ( rm<=0 ) then do; error=1; print “error” “Expetted ratio must be >0”; end;

if ( power>100 | power<1 ) then do; error=1;

print “error” “Power(％)must be in 1-100”; end;

if ( error = 1 )then stop;

if (error = 0) then do;

sd = ca/sqrt(2); SDD=2*sd; n=2;

do until(pwl+pw2-1 >=power/100);

L_ncp = abs(log(rm)-log(Lelr))#sqrt(n)/(sd#sqrt(2));

U_ncp = abs(log(rm)-log(Uelr))#sqrt(n)/(sd#sqn(2));

df=2*(n-1); t=tinv(1-a,df);

pw1 = 1-probt(t, df, L_ncp);

pw2 = 1-probt(t, df, U_ncp);

n = n+0.01; end;

n = ceil(n-0.01);

print a[1abel = “Test significance level”]

Lelr [ label =“Lower equivalence limitratio”]

Uelr [ 1abel =“Upper equivalence limitratio”]

rm [ label =“Expected ratio”]

ca [ 1abel = “Crossover ANOVA,sqrt(MSE)” ]

SDD [ label = “SD differences” ]

Power [ 1abel= “Power(％)” ]

N [ label=“n per sequence group”]; end;

finish MTE2co;

run MTE2co(0.025, 0.8, 1.25, 1, 0.198, 90)

run MTE2co(0.025, 0.8, 1.25, 1, 0.294, 90)

quit;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：sas代码 哪里出错了 significance Equivalence differences

   @lovexialulu ，请求大神帮忙啊！

金树2017 发表于 2017-12-14 09:49
@lovexialulu ，请求大神帮忙啊！

对proc iml; 不熟，对你的程序调了下，但是log还是有error：
NOTE: Module MTE2CO defined.
216        /*run MTE2co;*/
217        run MTE2co(0.025, 0.8, 1.25, 1, 0.294, 90);
ERROR: (execution) Matrix has not been set to a value.

1. PROC IML;
   
2. /*
   
3. a=0.025;
   
4. Lelr=0.8;
   
5. Uelr=1.25;
   
6. rm=1;
   
7. ca=0.294;
   
8. power=90;
   
9. 
   
10. start MTE2co;
    
11. */
    
12. start MTE2co (a, Lelr, Uelr, rm, ca, power );
    
13. error = 0;
    
14. if ( a>1 | a<0 ) then do; error =1;
    
15. print 'error' 'Test significance level must be in 0-1'; end;
    
16. If ( Lelr>Uelr ) then do; error=1;
    
17. print 'error' 'Lower equivalence limit ratio must be <=Upper equivalence limit ratio'; end;
    
18. if ( Uelr < Lelr ) then do; error=1;
    
19. print 'error' 'Upper equivalence limit ratio must be >=Lower equivalencelimit ratio'; end;
    
20. if ( rm <= Lelr | rm >=Uelr) then do; error=1;
    
21. print 'error' 'Expected ratio must be inLeld-Ueld';end;
    
22. if ( rm<=0 ) then do; error=1;
    
23. print 'error' 'Expetted ratio must be >0'; end;
    
24. if ( power>100 | power<1 ) then do; error=1;print 'error' 'Power(%)must be in 1-100'; end;
    
25. if ( error = 1 )then stop;
    
26. if (error = 0) then do;
    
27. 
    
28. sd = ca/sqrt(2);
    
29. SDD=2*sd;
    
30. n=2;
    
31.         do until(pwl+pw2-1 >=power/100);
    
32. L_ncp = abs(log(rm)-log(Lelr))#sqrt(n)/(sd#sqrt(2));
    
33. U_ncp = abs(log(rm)-log(Uelr))#sqrt(n)/(sd#sqrt(2));
    
34. df=2*(n-1); t=tinv(1-a,df);
    
35. pw1 = 1-probt(t, df, L_ncp);
    
36. pw2 = 1-probt(t, df, U_ncp);
    
37. n = n+0.01;
    
38.         end;
    
39. n = ceil(n-0.01);
    
40. 
    
41. print a [ label = 'Test significance level'] Lelr [ label ='Lower equivalence limitratio']
    
42. Uelr [ label ='Upper equivalence limitratio'] rm [ label ='Expected ratio']
    
43. ca [ label = 'Crossover ANOVA,sqrt(MSE)' ] SDD [ label = 'SD differences' ]
    
44. Power [ label= 'Power(%)' ] n [ label='n per sequence group'] ;
    
45. end;
    
46. 
    
47. finish;
    
48. /*run MTE2co;*/
    
49. run MTE2co(0.025, 0.8, 1.25, 1, 0.294, 90);
    
50. quit;

复制代码

lovexialulu 发表于 2017-12-14 13:07
对proc iml; 不熟，对你的程序调了下，但是log还是有error：
NOTE: Module MTE2CO defined.
216        ...

谢谢，感谢您的帮助！我是菜鸟啊，刚接触SAS，好多的都不懂，您能推荐几本比较好的学习SAS的书吗？

金树2017 发表于 2017-12-14 13:50
谢谢，感谢您的帮助！我是菜鸟啊，刚接触SAS，好多的都不懂，您能推荐几本比较好的学习SAS的书吗？

入门书 the little sas book