BookmarkSubscribeRSS Feed
🔒 This topic is solved and locked. Need further help from the community? Please sign in and ask a new question.
desireatem
Pyrite | Level 9

Please, I wish to ask if any could help correct this error:

NOTE: Module FUNC defined.

47   interval = {A, B};

48   call quad(answer, "Func", interval);

param must be numeric

ERROR: (execution) Character argument should be numeric

This is data set and code:

Data sim;

Input A B;

cards;

1   3

2   7

3   10

4   12

;

Run;

proc iml;

use sim;

read all  var{A B} into DM;

close;

A = DM[,1]; B = DM[,2];

n = nrow(DM);

start Func(x);

   return(  exp(2-x)#cdf("Normal", x) );

finish;

interval = {A, B};

call quad(answer, "Func", interval);

create kaplan1n var{A  B Answer };

append;

quit;

1 ACCEPTED SOLUTION

Accepted Solutions
Rick_SAS
SAS Super FREQ

You want to allocate a vector outside of the loop and fill it within the loop, as described in this article: Pre-allocate arrays to improve efficiency - The DO Loop

This and other efficiency tips that you've asked about are found in Chapter 2 of Statistical Programming with SAS/IML Software

answer = j(nrow(A),1);

do i = 1 to nrow(A);

   call quad(result, "Func", A||B);

   answer = result;

end;

View solution in original post

6 REPLIES 6
Rick_SAS
SAS Super FREQ

If you are trying to compute four integrals, use the i_th row of DM as the interval for the i_th computation:

do i = 1 to nrow(DM);

   call quad(answer, func, DM[i,]);

   ...

end;

desireatem
Pyrite | Level 9

Hello Rick, My interval is A, B. The lower limit is A and upper limit is B, , the code does not specified the upper and lower limit. I have to pick the correct limit .

do i = 1 to nrow(DM);

   call quad(answer, func, DM[i,]);

   ...

end;

I modified my code as such:

Data sim;

Input A B;

cards;

1 3

2 7

3 10

4 12

;

Run;

proc iml;

use sim;

read all  var{A B} into DM;

close;

A = DM[,1]; B = DM[,2];

n = nrow(DM);

start Func(x);

   return(  exp(2-x)#cdf("Normal", x) );

finish;

do i = 1 to nrow(DM);

call quad(answer, "Func", DM[i,]);

end;

create kaplan1n var{A  B Answer };

append;

quit;

statement : CALL at line 8281 column 1

ERROR: (execution) Matrices do not conform to the operation.

operation : QUAD at line 8307 column 1

operands  : *LIT1007, _TEM1001

*LIT1007      1 row       1 col     (character, size 4)

Func

_TEM1001      1 row       1 col     (numeric)

         3

statement : CALL at line 8307 column 1

Rick_SAS
SAS Super FREQ

Your error is that DM[1,2] is a scalar, whereas QUAD is expecting a vector with two elements, as in my code.  You could also use A || B, if you prefer.

Of course, you also will want to accumulate the results inside the loop. Currently you are just overwriting ANSWER four times, and the data set is only writing the answer to the last integration.

desireatem
Pyrite | Level 9

Thanks a lot. Yes you are right,

This just give me one answer and I need a result for each row.

This code below gives me and answer. How ever, just as you said , it is not giving me the answer for each row.  How do I loop it from 1 to 4 so that each row has an answer. ?

Data sim;

Input A B;

cards;

1 3

2 7

3 10

4 12

;

Run;

proc iml;

use sim;

read all  var{A B} into DM;

close;

A = DM[,1]; B = DM[,2];

n = nrow(DM);

start Func(x);

   return(  exp(2-x)#cdf("Normal", x) );

finish;

do i = 1 to nrow(DM);

call quad(answer, "Func", A || B,);

end;

create kaplan1n var{A  B Answer };

append;

quit;

Rick_SAS
SAS Super FREQ

You want to allocate a vector outside of the loop and fill it within the loop, as described in this article: Pre-allocate arrays to improve efficiency - The DO Loop

This and other efficiency tips that you've asked about are found in Chapter 2 of Statistical Programming with SAS/IML Software

answer = j(nrow(A),1);

do i = 1 to nrow(A);

   call quad(result, "Func", A||B);

   answer = result;

end;

desireatem
Pyrite | Level 9

Data sim;

Input A B;

cards;

1 3

2 7

3 10

4 12

;

Run;

proc iml;

use sim;

read all  var{A B} into DM;

close;

A = DM[,1]; B = DM[,2];

n = nrow(DM);

start Func(x);

   return(  exp(2-x)#cdf("Normal", x) );

finish;

answer = j(nrow(DM),1);

do i = 1 to nrow(DM);

call quad(result, "Func", A || B,);

answer = result;

end;

create kaplan1n var{A  B Answer };

append;

quit;

sas-innovate-2024.png

Today is the last day to save with the early bird rate! Register today for just $695 - $100 off the standard rate.

 

Plus, pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

Multiple Linear Regression in SAS

Learn how to run multiple linear regression models with and without interactions, presented by SAS user Alex Chaplin.

Find more tutorials on the SAS Users YouTube channel.

From The DO Loop
Want more? Visit our blog for more articles like these.
Discussion stats
  • 6 replies
  • 2219 views
  • 0 likes
  • 2 in conversation