jfb-h

data {
  int N; // number of rows
  int E; // number of events
  int L; // number of transaction categories
  int K; // number of covariates
  
  vector<lower = 0, upper = 1>[N] choice;
  matrix[N, K] X;