BookmarkSubscribeRSS Feed
desireatem
Pyrite | Level 9

The following runs without error but I do not know if the results id correct because I wish to perform the integral using "ni" as a by statement or class statement? I think my code might be treating the data set as a date with 12 observation. However it is a data set with 4 observation repeated  3 times.

Data sim;

Input ni A B;

cards;

1  1 3

1  2 7

1  3 10

1  4 12

2  1 3

2  2 9

2  3 10

2  3 10

3  1 3

3  2 9

3  3 10

3  3 10

;

Run;

proc iml;

use sim;

read all  var{ni A B} into DM;

close;

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

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{ni A  B Answer };

append;

quit;

1 REPLY 1
desireatem
Pyrite | Level 9

I think there is no difference since the integration is performed independently for each row

Catch up on SAS Innovate 2026

Nearly 200 sessions are now available on demand with the SAS Innovate Digital Pass.

Explore Now →
From The DO Loop
Want more? Visit our blog for more articles like these.
Discussion stats
  • 1 reply
  • 1429 views
  • 0 likes
  • 1 in conversation