# Not run
data(sensory)
Consumer <- factor(sensory$consumer)
Sacarose <- factor(sensory$sacarose)
#### Model
# Not run
dex1 <- Bayesthresh(cor ~ (1|Consumer) + Sacarose, burn = 0,
jump = 1, ef.iter = 10, data=sensory, Write=TRUE)
coef(dex1)
coef(dex1, HPDinterval=TRUE)
# Not run
dex2 <- Bayesthresh(cor ~ (1|Consumer) + Sacarose, burn = 0,
jump = 1, ef.iter = 10, data=sensory)
coef(dex2)
Run the code above in your browser using DataLab