复制粘贴出了点问题，漏掉了[i], model里应该是这样
logit(p[i])<- alpha+b.sex*sex[i]+b.aeti*aeti[i]+b.afp*AFP_value[i]+b.diff*difference[i]+b.age*age[i]        #Linear regression on logit
whether_hcc[i]~dbern(p[i])

有人给我说，要把data standardised 一下再来fit, 不知道是否有用

试了一下好像是可以的。

摘自用户手册：
a) 'undefined real result' indicates numerical overflow. Possible reasons include:
- initial values generated from a 'vague' prior distribution may be numerically extreme - specify
appropriate initial values;
- numerically impossible values such as log of a non-positive number - check, for example, that no
zero expectations have been given when Poisson modelling;
- numerical difficulties in sampling. Possible solutions include:
- better initial values;
- more informative priors - uniform priors might still be used but with their range restricted to
plausible values;
- better parameterisation to improve orthogonality;
- standardisation of covariates to have mean 0 and standard deviation 1.
- can happen if all initial values are equal.