High-order moment methods for LRFD including random variables with unknown probability distributions