I am trying to run logistic regression using informative prior assumptions with PROC GENMOD. I found an example of how to assign prior distributions explicitly for dichotomous variables. However, the variable I want to use has 4 levels(called VAR below). I just want my priors to reflect the VAR distributions. Could someone please help me with naming the variables in the prior dataset? Snapshot of data - VAR VARBIG EVENT 1 to 4 1 to 999 0/1 *Input prior distributions; data Prior; input _type_ $ Intercept ???????; *What do I call them to indicate VAR=1, VAR=2...?; datalines; Var 1 1 1 1 Mean 0 .2 .3 .5 ; proc genmod data=model_out; class VAR (ref='1') VARBIG / param=ref; model EVENT(event='1') = VAR VARBIG VAR*VARBIG / DIST=bin LINK=logit; bayes seed=1 coeffprior=normal(input=Prior); run; Thanks!!!
... View more
