运行IML程序时的问题

iamssj

1851

收藏 2012-12-16

悬赏 10 个论坛币已解决

我想对一个数据集作如下处理：给满足条件的某些观测前后一定范围的观测作出标记。
这里我用了IML处理，主要代码如下：
proc iml;
use b;
read all var _all_ into x;
do i=1 to nrow(x);
  if x[i,3] = 1
then
do while(i-50<=j<=i+50);
   x[j,6] = 1;
end;
end;
create new from x;
append from x;
quit;
但是日志却显示有错误：
ERROR: (execution) Matrix has not been set to a value.

operation : <= at line 336 column 18
operands  : _TEM1003, J

_TEM1003    1 row    1 col    (numeric)

   140

J    0 row    0 col    (type ?, size 0)

statement : DO at line 336 column 5
ERROR: Error evaluating DO WHILE expression; DO loop exited.

请问这是为什么？应该如何修改？？？
谢谢。

最佳答案

bbear 查看完整内容

proc iml; use b; read all var _all_ into x; do i=1 to nrow(x); if x = 1 then do; do j=-50 to 50; if ((j-50>=1) & (j+50 while(i-50

扫码加我拉你入群

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

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

全部回复

bbear

2012-12-16 17:36:20

proc iml;
use b;
read all var _all_ into x;
do i=1 to nrow(x);
if x[i,3] = 1  then  do;
      do j=-50 to 50;
         if ((j-50>=1) & (j+50<=nrow(x) )) then do;
            x[j,6] = 1;
         end;
      end; *do j;
   end; *if ;
end; *i;
create new from x;
append from x;
quit;

--> while(i-50<=j<=i+50); 這不是程序, J 沒有定義, 且二個條件寫在一起是不對的語法
--> 靠do j 定義J, 條件改為if , check 不要out of range [1 to row(x)].

扫码加我拉你入群

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

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

ziyenano

2012-12-16 18:13:39

proc iml;
use b;
read all var _all_ into x;
do i=1 to nrow(x);
  if x[i,3] = 1
then
j=max(1,i-50);
do while(i-50<=j<=min(i+50,nrow(x)));
   x[j,6] = 1;
   j=j+1;
end;
end;
create new from x;
append from x;
quit;

扫码加我拉你入群

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

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

ziyenano

2012-12-16 18:31:08

不过用矩阵的指针操作更方便：
proc iml;
use b;
read all var _all_ into x;
do i=1 to nrow(x);
if x[i,3] = 1
then
x[max(i-50,1):min(i+50,nrow(x)),6] = 1;
end;
create new from x;
append from x;
quit;

扫码加我拉你入群

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

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

最佳答案

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群