【问题标题】:Logistic regression Macro逻辑回归宏
【发布时间】:2016-04-22 05:10:15
【问题描述】:
    %macro intercept(i1= ,i2= );
      %let n = %sysfunc(countw(&i1));

       %do i = 1 %to &n;
          %let val_i1 = %scan(&i1,&i,'');
          %let val_i2 = %scan(&i2,&i,'');
          data scores;
            set repeat_score2;
            /* Segment 1 probablity score  */
            p1 = 0;
            z1 =   &val_i1      +
                    a         *      0.03            +
                    r         *    -0.0047841        +
                    p         *   -0.000916081      ;

            p1 = 1/(1+2.71828**-z1);

            /* Segment 2 probablity score  */
           p2 = 0;
           z2 =   &val_i2           +
                   r           *       0.09           +
                   m           *     0.012786245      +
                   c           *    -0.00179618        +
           p2 = 1/(1+2.71828**-z2);

           logit_score = 0;
           if max(p1,p2) =      p1 then logit_score = 1;
           else if max(p1,p2) = p2 then logit_score = 2;
         run;

         proc freq data = scores;
            table logit_score * clu_ /nocol norow nopercent;
         run;
    %end;
  %mend;
%intercept (i1=-0.456491042, i2=-3.207379842, i3=-1.380627318 , i4=0.035684096, i5=-0.855283373);
%intercept (i1=-0.456491042 0, i2=-3.207379842 -3.207379842, i3=-1.380627318 -1.380627318, i4=0.035684096 0.035684096,
            i5=-0.855283373 -0.855283373);

我有上面的宏,它获取上述两个模型的截距,然后计算概率分数,然后根据该概率分数为段分配一个值。

上述宏的第一个问题是,当我使用一个参数执行宏时,它会将宏变量“n”解析为 2 并执行两次。第一次迭代,它给出了正确的结果,而第二次它是错误的。

对于第二个实现(每个宏有两个参数),n 解析为 3,scan 一次解析为这两个值(例如,迭代本身的 i1 为 -0.45 和 0),如果我删除空间,然后它需要'。作为分隔符并将其解析为( 0,45,0 - 每次迭代一个)。我没有得到这个案子的任何结果。

我怎样才能让它以正确的方式工作?

谢谢!!!

【问题讨论】:

  • 我还要提供一个建议——将参数放入数据集中并将它们作为数据集传递。让您的宏循环遍历数据集的行,而不是多个宏变量。

标签: macros sas logistic-regression


【解决方案1】:

%SCANCOUNTW 函数默认将标点符号和空格视为分隔符。由于您的参数包含小数点,因此您需要明确声明 COUNTW%SCAN 的分隔符应为空白。你为%SCAN 做的,但不是为COUNTW 做的。 所以第二行代码应该是:

 %let n = %sysfunc(countw(&i1,' '))

我不确定这是拼写错误还是只是格式化的东西,但在您的 %SCAN 函数中,第三个参数看起来像两个引号一起 '',而不是应该是的 quote-blank-quote ' '

【讨论】:

    猜你喜欢
    • 2020-01-31
    • 2021-11-11
    • 2018-01-26
    • 2020-03-22
    • 2016-12-01
    • 2016-03-24
    • 1970-01-01
    • 1970-01-01
    • 2021-03-11
    相关资源
    最近更新 更多