I recently came across the following method for portfolio optimization: Let $Y$ be a random variable that describes the returns of $n$ assets. Fix a constraint matrix $A \in \mathbb{R}^{m \times n}$ and $b \in \mathbb{R}^m$. Then, we sample from $Y$ and calculate the weights $x$ that maximize returns under the constraints $Ax \ge b$. Finally, we average over all optimal allocations to compute the final weights $\bar x$.
Since the set $\{ Ax \ge b\} $ is convex, we can be sure that $\bar x$ also satisfies the constraints. Thus, I think the practical idea is that $\bar x$ will be a good feasible compromise between extreme cases of $Y$ while still somewhat maximizing the objective function.
But from a statistical point of view, it is unclear to be why this would be a good procedure. Why don't we use the mean of $Y$ directly? Is this some kind of Bayesian approach?
Update
I think my statistical concerns were essentially adressed here and the method was discussed in the paper mentioned there (though with a different objective function including a quadratic term that reflects risk)