There are a number of directories (Construct, Hidden, Inverse, Order and Symmetric). Each directory has subdirectories (IStudent, Normal, Student). The code for each sampler is called using: [r_burn,r_inf,model,details]=main(Y,X,n_burn,t_burn,n_inf,sampling,seed) where: Y is (n*p), X is (n*m), n_burn is no of samples used for burnin, t_burn is no of times used to select burn_in parameters, n_inf is no of samples in the inference phase, sampling is the spacing used to collect samples and seed is the seed to the random no generator.