I've run a probit regression in R with a random effect and can find no way to get the marginal effects with s.e. and p values. I have therefore tried to calculate the marginal effects 'by hand' by using the probit scalars and regression coefficients. However, I do not know how to get p values or standard errors and as far as I have found there is no easy way to do this for a mixed effects probit regression.
My model m1 is