Morphological integration and modularity are closely related concepts about how different traits of an organism are correlated. Integration is the overall pattern of intercorrelation; modularity is the partitioning of integration into evolutionarily or developmentally independent blocks of traits. Modularity and integration are usually studied using quantitative phenotypic data, which can be obtained either from extant or fossil organisms. Many methods are now available to study integration and modularity, all of which involve the analysis of patterns found in trait correlation or covariance matrices. We review matrix correlation, random skewers, fluctuating asymmetry, cluster analysis, Euclidean distance matrix analysis (EDMA), graphical modelling, two-block partial least squares, RV coefficients, and theoretical matrix modelling and discuss their similarities and differences. We also review different coefficients that are used to measure correlations. We apply all the methods to cranial landmark data from and ontogenetic series of Japanese macaques,Macaca fuscatato illustrate the methods and their individual strengths and weaknesses. We conclude that the exploratory approaches (cluster analyses of various sorts) were less informative and less consistent with one another than were the results of model testing or comparative approaches. Nevertheless, we found that competing models of modularity and integration are often similar enough that they are not statistically distinguishable; we expect, therefore, that several models will often be significantly correlated with observed data.