r - 从 anova 到锯齿状的 ancova

Question

我有一个基于分类变量 a[x] 预测 a0 的模型。现在，我需要添加另一个连续变量，从而将模型转换为 ANCOVA。很高兴听到任何建议。

    modelstring = "
    model {
    for ( i in 1:Ntotal ) {
      y[i] ~ dnorm( a0 + a[x[i]] , 1/ySigma^2 )
    }
    ySigma ~ dunif( ySD/100 , ySD*10 )
    a0 ~ dnorm( yMean , 1/(ySD*5)^2 ) 
    #
    for ( j in 1:NxLvl ) { a[j] ~ dnorm( 0.0 , 1/aSigma^2 ) } {}
    aSigma ~ dgamma( aGammaShRa[1] , aGammaShRa[2] ) 
    # Convert a0,a[] to sum-to-zero b0,b[] :
    for ( j in 1:NxLvl ) { m[j] <- a0 + a[j] } # cell means 
    b0 <- mean( m[1:NxLvl] )
    for ( j in 1:NxLvl ) { b[j] <- m[j] - b0 }
   }
   "