We focus on the Fisher information matrix used for design evaluation and optimization in nonlinear mixed effects multiple response models. We evaluate the appropriateness of its expression computed by linearization as proposed for a single response model. Using a pharmacokinetic-pharmacodynamic (PKPD) example, we first compare the computation of the Fisher information matrix with approximation to one derived from the observed matrix on a large simulation using the stochastic approximation expectation-maximization algorithm (SAEM). The expression of the Fisher information matrix for multiple responses is also evaluated by comparison with the empirical information obtained through a replicated simulation study using the first-order linearization estimation methods implemented in the NONMEM software (first-order (FO), first-order conditional estimate (FOCE)) and the SAEM algorithm in the MONOLIX software. The predicted errors given by the approximated information matrix are close to those given by the information matrix obtained without linearization using SAEM and to the empirical ones obtained with FOCE and SAEM. The simulation study also illustrates the accuracy of both FOCE and SAEM estimation algorithms when jointly modelling multiple responses and the major limitations of the FO method. This study highlights the appropriateness of the approximated Fisher information matrix for multiple responses, which is implemented in PFIM 3.0, an extension of the R function PFIM dedicated to design evaluation and optimization. It also emphasizes the use of this computing tool for designing population multiple response studies, as for instance in PKPD studies or in PK studies including the modelling of the PK of a drug and its active metabolite.