Hello Everyone,

I appreciate for your help on the following stan model. I ran into a computational issue when the data set gets larger to 350,000 observations and 40 variables. I wonder any guidance you can provide me on optimizing it.

Thank you very much!

data {

int<lower=1> n;

int<lower=1> p;

int<lower=0> V[n, p];

int<lower=0> I[n]; /

}

parameters{

vector[p] beta;

vector<lower=0>[p] alpha;

vector[n] theta;

}

model{

theta ~ normal(0, 1);

alpha ~ lognormal(1, 1);

beta ~ normal(0, 3);

for(k in 1:p){

V[,k] ~ binomial_logit(I, beta[k] + alpha[k] * theta); //model

}

}