这段程序，请教大牛们看看啊！

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

这段代码，运行之后出现一点问题，不太懂，求解！

datass_for_be_hvd;

keepalpha beta CV l u gmr power n_ema_3 n_ema_4;

inputalpha beta cv l U gmr;

z_alpha = probit (1-alpha);

z_beta = probit (1-beta);

SW = sqrt(log(cv**2+1));

*CV is defined on the original(not logarithmic) scale．;

*SW is defined as standard deviation on the logaritlunicscale．：

*For log-normal data, the following relationship exists：;

*CV=SQRT(Exp(SW*SW)-1), Sw=SQRT(Log(CV*CV+1));

*---------------------------------------------------------------;

*EMA highly variable drugs;

*---------------------------------------------------------------;

sw0_ema = 0.294;

lvema = log(l)*sw / sw0_ema;

uvema = log(U)*sw / sw0_ema;

l_ema = exp(lvema);

u_ema = exp(uvema);

ifCV >=0.5 then do;

l_ema = 0.6984：

u_ema = 1.4319; end;

n0 = ceil(2*((z_alpha + z_beta)*SW)**2 / (log(gmr)-log(u_ema))**2);

do until ( power>=(1-beta));

nc1 = sqrt(n0)*(log(gmr)-log(l_ema)) / sqrt(2*SW**2);

ifnc1>20 thennc1=20; else if nc1< -20 thennc1 =-20：

nc2 = sqrt(n0)*(log(gmr)-log(u_ema)) / sqrt(2*SW**2);

ifnc2>20 thennc2=20; else if nc2<-20 thennc2=-20;

beta1 = probt (-tinv(1-alpha，n0-2), n0-2,nc2);

beta2 = probt ( tinv(1-alpha，n0-2), n0-2,nc1);

power = beta1 - beta2;

n0 = n0 + 1;

end;

n_EMA = n0 - 1：

n_EMA_3 = ceil(n_EMA*3/4);

n_EMA_4 = ceil(n_EMA/2);

*--------------------------------------------------------------------------;

output;

cards;

0.05 0.20 0.60 0.8 1.25 0.85

0.05 0.20 0.60 0.8 1.25 0.9

0.05 0.20 0.60 0.8 1.25 0.95

0.05 0.20 0.60 0.8 1.25 1

0.05 0.20 0.60 0.8 1.25 1.05

0.05 0.20 0.60 0.8 1.25 1.1

0.05 0.20 0.60 0.8 1.25 1.15

;

Run;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：relationship logarithmic Deviation following relations

3点吧：1）， if  then  后面至少要空一格 比如‘ifCV >=0.5 then do;’ 要换成 ‘if CV >=0.5 then do;’ 2）sas语句 结束是英文‘;’ 不是中文‘；’或者‘：’ 3）代码要在英文状态下编码,别切到搜狗中文输入法了；

1. data a;
   
2. input alpha beta cv l U gmr;
   
3. cards;
   
4. 0.05 0.20 0.60 0.8 1.25 0.85
   
5. 0.05 0.20 0.60 0.8 1.25 0.9
   
6. 0.05 0.20 0.60 0.8 1.25 0.95
   
7. 0.05 0.20 0.60 0.8 1.25 1
   
8. 0.05 0.20 0.60 0.8 1.25 1.05
   
9. 0.05 0.20 0.60 0.8 1.25 1.1
   
10. 0.05 0.20 0.60 0.8 1.25 1.15
    
11. ;
    
12. 
    
13. Run;
    
14. 
    
15. data ss_for_be_hvd;
    
16. set a;
    
17. z_alpha = probit (1-alpha);
    
18. z_beta = probit (1-beta);
    
19. SW = sqrt(log(cv**2+1));
    
20. sw0_ema = 0.294;
    
21. lvema = log(l)*sw / sw0_ema;
    
22. uvema = log(U)*sw / sw0_ema;
    
23. l_ema = exp(lvema);
    
24. u_ema = exp(uvema);
    
25. if CV >=0.5 then do;
    
26. l_ema = 0.6984;
    
27. u_ema = 1.4319; end;
    
28. n0 = ceil(2*((z_alpha + z_beta)*SW)**2 / (log(gmr)-log(u_ema))**2);
    
29. 
    
30. 
    
31. do until ( power>=(1-beta));
    
32. nc1 = sqrt(n0)*(log(gmr)-log(l_ema)) / sqrt(2*SW**2);
    
33. if nc1>20 then nc1=20; else if nc1< -20 then nc1 =-20;
    
34. nc2 = sqrt(n0)*(log(gmr)-log(u_ema)) / sqrt(2*SW**2);
    
35. if nc2>20 then nc2=20; else if nc2<-20 then nc2=-20;
    
36. beta1 = probt (-tinv(1-alpha,n0-2), n0-2,nc2);
    
37. beta2 = probt ( tinv(1-alpha,n0-2), n0-2,nc1);
    
38. power = beta1 - beta2;
    
39. n0 = n0 + 1;
    
40. end;
    
41. n_EMA = n0 - 1;
    
42. 
    
43. n_EMA_3 = ceil(n_EMA*3/4);
    
44. 
    
45. n_EMA_4 = ceil(n_EMA/2);
    
46. keep alpha beta CV l u gmr power n_ema_3 n_ema_4;
    
47. run;

复制代码

嗯，您说的这些都有注意到了，复制到论坛上的时候会有点偏差，您可以帮忙运行一下吗？按您自己修改的运行一下，我真的是找不到别的错误了！

金树2017 发表于 2017-12-13 15:52
嗯，您说的这些都有注意到了，复制到论坛上的时候会有点偏差，您可以帮忙运行一下吗？按您自己修改的运行一 ...

我修改后run了下 log干净的，就是简单的计算；是不是你的计算rule 写的不对

lovexialulu 发表于 2017-12-13 15:47
3点吧：1）， if  then  后面至少要空一格 比如‘ifCV >=0.5 then do;’ 要换成 ‘if CV >=0.5 then do;’  ...

您好！您的这段代码，我也放到SAS上试过了，出现问题啊！

lovexialulu 发表于 2017-12-13 15:59
我修改后run了下 log干净的，就是简单的计算；是不是你的计算rule 写的不对

计算规则是FDA制定的，应该不会出问题吧！

金树2017 发表于 2017-12-13 16:25
您好！您的这段代码，我也放到SAS上试过了，出现问题啊！

我的sas 出来的 log是clean的，如下：

168        
169        data a;
170        input alpha beta cv l U gmr;
171        cards;

NOTE: The data set WORK.A has 7 observations and 6 variables.
NOTE: DATA statement used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds
      
179        ;

180        
181        Run;
182        
183        data ss_for_be_hvd;
184        set a;
185        z_alpha = probit (1-alpha);
186        z_beta = probit (1-beta);
187        SW = sqrt(log(cv**2+1));
188        sw0_ema = 0.294;
189        lvema = log(l)*sw / sw0_ema;
190        uvema = log(U)*sw / sw0_ema;
191        l_ema = exp(lvema);
192        u_ema = exp(uvema);
193        if CV >=0.5 then do;
194        l_ema = 0.6984;
195        u_ema = 1.4319; end;
196        n0 = ceil(2*((z_alpha + z_beta)*SW)**2 / (log(gmr)-log(u_ema))**2);
197        
198        
199        do until ( power>=(1-beta));
200        nc1 = sqrt(n0)*(log(gmr)-log(l_ema)) / sqrt(2*SW**2);
201        if nc1>20 then nc1=20; else if nc1< -20 then nc1 =-20;
202        nc2 = sqrt(n0)*(log(gmr)-log(u_ema)) / sqrt(2*SW**2);
203        if nc2>20 then nc2=20; else if nc2<-20 then nc2=-20;
204        beta1 = probt (-tinv(1-alpha,n0-2), n0-2,nc2);
205        beta2 = probt ( tinv(1-alpha,n0-2), n0-2,nc1);
206        power = beta1 - beta2;
207        n0 = n0 + 1;
208        end;
209        n_EMA = n0 - 1;
210        
211        n_EMA_3 = ceil(n_EMA*3/4);
212        
213        n_EMA_4 = ceil(n_EMA/2);
214        keep alpha beta CV l u gmr power n_ema_3 n_ema_4;
215        run;

NOTE: There were 7 observations read from the data set WORK.A.
NOTE: The data set WORK.SS_FOR_BE_HVD has 7 observations and 9 variables.
NOTE: DATA statement used (Total process time):
      real time           0.00 seconds
      cpu time            0.01 seconds

lovexialulu 发表于 2017-12-13 16:28
我的sas 出来的 log是clean的，如下：

168

是啊！我是直接复制的您这段代码，运行了一下还是出现问题，是有点奇怪了！

金树2017 发表于 2017-12-13 16:29
是啊！我是直接复制的您这段代码，运行了一下还是出现问题，是有点奇怪了！

我运行了“lovexiaolulu”的那段代码，
在log窗口也没有提示有问题.

log结果如下：
8        data a;
9        input alpha beta cv l U gmr;
10       cards;

NOTE: 数据集 WORK.A 有 7 个观测和 6 个变量。
NOTE: “DATA 语句”所用时间（总处理时间）:
      实际时间          0.93 秒
      CPU 时间          0.04 秒


18       ;
19
20       Run;
21
22       data ss_for_be_hvd;
23       set a;
24       z_alpha = probit (1-alpha);
25       z_beta = probit (1-beta);
26       SW = sqrt(log(cv**2+1));
27       sw0_ema = 0.294;
28       lvema = log(l)*sw / sw0_ema;
29       uvema = log(U)*sw / sw0_ema;
30       l_ema = exp(lvema);
31       u_ema = exp(uvema);
32       if CV >=0.5 then do;
33       l_ema = 0.6984;
34       u_ema = 1.4319; end;
35       n0 = ceil(2*((z_alpha + z_beta)*SW)**2 / (log(gmr)-log(u_ema))**2);
36
37
38       do until ( power>=(1-beta));
39       nc1 = sqrt(n0)*(log(gmr)-log(l_ema)) / sqrt(2*SW**2);
40       if nc1>20 then nc1=20; else if nc1< -20 then nc1 =-20;
41       nc2 = sqrt(n0)*(log(gmr)-log(u_ema)) / sqrt(2*SW**2);
42       if nc2>20 then nc2=20; else if nc2<-20 then nc2=-20;
43       beta1 = probt (-tinv(1-alpha,n0-2), n0-2,nc2);
44       beta2 = probt ( tinv(1-alpha,n0-2), n0-2,nc1);
45       power = beta1 - beta2;
46       n0 = n0 + 1;
47       end;
48       n_EMA = n0 - 1;
49
50       n_EMA_3 = ceil(n_EMA*3/4);
51
52       n_EMA_4 = ceil(n_EMA/2);
53       keep alpha beta CV l u gmr power n_ema_3 n_ema_4;
54       run;

NOTE: 从数据集 WORK.A. 读取了 7 个观测
NOTE: 数据集 WORK.SS_FOR_BE_HVD 有 7 个观测和 9 个变量。
NOTE: “DATA 语句”所用时间（总处理时间）:
      实际时间          0.83 秒
      CPU 时间          0.06 秒

priss111 发表于 2017-12-13 16:35
我运行了“lovexiaolulu”的那段代码，
在log窗口也没有提示有问题.

能将结果发给我看一下吗？谢谢！