Dear ballardw Thank you for your response. {proc nlmixed data=new qpoints=50 maxiter=1000; parms a1=-1.35 a2=-0.74 Time_a=-0.02 group_a=0.09 cs_a=0.97 LAS_a=-0.29 b1=-6.67 b2=-2.58 Time_b=-0.03 group_b=0.76 cs_b=0.21 LAS_b=0.03 s1=0.35 s2=0.5 s12=0.1; *y1 = HH; if j in (1,2,3,4,5,6,7)then do; eta1 = group_a*Time_a + cs_a*cs_a + LAS_a*LAS_a; delta11= sqrt(1+s1**2) * PROBIT( exp(a1+eta1)/(1+exp(a1+eta1))); delta12= sqrt(1+s1**2) * PROBIT( exp(a2+eta1)/(1+exp(a2+eta1))); if ordinal=1 then do; lik= probnorm (delta11+u1); end; if ordinal=2 then do; lik= 1- probnorm (delta12+u1); end; if (lik > 1e-10) then loglik = log(lik); else loglik = 1e100; end; *y2 = LL; if j in (1,2,3,4,5,6,7)then do; eta2 = group_b*Time_b + cs_b*cs_b + LAS_b*LAS_b; delta21 = sqrt(1+s2**2) * PROBIT( exp(b1+eta2)/(1+exp(b1+eta2))); delta22 = sqrt(1+s2**2) * PROBIT( exp(b2+eta2)/(1+exp(b2+eta2))); if ordinal=1 then do; lik = PROBNORM(delta21+u2); end; if ordinal=1 then do; lik= 1- probnorm (delta22+u2); end; if (lik > 1e-10) then loglik = log(lik); else loglik = 1e100; end; model ordinal ~ general(loglik); random u1 u2 ~ normal ([0,0] , [s1**2,s12,s2**2]) subject=id; estimate "Rho" s12/(s1*s2); run;} {NOTE: The parameter j is assigned the default starting value of 1.0, because it is not assigned an initial value with the PARMS statement. WARNING: Parameter Time_a is in the input data set. WARNING: Parameter group_a is in the input data set. WARNING: Parameter cs_a is in the input data set. WARNING: Parameter LAS_a is in the input data set. WARNING: Parameter Time_b is in the input data set. WARNING: Parameter group_b is in the input data set. WARNING: Parameter cs_b is in the input data set. WARNING: Parameter LAS_b is in the input data set. ERROR: Dependent variable is not in the input data set. NOTE: The SAS System stopped processing this step because of errors. NOTE: PROCEDURE NLMIXED used (Total process time): real time 1.10 seconds cpu time 0.12 seconds} Alphabetic List of Variables and Attributes# Variable Type Len Format Label9238101145671 HH Num 8 HH. HH LAS_a Num 8 F12. LAS_a LAS_b Num 8 F8.2 LAS_b LL Num 8 LL. LL Time_a Num 8 F8.2 Time_a Time_b Num 8 F8.2 Time_b cs_a Num 8 F12. cs_a cs_b Num 8 F8.2 cs_b group_a Num 8 GROUP_AA. group_a group_b Num 8 GROUP_BA. group_b id Num 8 F8.2 id I look forward to your reply Thanks
... View more