logitregressionmodel.txt.svn-base

matlab中实现openbugs或winbugs的功能调用

model {
  beta[1:p] ~ dmnorm(mu[], precMat[,])
  for (i in 1:n) {
    s[i] <- inprod(beta[], X[i,])
    logit(pr[i]) <- s[i]
    y[i] ~ dbern(pr[i])
  }
}