Overview of the method
We assume the non-Gaussian trait of interest follows a generalized linear mixed model of the form:
$$\begin{aligned} g(E(\pmb {y}|\pmb {\alpha }_1,\pmb {\alpha }_2))= & X_S \pmb {\beta }_S + X_c \pmb {\beta }_c + Z_1 \pmb…
We assume the non-Gaussian trait of interest follows a generalized linear mixed model of the form:
$$\begin{aligned} g(E(\pmb {y}|\pmb {\alpha }_1,\pmb {\alpha }_2))= & X_S \pmb {\beta }_S + X_c \pmb {\beta }_c + Z_1 \pmb…