Although not a major issue, we did find a coding error in lm.rrpp.ws, post-production of the article. The coding error did not change the overall conclusion but did change the likelihood ratio statistics reported in Table 1 of the article. Below is an R demonstration of the corrected table.
library(geomorph)data(pupfish.ws)# With phylogenyfit <-extended.pgls(f1 = coords~Species * Sex + Population, data = pupfish.ws, species ="Species",phy = pupfish.ws$phy,delta =0.1, print.progress =FALSE)lr_test(fit, print.progress =FALSE)
Wilks LR Z Pr(>LR)
Species 0.47016 731.3 6.1434 0.001 ***
Sex 0.89465 107.9 2.7601 0.005 **
Population 0.02250 3676.8 27.4201 0.001 ***
Species:Sex 0.42099 838.3 10.5811 0.001 ***
Full.Model 0.00517 5102.1 24.5535 0.001 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
References
Adams, Dean C, and Michael L Collyer. 2024. “Extending Phylogenetic Regression Models for Comparing Within-Species Patterns Across the Tree of Life.”Methods in Ecology and Evolution 15 (12): 2234–46. https://doi.org/10.1111/2041-210X.14438.