Hi,

I am trying to do modeling based on SDT framework. Therefore, I am following guidelines in Matti Vuorre’s posting (Matti’s homepage - Bayesian Estimation of Signal Detection Models). This has been very useful, but still, I am a bit confused in applying to my data analysis purpose. What I am trying to do is to find out whether criterion would differ depending on “sceneCond” (which has three conditions) factor in the code ~~below~~.

```
## model formula
glmm2 <- bf(samediffresp ~ Phi(dprime * changeCond - c),
dprime ~ sceneCond + (sceneCond | subindex),
c ~ sceneCond + (sceneCond | subindex),
nl = TRUE
)
prior <- get_prior(data = totaldata, glmm2, family = bernoulli(link="identity"))
prior <- c(
prior(normal(.5, 3), nlpar = "dprime"),
prior(normal(0, 1.5), nlpar = "c")
)
evsdt_glmm2 <- brm(glmm2,
family = bernoulli(link = "identity"),
data = totaldata,
prior = prior,
control = list(adapt_delta = .99),
cores = 4, init = 0,
# file = "sdtmodel2-2"
)
```

I want to ask you guys 1) whether I wrote the formula correctly, and if it is right, 2) how I can compute bayes factor of this model that incorporates the “sceneCond”, and 3) how to do post-hoc analysis between two conditions?

Regarding the second question, I want to check whether it is ~~to ~~appropriate to construct a null model like the below.

```
## model formula
glmm2 <- bf(samediffresp ~ Phi(dprime * changeCond - c),
dprime ~ (sceneCond | subindex),
c ~ (sceneCond | subindex),
nl = TRUE
)
prior <- get_prior(data = totaldata, glmm2, family = bernoulli(link="identity"))
prior <- c(
prior(normal(.5, 3), nlpar = "dprime"),
prior(normal(0, 1.5), nlpar = "c")
)
null_model <- brm(glmm2,
family = bernoulli(link = "identity"),
data = totaldata,
prior = prior,
control = list(adapt_delta = .99),
cores = 4, init = 0,
# file = "sdtmodel2-2"
)
```

This question is quite lengthy, and I would greatly appreciate receiving an answer to any of the questions I’ve posed!

Thank you!