Hey everyone - there is a way to use PROC MIXED to fit multiple response variables. A good example can be found in the documentation for the TYPE= option for the REPEATED statement. Look at TYPE=UN@AR(1), where Height and Weight for children measured over time is given. The code suggested is:
proc mixed;
class Var Year Child;
model Y = Var Year Var*Year;
repeated Var Year / type=un@ar(1)
subject=Child;
run;
A similar situation is included in the PROC GLIMMIX documentation, but is much more difficult to implement. Start with the example: Joint Modeling of Binary and Count Data. Note that this approach requires different distribution types for each of the dependent variables, which is not required for the PROC MIXED method.
SteveDenham
... View more