有谁知道下面程序中MTM1Delta怎么定义吗？因为缺少MTM1Delta这个函数，程序无法运行[s:15][s:15]
for (j in 2:nobs)
{
y[i,j]~dbern(p[i,j])
logit(p[i,j])<-delta[i,j]+gamma[i,j]*y[i,j-1]
etaM[i,j]<-beta0all+beta1all*t[(sumn-nobs+j)]
# MTM1Delta is a function written in the WinBUGS development
# interface for computing the intercept delta[i,j] to satisfy
# marginal mean and serial dependence simultaneously
delta[i,j]<-MTM1Delta(etaM[i,j],etaM[i,j-1],gamma[i,j])
gamma[i,j]<-cor }}
