given the types of intermittent maximum effort needed during a team sport such as handball , it is important to evaluate performance by measuring variables such as muscle strength , velocity and resistance ( marques , 2010 ; michalsik et al . , 2013 ; pvoas et al . , the factors that determine performance in handball are explosive actions such as jumping ( i.e. , vertical or horizontal ) , speed of reduced - space displacements ( i.e. , 10 and 20 m sprint time ) and ball release velocity ( brv ) . to assess physical performance of a handball player , several tests are commonly used encompassing sprinting , jumping , agility and maximal strength as often measured in the one repetition maximum ( 1rm ) , or ball release velocity tests ( marques and gonzlez  badillo , 2006 ) . thus , according to the scores obtained in these tests , athletes can be classified in terms of a particular skill ( e.g. , the player who is able to throw the ball faster or lift the heaviest weight in the bench press ) ( ingebrigtsen et al . , 2013 ; vila et al . , 2012 ) the problem with this type of classification is that , despite an existing correlation between the maximum weight lifted in the bench press and ball throwing velocity ( gorostiaga et al . , 2005 ; granados et al . , 2007 ; marques et al . , 2007 ) , the player who throws the ball more rapidly will not necessarily be the individual attaining the highest 1rm in the bench press . the second objective of such tests is to compare the effects of training programs on performance during a given period of time ( gorostiaga et al . , 1999 ) or during a competition season ( bresciani et al . , 2010 ; ; marques et al . , 2007 ) , the player who throws the ball more rapidly will not necessarily be the individual attaining the highest 1rm in the bench press . the second objective of such tests is to compare the effects of training programs on performance during a given period of time ( gorostiaga et al . , 1999 ) or during a competition season ( bresciani et al . , 2010 ; a search of the scientific literature indicates that such measures of performance are commonly used to determine the response to a given training intervention ( chelly et al . , 2009 , 2013 ; the researcher looks for differences between pre- and post test values of a certain variable . moreover , test scores may be used to measure the effect size produced by a training intervention on a study population ( peterson et al . , 2004 ; rhea et al . , 2003 when used in this manner , the cross - sectional perspective of these tests is completely lost . in other words , this information could be of vital importance to measure the physical state of a player before a training program compared with the rest of the players . this line of research has generated many studies in different sports such as soccer ( unnithan et al . , 2012 ; votteler and hner , 2013 ; william and reilly , 2000 ) , handball ( baker et al . , 2013 ; lidor et al . , 2005 ; mohamed et al . , 2009 ) , basketball ( hoare , 2000 ) or rugby ( gabbett , 2006 , 2011 ) in which the principle objective was to differentiate elite from amateur players based on variables related to physical performance , anthropometric measures or decision making ability ( 2011 ) developed a multidisciplinary model to identify talented young handball players ( under-14 [ u14 ] , under-16 [ u16 ] and under-18 [ u18 ] ) . results indicated that elite players could be distinguished from non - elite ones according to physical performance characteristics , although the discrimination capacity of this method varied according to the age group tested . ( 2011 ) clustering of handball players occurred naturally depending upon their respective team ( elite vs. non - elite team ) and age ( u14 , u16 and u18 ) . using multivariate statistics ( i.e. , cluster analysis ) and theoretical knowledge about performance , it is possible to create an objective physical performance classification system of the handball players . therefore , the aim of the present study was to identify different cluster groups of handball players according to their physical performance level evaluated in a series of physical assessments ( a countermovement jump [ cmj ] , a squat jump [ sj ] , 10 m sprint time [ st10 ] , 20 m sprint time [ st20 ] , the one repetition maximum [ 1rmest ] , peak power at the propulsive phase [ ppeakmp ] and mean power at the propulsive phase [ pmppmp ] in the bench press and ball release velocity [ brv ] ) , which could then be used to design a training program based on individual strengths and weaknesses , and to determine which of these variables best identified elite performance in a group of under-19 [ u19 ] national level handball players . participants age , body mass and height ( mean  sd ) were 18  0.4 years , 87.38  9.97 kg and 189  6.72 cm , respectively . written informed consent was obtained from the participants after the nature of the study and its objectives had been explained . the study protocol received university of granada review board approval and was performed in accordance with the tenets of the declaration of helsinki . a cross - sectional design was performed to measure physical performance variables in young handball players . this investigation was conducted one month prior to the european junior handball championship in 2013 . using a multivariate approach , a hierarchical cluster analysis was performed to create a grouping of training based on different important training variables in handball players . after that , a descriptive discriminant analysis was performed to indicate the most important variable to generate this grouping classification . the measured variables included in the analysis were in vertical dimension ( cmj and sj ) , in horizontal dimension ( st10 and st20 ) and cross - sectional dimension ( 1rmest , ppmppmp and pmppmp in the bench press and brv ) . the physical tests were performed during the training camp for the european junior handball championship in 2013 . sprint time , brv , cmj , sj and incremental - load bench press tests were applied since these are regularly used in testing physical performance in handball and measure variables such as muscle strength , power and velocity in the upper and lower body very well ( granados et al . , 2008 ; ingebrigtsen et al . , 2013 ; krger et al . , 2014 ; marques and gonzlez badillo , 2006 ) . these physical performance tests were conducted at the team s playing court and training facilities . before the tests , participants attended a familiarization session in which each test protocol and its key points were explained in detail . all participants undertook their standard 15 min warm up protocol involving sprints with direction changes and ball passes . after a 5 min rest period in which the subjects could drink water ad libitum , they completed the physical performance assessment tests established for the session . the tests were performed over 3 days as follows : day 1 , st10 and st20 and brv ; day 2 , cmj and sj ; and day 3 , an incremental - load bench press test . at least 24 h elapsed between each testing session . 
 sprint tests : a photocell system ( kit racetime 2 light radio , microgate , italy ) was used to measure the time ( s ) needed to cover 10 m and 20 m distances . the starting point was 1 m before the position of the photocell so that each subject selected the moment of starting . each subject was allowed three attempts at the test and players received feedback of the score obtained in every sprint . descriptive variables recorded ( provided as their mean , standard deviation , intraclass correlation coefficient [ icc2.k ] , mean 95% ci icc2.k and repeated measures analysis of variance [ rm anova ] ) . value used was the mean of repetitions 2 and 3 ( see procedures for details ) ppeakmp = peak power in the maximum power load . = intraclass correlation coefficient ( icc2.k ) = only one measure of this variable used . vertical jump tests : the vertical jump height tests ( i.e. , cmj and sj ) were performed using a jumping mat ( optojump , microgate , italy ) , which measures the time the feet are off the mat and translates the result into a jump height in cm . the technique was explained before each test and each athlete had 3 attempts at each jump . the sj was started from a static semisquatting position with the knee angle of 90 degrees of flexion , followed by subsequent action , during which the leg and hip extensor muscles contracted concentrically . at the beginning of the cmj test , each participant stood erect on the jump mat . the test consisted of preparatory movement down to approximately 90 degrees of knee flexion , stretching the leg extensor muscles followed by explosive maximal extension in the opposite direction . for the cmj , the mean height of the last 2 attempts was included in the analysis while for the sj , the mean of the last 3 attempts ( i.e. , repetitions 2 and 3 ) were used ( table 1 ) . ball release velocity test : brv was determined with a three - pace run up . peak velocity was measured using a radar stalker ats gun ( radar sales , minneapolis , mmn , usa ) positioned behind the goal . just after each throw the mean of repetitions 2 and 3 was entered in the analysis ( table 1 ) . upper body force - velocity test : in the last testing session , athletes performed an incremental - load free bench press test ( salter ) . the bench position , barbell grip and displacement were previously established and controlled in each attempt or repetition . participants performed a 15 min warm up and 2 sets of presses with a light load ( i.e. , barbell weight of 20 kg ) before starting the test protocol . the starting weight for the test was 25 kg and was increased by 10 kg up to a load of 85 kg . three repetitions were performed at each intensity . during the incremental test , the descent phase of the bar was controlled by verbal instructions from the researcher as follows :  down  ( 2 s ) ,  chest ( 1 s ) ,  go  . to avoid the rebound effect , the each subject was encouraged to lift the barbell as quickly as possible during the concentric phase ( i.e. , explosively ) until the elbows were completely extended . a linear displacement system was employed ( t - force system , ergotech , murcia , spain ) to record the mean velocity of the propulsive phase , mean power of the propulsive phase and the estimated 1 repetition maximum ( 1rmest ) . it had been shown that mean velocity could be used to precisely estimate loading intensity ( gonzlez - badillo and snchez - medina , 2010 ) . information of the movement velocity of the bar is transmitted and transformed from analogue to digital by a data acquisition board . vertical instantaneous velocity was directly sampled by the device at that the frequency of 1000 hz . the propulsive phase was defined as the portion of the concentric phase during which the measured acceleration ( a ) was greater than acceleration due to gravity ( i.e. , a > -9.8 ms ) ( gonzlez - badillo and snchez - medina , 2010 ) . once the test was completed and the data processed , the load corresponding to the mean maximum power ( loadmp ) was determined for each individual . then , for this load , we calculated the variables of mean power of the propulsive phase ( pmppmp ) and the peak power ( ppeakmp ) . the variables were transformed into the z - score in order to analyse the presence of outliers . the intrasession intraclass correlation coefficient was calculated for all tests ( except for the 1rmest variable ) using the equation 2.k ( icc2.k ) ( weir , 2005 ) . the coefficient of variation ( cv ) was calculated as 100(e  1 ) ( hopkins , 2011 ) . the strength of a correlation was defined as r = 0  0.1 ( trivial ) , r = 0.1  0.3 ( small ) , r = 0.3  0.5 ( moderate ) , r = 0.5  0.7 ( large ) , r = 0.7  0.9 ( very large ) and r = 0.9  1 ( nearly complete ) ( hopkins , 2003 ) . a repeated measures analysis of variance ( rm anova ) was used to assess differences in the means recorded for the three repetitions of each dependent variable . when the sphericity test was violated , a greenhouse - geisser estimation was used to adjust the degree of freedom . we used hierarchical cluster analysis to create several groupings according to the results of the performance tests . each variable was transformed into the z - score to avoid effects of different measurement scales . as a measure of similarity , the clustering ward s method was used . a two between group multivariate analysis of variance ( manova ) was used for upper trunk ( 1rmest , ppeakmp , pmppmp and brv ) and lower trunk ( st10 , st20 , cmj and sj ) variables , where the predictor variable was the cluster solution ( i.e. , cluster 1 , 2 , 3 and 4 ) . the players were classified into different clusters through a descriptive discriminant analysis . as the criterion variable , we used the 4-cluster classification generated by the cluster analysis , while lower trunk ( st10 , st20 , cmj and sj ) and upper trunk ( 1rmest , ppeakmp , pmppmp and brv ) variables were entered as the predictor variables . all statistical analyses were performed using spss software ( spss inc , illinois , usa ) . participants age , body mass and height ( mean  sd ) were 18  0.4 years , 87.38  9.97 kg and 189  6.72 cm , respectively . written informed consent was obtained from the participants after the nature of the study and its objectives had been explained . the study protocol received university of granada review board approval and was performed in accordance with the tenets of the declaration of helsinki . a cross - sectional design was performed to measure physical performance variables in young handball players . this investigation was conducted one month prior to the european junior handball championship in 2013 . using a multivariate approach , a hierarchical cluster analysis was performed to create a grouping of training based on different important training variables in handball players . after that , a descriptive discriminant analysis was performed to indicate the most important variable to generate this grouping classification . the measured variables included in the analysis were in vertical dimension ( cmj and sj ) , in horizontal dimension ( st10 and st20 ) and cross - sectional dimension ( 1rmest , ppmppmp and pmppmp in the bench press and brv ) . the physical tests were performed during the training camp for the european junior handball championship in 2013 . sprint time , brv , cmj , sj and incremental - load bench press tests were applied since these are regularly used in testing physical performance in handball and measure variables such as muscle strength , power and velocity in the upper and lower body very well ( granados et al . , 2008 ; these physical performance tests were conducted at the team s playing court and training facilities . before the tests , participants attended a familiarization session in which each test protocol and its key points were explained in detail . all participants undertook their standard 15 min warm up protocol involving sprints with direction changes and ball passes . after a 5 min rest period in which the subjects could drink water ad libitum , they completed the physical performance assessment tests established for the session . the tests were performed over 3 days as follows : day 1 , st10 and st20 and brv ; day 2 , cmj and sj ; and day 3 , an incremental - load bench press test . sprint tests : a photocell system ( kit racetime 2 light radio , microgate , italy ) was used to measure the time ( s ) needed to cover 10 m and 20 m distances . the starting point was 1 m before the position of the photocell so that each subject selected the moment of starting . each subject was allowed three attempts at the test and players received feedback of the score obtained in every sprint . descriptive variables recorded ( provided as their mean , standard deviation , intraclass correlation coefficient [ icc2.k ] , mean 95% ci icc2.k and repeated measures analysis of variance [ rm anova ] ) . value used was the mean of repetitions 2 and 3 ( see procedures for details ) ppeakmp = peak power in the maximum power load . = intraclass correlation coefficient ( icc2.k ) = only one measure of this variable used . vertical jump tests : the vertical jump height tests ( i.e. , cmj and sj ) were performed using a jumping mat ( optojump , microgate , italy ) , which measures the time the feet are off the mat and translates the result into a jump height in cm . the technique was explained before each test and each athlete had 3 attempts at each jump . the sj was started from a static semisquatting position with the knee angle of 90 degrees of flexion , followed by subsequent action , during which the leg and hip extensor muscles contracted concentrically . at the beginning of the cmj test , each participant stood erect on the jump mat . the test consisted of preparatory movement down to approximately 90 degrees of knee flexion , stretching the leg extensor muscles followed by explosive maximal extension in the opposite direction . for the cmj , the mean height of the last 2 attempts was included in the analysis while for the sj , the mean of the last 3 attempts ( i.e. , repetitions 2 and 3 ) were used ( table 1 ) . ball release velocity test : brv was determined with a three - pace run up . peak velocity was measured using a radar stalker ats gun ( radar sales , minneapolis , mmn , usa ) positioned behind the goal . just after each throw the mean of repetitions 2 and 3 was entered in the analysis ( table 1 ) . upper body force - velocity test : in the last testing session , athletes performed an incremental - load free bench press test ( salter ) . the bench position , barbell grip and displacement were previously established and controlled in each attempt or repetition . participants performed a 15 min warm up and 2 sets of presses with a light load ( i.e. , barbell weight of 20 kg ) before starting the test protocol . the starting weight for the test was 25 kg and was increased by 10 kg up to a load of 85 kg . three repetitions were performed at each intensity . during the incremental test , the descent phase of the bar was controlled by verbal instructions from the researcher as follows :  down  ( 2 s ) ,  chest  ( 1 s ) ,  go  . to avoid the rebound effect , the each subject was encouraged to lift the barbell as quickly as possible during the concentric phase ( i.e. , explosively ) until the elbows were completely extended . a linear displacement system was employed ( t - force system , ergotech , murcia , spain ) to record the mean velocity of the propulsive phase , mean power of the propulsive phase and the estimated 1 repetition maximum ( 1rmest ) . it had been shown that mean velocity could be used to precisely estimate loading intensity ( gonzlez - badillo and snchez - medina , 2010 ) . information of the movement velocity of the bar is transmitted and transformed from analogue to digital by a data acquisition board . vertical instantaneous velocity was directly sampled by the device at that the frequency of 1000 hz . the propulsive phase was defined as the portion of the concentric phase during which the measured acceleration ( a ) was greater than acceleration due to gravity ( i.e. , a > -9.8 ms ) ( gonzlez - badillo and snchez - medina , 2010 ) . once the test was completed and the data processed , the load corresponding to the mean maximum power ( loadmp ) was determined for each individual . then , for this load , we calculated the variables of mean power of the propulsive phase ( pmppmp ) and the peak power ( ppeakmp ) . the variables were transformed into the z - score in order to analyse the presence of outliers . the intrasession intraclass correlation coefficient was calculated for all tests ( except for the 1rmest variable ) using the equation 2.k ( icc2.k ) ( weir , 2005 ) . the coefficient of variation ( cv ) was calculated as 100(e  1 ) ( hopkins , 2011 ) . pearson correlation coefficients were used to determine interrelationships between variables . the strength of a correlation was defined as r = 0  0.1 ( trivial ) , r = 0.1  0.3 ( small ) , r = 0.3  0.5 ( moderate ) , r = 0.5  0.7 ( large ) , r = 0.7  0.9 ( very large ) and r = 0.9  1 ( nearly complete ) ( hopkins , 2003 ) . a repeated measures analysis of variance ( rm anova ) was used to assess differences in the means recorded for the three repetitions of each dependent variable . when the sphericity test was violated , a greenhouse - geisser estimation was used to adjust the degree of freedom . we used hierarchical cluster analysis to create several groupings according to the results of the performance tests . each variable was transformed into the z - score to avoid effects of different measurement scales . as a measure of similarity , a two between group multivariate analysis of variance ( manova ) was used for upper trunk ( 1rmest , ppeakmp , pmppmp and brv ) and lower trunk ( st10 , st20 , cmj and sj ) variables , where the predictor variable was the cluster solution ( i.e. , cluster 1 , 2 , 3 and 4 ) . the players were classified into different clusters through a descriptive discriminant analysis . as the criterion variable , we used the 4-cluster classification generated by the cluster analysis , while lower trunk ( st10 , st20 , cmj and sj ) and upper trunk ( 1rmest , ppeakmp , pmppmp and brv ) variables were entered as the predictor variables . all statistical analyses were performed using spss software ( spss inc , illinois , usa ) . table 1 provides the means (  sd ) , the 95% confident interval ( ci ) of the mean , icc 2.k and 95% ci icc2.k of the variables recorded ( lower trunk [ st10 , st20 , cmj and sj ] and upper trunk [ 1rmest , ppeakmp , pmppmp and brv ] ) . in parallel , the 3 repetitions of each test and the variables examined were compared by rm anova . significant differences were detected for the cmj ( f [ 1.43 , 21.41 ] = 2.2 , p = .021 , p = .13 ) and brv ( f [ 2 , 30 ] = 8 , p = .004 , p = .35 ) ; post hoc bonferroni correction indicated differences between set 1 with 2 ( p = .01 ) for the cmj variable and set 1 vs set 2 ( p = .020 ) and set 1 vs. set 3 ( p = .023 ) in brv . thus , we used the mean of set 2 and 3 for the final values of brv and the cmj . for the rest of the variables the mean of the three attempts pearson correlation coefficients between brv and bench press variables ( 1rmest , ppeakmp and pmppmp ) were respectively .68 ( p < .01 ) large , .55 ( p < .05 ) large and 0.54 ( p < .05 ) large . the hierarchical cluster analysis dendrogram constructed before examining any groupings is provided in figure 1 . in the analysis , each subject was first considered as an individual group i.e. , 16 subjects = 16 groups , in which each of the nine variables was analysed . then , differences between all subjects were analysed with each subject assigned to the same group in a cluster matrix . dendrogram generated by the cluster analysis . after creating a grouping classification using a cluster analysis , the players were classified in 4 clusters . the means  sds of all the variables of the grouped solution generated by the cluster analysis are provided in table 2 . mean sd of the performance indicators according to the clusters generated by the hierarchical cluster analysis . first multivariate statistical technique used ( i.e. , cluster analysis ) provided a grouping classification of the players based on a measure of similarity ( i.e. , squared euclidean distance ) and the agglomeration method ( i.e. , wards method ) . , a manova showed significant differences in the lower (  = .07 ; p = .004 ; p = .59 ) and upper trunk (  = .03 ; p = .0001 ; p = .69 ) . figure 2 shows z - score values of each variable ( i.e. , the upper and lower trunk ) at each cluster solution obtained . z - scores of each variable at both the upper ( a ) and lower ( b ) trunk at each cluster . the following values were obtained from discriminant analysis : function 1 (  = .001 ; p = .0001 ) , function 2 (  = .044 ; p = .013 ) and function 3 (  = .24 ; p = .046 ) . the percentage of variance explained by each of these functions was 75% , 15% and 10% , respectively . combined intra - group correlations between the discriminant variables and the typified canonical discriminant functions were r = .44 and -.39 ( brv and st20 for function 1 , respectively ) , r = .50 ( 1rmest for function 2 ) and r = .48 , r = .47 , r = -.42 , -.36 and .22 ( pmpfmp , ppeakmp , cmj , sj and st10 for function 3 , respectively ) . these variables at each canonical function reflected the most important variable to classify players in the cluster solutions . discriminant analysis provided a probability of each player and their percentage probabilities of belonging to each cluster . this analysis revealed that 100% of cases were correctly classified in the 4-group solution generated by the cluster analysis using the performance variables assessed . figure 3 shows the mean of the centroids ( square ) and the groupings of different players in their corresponding cluster ( circle ) . graphical representation of each player ( circle ) and the centroid mean for each group ( square ) . a qualitative analysis showed that players were distributed according to their playing positions as follows : cluster 1 ( two center and two right back players ) ; cluster 2 ( one center , one left wing , two left backs , one right back and one pivot ) ; cluster 3 ( one goalkeeper , one right back and one left wing ) and cluster 4 ( two pivots and one goalkeeper ) . the aim of the present study was to identify different cluster groups of handball players according to their physical performance level assessed in a series of tests , which could then be used to design a training program based on the players individual needs . moreover , our study sought to determine the variables that best reflected elite performance of a group of male handball players in the ui9 national team . to the best of our knowledge , this is the first attempt to create a classification scheme that analyses individual strengths and weaknesses in elite handball players based on their physical performance using a multivariate statistical approach . several studies have compared the performance of individual athletes considering a given skill or variable ( e.g. , brv , 1rm in the bench press , rfd , repeated sprint ability ) ( krger et al . , 2014 ; marques et al . , 2011 ; michalsik et al . , physiologic variables interact and can obscure which factors are the real determinants of performance . in handball , factors such as brv ( marques et al . , 2007 ) , training with external resistance ( hermassi et al . , 2011 ) or combining different training modalities ( van den tillaar , 2004 ) seem to play an important role in the physical performance of players . in the present study , physical performance was assessed at two levels : the upper trunk ( 1rmest , ppeakmp , pmppmp and brv ) and the lower trunk ( cmj , sj , st10 and st20 ) . before classifying players according to their physical performance , a cluster analysis was performed to identify individuals obtaining similar scores in the performance tests used , i.e. , according to their strengths and weaknesses . as it may be seen in the dendrogram generated by the cluster analysis ( figure 1 ) , we opted for a 4-cluster solution . table 2 provides the means and sds of each performance variable analysed according to the cluster each player was assigned to . this information could be useful in order to create different training groups based on physical weakness points . in other words , using this statistical approach makes it possible to provide specific training to each group of players and enhance the physical performance of the players based on their specific needs . this allows for the identification of strengths and weaknesses in the three physical assessment dimensions ( e.g. , cluster 3 should focus on improving cross - sectional and horizontal dimension performance and cluster 2 should improve vertical dimension performance ) ( figure 2 ) . therefore , once the identification of strengths and weaknesses was performed , it was possible to design specific training sessions for each cluster solution . it should be mentioned that the cluster analysis allocates the same weight to each performance variable examined . thus , once all the players were grouped according to their physical performance , a discriminant analysis was used to determine whether the players could be classified in the established groups according to the performance variables tested and to check the specific weight of each variable for this purpose . this classification is graphically shown in figure 3 , which represents the organization of the 16 players assessed . in figure 3 it may be observed how the cluster analysis generated four extremely homogeneous groups in themselves , but these were heterogeneous compared to each other ( each square represents the mean of the centroids of each group ) . in addition , discriminant analysis generated three discriminant canonical functions . analysis of these functions demonstrated statistical significance in all cases when classifying clusters of players . moreover , the first and second functions have greater discriminative power between players classification ( 75% and 15% , respectively ) . as it may be seen in the results section , in the first canonical function ( which was the function that explained more variance ) , brv was the variable with most discriminatory power for the upper trunk , and st20 for the lower trunk . our results support these of previous literature suggesting that brv is often considered as a key performance variable in handball players ( chelly et al . , 2010 ; hermassi et al . , 2011 ; marques and gonzlez  badillo , 2006 ) . in the second canonical discriminating function , this variable has a lot of significance as several studies have shown the importance of maximum dynamic strength to improve 1rmest in handball players . a large correlation coefficient was obtained between brv and 1rmest ( r = .68 , p < .01 ) . as described before , it seems that these variables usually can be described by a positive linear function ( chelly et al . , 2010 ; marques , 2006 ; marques and gonzlez - badillo , 2006 ) . finally , in the last canonical discriminatory function , the best performance indicators to correctly classify the players were ppeakmp , pmppmp , cmj , sj and st10 . as mentioned earlier , no prior study had examined ppeakmp and pmppmp variable in handball players , although marques et al . ( 2007 ) correlated the peak power attained during the bench press ( 26 and 36 kg ) with brv ( r = .63 , p = .017 ) . however , in our study these variables ( pmppmp and ppeakmp ) returned a large coefficient with brv ( r = .55 , p < .05 and r = .54 , p < .05 , respectively ) . the main limitation of our study was the sample size ( n = 16 ) , although due to the selective nature of the elite population used , this was the maximum sample size available . however , our intention was not to generalize the results to the population , but to create a multivariate systematic tool to assess intra - squad physical performance . it can be therefore applied in order to acknowledge which variables require a training focus within the squad for certain individuals . it is interesting to note that this is the first study that uses a multivariate approach to analyse physical performance with a practical application . information obtained in other age groups will help generate a broader database to classify potential players . this will help identify new performance measures that will be able to better discriminate between players whose physical performance is optimal in a given dimension . also future research needs to be conducted in order to make a cluster and discriminant analysis in the specific handball playing positions and moreover , including other types of variables ( i.e. , decision - making and endurance tests ) . this would allow identification of the most relevant strength training variables at each specific playing position . the present study provides new information on specific interactions between physical performance and cluster created solutions . according to performance measured in upper body variables ( 1rmest , ppeakmp , pmppmp and brv ) and others such as 10 m and 20 m sprint velocity , as well as cmj and sj height , we were able to classify 16 young male elite handball players into four groups through cluster analysis . the results revealed that within a young elite handball team there were various physical performance levels and therefore , these athletes required different training focuses based on their individual strengths and weaknesses . finally , a discriminant analysis identified brv , st20 , 1rmest , ppeakmp and pmppmp as the variables that best reflected elite physical performance in the handball players . the information emerging from this study can help researchers , coaches and strength and conditioning coaches objectively evaluate the results of physical fitness tests conducted on the players in their teams . notably , the classification system obtained through cluster and discriminant analyses enabled us to discriminate four groups of athletes according to their performance considering several physical variables rather than just one ( i.e. , a multivariate approach ) . using this method , players are assigned to clusters by physical variables that identify their strengths and weaknesses . as it can be seen in figure 2 , players assigned to cluster 3 were less physically fit and required improvement in the horizontal ( e.g. , st10 and st20 ) and cross - sectional dimensions ( e.g. , 1rmest and brv ) . to train the specific needs of this group , technique and maximum strength sessions could be implemented in their training schedule . whereas , players assigned to cluster 1 needed enhancement in the cross - sectional dimension of the upper trunk only ( i.e. , ppeakmp and pmppmp ) . training sessions based on high execution velocity would be a good choice for this cluster group to improve the variables related to muscular power . therefore , the cluster analysis allows coaches to group training sessions ensuring players will develop specific skills and abilities to enhance their physical performance parameters . finally , our discriminant analysis assigned greater relative weight to some of the performance variables assessed compared to others . based on this type of objective quantifiable information , training sessions could be designed to improve performance in these specific physical measures .