bioRxivEvolutionary Biology
JNplots: an R package to visualize outputs from the Johnson-Neyman technique for categorical and continuous moderators, including options for phylogenetic regressions.
Abstract
The analysis of two-way interactions in linear models is common in the fields of ecology and evolution, being often present in allometric, macroevolutionary, and experimental studies, among others. However, the interpretation of significant interactions can be incomplete when limited to the examination of model coefficients and significance tests. The Johnson-Neyman technique represents a step forward in the interpretation of significant two-way interactions, allowing the user to examine how changes in the moderator variable, it being categorical or continuous, affect the significance of the relationship between the dependent variable and the predictor. Despite its implementation in several software since its initial development, the available options to perform the method lack certain functionality aspects, including the visualization of regions of non-significance when the moderator is categorical, the implementation of phylogenetic corrections, and more intuitive graphical outputs. Here I present the R package JNplots, which aims to fill gaps left by previous software regarding the calculation and visualization of regions of non-significance when fitting two-way interaction models. JNplots includes two basic functions which allow the user to investigate different types of interaction models, including cases where the moderator variable is categorical or continuous. The user can also specify whether the model to explore should be phylogenetically informed and choose a particular phylogenetic correlation structure to be used. Finally, the functions of JNplots produce plots that are largely customizable and allow a more intuitive interpretation of the interaction term. Here I provide a walkthrough on the use of JNplots using three different examples based on empirical data, each representing a different common scenario in which the package can be useful. Additionally, I present the different customization options for the graphical outputs of JNplots.
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- BiMultiNetPlot: An R package for visualizing ecological bipartite multilayer networks
Hai‐Dong Li·bioRxiv·Sep 20, 2024
- hudson: A User-Friendly R Package to Extend Manhattan Plots
Anastasia Lucas, Anurag Verma, Marylyn DeRiggi Ritchie·bioRxiv·Jan 25, 2022
- When can a joint species distribution model capture biotic interactions? A simulation study across interaction mechanisms and levels of model integration
Jeanne Tous, Julien Chiquet·bioRxiv·May 22, 2026
- Biomod2EZ - An R script suite for visualizing projected niche model ensembles and reporting statistical results.
Paul R. Sesink Clee, Stephen Woloszynek, Mary Katherine Gonder·bioRxiv·May 24, 2017
- NLMR and landscapetools: An integrated environment for simulating and modifying neutral landscape models in R
M. Sciaini, M. Fritsch, … C. E. Simpkins·bioRxiv·Apr 27, 2018
- The Poisson-lognormal model as a versatile framework for the joint analysis of species abundances
J. Chiquet, M. Mariadassou, S. Robin·bioRxiv·Oct 7, 2020