Uncertainty quantification for hybrid random logistic models with harvesting via density functions