I modeled the declaration of tausq and tausq burnin after sigmasq

…eady

Once again for declaration of statements for tau, I copied those
statements for sigma closely.

This was modeled after get gibbs sample for sigsq closely.

Once again, I closely modified from sigsq

This is the main modification to get the random intercept parameters of
b_i and tausq. The estimated parameters are likely not consistent with
the truth since I have not been able to test the codes in R. However, we
have a developed a method that produces consistent estimates by using
bart from BayesTree to draw the regression trees and then drawing the
sigsqs, tausqs, and bis outside of bart. Hence, I should be able to make
the correct changes once I know how to get modifications on bartMachine
to work.

Similarly, copied closely from sigsq

Added tausq stuff by modifying from sigsq codes and for grouping
variables, copied cov_split_prior.

Updates so that R can recognize the addition of grouping variable and
initial values for tau.

Similarly to update R files so that groupings variable can be read into
the java files.

Random Intercept BART from Vincent