We consider the nonparametric estimation of the multivariate probability density function and its partial derivative with a support on [0, ∞) by dependent data. To this end we use the class of kernel estimators with asymmetric gamma kernel functions. The gamma kernels are nonnegative. They change their shape depending on the position on the semi-axis and possess good boundary properties for a wide class of densities. The theoretical asymptotic properties of the multivariate density and its partial derivative estimates like biases, variances and covariances are derived. We obtain the optimal bandwidth selection for both estimates as a minimum of the mean integrated squared error (MISE) assuming dependent data with a strong mixing. Optimal rates of convergence of the MISE both for the density and its derivative are found.