最初的实现将Bernoulli数的计算放在主要方法中。我做了一个新的类来返回单个Bernoulli数的计算。我做错了什么?= new BigRational[N+1]; bernoulli[1] = new BigRational[k] = bernoulli[k].minus(coef.times(bernoulli[i]));
我试图使用cmdstan2.17.1运行bernoulli示例,即CmdStanDir="/CMDSTAN_HOME"
对于我创建.bashPlease check your Stan program in variable 'bernoulli' and the contents of /Users/panky/Desktop/Uni/IO/T2/tmp
Bernoulli <- rbinom(1000, 1, 0.5) for (i in 1:1000) {#for each Bernoulli realization
if (Bernoulli[i] == 1){ #sample corresponding normal componentBernorm <- ifelse(Bernoulli == 1, rnorm(1, mean=10, sd=1),