Dynamical modeling for non-Gaussian data with high-dimensional sparse ordinary differential equations