the shape of the mandibular dental arch is an important element in diagnosis and treatment planning . as arch form configuration is a key requirement for stable treatment results , correct identification of arch form is essential during the diagnostic phase , especially given the popularity of superelastic archwires.1,2,3,4 there is no single ideal arch form that meets all treatment goals due to the inevitable variations between individuals.5,6,7 mandibular arch forms have been classified by various mathematical and statistical methods.8,9,10,11,12,13,14,15,16 for example , raberin et al.17 classified the mandibular arch into 5 ideal forms by k - means cluster analysis . however , most of the researchers selected landmarks on incisal edges and cusp tips , which may not appropriately represent arch forms . the validity of three - dimensional ( 3d ) virtual models has been confirmed in several studies.18,19,20,21 nevertheless , few studies have used this technology to evaluate dental arch forms.20,22,23 given the substantial overlap in the identification of mandibular arch forms , a revised template based on facial axis ( fa ) points has been suggested.22 recently , bayome et al.24 proposed a clinical classification of dental arch forms by using fa points derived from 3d models . a previous study evaluated the morphologic characteristics of lingual arch form in a limited sample size.25 although the average arch form of a normal occlusion sample was provided , a classification of the lingual arch forms according to shape or size has not been attempted . lombardo et al.26 classified lingual arch forms into three groups according to the x- and y - values of each tooth . such two - dimensional evaluation has inherent errors in landmark identification and addresses size but not shape . the purposes of this study were 1 ) to classify lingual dental arch form types based on the lingual bracket points and 2 ) to provide a new lingual arch form template based on this classification for clinical application through the analysis of 3d virtual models of normal occlusion sample . the sample consisted of 115 young adults with normal occlusion ( 64 males and 51 females ) selected from a database of 15,621 korean subjects . the subjects ' age ranged from 20 to 30 years , with a mean of 23.5 years . the inclusion criteria were as follows : 
 angle 's class i molar and canine relationships0 < anb < 40 mm < overbite and overjet < 4 mmminor tooth size - arch length discrepancy ( < 3 mm crowding , < 1 mm spacing)flat or slight curve of spee ( < 2 mm)absence of crossbite or dental midline deviationpermanent dentition with normal tooth size and shape , except third molarsno previous orthodontic treatmentno restorations involving contact areas , cusp tips , incisal edges , or facial surfaces 
  angle 's class i molar and canine relationships 0 mm < overbite and overjet < 4 mm minor tooth size - arch length discrepancy ( < 3 mm crowding , < 1 mm spacing ) flat or slight curve of spee ( < 2 mm ) absence of crossbite or dental midline deviation permanent dentition with normal tooth size and shape , except third molars no previous orthodontic treatment no restorations involving contact areas , cusp tips , incisal edges , or facial surfaces the institutional review board of the catholic university of korea ( mc12sssi0083 ) approved the study . mandibular and maxillary casts were placed in the occluded relationship and scanned with a 3d laser scanner ( orapix kod-300 ; orapix co. , ltd , seoul , korea ) , set at 20-m resolution . lingual bracket points , which are the points on the lingual surface of the teeth on which the lingual brackets are centered , were digitized on each tooth of the upper and lower virtual models using rapidform 2006 software ( inus technology , inc . , the transverse and anteroposterior directions were set as the x - axis and y - axis , respectively ; the plane perpendicular to these axes was set as the z - axis . for reorientation of casts , the original x - axis of each cast was adjusted so that the mean inclination of line a ( line connecting the right and left contact points between the first and second premolars ) and line b ( straight line connecting the right and left contact point between the second premolars and the first molars ) became parallel to the original x - axis . the difference in z - axis values between the teeth was very small in each case so they were nullified for simplification of the mathematical calculations . six linear measurements and three ratios were evaluated per each upper and lower arch form ( table 1 and figure 1 ) . intraclass correlation coefficient ( icc ) test revealed high reliability between the two assessments ( icc > 0.8 ) . the 115 casts were divided randomly into 2 datasets : dataset 1 ( 68 cases ) was used to classify arch forms , and dataset 2 ( 47 cases ) was used to examine the differences between clusters for validation purposes . k - means cluster analysis was performed with dataset 1 to classify the samples into arch form types using spss ( spss ver . clustering was based on mandibular intercanine , interpremolar , and intermolar widths and width / depth ratios . potential classifications into three to nine groups were run through the software . after inspecting the distributions in each classification and homogeneity of group sizes , the classification that demonstrated highly significant intergroup differences in the six parameters was chosen . chi - square test was applied to evaluate the frequency distribution of the clusters and association between clusters and gender . then , a multivariate test with bonferroni adjustment for multiple comparisons was performed to evaluate the properties of the clusters . the mean arch form of each cluster was calculated from the x- and y - coordinates of the lingual bracket points of each subject in the group . the mean coordinates of each point were inputted into a mathematical program ( matlab 7.5 [ r2007b ] ; mathworks , inc . , natick , ma , usa ) to draw the best - fit curve that represented the mean arch form , using different types of polynomial equations . subsequently , dataset 2 was mapped into the same space and assigned to the 4 clusters already constructed . a multivariate test with bonferroni adjustment for multiple comparisons then , members of both datasets in each cluster were combined and matlab was used to draw the arch forms . the mandibular arch forms were superimposed on the origin point and differences between arch forms were evaluated for clinical significance and application . pearson correlation analysis was applied to assess the relationship between maxillary and mandibular arch form variables . k - means cluster analysis was performed seven times to classify dataset 1 into three to nine clusters . the classification into 4 clusters showed the maximum intercluster distance in the 6 parameters and had the most homogeneous cluster size . chi - square test revealed an even frequency distribution of the subjects between the clusters ( p = 0.572 ) , and no significant association was noted between gender and clusters ( p = 0.304 ) . table 3 compares the arch dimensions of the 4 cluster groups using dataset 2 . in the mandible , the intercanine and intermolar depths showed significant differences ( p = 0.025 and 0.015 , respectively ) . however , there was no significant difference in the premolar depth . in the maxilla , figure 2 depicts the superimposition of the best - fit curves of the mandibular arches of the 4 clusters at the origin point . cluster 2 ( wide type ) has broad intercanine and intermolar widths ; cluster 3 ( narrow type ) illustrates narrow intercanine and intermolar widths ; cluster 1 ( tapering type ) demonstrates broad intermolar width and narrow intercanine width ; cluster 4 ( ovoid type ) shows relatively narrow intermolar width and broad intercanine width . figure 3 presents the suggested mandibular template of the 4 arch forms in 1:1 size . however , the best - fit curves of the maxillary arch forms had a high sum of square error ( sse ) due to the position of the canines ( figure 4 ) . currently , there is an enormous and growing interest in invisible braces due to the increased numbers of adult patients and social pressures in today 's society . despite improvements to the lingual orthodontic protocol , some patients experience discomfort and dysfunction.27,28 these problems could be minimized by selecting the most appropriate arch form . therefore , it is necessary to have a classification for lingual arch forms that corresponds with the different groups of patients . the existing arch wire templates lack good fit to the average dental arches and have substantial overlap among their classifications , thus several arch form classifications derived from fa points have been reported.2,22,24,29 therefore , the aim of this study was to introduce a new improved classification system of arch form based on the lingual bracket points to improve clinical application of lingual preformed arch wires . recently , customized orthodontic treatment has recently gained recognition because of the wide variation between individuals . as a result , cluster analysis has gained popularity in biomedical science as it considers both individual variation and better control of the outliers . in our study , a new , four - type classification of lingual arch form was found that considers size and shape . the classification was based on mandibular widths and width / depth ratios which were reported to be important for the stability of the treatment outcomes and represents the curvature of each segment of the arch.1,2,3 the validity of using the six parameters as determinants in k - means cluster analysis was tested . several variables were added as determinants in different combinations , but the clustering membership and homogeneity showed no significant changes . meanwhile , removing any of those six variables led to different results . also , the significant strong and moderate correlations found between the maxillary and mandibular parameters supported the inclusion of the mandibular parameters as the only clustering factors . hence , this classification considered both the shape and the size of the dental arch . our results support previous findings that there is no single , ideal arch form appropriate for all patients.2,17,24,29 the 4-cluster classification developed in our study might be more practical due to the difficulty in customizing superelastic archwires . the width / depth ratios of this group were somewhat similar to those of the narrow cluster of lombardo et al.26 the remaining three groups in our results had larger arch widths and ratios but smaller intercanine depth and similar intermolar depth to all the forms proposed by lombardo et al.26 these differences may be explained by differences in ethnicity and research methods . this gender homogeneity is consistent with the findings of raberin et al.17 this study could have been improved by using a larger sample , which might have also revealed rare arch forms . several studies recommended the use of fourth- or sixth - degree polynomial curves to represent the labial dental arch.2,8,13,14,15,30 in our study , several levels of polynomial equations were applied : the ninth - degree polynomial equation showed the least sse and , hence , better fit ; the resultant curves are smooth enough for clinical purposes . the description of the maxillary arch shapes might not be sufficiently accurate for clinical use because of the high sse even with the ninth - degree polynomial equation . further studies are recommended to evaluate the maxillary arch by combining two mathematical equations : one for the anterior segment and the other for the posterior area . clinically , the maxillary archwire can be assessed by excluding the canines from the arch curve equation and then adding them as an offset bend to the archwire . studies of the long - term stability of treatment outcomes based on the proposed classification are also recommended . a new classification of lingual arch form has been identified that considers arch size and shape . wide , narrow , tapering , and ovoid types were described according to the intercanine and intermolar widths . a new template of these arch forms was introduced for more efficient identification of patient arch forms . this 4-type template might be more practical for clinical choice of mandibular lingual superelastic wires .
objectivethe purposes of this study were 1 ) to classify lingual dental arch form types based on the lingual bracket points and 2 ) to provide a new lingual arch form template based on this classification for clinical application through the analysis of three - dimensional virtual models of normal occlusion sample.methodsmaxillary and mandibular casts of 115 young adults with normal occlusion were scanned in their occluded positions and lingual bracket points were digitized on the virtual models by using rapidform 2006 software . 
 sixty - eight cases ( dataset 1 ) were used in k - means cluster analysis to classify arch forms with intercanine , interpremolar and intermolar widths and width / depth ratios as determinants . 
 the best - fit curves of the mean arch forms were generated . 
 the remaining cases ( dataset 2 ) were mapped into the obtained clusters and a multivariate test was performed to assess the differences between the clusters.resultsfour-cluster classification demonstrated maximum intercluster distance . 
 wide , narrow , tapering , and ovoid types were described according to the intercanine and intermolar widths and their best - fit curves were depicted . 
 no significant differences in arch depths existed among the clusters . 
 strong to moderate correlations were found between maxillary and mandibular arch widths.conclusionslingual arch forms have been classified into 4 types based on their anterior and posterior dimensions . 
 a template of the 4 arch forms has been depicted . 
 three - dimensional analysis of the lingual bracket points provides more accurate identification of arch form and , consequently , archwire selection .