cc_byEvers, Michael P.Kontny, Markus2025-11-242025-11-242025https://doi.org/10.1002/mma.11226https://hohpublica.uni-hohenheim.de/handle/123456789/18282We provide a generalization of Faà di Bruno’s formula to represent the 𝑛-th total derivative of the multivariate and vector-valued composite 𝑓 ∘𝑔. To this end, we make use of properties of the Kronecker product and the 𝑛-th derivative of the left-composite 𝑓 , which allow the use of a multivariate and matrix-valued form of partial Bell polynomials to represent the generalized Faà di Bruno’s formula. We further show that standard recurrence relations that hold for the univariate partial Bell polynomial also hold for the multivariate partial Bell polynomial under a simple transformation. We apply this generalization of Faà di Bruno’s formula to the computation of multivariate moments of the normal distribution.engCommutation matricesGeneralized Faà di Bruno's formulaMultivariate and matrix‐valued Bell polynomials𝑛-th total derivative of multivariate and vector-valued composite functions510Article2025-11-04