add trace plot and extract parameters
fix stan model formula
add random effect model