powered by
This function will be used to calculate the P(coeficient variable cluster <=0)
computeprob(func, k)
is the inla marginals of the model parameter
is the cutoff
probability model coefficient <=k