Пробит регрессии с увеличением данных в стан

Я пытаюсь сделать пробную модель с дополнением данных с использованием stan. Вот где у нас есть результатыy либо 0/1, которые говорят нам знак скрытой переменнойystar, Это то, что я имею до сих пор, но я не уверен, как добавить информацию вmodel раздел оy, Какие-нибудь мысли?

data {
  int<lower=0> N; // number of obs
  int<lower=0> K; // number of predictors
  int<lower=0,upper=1> y[N]; // outcomes
  matrix[N, K] x; // predictor variables
}
parameters {
  vector[K] beta; // beta coefficients
  vector[N] ystar; // latent variable
}
model {
  vector[N] mu; 
  beta ~ normal(0, 100);
  mu <- x*beta;
  ystar ~ normal(mu, 1);
}

Ответы на вопрос(1)

Ваш ответ на вопрос