This is cmdstanr version 0.0.0.9000 - Online documentation and vignettes at mc-stan.org/cmdstanr - CmdStan path set to: /home/saumyap/.cmdstanr/cmdstan - Use set_cmdstan_path() to change the path Compiling Stan program... functions { real partial_sum_loc(int[] slice_location,int start, int end, vector attitude, row_vector eta_attitude_rating, matrix s, matrix eta_s_rating){ return categorical_lpmf(slice_location | softmax(to_vector(attitude* eta_attitude_rating + s *eta_s_rating))); } real partial_sum_rat(int[] slice_rating,int start, int end, vector attitude, row_vector eta_attitude_rating, matrix s, matrix eta_s_rating){ return categorical_lpmf(slice_rating | softmax(to_vector(attitude* eta_attitude_rating + s *eta_s_rating))); } real partial_sum_dom(int[] slice_domain,int start, int end, vector attitude, row_vector eta_attitude_rating, matrix s, matrix eta_s_rating){ return categorical_lpmf(slice_domain | softmax(to_vector(attitude* eta_attitude_rating + s *eta_s_rating))); } } data{ int N; // number of instances in the data int K; // number of covariates int R; // No. of possible ratings int L; // No. of Locations int D; // No. of Unique domains matrix[N, K] s ;// matrix of sensitive attributes int location[N]; // location int domain[N]; //domain int rating[N]; // rating } transformed data{ vector[K] zero_K; vector[K] one_K; vector[R] zero_R; vector[R] one_R; vector[L] zero_L; vector[L] one_L; vector[D] zero_D; vector[D] one_D; zero_K = rep_vector(0,K); one_K = rep_vector(1,K); zero_R = rep_vector(0,R); one_R = rep_vector(1,R); zero_L = rep_vector(0,L); one_L = rep_vector(1,L); zero_D = rep_vector(0,D); one_D = rep_vector(1,D); } parameters{ vector[N] attitude; //latent Variables matrix [K,R] eta_s_rating; row_vector[R] eta_attitude_rating; matrix [K,L] eta_s_location; row_vector[L] eta_attitude_location; matrix [K,D] eta_s_domain; row_vector[D] eta_attitude_domain; real sigma_g_Sq; } transformed parameters{ real sigma_g; // matrix[N,R] theta_rat_mu; // matrix[N,L] theta_loc_mu; // matrix[N,D] theta_dom_mu; sigma_g = sqrt(sigma_g_Sq); // SD // theta_rat_mu = attitude* eta_attitude_rating + s *eta_s_rating; // theta_loc_mu = attitude* eta_attitude_location + s * eta_s_location; // theta_dom_mu = attitude* eta_attitude_domain + s * eta_s_domain; } model{ int grainsize =1; attitude ~ normal(0,1); for (i in 1:K) eta_s_rating[i] ~ normal(zero_R, one_R); for (i in 1:1) eta_attitude_rating[i] ~ normal(zero_R, one_R); for (i in 1:K) eta_s_location[i] ~ normal(zero_L, one_L); for (i in 1:1) eta_attitude_location[i] ~ normal(zero_L, one_L); for (i in 1:K) eta_s_domain[i] ~ normal(zero_D, one_D); for (i in 1:1) eta_attitude_domain[i] ~ normal(zero_D, one_D); sigma_g_Sq ~ inv_gamma(1, 1); // Driven by normal prior and categorical distribution for (i in 1:N){ target += reduce_sum(partial_sum_loc, location, grainsize, attitude, eta_attitude_rating, s, eta_s_rating); } for (i in 1:N){ target += reduce_sum(partial_sum_rat, rating, grainsize, attitude, eta_attitude_rating, s, eta_s_rating); } for (i in 1:N){ target += reduce_sum(partial_sum_dom, domain, grainsize, attitude, eta_attitude_rating, s, eta_s_rating); } } Warning message: In readLines(self$stan_file()) : incomplete final line found on '~/.cmdstanr/cmdstan-2.23-rc2/examples/parallel_test/tptrain_parallelV2.stan' Running MCMC with 2 chain(s) on 1 core(s)... Running ./tptrain_parallelV2 'id=1' random 'seed=1220086105' data \ 'file=/tmp/Rtmp2iZMR1/standata-28e71a72fc21.json' output \ 'file=/tmp/Rtmp2iZMR1/tptrain_parallelV2-202004202244-1-b542d4.csv' \ 'refresh=100' 'method=sample' 'num_samples=100' 'num_warmup=100' \ 'save_warmup=0' 'algorithm=hmc' 'engine=nuts' adapt 'engaged=1'