### Ana Bucchi; Andrés Neumann; Fernanda Quevedo-Díaz & Gabriel M. Fonseca

#### The study of the shape variation in geometric morphometrics has an important limitation known as the Pinocchio effect. The Pinocchio effect produces artifactual variances of the landmarks and implies that it is not possible to know the morphological change structure of an object, other than by dividing the landmark sets and then comparing them. This, however, involves making prior assumptions about the pattern of variation of an object. In this study, we provide a code in R to iterate over a complete set of landmarks and test all possible combinations of landmarks until deliver those landmarks associated with the largest to the smallest morphological changes. We tested this on a sample of 28 landmarks in 143 3D models of human skulls. The results indicated that this process can result in a pooled variance of a subset of landmarks that is an order of magnitude larger than that of several other regions of the skull. This method makes it possible to describe the pattern of variation of any 2D or 3D object represented by fixed landmarks, to distinguish the shape features that have more morphological dispersion, and to avoid any aprioristic assumptions about how the morphological changes of an object behave.

#### KEY WORDS: Geometric morphometrics; Local shape variation; Pinocchio effect; Procrustes superimposition; Variance.

#### BUCCHI, A.; NEUMANN, A.; QUEVEDO-DÍAZ, F. & FONSECA, G. M. An R script to estimate which and how many landmarks represent the most variable parts of an object. *Int. J. Morphol., 42(2):*289-293, 2024.