I am using mri_segstats to obtain intensity values from each of the 34 cortical parcellations in the Desikan-Killiany Atlas; this provides metrics such as mean, min, max, range, and standard deviation. 

Is there a way to extract information such as the mode, or information to create a histogram for each of the particular regions? 

Thank you, 

Yash