The design of a sampling unit, whether a simple plot or a subplot within a clustered structure, including shape and size, has received little attention in inferential forestry research. The use of auxiliary variables from remote sensing impacts the precision of estimators from both model-assisted and model-based inference perspectives. In both cases, model parameters are estimated from a sample of field plots and information from pixels corresponding to these units. In studies assisted by remote sensing, the shape of the plot used to fit regression models (typically circular) often differs from the shape of the population elements for prediction, where the area of interest is divided into equal tessellated parts. This raises interest in understanding the effect of the sampling unit shape on the mean of variables in forest stands of interest. Therefore, the objective of this study was to evaluate the effect of circular and square subplots, concentrically overlapped and arranged in an inverted Y cluster structure, over tree density, basal area, and aboveground biomass in a managed temperate forest in central Mexico. We used a Multivariate Generalised Linear Mixed Model, which considers the Gamma distribution of the variables and accounts for spatial correlation between Secondary Sampling Units nested within the Primary Sampling Unit. The main findings of this study indicate that the type of secondary sampling unit of the same area and centroid, whether circular or square, does not significantly affect the mean tree density (trees), basal area (m2), and aerial biomass.