我的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

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

我运行了“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 秒

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

如图所示   ：

嗯，这次又可以运行不报错了！可能是我复制的有问题，但是这次是output里面不出结果啊！我要崩溃了！

Obs alpha beta cv l U gmr power n_EMA_3 n_EMA_4
1 0.05 0.2 0.6 0.8 1.25 0.85 0.80034 75 50
2 0.05 0.2 0.6 0.8 1.25 0.90 0.80111 46 31
3 0.05 0.2 0.6 0.8 1.25 0.95 0.80290 35 23
4 0.05 0.2 0.6 0.8 1.25 1.00 0.80873 33 22
5 0.05 0.2 0.6 0.8 1.25 1.05 0.80661 35 23
6 0.05 0.2 0.6 0.8 1.25 1.10 0.80202 43 29
7 0.05 0.2 0.6 0.8 1.25 1.15 0.80221 61 41

出现的就是这种结果！

在libraies里面的 work folder下面吧，产生的是临时数据集 WORK.SS_FOR_BE_HVD

错误1： l_ema = 0.6984：  : 改为;
错误2： nc1 =-20：  : 改为;
错误3：probt (-tinv(1-alpha，n0-2), n0-2,nc1);  
           probt (-tinv(1-alpha，n0-2), n0-2,nc2);         alpha， 使用了 中文，改为  英文,
错误4：n_EMA = n0 - 1：  : 改为;

大神，谢谢！牛啊！感激感激！

大神，求您也帮忙看看我的另外一个帖子，也是出现这样的问题啊！求救！

大神，求您也帮忙看看我的另外一个帖子，也是出现这样的问题啊！求救！

人生苦短，快用python