Thank you for your reply. Is there a way to write the following code without using proc iml? start pdfTN(x, mu, sigma, a, b); /* Support one-sided truncation. Define Phi(.M)=0 and Phi(.P)=1 */ cdf_a = choose(a=., 0, cdf("Normal", a, mu, sigma)); /* Phi(a) */ cdf_b = choose(b=., 1, cdf("Normal", b, mu, sigma)); /* Phi(b) */ return( pdf("Normal", x, mu, sigma) / (cdf_b - cdf_a) ); finish; Or how can I access the function pdfTN in my code? data pdf; set example; do x = 0 to 13 by 0.01; f1=Prior_1*pdf("Normal",x,Pred_1,Var_1); f2=Prior_2*pdfTN(x,Pred_2,Var_2,3,8); output pdf; keep f1 f2; end; run; When I run the code, the error message is "The function PDFTN is unknown, or cannot be accessed."
... View more
