All,
I understand that proc nlmixed do not have the class statement option. I need help to code a categorical variable Var5 with six levels: Ind_A, Ind_B, Ind_C, Ind_D, Ind_E, Ind_F.
The following is the code which uses 4 explanatory variables:
proc nlmixed data=data;
parms b0=x b1=x b2=x b3=x b4=x
gm1=x gm=x;
w=b0+b1*Var1+b2*Var2+b3*Var3+b4*Var;
ww=w+gm1*t;
consigma=gm2;
model Y~normal(ww,consigma);
random t~normal(0,1) subject=group;
predict ww out= output (keep=Y pred );
run;
The model converges but I need to incorporate the categorical variable.
Any suggestions?