The natural origin of the fibers combined with random production flaws results in significant uncertainties in the properties of natural fiber reinforced composites. A probabilistic assessment can help to characterize the uncertainties and evaluate the reliability of natural fiber composites, enabling their use in engineering designs. Toward this end, this study aims to quantify the uncertainties in the tensile strength and frequency response of a unidirectional flax/epoxy composite due to the variability of various input parameters, including the fiber material properties and manufacturing flaws. Based on the available data in the literature, the non-deterministic input variables were divided into normal and uniform variables using a statistical test. A computationally efficient response surface approach based on the polynomial chaos expansion was adopted to conduct the uncertainty analysis with multi-type uncertain variables. Moreover, the results were validated by the direct Monte Carlo simulation to demonstrate the accuracy and efficiency of the surrogate model.