如何为变更点逻辑模型的mcmc编写R代码

Question

我想为变更点逻辑模型编写mcmc代码， 我的模型是这样的：

@Value
public class Request {
    private String id;
    private ImmutableList<String> list;
}

其中是要估算的参数：logit{p(x)}=beta0+beta1*doseA+(beta2-beta1)*(doseA-w)*I(doseA>w) +beta3*doseB+(beta4-beta3)*(doseB-v)*I(doseB>v)

beta0-beta4,w,v;是提供的数据

我有一个用于简化模型的Metroplis Hasting代码：

x,doseA,doseB

但是我不知道如何处理模型中的变更点参数，我想知道是否应该将它们与其他参数一起采样或分离。另外，更改点之前和之后的模型具有不同数量的参数，因此它们使用不同的先验，但是我想知道是否应该为此分离数据。

我尝试将logit{p(x)}=beta0+beta1*doseA+beta2*doseB 用于mcmc，但是速度非常慢。我希望我可以自己编写mcmc代码，或使用“ mcmc”或“ HI”之类的软件包

rstan