<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def kde_statsmodels_m(data, grid, **kwargs):
 """ Multivariate Kernel Density Estimation with Statsmodels Parameters data : numpy.array Data points used to compute a density estimator. It has `n x p` dimensions, representing n points and p variables. grid : numpy.array Data points at which the desity will be estimated. It has `m x p` dimensions, representing m points and p variables. Returns ------- out : numpy.array Density estimate. Has `m x 1` dimensions """

    kde = KDEMultivariate(data, **kwargs)
    return kde.pdf(grid)