A closure approximation for liquid-crystalline polymer models based on parametric density estimation