SAS Proc IML: цикл для заполнения матрицы

Question

У меня есть следующий код, который прекрасно работает в MATLAB и я, чтобы транспонировать в SAS/PROC IML:

[row col] = size(coeff);

A_temp    = zeros(row,col);
for i = 1: row/6            
    A_temp(6*(i-1)+1:6*i,:) = coeff(6*(i-1)+1:6*i,4:col);end;

В Proc IML я делаю следующее:

proc iml;
  use i.coeff;
  read all var {...} into coeff;
  print coeff;

row=NROW(coeff);
print row;
col=NCOL(coeff);
print col;
A_temp=J(row,col,0); *create zero matrix;
print A_temp;

Do i=1 TO row/6;
A_temp[(6*(i-1)+1):(6*i),]=coeff[(6*(i-1)+1):(6*i),(4:col)];
END;
quit;

Код ломается в цикле DO "(исполнение) Матрицы не соответствуют операции."... почему? Если я правильно понимаю в PROC IML, если я хочу выбрать все столбцы (в MATLAB это будет ":"), но в SAS IML я просто оставляю это поле пустым

Answer 1

Вы должны указать это правильно. [Строки,] означает ВСЕ столбцы А, а не только их количество. Посмотрите этот упрощенный пример:

proc iml;
/*  use i.coeff;
  read all var {...} into coeff;
  print coeff;
*/
coeff = J(15,10,3);
row=NROW(coeff);
print row;
col=NCOL(coeff);
print col;
A_temp=J(row,col,0); *create zero matrix;
print A_temp;

Do i=1 TO row;
* does not work; *A_temp[i,]=coeff[i,(4:col)];
A_temp[i,1:col-3]=coeff[i,(4:col)];
END;
quit;