**Integrate GP Classification from GaussianProcess.jl with CRRao** **Possible code for user** ```Julia container = fit(y~x1+x2, dataframe, index_var = [t1, t2] , GP_claaification(), matern_kernel()) ```