Dear statlists,
I have a data set of 154 dental implants and 71 patients. i have recorded bone change in 4 times (year1,2,3,4) .
frist i have reshape my data set in long: reshape long level change, i(CODICEIMP) j(time)
then, i declare a panel data set
:
xtset CODICEIMP time
I have predictors at patient level (sex, age, smoking) and implant level (type of surgery, depth, porsthesis,diameter).
i want to predict my outcome (bone change) during the time on the basis of my predictors at both leve (implant and patients) so i need a random slope model.
I set my model null model :
mixed change i.time ||CODICEPZ: time, cov(unstruc) || CODICEIMP: time, cov(unstruc)mle
but if i add my first predictor at implant level
: mixed change i.time, cov(unstruc) ||CODICEPZ: time, cov(unstruc) || CODICEIMP: time DEPTH, cov(unstruc)mle
then STATA gives to me:
Hessian is not negative semidefinite
conformability error
conformability error
I need a nested model beacuse of my outcome of interest (bone change related to implant and patients variables) . how can i solve my problem?
Thankyou!
Comment