Estoy intentando hacer un modelo probit con aumento de datos usando stan. Aquí es donde tenemos resultadosy 0/1 que nos dice el signo de la variable latenteystar. Esto es lo que tengo hasta ahora, pero no estoy seguro de cómo agregar información en elmodel sección sobrey. ¿Alguna idea?

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);
}