Thanks Rick. Q is generated before proc nlmixed, see below code. My goal is to generate confidence bounds for the KM parameter using the likelihood profile method. How does what you propose doing accomplish that? Thank you, Rebekha /*calculate Q*/ proc print data='\\Esc-server1\home\rshaw\EH027\seidmanetal1986.sas7bdat'; run; PROC SQL; CREATE TABLE WORK.SEIDMANETAL1986_Q AS SELECT t1.TSFE_yrs_Range, t1.TSFE_yrs_Avg, t1.duration, t1.conc, t1.PY, t1.'Obs'n, /* Q */ (case when t1.TSFE_yrs_Avg=. then . when t1.TSFE_yrs_Avg<10 then 0 when t1.TSFE_yrs_Avg>(10+t1.duration) then (t1.TSFE_yrs_Avg-10)**3-(t1.TSFE_yrs_Avg-10-t1.duration)**3 else (t1.TSFE_yrs_Avg-10)**3 end) AS Q FROM '\\Esc-server1\home\rshaw\EH027\seidmanetal1986.sas7bdat' t1; run; proc nlmixed data=WORK.SEIDMANETAL1986_Q; parms KM 1e-8; pred = conc*KM*Q*PY; ll=LogPDF("POISSON",Obs,Pred); model Obs ~ general(ll); predict pred out=model1 alpha=0.1 /*df=5*/; Proc print data=model1; run;
... View more