“…Based on the joint GEE model formulation [ 40 ], we construct the multivariate linear model for describing the association relationship between K correlated traits and genetic variants, which is given as follows: where g −1 (•) is the inverse function of g (•) and is a response-specific link function [ 40 ], μ i = ( μ i 1 T , μ i 2 T ,⋯, μ iK T ) T is the ( n i × K ) × 1 vector of the expected mean of the multivariate traits y i = ( y i 1 T , y i 2 T ,⋯, y iK T ) T , α = ( α 1 T , α 2 T ,⋯, α K T ) T is the (( q + 1) × K ) × 1 vector of regression coefficients of the ( q + 1) nongenetic covariates for the K correlated traits, and β = ( β 1 T , β 2 T ,⋯, β K T ) T is the ( p × K ) × 1 vector of regression coefficients of the p genetic variants for the K correlated traits.…”