For model selection within an ABC framework when the models are nested, say model 1 is equal to model 2 on some subset of the parameter space, is it better to try and do parameter inference or use a model selection framework (for example abcrf)?
In more detail; my specific case is that I only have two models $M_1$ and $M_2$. $M_1$ has an associated parameter space $\Theta_1\subset \mathbb{R}_{\geq 0}$, while for $M_2$ we have an associated parameter space $\Theta_2\subset \mathbb{R}_{\geq 0}^2$. $M_2$ is equal to $M_1$ along $(\theta_1,0)$. I have some data $D$ and wish to know which model is more appropriate. From here
1) Is it better to estimate $Pr(M_2|D)/Pr(M_1|D)$ or $Pr(\theta_2>0|D,M_2)/Pr(\theta_2=0|D,M_2)$?
2) If I cannot analytically obtain the likelihood and so use abc with some chosen summary statistics in place of $D$, does the answer to 1 change?
Regarding 1, I'll compare the options. Assuming $Pr(M_1|D)=Pr(M_2|D)$ we have \begin{align} \frac{Pr(M_2|D)}{Pr(M_1|D)}=\frac{\int_0^\infty \int_0^\infty Pr(D|M_2,(\theta_1,\theta_2))Pr((\theta_1,\theta_2|M_2))d\theta_1 d\theta_2}{ \int_0^\infty Pr(D|M_1,\theta_1)Pr(\theta_1|M_1)d\theta_1 } \end{align} For the numerator we notice \begin{align} \int_0^{\infty}Pr(D|M_2,(\theta_1,0))Pr((\theta_1,0)|M_2))d\theta_1 \\+ \int_0^\infty \int_{0^+}^\infty Pr(D|M_2,(\theta_1,\theta_2))Pr((\theta_1,\theta_2|M_2))d\theta_1 d\theta_2 \\ = Pr(D|\theta_2=0,M_2)Pr(\theta_2=0|M_2) + Pr(D|\theta_2>0,M_2)Pr(\theta_2>0|M_2). \end{align} Concerning the denominator \begin{align} \int_0^\infty Pr(D|M_1,\theta_1)Pr(\theta_1|M_1)d\theta_1 \\ =Pr(D|\theta_2=0,M_2)Pr(\theta_2=0|M_2) \end{align} Which results in \begin{align} \frac{Pr(M_2|D)}{Pr(M_1|D)}=\frac{ Pr(D|\theta_2=0,M_2)Pr(\theta_2=0|M_2) + Pr(D|\theta_2>0,M_2)Pr(\theta_2>0|M_2)}{ Pr(D|\theta_2=0,M_2)Pr(\theta_2=0|M_2) }. \end{align} However the alternative is \begin{align} \frac{Pr(\theta_2>0|D,M_2)}{Pr(\theta_2=0|D,M_2)}=\frac{Pr(D|\theta_2>0,M_2)Pr(\theta_2>0|M_2)}{Pr(D|\theta_2=0,M_2)Pr(\theta_2=0|M_2)} \end{align} Have I made some error or is there some implicit assumption I've made that leads to this?