I am modeling a multivariate multilevel model (level 1: response (4 different responses; cross-sectional design); level 2: pupils (N = 3997); level 3: schools). For each predictor I have added separate coefficients for each response (interactions).
Some of my higher-level predictor variables (school characteristics) do not have a significant effect on all the response variables. Is it best to keep the non-significant interactions in the model or to exclude them?