Logical. If TRUE then scale the RW1 and RW2 and BESAG and BYM and BESAG2 and RW2D models so the their (generlized) variance is 1. Default value is inla.getOption("scale.model.default")
Now, say if I use my model,
Y`X+f(ID,X,model="besag")
In that case If I use scale.model=T, will it resolve identifiability issue? If not how can we do so? (as it is a spatially varying coefficient part included, I keep const=F)
P.S.: Identifiability issue comes in my opinion as the model seems to be Y=intercept+X \beta_1 + X \beta_ID,X = intercept+X(\beta1+\beta_ID,X)
Alokesh.
*Message sent from a system outside of UConn.*