From your error message, in the statement; 857     timeS[locC ] = time+ (time[i+1]-time)*((Col[i,j]-A)/(Col[i,j]+col[i+1,j])); the subscript in time[i+1] exceeds the dimension of time vector. This would happend if col[i,j] > A for some i >= the number of rows in time. You can add some print statements in the loop to see exactly what is happening, and then you will need to figure out how to handle that case, if it is legitimate. ... View more

I fixed it; proc means data=para noprint; class Variable; var estimate; output out=two mean(estimate)=col1 std(estimate)=col2 mean(StdErr)=col3; run; The SAS System           0 50 6.7 4.45480 1.11542 2 Intercept 1 25 2.44651 0.09489 0.30363 3 TIMES 1 25 10.99305 1.57000 1.92721           1 0 50 6.71978 4.45480 1.11542 2 Intercept 1 25 2.44651 0.09489 0.30363 3 TIMES 1 25 10.99305 1.57000 1.92721 ... View more

s =InverseC ; print s; ... View more

You can use the REMOVE function or use LOC and indexing. See Remove or keep: Which is faster? - The DO Loop for examples and comparison. ... View more