Tried using the following code and obtained the error message below. Would greatly appreciate any insight into how to correct code/what I am doing wrong. variables BMI and age are continuous, gender is binary. Code: proc mcmc data=pdata outpost=pdata seed=1234 nmc=20000; ods select PostSumInt; parms beta0 0 beta1 0 beta2 0 beta3 0 s2 1; prior s2 ~ igamma(0.01, s=0.01); prior beta: ~ general(0); w = beta0 + beta1*bmi + beta2*age + beta3*gender; random delta ~ normal(w, var=100) subject=childid; pi = logistic(delta); model sick ~ binomial(p=pi); run; ERROR: A hyperparameter of the random effect delta changed value in observation 2. The hyperparameter is a function of the data set variable bmi, which must remain constant within subjects.
... View more
