Data-Driven Modeling of Aggregate Flexibility Under Uncertain and Non-Convex Device Models