ObjectiveThis study aims to provide craniofacial features and incisor position prediction models among esthetics population after orthodontic treatment.MethodsA total of 338 patients with esthetic profiles were selected from 1 055 post-orthodontic adults based on inclusion and exclusion criteria and scored by five dental students. Mean values of cephalometric measurement were demonstrated with gender and sagittal skeletal pattern (ANB) differences assessed by generalized estimating equation. Correlations among age, skeletal, and soft tissue variables in different genders were determined. Multivariate analysis was used to assess the influence of skeletal and soft tissue variables on orthodontic incisor position.ResultsWithin the esthetic cohort after orthodontic treatment, gender differences existed in cranial and mandibular length, mandibular rotation, protrusion of pogonion, inclination of maxillary incisors, upper lip thickness, lower 1/3 face height, and nose protrusion. However, in class Ⅰ malocclusion, the development of the chin showed no significance, whereas the protrusion of the lower lip significantly differed between genders, with less protrusion in males. In males, the protrusion of lip and chin varied among groups, with thinner basic upper lip thickness of class Ⅱ than class Ⅰ and greater lower lip height of class Ⅲ than class Ⅰ. In females, the protrusion and height of the upper lip and chin varied in groups, with thicker basic upper lip thickness of class Ⅲ than class Ⅰ. Upper lip thickness, female basic upper lip thickness, and female lower lip height were correlated with age. The protrusion of the upper lip was mainly correlated with posterior cranial base, mandibular body length, ANB in male and SNA, ANB, sagittal maxillary length, mandibular rotation, and protrusion of pogonion in female. The protrusion of the lower lip in males was mainly correlated with ramus height, mandibular rotation, protrusion of pogonion in male and ANB, mandibular rotation, and protrusion of pogonion in females. After employing skeletal variables to fit dental parameters, optimal results were achieved in male IMPA (adjustedR2=0.712) and female G Vert-U1 (adjustedR2=0.795). After employing skeletal and soft tissue variables simultaneously, optimal results were achieved in male and female G Vert-U1, which was mainly influenced by the protrusion of the chin and nasal base in males (adjustedR2=0.836) and the protrusion of the nasal base as well as ANB in females (adjustedR2=0.842).ConclusionThis study demonstrated the craniofacial features of esthetics population after orthodontic treatment and revealed the correlation of gender, age, skeleton, and soft tissues to provide reliable prediction models for aesthetic orthodontic incisor position.