The ecological relationships of "pitayo" and "cardón" were studied in the Sayula, Jalisco Basin. A multivariate approach was used considering three sets of explanatory variables: geographical (FG), edaphic (E) and biotic (B). Three columnar cacti species Stenocereus queretaroensis, and S. dumortieri and Pachycereus grandis were found. According to a detrended correspondence analysis the length of the gradient was 3.59 S.D., which justifies the use of unimodal response model for the species data. Also the results of a correspondence analysis showed a sites and species pattern related to a supposed moisture gradient. Due to the presence of outliers in some variables, it was necessary to use certain passive sites in the canonical correspondence analysis. The selected variables were: a] geographical: latitude, longitude, elevation and exposure; b] edaphic: potassium content, cation exchange capacity and sand percentage. Of these components of species data variation, the first explains 28.94% and the second one 23.66%, while the interaction of both explains only 8.23%. At the biotic level, the accompanying shrubs and tree species data were related to "pitayo" and "cardón" attributes in two ways: a] as three separated populations, or b] as a complex. The better selected attributes were total cover and average height of individual in the complex, or else S. dumortieri and P. grandis, respectively, for the species strategy. However, there was a poor relation for the species (17.4%), with the variation of the associated shrub and tree species.