Probit Regression mit Datenerweiterung in Stan
Ich versuche, ein Probit-Modell mit Datenerweiterung unter Verwendung von Stan zu erstellen. Hier haben wir Ergebnissey
entweder 0/1, die uns das Vorzeichen der latenten Variablen @ sagystar
. Dies ist, was ich bisher habe, aber ich bin nicht sicher, wie ich Informationen in das @ hinzufügen somodel
Abschnitt übery
. Irgendwelche Gedanken?
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);
}